Skip to content

Commit 6e6dba1

Browse files
authored
[PWGEM/Dilepton] update matchingMFT.cxx (#10833)
1 parent fc1ff0a commit 6e6dba1

File tree

1 file changed

+12
-4
lines changed

1 file changed

+12
-4
lines changed

PWGEM/Dilepton/Tasks/matchingMFT.cxx

Lines changed: 12 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -60,7 +60,6 @@ struct matchingMFT {
6060
Configurable<float> maxEtaSA{"maxEtaSA", -2.5, "max. eta acceptance for MCH-MID"};
6161
Configurable<float> minEtaGL{"minEtaGL", -3.6, "min. eta acceptance for MFT-MCH-MID"};
6262
Configurable<float> maxEtaGL{"maxEtaGL", -2.5, "max. eta acceptance for MFT-MCH-MID"};
63-
Configurable<float> minRabsGL{"minRabsGL", 27.6, "min. R at absorber end for global muons (min. eta = -3.6)"}; // std::tan(2.f * std::atan(std::exp(- -3.6)) ) * -505.
6463
Configurable<float> minRabs{"minRabs", 17.6, "min. R at absorber end"};
6564
Configurable<float> midRabs{"midRabs", 26.5, "middle R at absorber end for pDCA cut"};
6665
Configurable<float> maxRabs{"maxRabs", 89.5, "max. R at absorber end"};
@@ -195,9 +194,6 @@ struct matchingMFT {
195194
if (chi2 < 0.f || maxChi2GL < chi2) {
196195
return false;
197196
}
198-
if (rAtAbsorberEnd < minRabsGL || maxRabs < rAtAbsorberEnd) {
199-
return false;
200-
}
201197
} else if (trackType == static_cast<uint8_t>(o2::aod::fwdtrack::ForwardTrackTypeEnum::MuonStandaloneTrack)) {
202198
if (eta < minEtaSA || maxEtaSA < eta) {
203199
return false;
@@ -228,6 +224,18 @@ struct matchingMFT {
228224
// LOGF(info, "mcParticle_MFTMCHMID.pdgCode() = %d, mcParticle_MCHMID.pdgCode() = %d, mcParticle_MFT.pdgCode() = %d", mcParticle_MFTMCHMID.pdgCode(), mcParticle_MCHMID.pdgCode(), mcParticle_MFT.pdgCode());
229225
// LOGF(info, "mcParticle_MFTMCHMID.globalIndex() = %d, mcParticle_MCHMID.globalIndex() = %d, mcParticle_MFT.globalIndex() = %d", mcParticle_MFTMCHMID.globalIndex(), mcParticle_MCHMID.globalIndex(), mcParticle_MFT.globalIndex());
230226

227+
if (fwdtrack.chi2MatchMCHMFT() > maxMatchingChi2MCHMFT) {
228+
return;
229+
}
230+
231+
if (fwdtrack.chi2() < 0.f || maxChi2GL < fwdtrack.chi2()) {
232+
return;
233+
}
234+
235+
if (fwdtrack.rAtAbsorberEnd() < minRabs || maxRabs < fwdtrack.rAtAbsorberEnd()) {
236+
return;
237+
}
238+
231239
if (std::abs(mcParticle_MCHMID.pdgCode()) != 13) { // select true muon
232240
return;
233241
}

0 commit comments

Comments
 (0)