Set reference time to CTA reference time in GCTAEventcube, GCTAEventList and GCT…

…ACubeExposure (#2295)

Set reference time to CTA reference time in GCTAEventcube, GCTAEventList and GCT…
…ACubeExposure (#2295)
Jürgen Knödlseder
1 parent e7028dbb
Showing 6 changed files with 21 additions and 10 deletions
ChangeLog
NEWS
doc/source/admin/release_history/1.5.rst
inst/cta/src/GCTACubeExposure.cpp
inst/cta/src/GCTAEventCube.cpp
inst/cta/src/GCTAEventList.cpp
@@ -3,6 +3,7 @@
         * Version 1.5.0 released
         =========================
+        Set reference time to CTA reference time in GCTAEventcube, GCTAEventList and GCTACubeExposure (#2295)
         Add GPha::counts_spectrum() and GPha::backscal_spectrum() methods (#2290)
         Add GNdarray sign(), operator*= and operator/= (#2290)
         Add GSkyMap::counts() and GSkyMap::flux() methods and clip() friend operator (#2290)
@@ -243,6 +243,9 @@ also position fitting is now possible for Fermi-LAT (#604, #2183).
 20. CTA interface
 -----------------
+Set the Good Time Interval time reference to the CTA time reference in
+GCTAEventList, GCTAEventCube and GCTACubeExposure (#2295).
+
 GCTAObservation and GCTAOnOffObservation handle statistic attribute in XML
 file (#2238).
@@ -30,6 +30,9 @@ Bug fixes
 Improvements
 ------------
+* [`2295 <https://cta-redmine.irap.omp.eu/issues/2295>`_] -
+  Set reference time to CTA reference time in ``GCTAEventcube``, ``GCTAEventList``
+  and ``GCTACubeExposure``
 * [`2290 <https://cta-redmine.irap.omp.eu/issues/2290>`_] -
   Add ``GSkyMap::counts()`` and ``GSkyMap::flux()`` methods and ``clip()``
   friend operator, ``GNdarray`` ``sign()``, ``operator*=`` and ``operator/=``,
 /***************************************************************************
  *          GCTACubeExposure.cpp - CTA cube analysis exposure class        *
  * ----------------------------------------------------------------------- *
- *  copyright (C) 2014-2017 by Chia-Chun Lu                                *
+ *  copyright (C) 2014-2018 by Chia-Chun Lu                                *
  * ----------------------------------------------------------------------- *
  *                                                                         *
  *  This program is free software: you can redistribute it and/or modify   *
@@ -33,11 +33,7 @@
 #include "GLog.hpp"
 #include "GObservations.hpp"
 #include "GSkyRegionCircle.hpp"
-#include "GCTACubeExposure.hpp"
-#include "GCTAObservation.hpp"
-#include "GCTAResponseIrf.hpp"
-#include "GCTAEventCube.hpp"
-#include "GCTAEventList.hpp"
+#include "GCTALib.hpp"
 /* __ Method name definitions ____________________________________________ */
 #define G_FILL_CUBE    "GCTACubeExposure::fill_cube(GCTAObservation&, GLog*)"
@@ -583,6 +579,9 @@ void GCTACubeExposure::init_members(void)
     m_wgt_left  = 0.0;
     m_wgt_right = 0.0;
+    // Set CTA time reference for GTIs
+    m_gti.reference(GTimeReference(G_CTA_MJDREF, "s", "TT", "LOCAL"));
+
     // Return
     return;
 }
 /***************************************************************************
  *           GCTAEventCube.cpp  -  CTA event bin container class           *
  * ----------------------------------------------------------------------- *
- *  copyright (C) 2010-2017 by Juergen Knoedlseder                         *
+ *  copyright (C) 2010-2018 by Juergen Knoedlseder                         *
  * ----------------------------------------------------------------------- *
  *                                                                         *
  *  This program is free software: you can redistribute it and/or modify   *
@@ -30,8 +30,7 @@
 #endif
 #include "GTools.hpp"
 #include "GFits.hpp"
-#include "GCTAException.hpp"
-#include "GCTAEventCube.hpp"
+#include "GCTALib.hpp"
 /* __ Method name definitions ____________________________________________ */
 #define G_NAXIS                                   "GCTAEventCube::naxis(int)"
@@ -710,6 +709,9 @@ void GCTAEventCube::init_members(void)
     // Prepare event bin
     init_bin();
+    // Set CTA time reference for GTIs
+    m_gti.reference(GTimeReference(G_CTA_MJDREF, "s", "TT", "LOCAL"));
+
     // Return
     return;
 }
@@ -43,7 +43,7 @@
 #include "GFitsTableULongCol.hpp"
 #include "GTime.hpp"
 #include "GTimeReference.hpp"
-#include "GCTAEventList.hpp"
+#include "GCTALib.hpp"
 #include "GCTASupport.hpp"
 /* __ Method name definitions ____________________________________________ */
@@ -974,6 +974,9 @@ void GCTAEventList::init_members(void)
     m_irf_names.clear();
     m_irf_values.clear();
+    // Set CTA time reference for GTIs
+    m_gti.reference(GTimeReference(G_CTA_MJDREF, "s", "TT", "LOCAL"));
+
     // Return
     return;
 }