Esempio n. 1
0
def get_tempo2_result(parfile, timfile, general2=None):
    """This function is to get the results from tempo and write them to a file.
    Parameter
    ---------
    parfile : str
        The file to read parameters.
    timfile : str
        The file to read toas.
    general2 : list/ None
        The values required from tempo2 general2 plugin.
    Return
    ----------
    A file named as parfile name ends with '.tempo2_test' stored residuals in the
    first column, binary delay in the second column and general2 results if
    general2 are provided.
    """
    if not has_tempo2_utils:
        raise ImportError(
            "To get tempo2 general2 results, tempo2_utils are"
            " required. See page"
            " https://github.com/demorest/tempo_utils"
        )
    residuals = t2u.general2(parfile, timfile, ["pre"])["pre"]
    outfile = parfile + ".tempo2_test"
    f = open(outfile, "w")
    outstr = "residuals "

    if general2 is not None and general2 != []:
        tempo2_vals = t2u.general2(parfile, timfile, general2)
        for keys in general2:
            outstr += keys + " "

    outstr += "\n"
    f.write(outstr)
    for ii in range(len(residuals)):
        outstr = longdouble2str(residuals[ii]) + " "
        if general2 is not None:
            for keys in general2:
                outstr += longdouble2str(tempo2_vals[keys][ii]) + " "
        outstr += "\n"
        f.write(outstr)
    f.close()
Esempio n. 2
0
def get_tempo2_result(parfile, timfile, general2=None):
    """This function is to get the results from tempo and write them to a file.
    Parameter
    ---------
    parfile : str
        The file to read parameters.
    timfile : str
        The file to read toas.
    general2 : list/ None
        The values required from tempo2 general2 plugin.
    Return
    ----------
    A file named as parfile name ends with '.tempo2_test' stored residuals in the
    first column, binary delay in the second column and general2 results if
    general2 are provided.
    """
    if not has_tempo2_utils:
        raise ImportError("To get tempo2 general2 results, tempo2_utils are"
                          " required. See page"
                          " https://github.com/demorest/tempo_utils")
    residuals = t2u.general2(parfile, timfile, ['pre'])['pre']
    outfile = parfile + '.tempo2_test'
    f = open(outfile, 'w')
    outstr = 'residuals '

    if general2 is not None and general2 != []:
        tempo2_vals = t2u.general2(parfile, timfile, general2)
        for keys in general2:
            outstr += keys+' '

    outstr += '\n'
    f.write(outstr)
    for ii in range(len(residuals)):
        outstr = longdouble2string(residuals[ii])  +' '
        if general2 is not None:
            for keys in general2:
                outstr += longdouble2string(tempo2_vals[keys][ii])+' '
        outstr += '\n'
        f.write(outstr)
    f.close()
Esempio n. 3
0
def get_tempo2_result(parfile, timfile, general2=None):
    """This function is to get the results from tempo and write them to a file.
    Parameter
    ---------
    parfile : str
        The file to read parameters.
    timfile : str
        The file to read toas.
    general2 : list/ None
        The values required from tempo2 general2 plugin.
    Return
    ----------
    A file named as parfile name ends with '.tempo2_test' stored residuals in the
    first column, binary delay in the second column and general2 results if
    general2 are provided.
    """
    psr = lt.tempopulsar(parfile, timfile)
    residuals = psr.residuals()
    binary_delay = psr.binarydelay()
    outfile = parfile + '.tempo2_test'
    f = open(outfile, 'w')
    outstr = 'residuals  BinaryDelay '

    if general2 is not None and general2 != []:
        if not has_tempo2_utils:
            raise ImportError("To get tempo2 general2 results, tempo2_utils are"
                              " required. See page"
                              " https://github.com/demorest/tempo_utils")
        else:
            tempo2_vals = t2u.general2(parfile, timfile, general2)
            for keys in general2:
                outstr += keys+' '

    outstr += '\n'
    f.write(outstr)
    for ii in range(len(residuals)):
        outstr = longdouble2string(residuals[ii]) + ' ' +longdouble2string(binary_delay[ii]) +' '
        if general2 is not None:
            for keys in general2:
                outstr += longdouble2string(tempo2_vals[keys][ii])+' '
        outstr += '\n'
        f.write(outstr)
    f.close()
