Skip to content

Commit 9144b82

Browse files
yhou-gityhou-git
committed
merged serval processes
1 parent 3591e45 commit 9144b82

File tree

1 file changed

+77
-77
lines changed

1 file changed

+77
-77
lines changed

‎PWGJE/Tasks/chargedJetHadron.cxx‎

Lines changed: 77 additions & 77 deletions
Original file line numberDiff line numberDiff line change
@@ -162,39 +162,37 @@ struct ChargedJetHadron {
162162
AxisSpec dphijetAxis = {160, -1.7, 6.3, "#Delta#varphi"};
163163
AxisSpec detaAxis = {160, -1.6, 1.6, "#Delta#eta"};
164164
AxisSpec drAxis = {200, 0.0, 5.0, "#Delta#it{R}"};
165+
AxisSpec axisBdtScore = {100, 0., 1., "Bdt score"};
165166

166-
if (doprocessCollisions || doprocessCollisionsWeighted) {
167+
if (doprocessCollisionsQA || doprocessCollisionsQAWeighted) {
167168
registry.add("h_collisions", "event status;event status; entries", {HistType::kTH1F, {{4, 0.0, 4.0}}});
168169
registry.add("h_fakecollisions", "event status;event status; entries", {HistType::kTH1F, {{4, 0.0, 4.0}}});
169170
registry.add("h2_centrality_occupancy", "centrality vs occupancy; centrality; occupancy", {HistType::kTH2F, {centralityAxis, {60, 0, 30000}}});
170171
registry.add("h_collisions_Zvertex", "position of collision; #it{Z} (cm)", {HistType::kTH1F, {{300, -15.0, 15.0}}});
171172
registry.add("h_collisions_multFT0", " multiplicity using multFT0; entries", {HistType::kTH1F, {{3000, 0, 10000}}});
172-
if (doprocessCollisionsWeighted) {
173+
if (doprocessCollisionsQAWeighted) {
173174
registry.add("h_collisions_weighted", "event status;event status;entries", {HistType::kTH1F, {{4, 0.0, 4.0}}});
174175
}
175-
}
176-
177-
if (doprocessTracksQC || doprocessTracksQCWeighted) {
176+
//if (doprocessTracksQC || doprocessTracksQCWeighted) {
178177
registry.add("h_track_pt", "track #it{p}_{T}; #it{p}_{T,track} (GeV/#it{c})", {HistType::kTH1F, {trackPtAxis}});
179178
registry.add("h2_track_eta_track_phi", "track #eta vs. track #phi; #eta; #phi; counts", {HistType::kTH2F, {etaAxis, phiAxis}});
180179
registry.add("h2_track_eta_pt", "track #eta vs. track #it{p}_{T}; #eta; #it{p}_{T,track} (GeV/#it{c}; counts", {HistType::kTH2F, {etaAxis, trackPtAxis}});
181180
registry.add("h2_track_phi_pt", "track #phi vs. track #it{p}_{T}; #phi; #it{p}_{T,track} (GeV/#it{c}; counts", {HistType::kTH2F, {phiAxis, trackPtAxis}});
182181
}
183182

184-
if (doprocessSpectraData || doprocessSpectraMCD || doprocessSpectraMCDWeighted) {
183+
184+
if (doprocessSpectraAreaSubData || doprocessSpectraAreaSubMCD || doprocessSpectraMCDWeighted) {
185+
if (doprocessSpectraMCDWeighted) {
186+
registry.add("h_jet_phat", "jet #hat{p};#hat{p} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0, 1000}}});
187+
registry.add("h_jet_phat_weighted", "jet #hat{p};#hat{p} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0, 1000}}});
188+
}
185189
registry.add("h_jet_pt", "jet pT; #it{p}_{T,jet} (GeV/#it{c}); counts", {HistType::kTH1F, {jetPtAxis}});
186190
registry.add("h_jet_eta", "jet eta; #eta_{jet}; counts", {HistType::kTH1F, {etaAxis}});
187191
registry.add("h_jet_phi", "jet phi; #phi_{jet}; counts", {HistType::kTH1F, {phiAxis}});
188192
registry.add("h_jet_area", "jet Area_{jet}; Area_{jet}; counts", {HistType::kTH1F, {{150, 0., 1.5}}});
189193
registry.add("h_jet_ntracks", "jet N_{jet tracks}; N_{jet, tracks; counts}", {HistType::kTH1F, {{200, -0.5, 199.5}}});
190194
registry.add("h2_jet_pt_track_pt", "jet #it{p}_{T,jet} vs. #it{p}_{T,track}; #it{p}_{T,jet} (GeV/#it{c}); #it{p}_{T,track} (GeV/#it{c})", {HistType::kTH2F, {jetPtAxis, trackPtAxis}});
191-
if (doprocessSpectraMCDWeighted) {
192-
registry.add("h_jet_phat", "jet #hat{p};#hat{p} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0, 1000}}});
193-
registry.add("h_jet_phat_weighted", "jet #hat{p};#hat{p} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0, 1000}}});
194-
}
195-
}
196195

197-
if (doprocessSpectraAreaSubData || doprocessSpectraAreaSubMCD) {
198196
registry.add("h_jet_pt_rhoareasubtracted", "jet pt; #it{p}_{T,jet} (GeV/#it{c}); counts", {HistType::kTH1F, {jetPtAxisRhoAreaSub}});
199197
registry.add("h_jet_eta_rhoareasubtracted", "jet eta; #eta_{jet}; counts", {HistType::kTH1F, {etaAxis}});
200198
registry.add("h_jet_phi_rhoareasubtracted", "jet phi; #phi_{jet}; counts", {HistType::kTH1F, {phiAxis}});
@@ -231,6 +229,10 @@ struct ChargedJetHadron {
231229
registry.add("h_d0_mass", ";m_{D^{0}} (GeV/c^{2});dN/dm_{D^{0}}", {HistType::kTH1F, {{1000, 0., 10.}}});
232230
registry.add("h_d0_eta", ";#eta_{D^{0}} (GeV/c^{2});dN/d#eta_{D^{0}}", {HistType::kTH1F, {{200, -5., 5.}}});
233231
registry.add("h_d0_phi", ";#varphi_{D^{0}} (GeV/c^{2});dN/d#varphi_{D^{0}}", {HistType::kTH1F, {{200, -10., 10.}}});
232+
registry.add("h_d0_bdtScorePrompt", "D0 BDT prompt score", {HistType::kTH1F, {axisBdtScore}});
233+
registry.add("h_d0_bdtScoreBkg", "D0 BDT bkg score", {HistType::kTH1F, {axisBdtScore}});
234+
registry.add("h_d0bar_bdtScorePrompt", "D0bar BDT prompt score", {HistType::kTH1F, {axisBdtScore}});
235+
registry.add("h_d0bar_bdtScoreBkg", "D0bar BDT bkg score", {HistType::kTH1F, {axisBdtScore}});
234236
registry.add("h2_d0jet_detadphi", "D^{0}-jets deta vs dphi; #Delta#eta; #Delta#phi", {HistType::kTH2F, {detaAxis, dphiAxis}});
235237
}
236238

@@ -305,8 +307,6 @@ struct ChargedJetHadron {
305307
registry.get<TH1>(HIST("h_mcColl_counts"))->GetXaxis()->SetBinLabel(4, "recoEvtSel");
306308
registry.get<TH1>(HIST("h_mcColl_counts"))->GetXaxis()->SetBinLabel(5, "centralitycut");
307309
registry.get<TH1>(HIST("h_mcColl_counts"))->GetXaxis()->SetBinLabel(6, "occupancycut");
308-
309-
registry.add("h_mc_zvertex", "position of collision ;#it{Z} (cm)", {HistType::kTH1F, {{300, -15.0, 15.0}}});
310310
registry.add("h_mc_mult", " multiplicity global tracks; entries", {HistType::kTH1F, {{3000, 0, 60000}}});
311311

312312
registry.add("h_jet_pt_part", "partvjet pT;#it{p}_{T,jet}^{part} (GeV/#it{c}); counts", {HistType::kTH1F, {jetPtAxis}});
@@ -349,7 +349,7 @@ struct ChargedJetHadron {
349349
registry.get<TH1>(HIST("h_mixjeth_event_stats_part"))->GetXaxis()->SetBinLabel(5, "Total j-h pairs with accepted");
350350
}
351351

352-
if (doprocessSpectraAreaSubMCP || doprocessMixLeadingJetHadronMCP) {
352+
if (doprocessLeadingJetHadronMCP || doprocessMixLeadingJetHadronMCP) {
353353
registry.add("h_mcColl_counts_areasub", " number of mc events; event status; entries", {HistType::kTH1F, {{10, 0, 10}}});
354354
registry.get<TH1>(HIST("h_mcColl_counts_areasub"))->GetXaxis()->SetBinLabel(1, "allMcColl");
355355
registry.get<TH1>(HIST("h_mcColl_counts_areasub"))->GetXaxis()->SetBinLabel(2, "vertexZ");
@@ -359,8 +359,10 @@ struct ChargedJetHadron {
359359
registry.get<TH1>(HIST("h_mcColl_counts_areasub"))->GetXaxis()->SetBinLabel(6, "centralitycut");
360360
registry.get<TH1>(HIST("h_mcColl_counts_areasub"))->GetXaxis()->SetBinLabel(7, "occupancycut");
361361

362+
registry.add("h_mcColl_zvertex", "position of collision ;#it{Z} (cm)", {HistType::kTH1F, {{300, -15.0, 15.0}}});
362363
registry.add("h_mcColl_rho", "mc collision rho;#rho (GeV/#it{c}); counts", {HistType::kTH1F, {{500, 0.0, 500.0}}});
363364
registry.add("h_mcColl_centrality", "mc collision centrality; centrality; counts", {HistType::kTH1F, {centralityAxis}});
365+
registry.add("h_mcColl_mult", " multiplicity global tracks; entries", {HistType::kTH1F, {{3000, 0, 60000}}});
364366

365367
registry.add("h_particle_pt", "particle #it{p}_{T}; #it{p}_{T,particle} (GeV/#it{c})", {HistType::kTH1F, {trackPtAxis}});
366368
registry.add("h2_particle_eta_phi", "particle #eta vs. particle #phi; #eta; #phi; counts", {HistType::kTH2F, {etaAxis, phiAxis}});
@@ -1227,7 +1229,8 @@ struct ChargedJetHadron {
12271229
// ==========================================================
12281230
//.............process staring...............................
12291231
// ==========================================================
1230-
void processCollisions(soa::Filtered<aod::JetCollisions>::iterator const& collision)
1232+
void processCollisionsQA(soa::Filtered<aod::JetCollisions>::iterator const& collision,
1233+
FilterJetTracks const& tracks)
12311234
{
12321235
registry.fill(HIST("h_collisions"), 0.5);
12331236
if (!jetderiveddatautilities::selectCollision(collision, eventSelectionBits, skipMBGapEvents)) {
@@ -1241,11 +1244,19 @@ struct ChargedJetHadron {
12411244
registry.fill(HIST("h2_centrality_occupancy"), getCentrality(collision), collision.trackOccupancyInTimeRange());
12421245
registry.fill(HIST("h_collisions_Zvertex"), collision.posZ());
12431246
registry.fill(HIST("h_collisions_multFT0"), getMultiplicity(collision)); // collision.MultFT0M()
1247+
1248+
for (auto const& track : tracks) {
1249+
if (!jetderiveddatautilities::selectTrack(track, trackSelection)) {
1250+
continue;
1251+
}
1252+
fillTrackHistograms(track);
1253+
}
12441254
}
1245-
PROCESS_SWITCH(ChargedJetHadron, processCollisions, "collisions Data and MCD", true);
1255+
PROCESS_SWITCH(ChargedJetHadron, processCollisionsQA, "collisions and tracks QC for Data and MCD", true);
12461256

12471257
void processCollisionsWeighted(soa::Join<aod::JetCollisions, aod::JMcCollisionLbs>::iterator const& collision,
1248-
aod::JetMcCollisions const&)
1258+
aod::JetMcCollisions const&,
1259+
FilterJetTracks const& tracks)
12491260
{
12501261
if (!collision.has_mcCollision()) {
12511262
registry.fill(HIST("h_fakecollisions"), 0.5);
@@ -1258,51 +1269,24 @@ struct ChargedJetHadron {
12581269
}
12591270
registry.fill(HIST("h_collisions"), 1.5);
12601271
registry.fill(HIST("h_collisions_weighted"), 1.5, eventWeight);
1261-
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
1272+
1273+
if (!isGoodCollision(collision, eventSelectionBits, skipMBGapEvents, trackOccupancyInTimeRangeMin, trackOccupancyInTimeRangeMax, centralityMin, centralityMax, vertexZCut)) {
12621274
return;
12631275
}
12641276
registry.fill(HIST("h_collisions"), 2.5);
12651277
registry.fill(HIST("h_collisions_weighted"), 2.5, eventWeight);
12661278
registry.fill(HIST("h2_centrality_occupancy"), getCentrality(collision), collision.trackOccupancyInTimeRange());
12671279
registry.fill(HIST("h_collisions_Zvertex"), collision.posZ(), eventWeight);
1268-
}
1269-
PROCESS_SWITCH(ChargedJetHadron, processCollisionsWeighted, "weighted collisions for MCD", false);
1270-
1271-
void processTracksQC(soa::Filtered<aod::JetCollisions>::iterator const& collision,
1272-
FilterJetTracks const& tracks)
1273-
{
1274-
if (!isGoodCollision(collision, eventSelectionBits, skipMBGapEvents, trackOccupancyInTimeRangeMin, trackOccupancyInTimeRangeMax, centralityMin, centralityMax, vertexZCut)) {
1275-
return;
1276-
}
1277-
for (auto const& track : tracks) {
1278-
if (!jetderiveddatautilities::selectTrack(track, trackSelection)) {
1279-
continue;
1280-
}
1281-
fillTrackHistograms(track);
1282-
}
1283-
}
1284-
PROCESS_SWITCH(ChargedJetHadron, processTracksQC, "collisions and track QC for Data and MCD", false);
12851280

1286-
void processTracksQCWeighted(soa::Join<aod::JetCollisions, aod::JMcCollisionLbs>::iterator const& collision,
1287-
aod::JetMcCollisions const&,
1288-
FilterJetTracks const& tracks)
1289-
{
1290-
float eventWeight = collision.weight();
1291-
if (!jetderiveddatautilities::selectCollision(collision, eventSelectionBits, skipMBGapEvents)) {
1292-
return;
1293-
}
1294-
if (std::abs(collision.posZ()) > vertexZCut) {
1295-
return;
1296-
}
12971281
for (auto const& track : tracks) {
12981282
if (!jetderiveddatautilities::selectTrack(track, trackSelection)) {
12991283
continue;
13001284
}
13011285
fillTrackHistograms(track, eventWeight);
13021286
}
13031287
}
1304-
PROCESS_SWITCH(ChargedJetHadron, processTracksQCWeighted, "weighted collisions and tracks QC for MCD", false);
1305-
1288+
PROCESS_SWITCH(ChargedJetHadron, processCollisionsQCWeighted, "weighted collisions and tracks for MCD", false);
1289+
/*
13061290
void processSpectraData(soa::Filtered<aod::JetCollisions>::iterator const& collision,
13071291
CorrChargedJets const& jets,
13081292
aod::JetTracks const&)
@@ -1321,26 +1305,7 @@ struct ChargedJetHadron {
13211305
}
13221306
}
13231307
PROCESS_SWITCH(ChargedJetHadron, processSpectraData, "jet spectra for Data", false);
1324-
1325-
void processSpectraMCD(soa::Filtered<aod::JetCollisions>::iterator const& collision,
1326-
CorrChargedMCDJets const& jets,
1327-
aod::JetTracks const&)
1328-
{
1329-
if (!isGoodCollision(collision, eventSelectionBits, skipMBGapEvents, trackOccupancyInTimeRangeMin, trackOccupancyInTimeRangeMax, centralityMin, centralityMax, vertexZCut)) {
1330-
return;
1331-
}
1332-
for (auto const& jet : jets) {
1333-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
1334-
continue;
1335-
}
1336-
if (!isAcceptedJet<aod::JetTracks>(jet)) {
1337-
continue;
1338-
}
1339-
fillJetHistograms(jet);
1340-
}
1341-
}
1342-
PROCESS_SWITCH(ChargedJetHadron, processSpectraMCD, "jet spectra for MCD", false);
1343-
1308+
*/
13441309
void processSpectraAreaSubData(FilterCollision const& collision,
13451310
CorrChargedJets const& jets,
13461311
aod::JetTracks const&)
@@ -1355,6 +1320,7 @@ struct ChargedJetHadron {
13551320
if (!isAcceptedJet<aod::JetTracks>(jet)) {
13561321
continue;
13571322
}
1323+
fillJetHistograms(jet);
13581324
fillJetAreaSubHistograms(jet, collision.rho());
13591325
}
13601326
}
@@ -1449,6 +1415,13 @@ struct ChargedJetHadron {
14491415
registry.fill(HIST("h_d0_pt"), candidate.pt());
14501416
registry.fill(HIST("h_d0_eta"), candidate.eta());
14511417
registry.fill(HIST("h_d0_phi"), candidate.phi());
1418+
/*
1419+
registry.fill(HIST("h_d0_bdtScorePrompt"), candidate.mlScorePromptD0());
1420+
registry.fill(HIST("h_d0_bdtScoreBkg"), candidate.mlScoreBkgD0());
1421+
registry.fill(HIST("h_d0bar_bdtScorePrompt"), candidate.mlScorePromptD0bar());
1422+
registry.fill(HIST("h_d0bar_bdtScoreBkg"), candidate.mlScoreBkgD0bar());
1423+
//registry.fill(HIST("h_d0_"), candidate.());
1424+
*/
14521425
for (const auto& jet : jets) {
14531426
double deltaeta = candidate.eta() - jet.eta();
14541427
double deltaphi = candidate.phi() - jet.phi();
@@ -1458,8 +1431,28 @@ struct ChargedJetHadron {
14581431
}
14591432
}
14601433
PROCESS_SWITCH(ChargedJetHadron, processHFJetCorrelation, "D0-jet for Data", false);
1461-
//.................................................................................
14621434

1435+
//........MCD.........................................................................
1436+
/*
1437+
void processSpectraMCD(soa::Filtered<aod::JetCollisions>::iterator const& collision,
1438+
CorrChargedMCDJets const& jets,
1439+
aod::JetTracks const&)
1440+
{
1441+
if (!isGoodCollision(collision, eventSelectionBits, skipMBGapEvents, trackOccupancyInTimeRangeMin, trackOccupancyInTimeRangeMax, centralityMin, centralityMax, vertexZCut)) {
1442+
return;
1443+
}
1444+
for (auto const& jet : jets) {
1445+
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
1446+
continue;
1447+
}
1448+
if (!isAcceptedJet<aod::JetTracks>(jet)) {
1449+
continue;
1450+
}
1451+
fillJetHistograms(jet);
1452+
}
1453+
}
1454+
PROCESS_SWITCH(ChargedJetHadron, processSpectraMCD, "jet spectra for MCD", false);
1455+
*/
14631456
void processSpectraAreaSubMCD(FilterCollision const& collision,
14641457
CorrChargedMCDJets const& jets,
14651458
aod::JetTracks const&)
@@ -1474,6 +1467,7 @@ struct ChargedJetHadron {
14741467
if (!isAcceptedJet<aod::JetTracks>(jet)) {
14751468
continue;
14761469
}
1470+
fillJetHistograms(jet);
14771471
fillJetAreaSubHistograms(jet, collision.rho());
14781472
}
14791473
}
@@ -1614,7 +1608,6 @@ struct ChargedJetHadron {
16141608
return;
16151609
}
16161610
registry.fill(HIST("h_mcColl_counts"), 5.5);
1617-
registry.fill(HIST("h_mc_zvertex"), mccollision.posZ());
16181611
// registry.fill(HIST("h_mc_mult"), getMultiplicity(mccollision));//problem for mccollision' has no member named 'multFT0M'; did you mean 'multFT0A'? // return coll.multFT0M();
16191612
for (auto const& collision : collisions) {
16201613
registry.fill(HIST("h_mc_mult"), getMultiplicity(collision));
@@ -1632,16 +1625,17 @@ struct ChargedJetHadron {
16321625
}
16331626
PROCESS_SWITCH(ChargedJetHadron, processSpectraMCP, "jet spectra for MCP", false);
16341627

1635-
void processSpectraAreaSubMCP(McParticleCollision const& mccollision,
1636-
soa::SmallGroups<aod::JetCollisionsMCD> const& collisions,
1637-
CorrChargedMCPJets const& jets,
1638-
soa::Filtered<aod::JetParticles> const& particles)
1628+
void processLeadingJetHadronMCP(McParticleCollision const& mccollision,
1629+
soa::SmallGroups<aod::JetCollisionsMCD> const& collisions,
1630+
CorrChargedMCPJets const& jets,
1631+
soa::Filtered<aod::JetParticles> const& particles)
16391632
{
16401633
bool mcLevelIsParticleLevel = true;
16411634
bool hasSel8Coll = false;
16421635
bool centralityIsGood = false;
16431636
bool occupancyIsGood = false;
16441637
float centrality = -999.0;
1638+
float multiplicity = -999.0;
16451639

16461640
registry.fill(HIST("h_mcColl_counts_areasub"), 0.5);
16471641
if (std::abs(mccollision.posZ()) > vertexZCut) {
@@ -1659,6 +1653,7 @@ struct ChargedJetHadron {
16591653

16601654
if (acceptSplitCollisions == SplitOkCheckFirstAssocCollOnly || acceptSplitCollisions == NonSplitOnly) {
16611655
centrality = getCentrality(collisions.begin());
1656+
multiplicity = getMultiplicity(collisions.begin());
16621657
if (jetderiveddatautilities::selectCollision(collisions.begin(), eventSelectionBits, skipMBGapEvents)) {
16631658
hasSel8Coll = true;
16641659
}
@@ -1671,6 +1666,7 @@ struct ChargedJetHadron {
16711666
} else {
16721667
for (auto const& collision : collisions) {
16731668
centrality = getCentrality(collision);
1669+
multiplicity = getMultiplicity(collision);
16741670
if (jetderiveddatautilities::selectCollision(collision, eventSelectionBits, skipMBGapEvents)) {
16751671
hasSel8Coll = true;
16761672
}
@@ -1686,8 +1682,11 @@ struct ChargedJetHadron {
16861682
return;
16871683
}
16881684
registry.fill(HIST("h_mcColl_counts_areasub"), 6.5);
1685+
registry.fill(HIST("h_mcColl_zvertex"), mccollision.posZ());
16891686
registry.fill(HIST("h_mcColl_rho"), mccollision.rho());
16901687
registry.fill(HIST("h_mcColl_centrality"), centrality);
1688+
registry.fill(HIST("h_mcColl_mult"), multiplicity);
1689+
//problem for mccollision' has no member named 'multFT0M'; did you mean 'multFT0A'? // return coll.multFT0M();
16911690

16921691
// particle QA...........
16931692
for (auto const& particle : particles) {
@@ -1704,7 +1703,7 @@ struct ChargedJetHadron {
17041703
}
17051704
fillMCPLeadingJetHadronHistograms(mccollision, jets, particles);
17061705
}
1707-
PROCESS_SWITCH(ChargedJetHadron, processSpectraAreaSubMCP, "jet spectra with area-based & SM leading jet-hadron for MCP", false);
1706+
PROCESS_SWITCH(ChargedJetHadron, processLeadingJetHadronMCP, "jet spectra with area-based & SM leading jet-hadron for MCP", false);
17081707

17091708
void processMixLeadingJetHadronMCP(McParticleCollisions const& mccollisions,
17101709
soa::SmallGroups<aod::JetCollisionsMCD> const& collisions,
@@ -1850,10 +1849,11 @@ struct ChargedJetHadron {
18501849
if (!hasSel8Coll || !centralityIsGood || !occupancyIsGood) {
18511850
return;
18521851
}
1853-
fillMCPMixJetHadronHistograms(mccollisions, jets, particles); // MCP mixed problem: error: no matching function for call to 'o2::aod::mult::MultNTracksGlobal::MultNTracksGlobal()
1852+
fillMCPMixJetHadronHistograms(mccollisions, jets, particles);
18541853
}
18551854
PROCESS_SWITCH(ChargedJetHadron, processMixJetHadronMCP, "mixed event jet-hadron for MCP", false);
18561855

1856+
//............weighted................
18571857
void processSpectraMCPWeighted(aod::JetMcCollision const& mccollision,
18581858
soa::SmallGroups<aod::JetCollisionsMCD> const& collisions,
18591859
soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents, aod::ChargedMCParticleLevelJetEventWeights> const& jets,

0 commit comments

Comments
 (0)