Skip to content

Commit 9761502

Browse files
authored
[PWGLF] Corrected MID cut (#14478)
1 parent 8982531 commit 9761502

File tree

1 file changed

+28
-10
lines changed

1 file changed

+28
-10
lines changed

PWGLF/Tasks/Resonances/kstarqa.cxx

Lines changed: 28 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -2116,14 +2116,14 @@ struct Kstarqa {
21162116
// if (selectionMID(track1, 2)) // Pion misidentified as proton
21172117
// continue;
21182118
}
2119-
rEventSelection.fill(HIST("recMCparticles"), 13.5);
21202119

21212120
if (selectionConfig.isApplyParticleMIDPtDep) {
2122-
if (selectionMIDPtDep(track1, 0)) // Kaon misidentified as pion
2121+
if (selectionMIDPtDep(track1, 1)) // Kaon misidentified as pion
21232122
continue;
2124-
if (selectionMIDPtDep(track2, 1)) // Pion misidentified as kaon
2123+
if (selectionMIDPtDep(track2, 0)) // Pion misidentified as kaon
21252124
continue;
21262125
}
2126+
rEventSelection.fill(HIST("recMCparticles"), 13.5);
21272127

21282128
if (std::abs(track1.rapidity(o2::track::PID::getMass(o2::track::PID::Pion))) > selectionConfig.ctrackRapidity)
21292129
continue;
@@ -2140,15 +2140,33 @@ struct Kstarqa {
21402140
}
21412141

21422142
} else if (selectionConfig.isPDGCheckMC && (track1PDG == PDG_t::kKPlus)) {
2143-
if (!selectionConfig.isapplypTdepPID && !selectionConfig.isapplypTdepPIDTOF && !(selectionPID(track1, 0) && selectionPID(track2, 1))) { // pion and kaon
2143+
if (!selectionConfig.isapplypTdepPID && !selectionConfig.isapplypTdepPIDTOF && !(selectionPID(track1, 1) && selectionPID(track2, 0))) { // pion and kaon
21442144
continue;
2145-
} else if (selectionConfig.isapplypTdepPID && !selectionConfig.isapplypTdepPIDTOF && !(selectionPIDPtDep(track1, 0) && selectionPIDPtDep(track2, 1))) { // pion and kaon
2145+
} else if (selectionConfig.isapplypTdepPID && !selectionConfig.isapplypTdepPIDTOF && !(selectionPIDPtDep(track1, 1) && selectionPIDPtDep(track2, 0))) { // pion and kaon
21462146
continue;
2147-
} else if (!selectionConfig.isapplypTdepPID && selectionConfig.isapplypTdepPIDTOF && !(selectionPIDptDepTOF(track1, 0) && selectionPIDptDepTOF(track2, 1))) { // pion and kaon
2147+
} else if (!selectionConfig.isapplypTdepPID && selectionConfig.isapplypTdepPIDTOF && !(selectionPIDptDepTOF(track1, 1) && selectionPIDptDepTOF(track2, 0))) { // pion and kaon
21482148
continue;
21492149
}
21502150
rEventSelection.fill(HIST("recMCparticles"), 12.5);
21512151

2152+
if (selectionConfig.isApplyParticleMID) {
2153+
if (selectionMID(track2, 1)) // Pion misidentified as kaon
2154+
continue;
2155+
// if (selectionMID(track2, 2)) // Kaon misidentified as proton
2156+
// continue;
2157+
if (selectionMID(track1, 0)) // Kaon misidentified as pion
2158+
continue;
2159+
// if (selectionMID(track1, 2)) // Pion misidentified as proton
2160+
// continue;
2161+
}
2162+
2163+
if (selectionConfig.isApplyParticleMIDPtDep) {
2164+
if (selectionMIDPtDep(track1, 0)) // Kaon misidentified as pion
2165+
continue;
2166+
if (selectionMIDPtDep(track2, 1)) // Pion misidentified as kaon
2167+
continue;
2168+
}
2169+
21522170
if (std::abs(track1.rapidity(o2::track::PID::getMass(o2::track::PID::Kaon))) > selectionConfig.ctrackRapidity)
21532171
continue;
21542172

@@ -2360,21 +2378,21 @@ struct Kstarqa {
23602378
rEventSelection.fill(HIST("recMCparticles"), 11.5);
23612379

23622380
if (selectionConfig.isApplyParticleMID) {
2363-
if (selectionMID(track1, 0)) // Kaon misidentified as pion
2381+
if (selectionMID(track1, 1)) // Pion misidentified as kaon
23642382
continue;
23652383
// if (selectionMID(track1, 2)) // Kaon misidentified as proton
23662384
// continue;
2367-
if (selectionMID(track2, 1)) // Pion misidentified as kaon
2385+
if (selectionMID(track2, 0)) // Kaon misidentified as pion
23682386
continue;
23692387
// if (selectionMID(track2, 2)) // Pion misidentified as proton
23702388
// continue;
23712389
}
23722390
rEventSelection.fill(HIST("recMCparticles"), 12.5);
23732391

23742392
if (selectionConfig.isApplyParticleMIDPtDep) {
2375-
if (selectionMIDPtDep(track1, 0)) // Kaon misidentified as pion
2393+
if (selectionMIDPtDep(track1, 1)) // Pion misidentified as kaon
23762394
continue;
2377-
if (selectionMIDPtDep(track2, 1)) // Pion misidentified as kaon
2395+
if (selectionMIDPtDep(track2, 0)) // Kaon misidentified as pion
23782396
continue;
23792397
}
23802398

0 commit comments

Comments
 (0)