Differences

This shows you the differences between two versions of the page.

--- pulsar_finding [2018/05/06 19:33] – jrseti
+++ pulsar_finding [2023/10/28 18:36] (current) – external edit 127.0.0.1
@@ Line 8: / Line 8: @@
   * [[https://www.cv.nrao.edu/~sransom/PRESTO_search_tutorial.pdf|Presto tutorial]]
   * [[http://www.atnf.csiro.au/research/pulsar/psrcat/|Pulsar values lookup]]
+  * [[http://www.atnf.csiro.au/people/pulsar/tempo/ref_man_sections/obsys.txt|Tempo observatory file format]]
 ==B0329+54 May 05, 2018==
@@ Line 142: / Line 143: @@
 readline reports 10252288 Spectra per file. I chose the closest power of 2 as the -numout 8388608
-prepdata -dm 26.7641 -o outfilename_DM -numout 8388608 -mask themask_rfifind.mask /data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018.fil
+prepdata -nobary -dm 26.7641 -o outfilename_DM -numout 8388608 -mask themask_rfifind.mask /data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018.fil
+NOTE! Had to edit ~/presto/tempo/obsys.dat and add the ATA, otherwise tempo does not work
+Replaced:
+.30     795056.36        893.7            GB 300FT            9  G3
+    with
+    40490264.0     1212818.48       1018.69           ATA                 9  AT
+NOTE: Not removing the G3 line first caused the following error:
+  Generating barycentric corrections...
+  Error in chkfopen(): No such file or directory
+     path = 'resid2.tmp'
+Now it runs (removed the G3 line):
+             Pulsar Data Preparation Routine
+      Type conversion, de-dispersion, barycentering.
+                   by Scott M. Ransom
+  Assuming the data are SIGPROC filterbank format...
+  Reading SIGPROC filterbank data from 1 file:
+    '/data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018.fil'
+      Number of files = 1
+         Num of polns = 2 (summed)
+    Center freq (MHz) = 5990
+       Num of channels = 2048
+      Sample time (s) = 5.8593748690328e-05
+       Spectra/subint = 2400
+     Total points (N) = 10252288
+       Total time (s) = 600.71998657286
+       Clipping sigma = 6.000
+     Invert the band? = False
+            Byteswap? = False
+       Remove zeroDM? = False
+  File  Start Spec   Samples     Padding        Start MJD
+  ----  ----------  ----------  ----------  --------------------
+              0    10252288           0  58243.71099537036935
+  This is the ATA!
+  Read mask information from 'themask_rfifind.mask'
+  Attempting to read the data statistics from 'themask_rfifind.stats'...
+  ...succeded.  Set the padding values equal to the mid-80% channel averages.
+  Writing output data to 'outfilename_DM.dat'.
+  Writing information to 'outfilename_DM.inf'.
+  Generating barycentric corrections...
+     Average topocentric velocity (c) = 2.931167e-05
+     Maximum topocentric velocity (c) = 2.931716e-05
+     Minimum topocentric velocity (c) = 2.930619e-05
+  Collecting and barycentering /data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018.fil...
+  The DM of 26.76 at the barycentric observing freq of 6042.580 MHz
+   causes a delay of 0.003042 seconds compared to infinite freq.
+   This delay is removed from the barycented times.
+  Topocentric epoch (at data start) is:
+.71099537037
+  Barycentric epoch (infinite obs freq at data start) is:
+.70731137398
+  Massaging the data...
+  Amount Complete = 100%
+  Done.
+  Simple statistics of the output data:
+               Data points written:  8388608
+    Bins removed for barycentering:  246
+             Maximum value of data:  1055630491648.00
+             Minimum value of data:  804496343040.00
+                Data average value:  930264209988.64
+           Data standard deviation:  39408852841.53
+prepfold -dm 26.7641 -p 0.714519699726 /data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018.fil -mask themask_rfifind.mask -nsub 2048
+          Pulsar Raw-Data Folding Search Routine
+   Used for DM, Period, and P-dot tweaking of PSR candidates.
+                   by Scott M. Ransom
+  Assuming the data are SIGPROC filterbank format...
+  Reading SIGPROC filterbank data from 1 file:
+    '/data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018.fil'
+      Number of files = 1
+         Num of polns = 2 (summed)
+    Center freq (MHz) = 5990
+      Num of channels = 2048
+      Sample time (s) = 5.8593748690328e-05
+       Spectra/subint = 2400
+     Total points (N) = 10252288
+       Total time (s) = 600.71998657286
+       Clipping sigma = 6.000
+     Invert the band? = False
+            Byteswap? = False
+       Remove zeroDM? = False
+  File  Start Spec   Samples     Padding        Start MJD
+  ----  ----------  ----------  ----------  --------------------
+              0    10252288           0  58243.71099537036935
+  Read mask information from 'themask_rfifind.mask'
+  Attempting to read the data statistics from 'themask_rfifind.stats'...
+  ...succeded.  Set the padding values equal to the mid-80% channel averages.
+  This is the ATA!
+  Folding a 714.52ms_Cand candidate.
+  Output data file is '/data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018_714.52ms_Cand.pfd'.
+  Output plot file is '/data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018_714.52ms_Cand.pfd.ps'.
+  Best profile is in  '/data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018_714.52ms_Cand.pfd.bestprof'.
+  Folding (topo) epoch  (MJD)  =  58243.710995370369
+  Folding (bary) epoch  (MJD)  =  58243.707311373975
+  Data pt duration (dt)   (s)  =  5.85937486903e-05
+  Total number of data points  =  10137600
+  Number of profile bins       =  64
+  Folding period          (s)  =  0.714519699726
+  Folding frequency      (hz)  =  1.39954153872
+  True starting fraction       =  0
+  Generating barycentric corrections...
+  The average topocentric velocity is 2.93118e-05 (units of c).
+  Barycentric folding frequency    (hz)  =  1.39954153872
+  Barycentric folding f-dot      (hz/s)  =  0
+  Barycentric folding f-dotdot (hz/s^2)  =  0
+  Topocentric folding frequency    (hz)  =  1.39950050823
+  Topocentric folding f-dot      (hz/s)  =  2.374989e-11
+  Topocentric folding f-dotdot (hz/s^2)  =  -1.8167136e-14
+  Will search 385 DMs from 0.000 to 73046.868 (ddm = 189.7321)
+  Starting work on 'psrB0329+54_5990mhz_58us_6ants_10min_may5_2018.fil'...
+    Folded 10137600 points of 10137600
+  Optimizing...
+    Searching 385 DMs, 257 periods, and 257 p-dots...
+    Amount Complete = 100%  Done searching.
+  Maximum reduced chi-squared found  =  13.49789
+  Best DM     (pc cm^-3)  =  53314.7269
+  Best topocentric period        (s)  =   0.7161020(34)
+  Best topocentric p-dot       (s/s)  =   1.81(44)x10^-7
+  Best topocentric p-dotdot  (s/s^2)  =  0.0(4.8)x10^-10
+  Best barycentric period        (s)  =   0.7160810(34)
+  Best barycentric p-dot       (s/s)  =   1.81(44)x10^-7
+  Best barycentric p-dotdot  (s/s^2)  =  0.0(4.8)x10^-10
+  Making plots.
+  Effective number of DOF = 60.48
+  Off-pulse Reduced chi^2 correction factor = 0.27
+   Type <RETURN> for next page:
+  Done.
+{{ :psrb0329_54_may05_2018_try1.png?600 |}}