adding the possibility to cut on the MFT/MCH chi^2

AliceO2Group · Aug 9, 2024 · 492e52b · 492e52b
1 parent 4c8dac4
commit 492e52b
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/Modules/MUON/Common/src/TracksTask.cxx b/Modules/MUON/Common/src/TracksTask.cxx
@@ -103,6 +103,7 @@ void TracksTask::createTrackHistos(const Activity& activity)
   double cutChi2Max = getConfigurationParameter<double>(mCustomParameters, "cutChi2Max", 1000, activity);
   double nSigmaPDCA = getConfigurationParameter<double>(mCustomParameters, "nSigmaPDCA", 6, activity);
   double matchScoreMaxMFT = getConfigurationParameter<double>(mCustomParameters, "matchScoreMaxMFT", 1000, activity);
+  double matchChi2MaxMFT = getConfigurationParameter<double>(mCustomParameters, "matchChi2MaxMFT", 1000, activity);
   double diMuonTimeCut = getConfigurationParameter<double>(mCustomParameters, "diMuonTimeCut", 100, activity) / 1000;
 
   int etaBins = getConfigurationParameter<int>(mCustomParameters, "etaBins", 200, activity);
@@ -165,6 +166,12 @@ void TracksTask::createTrackHistos(const Activity& activity)
         return false;
       return true;
     },
+    // MFT-MCH chi2 score
+    [matchChi2MaxMFT](const MuonTrack& t) {
+      if (t.hasMFT() && t.hasMCH() && t.getMatchInfoFwd().getMFTMCHMatchingChi2() > matchChi2MaxMFT)
+        return false;
+      return true;
+    },
     // MCH chi2 cut
     [cutChi2Min, cutChi2Max](const MuonTrack& t) { return ((t.getChi2OverNDFMCH() >= cutChi2Min) && (t.getChi2OverNDFMCH() <= cutChi2Max)); }
   };