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/bruce
to$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/bruce
on 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/')