diff --git a/Modules/MUON/Common/src/TracksTask.cxx b/Modules/MUON/Common/src/TracksTask.cxx index 626f52e84..d4788c3d8 100644 --- 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(mCustomParameters, "cutChi2Max", 1000, activity); double nSigmaPDCA = getConfigurationParameter(mCustomParameters, "nSigmaPDCA", 6, activity); double matchScoreMaxMFT = getConfigurationParameter(mCustomParameters, "matchScoreMaxMFT", 1000, activity); + double matchChi2MaxMFT = getConfigurationParameter(mCustomParameters, "matchChi2MaxMFT", 1000, activity); double diMuonTimeCut = getConfigurationParameter(mCustomParameters, "diMuonTimeCut", 100, activity) / 1000; int etaBins = getConfigurationParameter(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)); } };