Pulse Phase Calculation
(gtpphase Tutorial)

Prerequisites

• Event data file in FT1 format a.k.a. a photon data file. (See Extract LAT Data).

• Orbit file to use for the barycentric correction

• Ephemeris information of suspected pulsation provided in one of the following forms:

• Manually input the source (pulsar) location (for the barycentric correction), pulse frequency or period, and related information

• Automatically extracted from a pulsar ephemerides database, available online.

Sample Files. To try the examples in this section, you can download the following fake data files:

SciTools Reference Pages: If you want the SciTools Reference pages to open in a popup window, click on: Open Popup Reference Window.

How To Assign Pulse Phases

There are several ways to input pulsar ephmeris information to the gtpphase tool; it can:

• Extract a series of ephemerides from an existing pulsar ephemerides database file, or

• Accept a single ephemeris entered manually by the user.

In each case, the tool does not produce a separate output file, but operates directly on the input file in place.

Note: If execution is successful, the tool produces no textual output.

Using A Pulsar Ephemerides Database File

The simplest way to run the gtpphase tool, which may suffice for most applications, is for the tool to extract the information from a master pulsar database. The exact process the tool uses for choosing ephemerides will be described below.

The following plot of the modified events file was produced using the Ftools fhisto and fplot, and shows the distribution of counts in the PULSE_PHASE column in the events file following the phase assignment performed above.

To an observer familiar with this pulsar from other observations, the characteristic shape of the pulse is easily recognizable.

Inputing A Pulsar Ephemeris Manually

Another option is to supply the essential ephemeris information manually. For example, this may be useful when an ephemeris is not available in a database, but the necessary information is available in the literature. In the example shown below, the same values which were extracted from the database in the previous example were entered manually.

The result of this run is identical to that of the previous example.

Getting More Information on Computation

If gtpphase is invoked with chatter parameter set to 4 or larger, it will display brief summaries of loaded ephemerides, arrival time corrections, and time systems being used. This feature shows you additional information on computations being performed, which may help you understand the tool. Also, the extra information printed may contain a clue to diagnose certain types of problem. So, try a high chatter when you are unsure on the tool's behavior.