Update 3fhl data

gammapy · May 18, 2017 · 34848e4 · 34848e4
1 parent 43e9c77
commit 34848e4
Show file tree

Hide file tree

Showing 8 changed files with 15 additions and 14 deletions.
diff --git a/2fhl/fermi_2fhl_data_prepare.sh b/2fhl/fermi_2fhl_data_prepare.sh
@@ -35,9 +35,9 @@ PSF=fermi_2fhl_psf_gc.fits
 # gtmktime scfile=$SPACECRAFT filter="(DATA_QUAL>0)&&(LAT_CONFIG==1)" \
 #          roicut=no evfile=$EVENTS_SELECTED outfile=$EVENTS
 
-# # # Compute livetime cube
-# # gtltcube zmax=$ZMAX evfile=$EVENTS_SELECTED_GTI scfile=$SPACECRAFT \
-# #          outfile=$LIVETIME dcostheta=0.025 binsz=1
+# Compute livetime cube
+gtltcube zmax=$ZMAX evfile=$EVENTS scfile=$SPACECRAFT \
+         outfile=$LIVETIME dcostheta=0.025 binsz=1
 
 # Counts cube
 gtbin evfile=$EVENTS outfile=$COUNTS algorithm=HEALPIX SCFILE=$SPACECRAFT \

diff --git a/3fhl/background_norm/background_norm.py b/3fhl/background_norm/background_norm.py
@@ -45,4 +45,4 @@
 if SMOOTH:
   data = hp.smoothing(map_in=data, sigma= Angle('5deg').rad)
 
-hp.write_map('bkg_corr_3fhl.fits', m=data, nest=NEST, coord='G')
+hp.write_map('background_norm_hpx.fits', m=data, nest=NEST, coord='G')
diff --git a/3fhl/fermi_3fhl_counts_cube_hpx.fits.gz b/3fhl/fermi_3fhl_counts_cube_hpx.fits.gz
diff --git a/3fhl/fermi_3fhl_data_config.yaml b/3fhl/fermi_3fhl_data_config.yaml
@@ -9,3 +9,4 @@ filenames:
   livetime  : fermi_3fhl_livetime_cube.fits.gz
   counts    : fermi_3fhl_counts_cube_hpx.fits.gz
   isodiff   : ../isodiff/iso_P8R2_SOURCE_V6_v06.txt
+  bkg_norm  : bkg_norm/bkg_norm_image_hpx.fits.gz
diff --git a/3fhl/fermi_3fhl_data_prepare.sh b/3fhl/fermi_3fhl_data_prepare.sh
@@ -11,8 +11,8 @@ EBINALG=LOG
 NEBINS=17
 
 # Time cuts taken from 3FHL paper
-TMIN=239557417
-TMAX=460250000
+TMIN=239557417.
+TMAX=460250000.
 
 # Raw data files
 EVENTS_LIST=$FERMI_FT1_FILE # Fermi FT1 file environment variable
@@ -26,17 +26,17 @@ EXPOSURE=fermi_3fhl_exposure_cube_hpx.fits
 COUNTS=fermi_3fhl_counts_cube_hpx.fits
 PSF=fermi_3fhl_psf_gc.fits
 
-# Merge weekly photon files and apply energy cut
-gtselect infile=$EVENTS_LIST outfile=$EVENTS_SELECTED ra=0 dec=0 rad=180 \
-         tmin=$TMIN tmax=$TMAX emin=$EMIN emax=$EMAX zmax=$ZMAX \
-         evclass=$EVENT_CLASS evtype=$EVENT_TYPE
+## Merge weekly photon files and apply energy cut
+#gtselect infile=$EVENTS_LIST outfile=$EVENTS_SELECTED ra=0 dec=0 rad=180 \
+#         tmin=$TMIN tmax=$TMAX emin=$EMIN emax=$EMAX zmax=$ZMAX \
+#         evclass=$EVENT_CLASS evtype=$EVENT_TYPE
 
-# Update GTI list
-gtmktime scfile=$SPACECRAFT filter="(DATA_QUAL>0)&&(LAT_CONFIG==1)" \
-         roicut=no evfile=$EVENTS_SELECTED outfile=$EVENTS
+## Update GTI list
+#gtmktime scfile=$SPACECRAFT filter="(DATA_QUAL>0)&&(LAT_CONFIG==1)" \
+#         roicut=no evfile=$EVENTS_SELECTED outfile=$EVENTS
 
 # Compute livetime cube
-gtltcube zmax=$ZMAX evfile=$EVENTS_SELECTED_GTI scfile=$SPACECRAFT \
+gtltcube zmax=$ZMAX evfile=$EVENTS scfile=$SPACECRAFT \
          outfile=$LIVETIME dcostheta=0.025 binsz=1
 
 # Counts cube

diff --git a/3fhl/fermi_3fhl_exposure_cube_hpx.fits.gz b/3fhl/fermi_3fhl_exposure_cube_hpx.fits.gz
diff --git a/3fhl/fermi_3fhl_livetime.fits.gz b/3fhl/fermi_3fhl_livetime.fits.gz
diff --git a/3fhl/fermi_3fhl_psf_gc.fits.gz b/3fhl/fermi_3fhl_psf_gc.fits.gz