precompute the resolution map based on input and configured resolution value

Author: GiorgioAlbertoLucia

Detectors/Upgrades/ALICE3/IOTOF/simulation/include/IOTOFSimulation/Digitizer.h

@@ -95,6 +95,8 @@ class Digitizer : public TObject
   /// \param mapName Name of the histogram inside the ROOT file (ignored for CCDB)
   void loadMap(TH2D*& map, const std::string& path, const char* mapName);
 
+  void prepareScaledResolutionMap();
+
   /// Check if the hit passes efficiency cut
   /// \param x Detector local coordinate x in cm with respect to the center of the sensitive volume.
   /// \param z Detector local coordinate z in cm with respect to the center of the sensitive volume.
@@ -120,6 +122,7 @@ class Digitizer : public TObject
   const o2::iotof::GeometryTGeo* mGeometry = nullptr; ///< IOTOF geometry
   TH2D* mEfficiencyMap = nullptr;                     ///< Efficiency map for the detector
   TH2D* mResolutionMap = nullptr;                     ///< Resolution map for the detector
+  TH2D* mSaledResolutionMap = nullptr;                ///< Scaled resolution map for the detector
   TH2D* mTimeOfArrivalMap = nullptr;                  ///< Time of arrival map for the detector
 
   std::vector<o2::iotof::Chip> mChips;                                               //! Chips in the detector, indexed by chip ID

Detectors/Upgrades/ALICE3/IOTOF/simulation/src/Digitizer.cxx

@@ -81,6 +81,7 @@ void Digitizer::init()
   }
 
   loadMap(mResolutionMap, digitizerParams.resolutionMapPath, "hResolutionMap");
+  prepareScaledResolutionMap();
   if (!mResolutionMap) {
     LOG(info) << "No resolution map loaded, using uniform time resolution: " << digitizerParams.timeResolution * 1e3 << " ps";
   }
@@ -432,6 +433,33 @@ void Digitizer::loadMap(TH2D*& map, const std::string& path, const char* mapName
   }
 }
 
+void Digitizer::prepareScaledResolutionMap()
+{
+  if (!mResolutionMap) {
+    LOG(warn) << "No resolution map available to prepare scaled resolution map.";
+    return;
+  }
+
+  const auto& digitizerParams = o2::iotof::DPLDigitizerParam::Instance();
+  const float nominalTimeResolution = digitizerParams.timeResolution;
+  const float minimumResolution = mResolutionMap->GetMinimum();
+  if (minimumResolution <= 0) {
+    LOG(warn) << "Minimum resolution in the map is non-positive, cannot prepare scaled resolution map.";
+    return;
+  }
+
+  mResolutionScalingMap = dynamic_cast<TH2D*>(mResolutionMap->Clone("hScaledResolutionMap"));
+  mResolutionScalingMap->SetDirectory(nullptr); // Detach from file to avoid deletion when file is closed
+
+  for (int binX = 1; binX <= mResolutionScalingMap->GetNbinsX(); ++binX) {
+    for (int binY = 1; binY <= mResolutionScalingMap->GetNbinsY(); ++binY) {
+      float originalValue = mResolutionScalingMap->GetBinContent(binX, binY);
+      float scalingValue = originalValue / minimumResolution;
+      mResolutionScalingMap->SetBinContent(binX, binY, scalingValue * nominalTimeResolution);
+    }
+  }
+}
+
 //_______________________________________________________________________
 bool Digitizer::isEfficient(const float x, const float z) const
 {