Differences

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

--- pulsar_finding [2018/05/06 21:01] – jrseti
+++ pulsar_finding [2023/10/28 18:36] (current) – external edit 127.0.0.1
@@ Line 152: / Line 152: @@
     40490264.0     1212818.48       1018.69           ATA                 9  AT
-Had to use -nobary to avoid the following error:
+NOTE: Not removing the G3 line first caused the following error:
   Generating barycentric corrections...
@@ Line 159: / Line 159: @@
      path = 'resid2.tmp'
-Now it runs:
+Now it runs (removed the G3 line):
              Pulsar Data Preparation Routine
@@ Line 172: / Line 172: @@
          Num of polns = 2 (summed)
     Center freq (MHz) = 5990
-      Num of channels = 2048
+       Num of channels = 2048
       Sample time (s) = 5.8593748690328e-05
        Spectra/subint = 2400
@@ Line 194: / Line 194: @@
   Writing information to 'outfilename_DM.inf'.
-  Massaging the data ...
+  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%
@@ Line 201: / Line 219: @@
   Simple statistics of the output data:
-             Data points written:  8388608
+               Data points written:  8388608
-           Maximum value of data:  1055518818304.00
+    Bins removed for barycentering:  246
-           Minimum value of data:  804456235008.00
+             Maximum value of data:  1055630491648.00
-              Data average value:  930263638197.76
+             Minimum value of data:  804496343040.00
-         Data standard deviation:  39409907532.20
+                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
-Gives the error:
+          Pulsar Raw-Data Folding Search Routine
+   Used for DM, Period, and P-dot tweaking of PSR candidates.
+                   by Scott M. Ransom
-    Error in chkfopen(): No such file or directory
+  Assuming the data are SIGPROC filterbank format...
-       path = 'themask_rfifind.mask'
+  Reading SIGPROC filterbank data from 1 file:
+    '/data/psrB0329+54_5990mhz_58us_6ants_10min_may5_2018.fil'
-Just like prepdata did before I used the -nobary. Probably not best to use -nobary anyway. There is no -nobary for prefold.
+      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 |}}