TPC: make Kr-83m calibration SpecialPhysicsCuts threshold configurable

Author: Ankur Yadav

Detectors/TPC/base/include/TPCBase/ParameterDetector.h

@@ -34,6 +34,7 @@ struct ParameterDetector : public o2::conf::ConfigurableParamHelper<ParameterDet
   float DriftTimeOffset = 7.3;   ///< drift time offset in time bins (we observe ~2.4\mus before October 2023 and ~1.45 \mus after)
   bool ExcludeFCGap = true;      ///< exclude electrons created in the gap between the IFC vessel and OFC vessel and FC strips
   bool UseGeant4Edep = false;    ///< use Geant4 energy deposit directly for ionisation (for Kr-83m calibration runs)
+  float SpecialCutsGeV = 1e-6f;  ///< CUTELE/CUTGAM/DCUTE/BCUTE applied to TPC drift-gas media by SetSpecialPhysicsCuts() when UseGeant4Edep is set (Kr-83m calibration isolation test); default 1e-6 GeV (1 keV) matches the prior hardcoded value
 
   O2ParamDef(ParameterDetector, "TPCDetParam");
 };

Detectors/TPC/simulation/src/Detector.cxx

@@ -3258,11 +3258,12 @@ void Detector::SetSpecialPhysicsCuts()
   LOG(info) << "TPC SetSpecialPhysicsCuts: UseGeant4Edep=" << detParam.UseGeant4Edep;
   if (detParam.UseGeant4Edep) {
     auto& matmgr = o2::base::MaterialManager::Instance();
+    const float specialCut = detParam.SpecialCutsGeV;
     for (int med : {(int)kDriftGas1, (int)kDriftGas2, (int)kCO2}) {
-      matmgr.SpecialCut(GetName(), med, o2::base::ECut::kCUTELE, 1e-6f);
-      matmgr.SpecialCut(GetName(), med, o2::base::ECut::kCUTGAM, 1e-6f);
-      matmgr.SpecialCut(GetName(), med, o2::base::ECut::kDCUTE,  1e-6f);
-      matmgr.SpecialCut(GetName(), med, o2::base::ECut::kBCUTE,  1e-6f);
+      matmgr.SpecialCut(GetName(), med, o2::base::ECut::kCUTELE, specialCut);
+      matmgr.SpecialCut(GetName(), med, o2::base::ECut::kCUTGAM, specialCut);
+      matmgr.SpecialCut(GetName(), med, o2::base::ECut::kDCUTE,  specialCut);
+      matmgr.SpecialCut(GetName(), med, o2::base::ECut::kBCUTE,  specialCut);
     }
   }
   o2::base::Detector::SetSpecialPhysicsCuts();