Get MID HV values within the run limits

Author: dstocco

Detectors/MUON/MID/Simulation/include/MIDSimulation/ChamberHV.h

@@ -17,6 +17,7 @@
 #define O2_MID_CHAMBERHV_H
 
 #include <array>
+#include <limits>
 #include <vector>
 #include <unordered_map>
 
@@ -43,7 +44,7 @@ class ChamberHV
 
   /// \brief Sets the HV from the DCS data points
   /// \param dpMap Map with DCS data points
-  void setHV(const std::unordered_map<o2::dcs::DataPointIdentifier, std::vector<o2::dcs::DataPointValue>>& dpMap);
+  void setHV(const std::unordered_map<o2::dcs::DataPointIdentifier, std::vector<o2::dcs::DataPointValue>>& dpMap, uint64_t startTS = 0, uint64_t endTS = std::numeric_limits<uint64_t>::max());
 
  private:
   std::array<double, detparams::NDetectionElements> mHV; ///< High voltage values

Detectors/MUON/MID/Simulation/include/MIDSimulation/ChamberResponse.h

@@ -17,6 +17,7 @@
 #ifndef O2_MID_CHAMBERRESPONSE_H
 #define O2_MID_CHAMBERRESPONSE_H
 
+#include <limits>
 #include <unordered_map>
 
 #include "DetectorsDCS/DataPointIdentifier.h"
@@ -70,7 +71,7 @@ class ChamberResponse
 
   /// @brief Sets the HV from the DCS data points
   /// @param dpMap Map with DCS data points
-  inline void setHV(const std::unordered_map<o2::dcs::DataPointIdentifier, std::vector<o2::dcs::DataPointValue>>& dpMap) { mHV.setHV(dpMap); }
+  inline void setHV(const std::unordered_map<o2::dcs::DataPointIdentifier, std::vector<o2::dcs::DataPointValue>>& dpMap, uint64_t startTS = 0, uint64_t endTS = std::numeric_limits<uint64_t>::max()) { mHV.setHV(dpMap, startTS, endTS); }
 
  private:
   ChamberResponseParams mParams; ///< Chamber response parameters

Detectors/MUON/MID/Simulation/src/ChamberHV.cxx

@@ -41,23 +41,27 @@ std::vector<std::pair<uint64_t, double>> getValues(const std::unordered_map<o2::
   return values;
 }
 
-double getAverage(std::vector<std::pair<uint64_t, double>>& values)
+double getAverage(std::vector<std::pair<uint64_t, double>>& values, uint64_t startTS, uint64_t endTS)
 {
   double num = 0., den = 0.;
   for (size_t ival = 1; ival < values.size(); ++ival) {
-    double delta = values[ival].first - values[ival - 1].first;
+    auto ts = values[ival - 1].first;
+    if (ts < startTS || ts > endTS) {
+      continue;
+    }
+    double delta = values[ival].first - ts;
     num += values[ival - 1].second * delta;
     den += delta;
   }
   return (den > 0.) ? num / den : num;
 }
 
-void ChamberHV::setHV(const std::unordered_map<o2::dcs::DataPointIdentifier, std::vector<o2::dcs::DataPointValue>>& dpMap)
+void ChamberHV::setHV(const std::unordered_map<o2::dcs::DataPointIdentifier, std::vector<o2::dcs::DataPointValue>>& dpMap, uint64_t startTS, uint64_t endTS)
 {
   for (int deId = 0; deId < detparams::NDetectionElements; ++deId) {
     auto alias = detElemId2DCSAlias(deId, dcs::MeasurementType::HV_V);
     auto values = getValues(dpMap, alias);
-    auto hv = getAverage(values);
+    auto hv = getAverage(values, startTS, endTS);
     setHV(deId, hv);
   }
 }

Steer/DigitizerWorkflow/src/MIDDigitizerSpec.cxx

@@ -10,6 +10,8 @@
 // or submit itself to any jurisdiction.
 
 #include "MIDDigitizerSpec.h"
+
+#include <limits>
 #include "TChain.h"
 #include "Framework/CCDBParamSpec.h"
 #include "Framework/ConfigParamRegistry.h"
@@ -64,7 +66,8 @@ class MIDDPLDigitizerTask : public o2::base::BaseDPLDigitizer
     }
     if (matcher == ConcreteDataMatcher(header::gDataOriginMID, "CHAMBER_HV", 0)) {
       auto* dpMap = static_cast<DPMAP*>(obj);
-      mDigitizer->getChamberResponse().setHV(*dpMap);
+      LOG(debug) << "Setting HV with startTS: " << mStartTS << " endTS: " << mEndTS;
+      mDigitizer->getChamberResponse().setHV(*dpMap, mStartTS, mEndTS);
       return;
     }
   }
@@ -79,6 +82,8 @@ class MIDDPLDigitizerTask : public o2::base::BaseDPLDigitizer
 
     // read collision context from input
     auto context = pc.inputs().get<o2::steer::DigitizationContext*>("collisioncontext");
+    mStartTS = context->getGRP().getTimeStart();
+    mEndTS = context->getGRP().getTimeEnd();
     // Triggers reading from CCDB
     pc.inputs().get<std::vector<ChEffCounter>*>("mid_ch_eff");
     pc.inputs().get<DPMAP*>("mid_ch_hv");
@@ -148,6 +153,8 @@ class MIDDPLDigitizerTask : public o2::base::BaseDPLDigitizer
   // RS: at the moment using hardcoded flag for continuos readout
   o2::parameters::GRPObject::ROMode mROMode = o2::parameters::GRPObject::CONTINUOUS; // readout mode
   ElectronicsDelay mElectronicsDelay;                                                // Electronics delay
+  uint64_t mStartTS = 0;                                                             // Start timestamp of run
+  uint64_t mEndTS = std::numeric_limits<uint64_t>::max();                            // End timestamp of run
 };
 
 o2::framework::DataProcessorSpec getMIDDigitizerSpec(int channel, bool mctruth)