Added support to handle GTIs in a separate file in ctselect. Added a test case for this purpose

Michael Mayer
1 parent 3af8c312
Showing 3 changed files with 46 additions and 0 deletions
src/ctselect/ctselect.cpp
test/data/obs_unbinned_gti.xml
test/test_ctools.sh
@@ -1284,6 +1284,17 @@ void ctselect::save_event_list(const GCTAObservation* obs,
             // Save observation into FITS file
             obs->save(outfile, clobber());
+            // If separate GTI file was present than create new filename
+            // and save to new file
+            if (!obs->gtifile().empty()) {
+
+                // Set event output file name
+                std::string outgti = set_outfile_name(obs->gtifile());
+
+                // Save GTI extension to different file
+                obs->save_gti(outgti, true);
+            }
+
             // Copy all extensions other than EVENTS and GTI from the input to
             // the output event list. The EVENTS and GTI extensions are written
             // by the save method, all others that may eventually be present
+<?xml version="1.0" encoding="UTF-8" standalone="no"?>
+<observation_list title="observation list">
+  <observation name="" id="00001" instrument="CTA">
+    <parameter name="EventList"           file="data/crab_events.fits" />
+    <parameter name="GoodTimeIntervals"   file="data/crab_events.fits" />
+    <parameter name="EffectiveArea"       file="irf/irf_file.fits" />
+    <parameter name="PointSpreadFunction" file="irf/irf_file.fits" />
+    <parameter name="EnergyDispersion"    file="irf/irf_file.fits" />
+    <parameter name="Background"          file="irf/irf_file.fits" />
+  </observation>
+</observation_list>
@@ -443,6 +443,30 @@ $ECHO &quot; ok&quot;
 #
+# Test ctselect (XML file)
+# ========================
+$ECHO -n "Test ctselect (observation XML file): "
+$ctselect inobs="data/obs_unbinned_gti.xml" \
+outobs="selected_obs.xml" \
+ra=83.63 \
+dec=22.01 \
+rad=10.0 \
+tmin=0.0 \
+tmax=1800.0 \
+emin=0.1 \
+emax=100.0
+$ECHO -n "."
+if [ -s "selected_obs.xml" ]
+then
+$ECHO -n "."
+else
+$ECHO " selected_obs.xml file is not found"
+exit 1
+fi
+$ECHO " ok"
+
+
+#
 # Test ctlike
 # ===========
 $ECHO -n "Test ctlike: "