Skip to content

Commit 0d2acd4

Browse files
Copilotfeisenhu
andcommitted
Fix cfgRequireTrueAssociation missing in TTCA weight calculation for DileptonMC.h, SingleTrackQCMC.h, and checkMCPairTemplate.cxx
Agent-Logs-Url: https://github.com/feisenhu/O2Physics/sessions/1d6d35c8-a846-4be0-b49e-97d5dba066f8 Co-authored-by: feisenhu <53603353+feisenhu@users.noreply.github.com>
1 parent 6533bee commit 0d2acd4

File tree

3 files changed

+24
-0
lines changed

3 files changed

+24
-0
lines changed

PWGEM/Dilepton/Core/DileptonMC.h

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -2222,6 +2222,9 @@ struct DileptonMC {
22222222
if (cfgEventGeneratorType >= 0 && mccollision_from_neg.getSubGeneratorId() != cfgEventGeneratorType) {
22232223
continue;
22242224
}
2225+
if (cfgRequireTrueAssociation && (mcpos.emmceventId() != collision.emmceventId() || mcneg.emmceventId() != collision.emmceventId())) {
2226+
continue;
2227+
}
22252228

22262229
if (isPairOK(pos, neg, cut, tracks)) {
22272230
passed_pairIds.emplace_back(std::make_pair(pos.globalIndex(), neg.globalIndex()));
@@ -2238,6 +2241,9 @@ struct DileptonMC {
22382241
if (cfgEventGeneratorType >= 0 && mccollision_from_pos2.getSubGeneratorId() != cfgEventGeneratorType) {
22392242
continue;
22402243
}
2244+
if (cfgRequireTrueAssociation && (mcpos1.emmceventId() != collision.emmceventId() || mcpos2.emmceventId() != collision.emmceventId())) {
2245+
continue;
2246+
}
22412247

22422248
if (isPairOK(pos1, pos2, cut, tracks)) {
22432249
passed_pairIds.emplace_back(std::make_pair(pos1.globalIndex(), pos2.globalIndex()));
@@ -2254,6 +2260,9 @@ struct DileptonMC {
22542260
if (cfgEventGeneratorType >= 0 && mccollision_from_neg2.getSubGeneratorId() != cfgEventGeneratorType) {
22552261
continue;
22562262
}
2263+
if (cfgRequireTrueAssociation && (mcneg1.emmceventId() != collision.emmceventId() || mcneg2.emmceventId() != collision.emmceventId())) {
2264+
continue;
2265+
}
22572266
if (isPairOK(neg1, neg2, cut, tracks)) {
22582267
passed_pairIds.emplace_back(std::make_pair(neg1.globalIndex(), neg2.globalIndex()));
22592268
}

PWGEM/Dilepton/Core/SingleTrackQCMC.h

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1055,6 +1055,9 @@ struct SingleTrackQCMC {
10551055
if (cfgEventGeneratorType >= 0 && mccollision_from_track.getSubGeneratorId() != cfgEventGeneratorType) {
10561056
continue;
10571057
}
1058+
if (cfgRequireTrueAssociation && (mctrack.emmceventId() != collision.emmceventId())) {
1059+
continue;
1060+
}
10581061

10591062
if (dielectroncuts.cfg_pid_scheme == static_cast<int>(DielectronCut::PIDSchemes::kPIDML)) {
10601063
if (!cut.template IsSelectedTrack<false>(track)) {
@@ -1074,6 +1077,9 @@ struct SingleTrackQCMC {
10741077
if (cfgEventGeneratorType >= 0 && mccollision_from_track.getSubGeneratorId() != cfgEventGeneratorType) {
10751078
continue;
10761079
}
1080+
if (cfgRequireTrueAssociation && (mctrack.emmceventId() != collision.emmceventId())) {
1081+
continue;
1082+
}
10771083
if (!cut.template IsSelectedTrack<false>(track)) {
10781084
continue;
10791085
}

PWGEM/Dilepton/Tasks/checkMCPairTemplate.cxx

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -2804,6 +2804,9 @@ struct checkMCPairTemplate {
28042804
if (cfgEventGeneratorType >= 0 && mccollision_from_neg.getSubGeneratorId() != cfgEventGeneratorType) {
28052805
continue;
28062806
}
2807+
if (cfgRequireTrueAssociation && (mcpos.emmceventId() != collision.emmceventId() || mcneg.emmceventId() != collision.emmceventId())) {
2808+
continue;
2809+
}
28072810

28082811
if (isPairOK(pos, neg, cut, tracks)) {
28092812
passed_pairIds.emplace_back(std::make_pair(pos.globalIndex(), neg.globalIndex()));
@@ -2820,6 +2823,9 @@ struct checkMCPairTemplate {
28202823
if (cfgEventGeneratorType >= 0 && mccollision_from_pos2.getSubGeneratorId() != cfgEventGeneratorType) {
28212824
continue;
28222825
}
2826+
if (cfgRequireTrueAssociation && (mcpos1.emmceventId() != collision.emmceventId() || mcpos2.emmceventId() != collision.emmceventId())) {
2827+
continue;
2828+
}
28232829

28242830
if (isPairOK(pos1, pos2, cut, tracks)) {
28252831
passed_pairIds.emplace_back(std::make_pair(pos1.globalIndex(), pos2.globalIndex()));
@@ -2836,6 +2842,9 @@ struct checkMCPairTemplate {
28362842
if (cfgEventGeneratorType >= 0 && mccollision_from_neg2.getSubGeneratorId() != cfgEventGeneratorType) {
28372843
continue;
28382844
}
2845+
if (cfgRequireTrueAssociation && (mcneg1.emmceventId() != collision.emmceventId() || mcneg2.emmceventId() != collision.emmceventId())) {
2846+
continue;
2847+
}
28392848
if (isPairOK(neg1, neg2, cut, tracks)) {
28402849
passed_pairIds.emplace_back(std::make_pair(neg1.globalIndex(), neg2.globalIndex()));
28412850
}

0 commit comments

Comments
 (0)