|Ask a Question||Search PSRCHIVE:|
PSRCHIVE user documentation: psrspa
1.0 PurposeThe PSR Single Pulse Analysis (psrspa) is designed to facilitate various analysis of single pulse data. It consists of three main parts:
2.1 Basicspsrspa supports several standard options such as preprocessing of the input archives using the "-j" or "-J" switches, meta file containing a list of input archives ("-M"), verbosity settings ("-v" and "-V") as well as "-h" which will print a summary of available options.
2.2 ModesThe mode in which psrspa is to be run needs to be chosen, i.e., at least one of the flags "-hf", "-hd", "-ha", "-hl", "-hc", "-fm" or "-sla" needs to be set. The first two options correspond to the first part of the program. Please see the spa documentation for more information. All of the "-hX" options are used to calculate phase resolved histograms. Finally, the "-fm" will find maximum amplitudes in chosen phase range.
2.3 Phase RangeThe phase range choice is also affecting the phase resolved histograms. Multiple phase ranges can be specified with the "-pr" option as a list of comma seperate beginnings and ends of a phase.
2.4 Output Files Naming and FormatThe resulting files with the chosen phase resolved histograms are stored in files named in a following way:
path/prefix_type_bin.extAll of the above with the exception of type can be specified by the user (see the section "Output of the histograms" in the builtin help). The value of type is chosen by psrspa as poldeg polang or flux for the polarisation degree, its angle or phase resolved flux, respectively. The format of the output files is following:
bin_start bin_end countThe histograms are not normalised and all the post-processing has to be performed by the user. Naming scheme of the files with maximum flux in the chosen phase range is similar to the previously specified:
path/prefix_type.extwith type set to maxamp. The format of these files is compatible with the GNU Scientific Library:
file subint bscrunch_factor phase_min phase_max max_bin max_amplitudewhere the bscrunch_factor specifies wheter the archive was bscrunched before the calculation of the following values when using automatic looped bscrunching with the "-bl" option.
3.0 AlgorithmsThe part of psrspa corresponding to the original spa is directly based on that code. The phase resolved histograms and search for maximal phase are based on the
4.0 Testing and examplesTo perform all the calculations with the standard settings simply run:
psrspa -hf -hd -ha -hc -hl -fm -sla input_file.arThe "input_file.ar" can be replaced by multiple files, regular expression or a metafile containing list of archives. To calculate only the polarisation degree for pulse profile phases falling between 0.45 and 0.50 as well as between 0.60 to 0.70 run:
psrspa -hd -pr 0.45,0.50,0.60,0.70 input_file.arNote the format of the phase range specification: each phase range has to be identified by two numbers, the beginning and end of the phase range. Multiple ranges can be specified by simply being separated by commas. When no phase range is chosen, the full range is used.
5.0 Known bugs and features that require implementation