Skip to content

Commit 1d7799f

Browse files
authored
[PWGJE] Fixing missing table bug in MCD task (#15121)
1 parent 07baa23 commit 1d7799f

File tree

1 file changed

+5
-14
lines changed

1 file changed

+5
-14
lines changed

PWGJE/Tasks/jetHadronRecoil.cxx

Lines changed: 5 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -188,8 +188,6 @@ struct JetHadronRecoil {
188188
registry.add("hDeltaRpTReference", "jet p_{T} vs #DeltaR;p_{T,jet};#DeltaR", {HistType::kTH2F, {{500, -100, 400}, dRAxis}}, doSumw);
189189
registry.add("hDeltaRpTDPhiReference", "jet p_{T} vs #DeltaR vs #Delta#phi;p_{T,jet};#Delta#phi;#DeltaR", {HistType::kTH3F, {{500, -100, 400}, {100, 0, o2::constants::math::TwoPI}, dRAxis}}, doSumw);
190190
registry.add("hDeltaRpTDPhiReferenceShifts", "testing shifts;p_{T,jet};#Delta#phi;#DeltaR;shifts", {HistType::kTHnSparseD, {{500, -100, 400}, {100, 0, o2::constants::math::TwoPI}, dRAxis, {20, 0.0, 2.0}}}, doSumw);
191-
registry.add("hPtTrackMatched", "Track p_{T};p_{T};entries", {HistType::kTH1F, {{200, 0, 200}}}, doSumw);
192-
registry.add("hPtTrackMatchedToCollisions", "Track p_{T};p_{T};entries", {HistType::kTH1F, {{200, 0, 200}}}, doSumw);
193191
}
194192

195193
if (doprocessMCP || doprocessMCPWeighted || doprocessMCPOnTheFly || doprocessMCPOnTheFlyWeighted) {
@@ -383,7 +381,7 @@ struct JetHadronRecoil {
383381
}
384382

385383
template <typename T, typename U>
386-
void fillHistogramsMCD(T const& jets, U const& tracks, float weight = 1.0, float rho = 0.0, float pTHat = 999.0, auto collisionID = 0)
384+
void fillHistogramsMCD(T const& jets, U const& tracks, float weight = 1.0, float rho = 0.0, float pTHat = 999.0)
387385
{
388386
bool isSigCol;
389387
std::vector<double> phiTTAr;
@@ -433,13 +431,6 @@ struct JetHadronRecoil {
433431
registry.fill(HIST("hPhiTrack"), track.phi(), weight);
434432
registry.fill(HIST("hTrack3D"), track.pt(), track.eta(), track.phi(), weight);
435433
registry.fill(HIST("hPtTrackPtHard"), track.pt() / pTHat, track.pt(), weight);
436-
if (track.has_mcParticle()) {
437-
registry.fill(HIST("hPtTrackMatched"), track.pt(), weight);
438-
auto mcParticle = track.mcParticle();
439-
if (mcParticle.mcCollisionId() == collisionID) {
440-
registry.fill(HIST("hPtTrackMatchedToCollisions"), track.pt(), weight);
441-
}
442-
}
443434
}
444435
if (nTT > 0) {
445436
int trigNumber = rand->Integer(nTT);
@@ -891,7 +882,7 @@ struct JetHadronRecoil {
891882
return;
892883
}
893884
registry.fill(HIST("hZvtxSelected"), collision.posZ());
894-
fillHistogramsMCD(jets, tracks, 1.0, 0.0, collision.mcCollision().ptHard(), collision.mcCollisionId());
885+
fillHistogramsMCD(jets, tracks, 1.0, 0.0, collision.mcCollision().ptHard());
895886
}
896887
PROCESS_SWITCH(JetHadronRecoil, processMCD, "process MC detector level", false);
897888

@@ -913,7 +904,7 @@ struct JetHadronRecoil {
913904
return;
914905
}
915906
registry.fill(HIST("hZvtxSelected"), collision.posZ());
916-
fillHistogramsMCD(jets, tracks, 1.0, collision.rho(), collision.mcCollision().ptHard(), collision.mcCollisionId());
907+
fillHistogramsMCD(jets, tracks, 1.0, collision.rho(), collision.mcCollision().ptHard());
917908
}
918909
PROCESS_SWITCH(JetHadronRecoil, processMCDWithRhoSubtraction, "process MC detector level with rho subtraction", false);
919910

@@ -935,7 +926,7 @@ struct JetHadronRecoil {
935926
return;
936927
}
937928
registry.fill(HIST("hZvtxSelected"), collision.posZ(), collision.mcCollision().weight());
938-
fillHistogramsMCD(jets, tracks, collision.mcCollision().weight(), 0.0, collision.mcCollision().ptHard(), collision.mcCollisionId());
929+
fillHistogramsMCD(jets, tracks, collision.mcCollision().weight(), 0.0, collision.mcCollision().ptHard());
939930
}
940931
PROCESS_SWITCH(JetHadronRecoil, processMCDWeighted, "process MC detector level with event weights", false);
941932

@@ -957,7 +948,7 @@ struct JetHadronRecoil {
957948
return;
958949
}
959950
registry.fill(HIST("hZvtxSelected"), collision.posZ(), collision.mcCollision().weight());
960-
fillHistogramsMCD(jets, tracks, collision.mcCollision().weight(), collision.rho(), collision.mcCollision().ptHard(), collision.mcCollisionId());
951+
fillHistogramsMCD(jets, tracks, collision.mcCollision().weight(), collision.rho(), collision.mcCollision().ptHard());
961952
}
962953
PROCESS_SWITCH(JetHadronRecoil, processMCDWeightedWithRhoSubtraction, "process MC detector level with event weights and rho subtraction", false);
963954

0 commit comments

Comments
 (0)