def METHOD_NAME():
    # locate the k-space raw data file
    input_file = existing_filepath(data_path, data_file)
    # acquisition data will be read from an HDF file input_file
    # AcquisitionData.set_storage_scheme('memory')
    acq_data = AcquisitionData(input_file)
    print('---\n acquisition data norm: %e' % acq_data.norm())
    # pre-process acquisition data only for cartesian readouts
    if trajtype != 'radial' and trajtype != 'goldenangle':
        print('---\n pre-processing acquisition data...')
        processed_data  = preprocess_acquisition_data(acq_data)
    else:
        processed_data = acq_data
    #set the trajectory
    print('---\n setting the trajectory...')
    if trajtype == 'cartesian':
        pass
    elif trajtype == 'grpe':
        processed_data = set_grpe_trajectory(processed_data)
    elif trajtype == 'radial':
        processed_data = set_radial2D_trajectory(processed_data)
    elif trajtype == 'goldenangle':
        processed_data = set_goldenangle2D_trajectory(processed_data)
    else:
        raise NameError('Please submit a trajectory name of the following list: (cartesian, grpe, radial). You gave {}'\
                        .format(trajtype))
    # sort processed acquisition data;
    print('---\n sorting acquisition data...')
    processed_data.sort()
    if run_recon:
        recon = ConjugateGradient(processed_data, num_iter = 20, stop_criterion = 1e-7)
        if show_plot:
            recon.show(title = 'Reconstructed images using CG() (magnitude)')
    else:
        print('---\n Skipping non-cartesian code...')