Adding new metadata variables to manage the global time reset

Author: juanan

inc/TRestGeant4Metadata.h

@@ -117,6 +117,14 @@ class TRestGeant4Metadata : public TRestMetadata {
     /// (fFullChain=true), or just a single decay (fFullChain=false).
     Bool_t fFullChain = true;
 
+    /// \brief Reset global time in case of a long radioactive decay, the global time of
+    /// the decay is stored in the variable offset defined on TRestGeant4Track
+    Bool_t fResetGlobalTime = true;
+
+    /// \brief Time precision, in us, to determine if the global time has to be reset
+    /// or not.
+    Double_t fResetTimePrecision = 0.001;
+
     /// \brief If defined, it will stop the full chain decay simulation when one of these isotope appears.
     std::set<std::string> fFullChainStopIsotopes;
 
@@ -216,6 +224,11 @@ class TRestGeant4Metadata : public TRestMetadata {
     /// \brief Returns true in case full decay chain simulation is enabled.
     inline Bool_t isFullChainActivated() const { return fFullChain; }
 
+    /// \brief Returns true in case global time is reset for long radioactive decays.
+    inline Bool_t isGlobalTimeReset() const { return fResetGlobalTime; }
+
+    inline Double_t GetResetTimePrecision() const { return fResetTimePrecision; }
+
     /// \brief Returns the isotopes that will stop the full chain decay simulation.
     inline std::set<std::string> GetFullChainStopIsotopes() const { return fFullChainStopIsotopes; }
 
@@ -417,7 +430,7 @@ class TRestGeant4Metadata : public TRestMetadata {
     TRestGeant4Metadata(const TRestGeant4Metadata& metadata);
     TRestGeant4Metadata& operator=(const TRestGeant4Metadata& metadata);
 
-    ClassDefOverride(TRestGeant4Metadata, 19);
+    ClassDefOverride(TRestGeant4Metadata, 20);
 
     // Allow modification of otherwise inaccessible / immutable members that shouldn't be modified by the user
     friend class SteppingAction;

src/TRestGeant4Metadata.cxx

@@ -850,6 +850,8 @@ void TRestGeant4Metadata::InitFromConfigFile() {
     Double_t defaultTime = 1. / REST_Units::s;
     fSubEventTimeDelay = GetDblParameterWithUnits("subEventTimeDelay", defaultTime);
 
+    fResetTimePrecision = GetDblParameterWithUnits("resetTimePrecision", defaultTime);
+
     auto nEventsString = GetParameter("nEvents");
     if (nEventsString == PARAMETER_NOT_FOUND_STR) {
         nEventsString = GetParameter("Nevents");  // old name
@@ -872,6 +874,9 @@ void TRestGeant4Metadata::InitFromConfigFile() {
     fRegisterEmptyTracks = ToUpper(GetParameter("registerEmptyTracks", "false")) == "TRUE" ||
                            ToUpper(GetParameter("registerEmptyTracks", "off")) == "ON";
 
+    fResetGlobalTime = ToUpper(GetParameter("resetGlobalTime", "true")) == "TRUE" ||
+                     ToUpper(GetParameter("resetGlobalTime", "on")) == "ON";
+
     ReadGenerator();
     // Detector (old storage) section is processed after initializing geometry info in Detector Construction
     // This allows to use regular expression to match logical or physical volumes etc.
@@ -1508,6 +1513,14 @@ void TRestGeant4Metadata::PrintMetadata() {
     RESTMetadata << "GDML geometry: " << GetGdmlReference() << RESTendl;
     RESTMetadata << "GDML materials reference: " << GetMaterialsReference() << RESTendl;
     RESTMetadata << "Sub-event time delay: " << GetSubEventTimeDelay() << " us" << RESTendl;
+
+      if(isGlobalTimeReset()){
+        RESTMetadata << "Reset global time: enabled" << RESTendl;
+        RESTMetadata << "Reset Time precision " << GetResetTimePrecision() << " us" << RESTendl;
+      } else {
+        RESTMetadata << "Reset global time: disabled" << RESTendl;
+      }
+
     Double_t mx = GetMagneticField().X();
     Double_t my = GetMagneticField().Y();
     Double_t mz = GetMagneticField().Z();
@@ -1703,6 +1716,8 @@ TRestGeant4Metadata& TRestGeant4Metadata::operator=(const TRestGeant4Metadata& m
     fBiasingVolumes = metadata.fBiasingVolumes;
     fMaxTargetStepSize = metadata.fMaxTargetStepSize;
     fSubEventTimeDelay = metadata.fSubEventTimeDelay;
+    fResetTimePrecision = metadata.fResetTimePrecision;
+    fResetGlobalTime = metadata.fResetGlobalTime;
     fFullChain = metadata.fFullChain;
     fSensitiveVolumes = metadata.fSensitiveVolumes;
     fNEvents = metadata.fNEvents;