Esempio n. 4
0
import matplotlib.pyplot as plt
import tempo2_utils
import astropy.units as u
from pint.residuals import resids
import numpy as np
import os, unittest
import tempo2_utils

# Using Nanograv data B1855
datadir = '../tests/datafile'
parfile = os.path.join(datadir, 'B1855+09_NANOGrav_dfg+12_TAI_FB90.par')
timfile = os.path.join(datadir, 'B1855+09_NANOGrav_dfg+12.tim')

# libstempo calculation
print "tempo2 calculation"
tempo2_vals = tempo2_utils.general2(parfile, timfile, ['pre'])
# Build PINT model
print "PINT calculation"
mdd = mb.get_model(parfile)
# Get toas to pint
toas = toa.get_TOAs(timfile, planets=False, ephem='DE405')
tt = toas.table
# Get residuals
t2_resids = tempo2_vals['pre']
presids_us = resids(toas, mdd).time_resids
# Plot residuals
plt.errorbar(toas.get_mjds(high_precision=False),
             presids_us.value,
             toas.get_errors(),
             fmt='x')
plt.title("%s Pre-Fit Timing Residuals" % mdd.PSR.value)
Esempio n. 5
0
try:
    import tempo2_utils
except ImportError:
    log.error(
        "This example requires tempo_utils, download from: https://github.com/demorest/tempo_utils and 'pip install .'"
    )
    raise

# Using Nanograv data J0623-0200
datadir = "../../tests/datafile"
parfile = os.path.join(datadir, "J0613-0200_NANOGrav_dfg+12_TAI_FB90.par")
timfile = os.path.join(datadir, "J0613-0200_NANOGrav_dfg+12.tim")

# libstempo calculation
print("tempo2 calculation")
tempo2_vals = tempo2_utils.general2(parfile, timfile, ["pre"])
# Build PINT model
print("PINT calculation")
m = mb.get_model(parfile)
# Get toas to pint
toas = toa.get_TOAs(timfile, planets=False, ephem="DE405", include_bipm=False)

t2_resids = tempo2_vals["pre"]
presids_us = resids(toas, m).time_resids.to(u.us)
# Plot residuals
plt.errorbar(toas.get_mjds().value,
             presids_us.value,
             toas.get_errors().value,
             fmt="x")
plt.title("%s Pre-Fit Timing Residuals" % m.PSR.value)
plt.xlabel("MJD")
Esempio n. 6
0
    for ii, tt in enumerate(t.table):
        p = m.phase(tt)
        resids[ii] = p.frac
        ss_roemer[ii] = m.solar_system_geometric_delay(tt)
        ss_shapiro[ii] = m.solar_system_shapiro_delay(tt)

    time_phase = time.time() - t0
    sys.stderr.write("Computed phases in %.3f sec\n" % time_phase)

    # resids in (approximate) us:
    resids_us = resids / float(m.F0.value) * 1e6
    sys.stderr.write("RMS PINT residuals are %.3f us\n" % resids_us.std())

    # Get some general2 stuff
    tempo2_vals = tempo2_utils.general2(parfile, timfile,
                                        ['tt2tb', 'roemer', 'post_phase',
                                         'shapiro', 'shapiroJ'])
    t2_resids = tempo2_vals['post_phase'] / float(m.F0.value) * 1e6
    diff_t2 = resids_us - t2_resids
    diff_t2 -= diff_t2.mean()

    # run tempo1 also, if the tempo_utils module is available
    try:
        import tempo_utils
        t1_toas = tempo_utils.read_toa_file(timfile)
        tempo_utils.run_tempo(t1_toas, t1_parfile)
        t1_resids = t1_toas.get_resids(units='phase') / float(m.F0.value) * 1e6
        diff_t1 = resids_us - t1_resids
        diff_t1 -= diff_t1.mean()

        diff_t2_t1 = t2_resids - t1_resids
