Skip to content

Commit d25a5e4

Browse files
feisenhufeisenhu
committed
[PWGDQ] add LS histograms to output of MC selected signals
1 parent aa69c6b commit d25a5e4

1 file changed

Lines changed: 25 additions & 11 deletions

File tree

‎PWGDQ/Tasks/dqEfficiency_withAssoc.cxx‎

Lines changed: 25 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -1208,6 +1208,8 @@ struct AnalysisSameEventPairing {
12081208
if (fConfigQA) {
12091209
names = {
12101210
Form("PairsBarrelSEPM_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
1211+
Form("PairsBarrelSEPP_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
1212+
Form("PairsBarrelSEMM_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
12111213
Form("PairsBarrelSEPMCorrectAssoc_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
12121214
Form("PairsBarrelSEPMIncorrectAssoc_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
12131215
Form("PairsBarrelSEPM_ambiguousInBunch_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
@@ -1216,11 +1218,13 @@ struct AnalysisSameEventPairing {
12161218
Form("PairsBarrelSEPM_ambiguousOutOfBunch_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
12171219
Form("PairsBarrelSEPM_ambiguousOutOfBunchCorrectAssoc_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
12181220
Form("PairsBarrelSEPM_ambiguousOutOfBunchIncorrectAssoc_%s_%s", objArray->At(icut)->GetName(), sig.GetName())};
1219-
histNames += Form("%s;%s;%s;%s;%s;%s;%s;%s;%s;", names[0].Data(), names[1].Data(), names[2].Data(), names[3].Data(), names[4].Data(), names[5].Data(), names[6].Data(), names[7].Data(), names[8].Data());
1221+
histNames += Form("%s;%s;%s;%s;%s;%s;%s;%s;%s;", names[0].Data(), names[1].Data(), names[2].Data(), names[3].Data(), names[4].Data(), names[5].Data(), names[6].Data(), names[7].Data(), names[8].Data(), names[9].Data(), names[10].Data());
12201222
} else {
12211223
names = {
1222-
Form("PairsBarrelSEPM_%s_%s", objArray->At(icut)->GetName(), sig.GetName())};
1223-
histNames += Form("%s;", names[0].Data());
1224+
Form("PairsBarrelSEPM_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
1225+
Form("PairsBarrelSEPP_%s_%s", objArray->At(icut)->GetName(), sig.GetName()),
1226+
Form("PairsBarrelSEMM_%s_%s", objArray->At(icut)->GetName(), sig.GetName())};
1227+
histNames += Form("%s;%s;%s;", names[0].Data(), names[1].Data(), names[2].Data());
12241228
}
12251229
fBarrelHistNamesMCmatched.try_emplace(icut * fRecMCSignals.size() + isig, names);
12261230
} // end loop over MC signals
@@ -1582,24 +1586,24 @@ struct AnalysisSameEventPairing {
15821586
if (fConfigQA) {
15831587
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size()].Data(), VarManager::fgValues); // matched signal
15841588
if (isCorrectAssoc_leg1 && isCorrectAssoc_leg2) { // correct track-collision association
1585-
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 1].Data(), VarManager::fgValues);
1589+
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 3].Data(), VarManager::fgValues);
15861590
} else { // incorrect track-collision association
1587-
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 2].Data(), VarManager::fgValues);
1591+
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 4].Data(), VarManager::fgValues);
15881592
}
15891593
if (isAmbiInBunch) { // ambiguous in bunch
1590-
fHistMan->FillHistClass(histNames[icut][isig * fRecMCSignals.size() + 3].Data(), VarManager::fgValues);
1594+
fHistMan->FillHistClass(histNames[icut][isig * fRecMCSignals.size() + 5].Data(), VarManager::fgValues);
15911595
if (isCorrectAssoc_leg1 && isCorrectAssoc_leg2) {
1592-
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 4].Data(), VarManager::fgValues);
1596+
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 6].Data(), VarManager::fgValues);
15931597
} else {
1594-
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 5].Data(), VarManager::fgValues);
1598+
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 7].Data(), VarManager::fgValues);
15951599
}
15961600
}
15971601
if (isAmbiOutOfBunch) { // ambiguous out of bunch
1598-
fHistMan->FillHistClass(histNames[icut][isig * fRecMCSignals.size() + 6].Data(), VarManager::fgValues);
1602+
fHistMan->FillHistClass(histNames[icut][isig * fRecMCSignals.size() + 8].Data(), VarManager::fgValues);
15991603
if (isCorrectAssoc_leg1 && isCorrectAssoc_leg2) {
1600-
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 7].Data(), VarManager::fgValues);
1604+
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 9].Data(), VarManager::fgValues);
16011605
} else {
1602-
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 8].Data(), VarManager::fgValues);
1606+
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 10].Data(), VarManager::fgValues);
16031607
}
16041608
}
16051609
} else {
@@ -1610,6 +1614,11 @@ struct AnalysisSameEventPairing {
16101614
} else {
16111615
if (sign1 > 0) { // ++ pairs
16121616
fHistMan->FillHistClass(histNames[icut][1].Data(), VarManager::fgValues);
1617+
for (unsigned int isig = 0; isig < fRecMCSignals.size(); isig++) { // loop over MC signals
1618+
if (mcDecision & (uint32_t(1) << isig)) {
1619+
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 1].Data(), VarManager::fgValues);
1620+
}
1621+
}
16131622
if (fConfigQA) {
16141623
if (isAmbiInBunch) {
16151624
fHistMan->FillHistClass(histNames[icut][4].Data(), VarManager::fgValues);
@@ -1620,6 +1629,11 @@ struct AnalysisSameEventPairing {
16201629
}
16211630
} else { // -- pairs
16221631
fHistMan->FillHistClass(histNames[icut][2].Data(), VarManager::fgValues);
1632+
for (unsigned int isig = 0; isig < fRecMCSignals.size(); isig++) { // loop over MC signals
1633+
if (mcDecision & (uint32_t(1) << isig)) {
1634+
fHistMan->FillHistClass(histNamesMC[icut][isig * fRecMCSignals.size() + 2].Data(), VarManager::fgValues);
1635+
}
1636+
}
16231637
if (fConfigQA) {
16241638
if (isAmbiInBunch) {
16251639
fHistMan->FillHistClass(histNames[icut][5].Data(), VarManager::fgValues);

0 commit comments

Comments
 (0)