Skip to content

Commit 24a0312

Browse files
[PWGJE] jetSpectraCharged: jetEventWeight -> collisionWeight (#13970)
1 parent 3ffc60f commit 24a0312

File tree

1 file changed

+14
-17
lines changed

1 file changed

+14
-17
lines changed

PWGJE/Tasks/jetSpectraCharged.cxx

Lines changed: 14 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -937,8 +937,7 @@ struct JetSpectraCharged {
937937
if (!isAcceptedJet<aod::JetTracks>(jet)) {
938938
continue;
939939
}
940-
float jetweight = jet.eventWeight();
941-
fillJetAreaSubHistograms(jet, centrality, collision.rho(), jetweight);
940+
fillJetAreaSubHistograms(jet, centrality, collision.rho(), eventWeight);
942941
}
943942
}
944943
PROCESS_SWITCH(JetSpectraCharged, processSpectraAreaSubMCDWeighted, "jet spectra with rho-area subtraction for MCD", false);
@@ -1091,45 +1090,45 @@ struct JetSpectraCharged {
10911090
}
10921091
}
10931092

1093+
float eventWeight = mccollision.weight();
10941094
for (auto const& jet : jets) {
10951095
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
10961096
continue;
10971097
}
10981098
if (!isAcceptedJet<aod::JetParticles>(jet, mcLevelIsParticleLevel)) {
10991099
continue;
11001100
}
1101-
float jetweight = jet.eventWeight();
1102-
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 1.0, jetweight); // INEL
1101+
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 1.0, eventWeight); // INEL
11031102

11041103
if (!passesZvtxCut) {
11051104
continue;
11061105
}
1107-
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 2.0, jetweight); // zvtx
1106+
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 2.0, eventWeight); // zvtx
11081107

11091108
if (!hasRecoColl) {
11101109
continue;
11111110
}
1112-
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 3.0, jetweight); // noRecoColl
1111+
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 3.0, eventWeight); // noRecoColl
11131112

11141113
if (!passesSplitCollCut) {
11151114
continue;
11161115
}
1117-
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 4.0, jetweight); // splitColl
1116+
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 4.0, eventWeight); // splitColl
11181117

11191118
if (!hasSel8Coll) {
11201119
continue;
11211120
}
1122-
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 5.0, jetweight); // recoEvtSel
1121+
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 5.0, eventWeight); // recoEvtSel
11231122

11241123
if (!centralityIsGood) {
11251124
continue;
11261125
}
1127-
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 6.0, jetweight); // centralitycut
1126+
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 6.0, eventWeight); // centralitycut
11281127

11291128
if (!occupancyIsGood) {
11301129
continue;
11311130
}
1132-
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 7.0, jetweight); // occupancycut
1131+
registry.fill(HIST("h2_jet_pt_part_eventselection"), jet.pt(), 7.0, eventWeight); // occupancycut
11331132
}
11341133
}
11351134
PROCESS_SWITCH(JetSpectraCharged, processCrossSectionEfficiencyWeighted, "jet spectra QC for MC particle level with step-by-step cuts (weighted)", false);
@@ -1179,15 +1178,14 @@ struct JetSpectraCharged {
11791178
if (!isAcceptedJet<aod::JetParticles>(jet, mcLevelIsParticleLevel)) {
11801179
continue;
11811180
}
1182-
float jetweight = jet.eventWeight();
1183-
double pTHat = 10. / (std::pow(jetweight, 1.0 / pTHatExponent));
1181+
double pTHat = 10. / (std::pow(eventWeight, 1.0 / pTHatExponent));
11841182
int Nmax = 21;
11851183
for (int N = 1; N < Nmax; N++) {
11861184
if (jet.pt() < N * 0.25 * pTHat && jet.r() == round(selectedJetsRadius * 100.0f)) {
1187-
registry.fill(HIST("h2_jet_ptcut_part"), jet.pt(), N * 0.25, jetweight);
1185+
registry.fill(HIST("h2_jet_ptcut_part"), jet.pt(), N * 0.25, eventWeight);
11881186
}
11891187
}
1190-
fillMCPHistograms(jet, jetweight);
1188+
fillMCPHistograms(jet, eventWeight);
11911189
}
11921190
}
11931191
PROCESS_SWITCH(JetSpectraCharged, processSpectraMCPWeighted, "jet spectra for MC particle level weighted", false);
@@ -1214,8 +1212,7 @@ struct JetSpectraCharged {
12141212
if (!isAcceptedJet<aod::JetParticles>(jet, mcLevelIsParticleLevel)) {
12151213
continue;
12161214
}
1217-
float jetweight = jet.eventWeight();
1218-
fillMCPAreaSubHistograms(jet, mccollision.rho(), jetweight);
1215+
fillMCPAreaSubHistograms(jet, mccollision.rho(), eventWeight);
12191216
}
12201217
}
12211218
PROCESS_SWITCH(JetSpectraCharged, processSpectraAreaSubMCPWeighted, "jet spectra with area-based subtraction for MC particle level", false);
@@ -1300,7 +1297,7 @@ struct JetSpectraCharged {
13001297
if (!isAcceptedJet<aod::JetTracks>(mcdjet)) {
13011298
continue;
13021299
}
1303-
fillMatchedHistograms<ChargedMCDMatchedJetsWeighted::iterator, ChargedMCPMatchedJetsWeighted>(mcdjet, mcdjet.eventWeight());
1300+
fillMatchedHistograms<ChargedMCDMatchedJetsWeighted::iterator, ChargedMCPMatchedJetsWeighted>(mcdjet, eventWeight);
13041301
}
13051302
}
13061303
PROCESS_SWITCH(JetSpectraCharged, processJetsMatchedWeighted, "matched mcp and mcd jets with weighted events", false);

0 commit comments

Comments
 (0)