Allow to switch off the cross-talk parameter in the MID response

Author: dstocco

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

@@ -44,20 +44,22 @@ class ChamberResponse
   /// @param distance Distance between the hit and the current strip
   /// @param cathode Anode or cathode
   /// @param deId Detection element ID
+  /// @param hasCrossTalk If false, allows the function to go to 0 at infinity (default: true)
   /// @param theta Particle impact angle
   /// @return true if the strip is fired
-  inline bool isFired(double prob, double distance, int cathode, int deId, double theta = 0.) const
+  inline bool isFired(double prob, double distance, int cathode, int deId, bool hasCrossTalk = true, double theta = 0.) const
   {
-    return (prob < getFiredProbability(distance, cathode, deId, theta));
+    return (prob < getFiredProbability(distance, cathode, deId, hasCrossTalk, theta));
   }
 
   /// @brief Returns the fired probability
   /// @param distance Distance between the hit and the current strip
   /// @param cathode Anode or cathode
   /// @param deId Detection element ID
+  /// @param hasCrossTalk If false, allows the function to go to 0 at infinity (default: true)
   /// @param theta Particle impact angle
   /// @return The probability that the strip is fired
-  double getFiredProbability(double distance, int cathode, int deId, double theta = 0.) const;
+  double getFiredProbability(double distance, int cathode, int deId, bool hasCrossTalk = true, double theta = 0.) const;
 
   /// @brief Fired probability distribution
   /// @param var Pointer with function variables

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

@@ -38,17 +38,22 @@ ChamberResponse::ChamberResponse(const ChamberResponseParams& params, const Cham
 }
 
 //______________________________________________________________________________
-double ChamberResponse::getFiredProbability(double distance, int cathode, int deId, double theta) const
+double ChamberResponse::getFiredProbability(double distance, int cathode, int deId, bool hasCrossTalk, double theta) const
 {
   /// Get fired probability
 
   // Need to convert the distance from cm to mm
   double distMM = distance * 10.;
-  double parA = mParams.getParA(mHV.getHV(deId));
+  double hv = mHV.getHV(deId);
+  double parA = mParams.getParA(hv);
   double parB = mParams.getParB(cathode, deId);
-  double parC = mParams.getParC(mHV.getHV(deId));
   double costheta = std::cos(theta);
-  return (parC + parA / (parA + costheta * std::pow(distMM, parB))) / (1 + parC);
+  double core = parA / (parA + costheta * std::pow(distMM, parB));
+  if (hasCrossTalk) {
+    double parC = mParams.getParC(hv);
+    return (parC + core) / (1 + parC);
+  }
+  return core;
 }
 
 //______________________________________________________________________________
@@ -76,7 +81,7 @@ double ChamberResponse::firedProbabilityFunction(double* var, double* par)
   mParams.setParA(par[4], par[5]);
   mParams.setParC(par[6], par[7]);
 
-  return getFiredProbability(var[0], cathode, deId, par[2]);
+  return getFiredProbability(var[0], cathode, deId, true, par[2]);
 }
 
 //______________________________________________________________________________