Esempio n. 7
0
import matplotlib.pyplot as plt
import tempo2_utils
import astropy.units as u
from pint.residuals import resids
import numpy as np
import os, unittest
import tempo2_utils

# Using Nanograv data B1855
datadir = '../tests/datafile'
parfile = os.path.join(datadir, 'B1855+09_NANOGrav_dfg+12_TAI_FB90.par')
timfile = os.path.join(datadir, 'B1855+09_NANOGrav_dfg+12.tim')

# libstempo calculation
print "tempo2 calculation"
tempo2_vals = tempo2_utils.general2(parfile, timfile,['pre'])
# Build PINT model
print "PINT calculation"
mdd = mb.get_model(parfile)
# Get toas to pint
toas = toa.get_TOAs(timfile, planets=False, ephem='DE405')
tt = toas.table
# Get residuals
t2_resids = tempo2_vals['pre']
presids_us = resids(toas, mdd).time_resids
# Plot residuals
plt.errorbar(toas.get_mjds(high_precision=False), presids_us.value,
            toas.get_errors(), fmt='x')
plt.title("%s Pre-Fit Timing Residuals" % mdd.PSR.value)
plt.xlabel('MJD')
plt.ylabel('Residual (us)')
Esempio n. 8
0
datapath = os.path.join(os.environ['PINT'],'tests','datafile')
# Using Nanograv data B1855
parfile = os.path.join(datapath, 'B1855+09_NANOGrav_dfg+12_modified.par')
timfile = os.path.join(datapath, 'B1855+09_NANOGrav_dfg+12.tim')
# libstempo calculation
print "libstempo calculation"
psr = lt.tempopulsar(parfile, timfile)
# Build PINT model
print "PINT calculation"
mdd = mb.get_model(parfile)
# Get toas to pint
toas = toa.get_TOAs(timfile, planets=True)
tt = toas.table
# Run tempo2 general2 pluging
tempo2_vals = tempo2_utils.general2(parfile, timfile,
                                    ['tt2tb', 'roemer', 'post_phase',
                                     'shapiro', 'shapiroJ','bat','clock0',
                                     'clock1','clock2','clock3','clock4','sat',
                                     'tropo'])
# compute residules
t2_resids = tempo2_vals['post_phase'] / float(mdd.F0.value) * 1e6 * u.us
presids_us = resids(toas, mdd).time_resids.to(u.us)
toas = psr.toas()
toas.sort()
plt.plot(toas,presids_us-t2_resids)
plt.xlabel('Mjd (DAY)')
plt.ylabel('residule (us)')
plt.title('Residule difference between PINT and tempo2')
plt.show()
Esempio n. 9
0
    for ii, tt in enumerate(t.table):
        p = m.phase(tt)
        resids[ii] = p.frac
        ss_roemer[ii] = m.solar_system_geometric_delay(tt)
        ss_shapiro[ii] = m.solar_system_shapiro_delay(tt)

    time_phase = time.time() - t0
    sys.stderr.write("Computed phases in %.3f sec\n" % time_phase)

    # resids in (approximate) us:
    resids_us = resids / float(m.F0.value) * 1e6
    sys.stderr.write("RMS PINT residuals are %.3f us\n" % resids_us.std())

    # Get some general2 stuff
    tempo2_vals = tempo2_utils.general2(
        parfile, timfile,
        ['tt2tb', 'roemer', 'post_phase', 'shapiro', 'shapiroJ'])
    t2_resids = tempo2_vals['post_phase'] / float(m.F0.value) * 1e6
    diff_t2 = resids_us - t2_resids
    diff_t2 -= diff_t2.mean()

    # run tempo1 also, if the tempo_utils module is available
    try:
        import tempo_utils
        t1_toas = tempo_utils.read_toa_file(timfile)
        tempo_utils.run_tempo(t1_toas, t1_parfile)
        t1_resids = t1_toas.get_resids(units='phase') / float(m.F0.value) * 1e6
        diff_t1 = resids_us - t1_resids
        diff_t1 -= diff_t1.mean()

        diff_t2_t1 = t2_resids - t1_resids