-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathGenFermiData_P8R2_CLEAN_V15_calore_.sh
executable file
·30 lines (19 loc) · 2 KB
/
GenFermiData_P8R2_CLEAN_V15_calore_.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
cd /data/GCE_sys/
echo "Get a beer. This will take a while..."
echo "running gtselect"
gtselect /data/fermi_data_6-26-15/phfile.txt photons_merged_P8R2_CLEAN_V15_calore.fits ra=INDEF dec=INDEF rad=INDEF tmin=INDEF tmax=INDEF zmax=90 emin=300 emax=500000.0 convtype=-1 evclass=256 clobber=True evtype=INDEF
echo "running gtmktime"
gtmktime scfile=/data/fermi_data_1-8-14/lat_spacecraft_merged.fits filter="DATA_QUAL>0 && LAT_CONFIG==1 && ABS(ROCK_ANGLE)<52" roicut=no evfile=photons_merged_P8R2_CLEAN_V15_calore.fits outfile=photons_merged_cut_P8R2_CLEAN_V15_calore.fits clobber=True
rm photons_merged_P8R2_CLEAN_V15_calore.fits
echo "running gtltcube"
gtltcube evfile=photons_merged_cut_P8R2_CLEAN_V15_calore.fits scfile=/data/fermi_data_1-8-14/lat_spacecraft_merged.fits outfile="cube_P8R2_CLEAN_V15_calore.fits" dcostheta=0.1 binsz=1 zmin=0 zmax=180 clobber=True
# make ebin file
echo "running gtbindef"
gtbindef bintype=E binfile=bin_edges_P8R2_CLEAN_V15_calore.dat outfile=ebins_P8R2_CLEAN_V15_calore.fits energyunits=MeV clobber=True
# make psf file
echo "running gtpsf"
gtpsf expcube="cube_P8R2_CLEAN_V15_calore.fits" outfile=gtpsf_P8R2_CLEAN_V15_calore.fits irfs=P8R2_CLEAN_V6 emin=20 emax=1e6 nenergies=50 clobber=True ra=0 dec=0 thetamax=10 ntheta=200
echo "running gtbin"
gtbin evfile=photons_merged_cut_P8R2_CLEAN_V15_calore.fits scfile=/data/fermi_data_1-8-14/lat_spacecraft_merged.fits outfile="gtbin_P8R2_CLEAN_V15_calore.fits" algorithm=CCUBE nxpix=721 nypix=361 ebinalg=FILE ebinfile=ebins_P8R2_CLEAN_V15_calore.fits coordsys=GAL proj=CAR xref=0 yref=0 axisrot=0 binsz=0.5 clobber=True
echo "running gtexpcube2"
gtexpcube2 infile="cube_P8R2_CLEAN_V15_calore.fits" cmap="gtbin_P8R2_CLEAN_V15_calore.fits" coordsys=GAL outfile="gtexpcube2_P8R2_CLEAN_V15_calore.fits" irf=P8R2_CLEAN_V6 ebinfile=ebins_P8R2_CLEAN_V15_calore.fits ebinalg=FILE clobber=True