Flux Monitoring Egret Data Reduction

After an observation:

 export FLUX=/Volumes/tpac/vlbi/jlovell/Flux_Density_Monitoring
 export PATH="${FLUX}/bin:${PATH}"
  1. Copy bruce logs from smerd:/scratch/bruce to $FLUX/bruce_logs
  2. cd $FLUX/bruce_logs
  3. Look at the log file to find out when the last observation occurred.
    brucelog_summary.pl bruce_07278.log
  4. Run merge2fits on smerd in /scratch/bruce on last ~day of data
    merge2fits -o bruce_07225.fit 22500 22505
  5. Copy FITS data to $FLUX/Raw_data
    > matlab

When installing scripts for the first time, after unpacking

 File -> Set Path
 select 'script', click [OK]
 [Save]
 [Yes]
 [Save]
 [close]
 >>Fits2Mat
  • choose multiple input files (if necessary)
  • save as yyyy_ddd.mat
 >>ScanPSFit('<filename>',<if1clip>,<if2clip>,<fitmode>,<debug>)
  • where ifNclip is the central width of the scan to keep in deg. 0 = keep all. For IF1 (x-band at Hobart), choose 0.4, IF2, choose 0
  • fitmode should be 0
  • if debug = 1, each scan is shown. press [space] to go to the next one, Ctrl-C to quit.
 usually ScanPSFit('<filename>',0.4,0,0,0)

 >>ScanPSCalibration('filenamefreefit.mat',<FWHM1>,<FWHM2>,<gaincurve>)
  • The FWHM are in arcmin, gaincurve = 1
 X-band 5.6172675
 S-band 20.982175
 e.g ScanPSCalibration('2007_206freefit.mat',5.6172675,20.982175,1)

 >>matsplit('filenamefreefitcal.mat')
 >>CalData=load('2007_186freefitcalflux.mat')
 >>[CalData.FreeFitParams421,CalData.flg]=StandardCal('1934-638',2.94,1,CalData.flg,CalData.FreeFitParams421,'HoBaRt','g')
 >>[CalData.FreeFitParams421,CalData.flg]=StandardCal('1934-638',12.02,2,CalData.flg,CalData.FreeFitParams421,'HoBaRt','g')
 >>ScanExport(CalData.FreeFitParams421,CalData.flg,'','g','data.mat','./out/yyyy_ddd/')