Skip to content

Commit 67d38bc

Browse files
bhagyasahoalibuild
andauthored
[PWGDQ] Fixing McDecision issue for two prong signal (#13142)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent 29e869a commit 67d38bc

File tree

1 file changed

+8
-7
lines changed

1 file changed

+8
-7
lines changed

PWGDQ/Tasks/dqEfficiency_withAssoc.cxx

Lines changed: 8 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -2257,17 +2257,11 @@ struct AnalysisSameEventPairing {
22572257
continue;
22582258
}
22592259
if (sig->CheckSignal(true, t1_raw, t2_raw)) {
2260-
mcDecision |= (static_cast<uint32_t>(1) << isig);
22612260
VarManager::FillPairMC<VarManager::kDecayToMuMu>(t1, t2); // NOTE: This feature will only work for muons
22622261
fHistMan->FillHistClass(Form("MCTruthGenPair_%s", sig->GetName()), VarManager::fgValues);
2263-
if (useMiniTree.fConfigMiniTree) {
2264-
// WARNING! To be checked
2265-
dileptonMiniTreeGen(mcDecision, -999, t1.pt(), t1.eta(), t1.phi(), t2.pt(), t2.eta(), t2.phi());
2266-
}
22672262
}
22682263
}
22692264
}
2270-
isig++;
22712265
}
22722266
}
22732267
// Fill Generated PAIR histograms taking into account selected collisions
@@ -2290,15 +2284,22 @@ struct AnalysisSameEventPairing {
22902284
auto t1_raw = mcTracks.rawIteratorAt(t1.globalIndex());
22912285
auto t2_raw = mcTracks.rawIteratorAt(t2.globalIndex());
22922286
if (t1_raw.reducedMCeventId() == t2_raw.reducedMCeventId()) {
2287+
mcDecision = 0;
2288+
isig = 0;
22932289
for (auto& sig : fGenMCSignals) {
22942290
if (sig->GetNProngs() != 2) { // NOTE: 2-prong signals required here
22952291
continue;
22962292
}
22972293
if (sig->CheckSignal(true, t1_raw, t2_raw)) {
2298-
// mcDecision |= (static_cast<uint32_t>(1) << isig);
2294+
mcDecision |= (static_cast<uint32_t>(1) << isig);
22992295
VarManager::FillPairMC<VarManager::kDecayToMuMu>(t1, t2); // NOTE: This feature will only work for muons
23002296
fHistMan->FillHistClass(Form("MCTruthGenPairSel_%s", sig->GetName()), VarManager::fgValues);
2297+
if (useMiniTree.fConfigMiniTree) {
2298+
// WARNING! To be checked
2299+
dileptonMiniTreeGen(mcDecision, -999, t1.pt(), t1.eta(), t1.phi(), t2.pt(), t2.eta(), t2.phi());
2300+
}
23012301
}
2302+
isig++;
23022303
}
23032304
}
23042305
}

0 commit comments

Comments
 (0)