Flux Monitoring Egret Data Reduction
After an observation:
export FLUX=/Volumes/tpac/vlbi/jlovell/Flux_Density_Monitoring
export PATH="${FLUX}/bin:${PATH}"
- Copy bruce logs from
smerd:/scratch/bruceto$FLUX/bruce_logs cd $FLUX/bruce_logs- Look at the log file to find out when the last observation occurred.
brucelog_summary.pl bruce_07278.log - Run merge2fits on smerd in
/scratch/bruceon last ~day of datamerge2fits -o bruce_07225.fit 22500 22505 - 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/')
