Skip to content

Commit b436134

Browse files
spucillospucillo
authored
[PWGLF] Update MC centrality (#18233)
1 parent d56fd1a commit b436134

1 file changed

Lines changed: 150 additions & 43 deletions

File tree

‎PWGLF/Tasks/Strangeness/cascadeAnalysisLightIonsDerivedData.cxx‎

Lines changed: 150 additions & 43 deletions
Original file line numberDiff line numberDiff line change
@@ -17,6 +17,7 @@
1717

1818
#include "PWGLF/DataModel/LFStrangenessPIDTables.h"
1919
#include "PWGLF/DataModel/LFStrangenessTables.h"
20+
#include "PWGLF/DataModel/mcCentrality.h"
2021

2122
#include "Common/CCDB/EventSelectionParams.h"
2223
#include "Common/CCDB/RCTSelectionFlags.h"
@@ -56,7 +57,7 @@ using SimCollisions = soa::Join<aod::StraCollisions, aod::StraEvSels, aod::StraC
5657
using CascadeCandidates = soa::Join<aod::CascCollRefs, aod::CascCores, aod::CascExtras, aod::CascTOFPIDs, aod::CascTOFNSigmas>;
5758
using CascadeMCCandidates = soa::Join<aod::CascCollRefs, aod::CascCores, aod::CascExtras, aod::CascTOFPIDs, aod::CascTOFNSigmas, aod::CascCoreMCLabels>;
5859
using DaughterTracks = soa::Join<aod::DauTrackExtras, aod::DauTrackTPCPIDs>;
59-
using CollisionMCTrueTable = soa::Join<aod::StraMCCollisions, aod::StraMCCollMults>;
60+
using CollisionMCTrueTable = soa::Join<aod::StraMCCollisions, aod::StraMCCollMults, aod::McCentFV0As, aod::McCentFT0Ms, aod::McCentFT0Cs, aod::McCentNGlobals>;
6061
using CascadeMCCores = soa::Join<aod::CascMCCores, aod::CascMCCollRefs>;
6162

6263
struct CascadeAnalysisLightIonsDerivedData {
@@ -151,6 +152,7 @@ struct CascadeAnalysisLightIonsDerivedData {
151152

152153
// Centrality estimator
153154
Configurable<int> centralityEstimator{"centralityEstimator", 0, "0 = FT0C, 1 = FTOM, 2 = FV0A, 3 = NGlobal"};
155+
Configurable<bool> useMcCentrality{"useMcCentrality", false, "Use MC centrality?"};
154156

155157
// List of estimators
156158
enum Option { kFT0C,
@@ -299,6 +301,8 @@ struct CascadeAnalysisLightIonsDerivedData {
299301
// Multiplicity Histograms
300302
registryMC.add("hCentEstimator_truerec", "hCentEstimator_truerec", HistType::kTH1D, {{101, 0.0f, 101.0f}});
301303
registryMC.add("hCentralityVsNch_truerec", "hCentralityVsNch_truerec", HistType::kTH2D, {{101, 0.0f, 101.0f}, axisNch});
304+
if (useMcCentrality)
305+
registryMC.add("hEventVsGenCentrality_truerec", "hEventVsGenCentrality_truerec", HistType::kTH2D, {{101, 0.0f, 101.0f}, {101, 0.0f, 101.0f}});
302306
if (centralityEstimator == Option::kNGlobal)
303307
registryMC.add("hCentralityVsMultiplicity_truerec", "hCentralityVsMultiplicity_truerec", HistType::kTH2D, {{101, 0.0f, 101.0f}, axisNch});
304308
else
@@ -322,27 +326,46 @@ struct CascadeAnalysisLightIonsDerivedData {
322326
// QC Histograms
323327
registryQC.add("hVertexZGen", "hVertexZGen", HistType::kTH1D, {{vertexZAxis}});
324328
// Histograms for xi (mc)
325-
registryMC.add("h2dGenXiMinusVsMultMC_RecoedEvt", "h2dGenXiMinusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
326-
registryMC.add("h2dGenXiPlusVsMultMC_RecoedEvt", "h2dGenXiPlusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
327-
registryMC.add("h2dGenXiMinusVsMultMC", "h2dGenXiMinusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
328-
registryMC.add("h2dGenXiPlusVsMultMC", "h2dGenXiPlusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
329329
registryMC.add("h2dGenXiMinus", "h2dGenXiMinus", HistType::kTH2D, {centAxis, ptAxis});
330330
registryMC.add("h2dGenXiPlus", "h2dGenXiPlus", HistType::kTH2D, {centAxis, ptAxis});
331331

332332
// Histograms for omega (mc)
333-
registryMC.add("h2dGenOmegaMinusVsMultMC_RecoedEvt", "h2dGenOmegaMinusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
334-
registryMC.add("h2dGenOmegaPlusVsMultMC_RecoedEvt", "h2dGenOmegaPlusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
335-
registryMC.add("h2dGenOmegaMinusVsMultMC", "h2dGenOmegaMinusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
336-
registryMC.add("h2dGenOmegaPlusVsMultMC", "h2dGenOmegaPlusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
337333
registryMC.add("h2dGenOmegaMinus", "h2dGenOmegaMinus", HistType::kTH2D, {centAxis, ptAxis});
338334
registryMC.add("h2dGenOmegaPlus", "h2dGenOmegaPlus", HistType::kTH2D, {centAxis, ptAxis});
339335

340-
// Histograms for event loss/splitting
341-
registryMC.add("hGenEvents", "hGenEvents", HistType::kTH2D, {{axisNch}, {4, -0.5f, +3.5f}});
342-
registryMC.get<TH2>(HIST("hGenEvents"))->GetYaxis()->SetBinLabel(1, "All gen. events");
343-
registryMC.get<TH2>(HIST("hGenEvents"))->GetYaxis()->SetBinLabel(2, "All gen. events in INEL > 0");
344-
registryMC.get<TH2>(HIST("hGenEvents"))->GetYaxis()->SetBinLabel(3, "Gen. with at least 1 rec. events");
345-
registryMC.get<TH2>(HIST("hGenEvents"))->GetYaxis()->SetBinLabel(4, "Gen. with at least 1 rec. events in INEL > 0");
336+
if (useMcCentrality) {
337+
registryMC.add("h2dGenXiMinusVsMcCentrality_RecoedEvt", "h2dGenXiMinusVsMcCentrality_RecoedEvt", HistType::kTH2D, {axisCentEstimator, ptAxis});
338+
registryMC.add("h2dGenXiPlusVsMcCentrality_RecoedEvt", "h2dGenXiPlusVsMcCentrality_RecoedEvt", HistType::kTH2D, {axisCentEstimator, ptAxis});
339+
registryMC.add("h2dGenXiMinusVsMcCentrality", "h2dGenXiMinusVsMcCentrality", HistType::kTH2D, {axisCentEstimator, ptAxis});
340+
registryMC.add("h2dGenXiPlusVsMcCentrality", "h2dGenXiPlusVsMcCentrality", HistType::kTH2D, {axisCentEstimator, ptAxis});
341+
342+
registryMC.add("h2dGenOmegaMinusVsMcCentrality_RecoedEvt", "h2dGenOmegaMinusVsMcCentrality_RecoedEvt", HistType::kTH2D, {axisCentEstimator, ptAxis});
343+
registryMC.add("h2dGenOmegaPlusVsMcCentrality_RecoedEvt", "h2dGenOmegaPlusVsMcCentrality_RecoedEvt", HistType::kTH2D, {axisCentEstimator, ptAxis});
344+
registryMC.add("h2dGenOmegaMinusVsMcCentrality", "h2dGenOmegaMinusVsMcCentrality", HistType::kTH2D, {axisCentEstimator, ptAxis});
345+
registryMC.add("h2dGenOmegaPlusVsMcCentrality", "h2dGenOmegaPlusVsMcCentrality", HistType::kTH2D, {axisCentEstimator, ptAxis});
346+
347+
registryMC.add("hGenEventsVsMcCentrality", "hGenEventsVsMcCentrality", HistType::kTH2D, {axisCentEstimator, {4, -0.5f, +3.5f}});
348+
registryMC.get<TH2>(HIST("hGenEventsVsMcCentrality"))->GetYaxis()->SetBinLabel(1, "All gen. events");
349+
registryMC.get<TH2>(HIST("hGenEventsVsMcCentrality"))->GetYaxis()->SetBinLabel(2, "All gen. events in INEL > 0");
350+
registryMC.get<TH2>(HIST("hGenEventsVsMcCentrality"))->GetYaxis()->SetBinLabel(3, "Gen. with at least 1 rec. events");
351+
registryMC.get<TH2>(HIST("hGenEventsVsMcCentrality"))->GetYaxis()->SetBinLabel(4, "Gen. with at least 1 rec. events in INEL > 0");
352+
} else {
353+
registryMC.add("h2dGenXiMinusVsMultMC_RecoedEvt", "h2dGenXiMinusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
354+
registryMC.add("h2dGenXiPlusVsMultMC_RecoedEvt", "h2dGenXiPlusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
355+
registryMC.add("h2dGenXiMinusVsMultMC", "h2dGenXiMinusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
356+
registryMC.add("h2dGenXiPlusVsMultMC", "h2dGenXiPlusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
357+
358+
registryMC.add("h2dGenOmegaMinusVsMultMC_RecoedEvt", "h2dGenOmegaMinusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
359+
registryMC.add("h2dGenOmegaPlusVsMultMC_RecoedEvt", "h2dGenOmegaPlusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
360+
registryMC.add("h2dGenOmegaMinusVsMultMC", "h2dGenOmegaMinusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
361+
registryMC.add("h2dGenOmegaPlusVsMultMC", "h2dGenOmegaPlusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
362+
363+
registryMC.add("hGenEvents", "hGenEvents", HistType::kTH2D, {axisNch, {4, -0.5f, +3.5f}});
364+
registryMC.get<TH2>(HIST("hGenEvents"))->GetYaxis()->SetBinLabel(1, "All gen. events");
365+
registryMC.get<TH2>(HIST("hGenEvents"))->GetYaxis()->SetBinLabel(2, "All gen. events in INEL > 0");
366+
registryMC.get<TH2>(HIST("hGenEvents"))->GetYaxis()->SetBinLabel(3, "Gen. with at least 1 rec. events");
367+
registryMC.get<TH2>(HIST("hGenEvents"))->GetYaxis()->SetBinLabel(4, "Gen. with at least 1 rec. events in INEL > 0");
368+
}
346369
registryMC.add("hGenEventCentrality", "hGenEventCentrality", kTH1D, {{101, 0.0f, 101.0f}});
347370

348371
registryMC.add("hCentralityVsNcoll_beforeEvSel", "hCentralityVsNcoll_beforeEvSel", HistType::kTH2D, {centAxis, {50, -0.5f, 49.5f}});
@@ -684,12 +707,29 @@ struct CascadeAnalysisLightIonsDerivedData {
684707
if (applyZVtxSelOnMCPV && std::fabs(mcCollision.posZ()) > zVtx)
685708
continue;
686709

687-
registryMC.fill(HIST("hGenEvents"), mcCollision.multMCNParticlesEta05(), 0 /* all gen. events*/);
710+
float centralityMC = -1.f;
711+
if (centralityEstimator == Option::kFT0C)
712+
centralityMC = mcCollision.centFT0C();
713+
if (centralityEstimator == Option::kFT0M)
714+
centralityMC = mcCollision.centFT0M();
715+
if (centralityEstimator == Option::kFV0A)
716+
centralityMC = mcCollision.centFV0A();
717+
if (centralityEstimator == Option::kNGlobal)
718+
centralityMC = mcCollision.centNGlobal();
719+
720+
if (useMcCentrality)
721+
registryMC.fill(HIST("hGenEventsVsMcCentrality"), centralityMC, 0 /* all gen. events*/);
722+
else
723+
registryMC.fill(HIST("hGenEvents"), mcCollision.multMCNParticlesEta05(), 0 /* all gen. events*/);
688724

689725
if (requireInel0OnMC && mcCollision.multMCNParticlesEta10() < 1) {
690726
continue;
691727
}
692-
registryMC.fill(HIST("hGenEvents"), mcCollision.multMCNParticlesEta05(), 1 /* all gen. events in INEL > 0*/);
728+
729+
if (useMcCentrality)
730+
registryMC.fill(HIST("hGenEventsVsMcCentrality"), centralityMC, 1 /* all gen. events in INEL > 0*/);
731+
else
732+
registryMC.fill(HIST("hGenEvents"), mcCollision.multMCNParticlesEta05(), 1 /* all gen. events in INEL > 0*/);
693733

694734
auto groupedCollisions = getGroupedCollisions(collisions, mcCollision.globalIndex());
695735
// Check if there is at least one of the reconstructed collisions associated to this MC collision
@@ -740,14 +780,18 @@ struct CascadeAnalysisLightIonsDerivedData {
740780

741781
if (biggestNContribs < collision.multPVTotalContributors()) {
742782
biggestNContribs = collision.multPVTotalContributors();
743-
if (centralityEstimator == Option::kFT0C)
744-
centralitydata = collision.centFT0C();
745-
if (centralityEstimator == Option::kFT0M)
746-
centralitydata = collision.centFT0M();
747-
if (centralityEstimator == Option::kFV0A)
748-
centralitydata = collision.centFV0A();
749-
if (centralityEstimator == Option::kNGlobal)
750-
centralitydata = collision.centNGlobal();
783+
if (useMcCentrality) {
784+
centralitydata = centralityMC;
785+
} else {
786+
if (centralityEstimator == Option::kFT0C)
787+
centralitydata = collision.centFT0C();
788+
if (centralityEstimator == Option::kFT0M)
789+
centralitydata = collision.centFT0M();
790+
if (centralityEstimator == Option::kFV0A)
791+
centralitydata = collision.centFV0A();
792+
if (centralityEstimator == Option::kNGlobal)
793+
centralitydata = collision.centNGlobal();
794+
}
751795
}
752796
nCollisions++;
753797

@@ -761,8 +805,13 @@ struct CascadeAnalysisLightIonsDerivedData {
761805
registryQC.fill(HIST("hVertexZGen"), mcCollision.posZ());
762806

763807
if (atLeastOne) {
764-
registryMC.fill(HIST("hGenEvents"), mcCollision.multMCNParticlesEta05(), 2 /* at least 1 rec. event*/);
765-
registryMC.fill(HIST("hGenEvents"), mcCollision.multMCNParticlesEta05(), 3 /* at least 1 rec. event in INEL > 0*/);
808+
if (useMcCentrality) {
809+
registryMC.fill(HIST("hGenEventsVsMcCentrality"), centralityMC, 2 /* at least 1 rec. event*/);
810+
registryMC.fill(HIST("hGenEventsVsMcCentrality"), centralityMC, 3 /* at least 1 rec. event in INEL > 0*/);
811+
} else {
812+
registryMC.fill(HIST("hGenEvents"), mcCollision.multMCNParticlesEta05(), 2 /* at least 1 rec. event*/);
813+
registryMC.fill(HIST("hGenEvents"), mcCollision.multMCNParticlesEta05(), 3 /* at least 1 rec. event in INEL > 0*/);
814+
}
766815
registryMC.fill(HIST("hGenEventCentrality"), centralitydata);
767816
}
768817
}
@@ -955,6 +1004,25 @@ struct CascadeAnalysisLightIonsDerivedData {
9551004
multiplicityMcRec = RecCol.multNTracksGlobal();
9561005
}
9571006

1007+
if (useMcCentrality) {
1008+
const float recoCentrality = centralityMcRec;
1009+
centralityMcRec = -1.f;
1010+
1011+
if (RecCol.has_straMCCollision()) {
1012+
auto mcCollision = RecCol.template straMCCollision_as<CollisionMCTrueTable>();
1013+
if (centralityEstimator == Option::kFT0C)
1014+
centralityMcRec = mcCollision.centFT0C();
1015+
if (centralityEstimator == Option::kFT0M)
1016+
centralityMcRec = mcCollision.centFT0M();
1017+
if (centralityEstimator == Option::kFV0A)
1018+
centralityMcRec = mcCollision.centFV0A();
1019+
if (centralityEstimator == Option::kNGlobal)
1020+
centralityMcRec = mcCollision.centNGlobal();
1021+
}
1022+
1023+
registryMC.fill(HIST("hEventVsGenCentrality_truerec"), recoCentrality, centralityMcRec);
1024+
}
1025+
9581026
registryMC.fill(HIST("number_of_events_mc_rec_vs_centrality"), 0, centralityMcRec);
9591027

9601028
// event selections
@@ -1148,47 +1216,86 @@ struct CascadeAnalysisLightIonsDerivedData {
11481216
}
11491217

11501218
float centralityMC = 100.5f;
1219+
float centralityMCtrue = -1.f;
1220+
1221+
if (centralityEstimator == Option::kFT0C)
1222+
centralityMCtrue = mcCollision.centFT0C();
1223+
if (centralityEstimator == Option::kFT0M)
1224+
centralityMCtrue = mcCollision.centFT0M();
1225+
if (centralityEstimator == Option::kFV0A)
1226+
centralityMCtrue = mcCollision.centFV0A();
1227+
if (centralityEstimator == Option::kNGlobal)
1228+
centralityMCtrue = mcCollision.centNGlobal();
11511229

11521230
if (listBestCollisionIdx[mcCollision.globalIndex()] > -1) {
11531231
auto collision = RecCols.iteratorAt(listBestCollisionIdx[mcCollision.globalIndex()]);
1154-
if (centralityEstimator == Option::kFT0C)
1155-
centralityMC = collision.centFT0C();
1156-
if (centralityEstimator == Option::kFT0M)
1157-
centralityMC = collision.centFT0M();
1158-
if (centralityEstimator == Option::kFV0A)
1159-
centralityMC = collision.centFV0A();
1160-
if (centralityEstimator == Option::kNGlobal)
1161-
centralityMC = collision.centNGlobal();
1232+
1233+
if (useMcCentrality) {
1234+
centralityMC = centralityMCtrue;
1235+
} else {
1236+
if (centralityEstimator == Option::kFT0C)
1237+
centralityMC = collision.centFT0C();
1238+
if (centralityEstimator == Option::kFT0M)
1239+
centralityMC = collision.centFT0M();
1240+
if (centralityEstimator == Option::kFV0A)
1241+
centralityMC = collision.centFV0A();
1242+
if (centralityEstimator == Option::kNGlobal)
1243+
centralityMC = collision.centNGlobal();
1244+
}
11621245

11631246
if (cascMC.pdgCode() == kXiMinus && std::abs(cascMC.rapidityMC(0)) < rapcut) {
1164-
registryMC.fill(HIST("h2dGenXiMinusVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), ptmc);
1247+
if (useMcCentrality)
1248+
registryMC.fill(HIST("h2dGenXiMinusVsMcCentrality_RecoedEvt"), centralityMCtrue, ptmc);
1249+
else
1250+
registryMC.fill(HIST("h2dGenXiMinusVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), ptmc);
11651251
}
11661252
if (cascMC.pdgCode() == kXiPlusBar && std::abs(cascMC.rapidityMC(0)) < rapcut) {
1167-
registryMC.fill(HIST("h2dGenXiPlusVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), ptmc);
1253+
if (useMcCentrality)
1254+
registryMC.fill(HIST("h2dGenXiPlusVsMcCentrality_RecoedEvt"), centralityMCtrue, ptmc);
1255+
else
1256+
registryMC.fill(HIST("h2dGenXiPlusVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), ptmc);
11681257
}
11691258
if (cascMC.pdgCode() == kOmegaMinus && std::abs(cascMC.rapidityMC(2)) < rapcut) {
1170-
registryMC.fill(HIST("h2dGenOmegaMinusVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), ptmc);
1259+
if (useMcCentrality)
1260+
registryMC.fill(HIST("h2dGenOmegaMinusVsMcCentrality_RecoedEvt"), centralityMCtrue, ptmc);
1261+
else
1262+
registryMC.fill(HIST("h2dGenOmegaMinusVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), ptmc);
11711263
}
11721264
if (cascMC.pdgCode() == kOmegaPlusBar && std::abs(cascMC.rapidityMC(2)) < rapcut) {
1173-
registryMC.fill(HIST("h2dGenOmegaPlusVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), ptmc);
1265+
if (useMcCentrality)
1266+
registryMC.fill(HIST("h2dGenOmegaPlusVsMcCentrality_RecoedEvt"), centralityMCtrue, ptmc);
1267+
else
1268+
registryMC.fill(HIST("h2dGenOmegaPlusVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), ptmc);
11741269
}
11751270
}
11761271

11771272
if (cascMC.pdgCode() == kXiMinus && std::abs(cascMC.rapidityMC(0)) < rapcut) {
11781273
registryMC.fill(HIST("h2dGenXiMinus"), centralityMC, ptmc);
1179-
registryMC.fill(HIST("h2dGenXiMinusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
1274+
if (useMcCentrality)
1275+
registryMC.fill(HIST("h2dGenXiMinusVsMcCentrality"), centralityMCtrue, ptmc);
1276+
else
1277+
registryMC.fill(HIST("h2dGenXiMinusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
11801278
}
11811279
if (cascMC.pdgCode() == kXiPlusBar && std::abs(cascMC.rapidityMC(0)) < rapcut) {
11821280
registryMC.fill(HIST("h2dGenXiPlus"), centralityMC, ptmc);
1183-
registryMC.fill(HIST("h2dGenXiPlusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
1281+
if (useMcCentrality)
1282+
registryMC.fill(HIST("h2dGenXiPlusVsMcCentrality"), centralityMCtrue, ptmc);
1283+
else
1284+
registryMC.fill(HIST("h2dGenXiPlusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
11841285
}
11851286
if (cascMC.pdgCode() == kOmegaMinus && std::abs(cascMC.rapidityMC(2)) < rapcut) {
11861287
registryMC.fill(HIST("h2dGenOmegaMinus"), centralityMC, ptmc);
1187-
registryMC.fill(HIST("h2dGenOmegaMinusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
1288+
if (useMcCentrality)
1289+
registryMC.fill(HIST("h2dGenOmegaMinusVsMcCentrality"), centralityMCtrue, ptmc);
1290+
else
1291+
registryMC.fill(HIST("h2dGenOmegaMinusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
11881292
}
11891293
if (cascMC.pdgCode() == kOmegaPlusBar && std::abs(cascMC.rapidityMC(2)) < rapcut) {
11901294
registryMC.fill(HIST("h2dGenOmegaPlus"), centralityMC, ptmc);
1191-
registryMC.fill(HIST("h2dGenOmegaPlusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
1295+
if (useMcCentrality)
1296+
registryMC.fill(HIST("h2dGenOmegaPlusVsMcCentrality"), centralityMCtrue, ptmc);
1297+
else
1298+
registryMC.fill(HIST("h2dGenOmegaPlusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
11921299
}
11931300
} // cascMC loop
11941301
}

0 commit comments

Comments
 (0)