Skip to content

Commit 6476673

Browse files
spucillospucillo
authored
Add 2D histograms for corrections + improve binning
1 parent dae1154 commit 6476673

File tree

1 file changed

+26
-14
lines changed

1 file changed

+26
-14
lines changed

‎PWGLF/Tasks/Strangeness/cascadeAnalysisLightIonsDerivedData.cxx‎

Lines changed: 26 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -185,7 +185,7 @@ struct CascadeAnalysisLightIonsDerivedData {
185185
registryData.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(9, "kIsVertexTRDmatched");
186186
registryData.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(10, "kNoSameBunchPileup");
187187

188-
registryData.add("number_of_events_data_vs_centrality", "number of events in data vs centrality", HistType::kTH2D, {{20, -0.5f, +19.5f}, axisCentEstimator});
188+
registryData.add("number_of_events_data_vs_centrality", "number of events in data vs centrality", HistType::kTH2D, {{20, -0.5f, +19.5f}, {101, 0.0f, 101.0f});
189189
registryData.get<TH2>(HIST("number_of_events_data_vs_centrality"))->GetXaxis()->SetBinLabel(1, "All collisions");
190190
registryData.get<TH2>(HIST("number_of_events_data_vs_centrality"))->GetXaxis()->SetBinLabel(2, "sel8 cut");
191191
registryData.get<TH2>(HIST("number_of_events_data_vs_centrality"))->GetXaxis()->SetBinLabel(3, "posZ cut");
@@ -212,9 +212,9 @@ struct CascadeAnalysisLightIonsDerivedData {
212212
registryQC.add("hdcav0topvdata", "hdcav0topvdata", HistType::kTH1D, {{nBins, 0.0f, 2.0f}});
213213

214214
// Multiplicity Histograms
215-
registryData.add("hCentEstimator", "hCentEstimator", HistType::kTH1D, {axisCentEstimator});
216-
registryData.add("hCentralityVsNch", "hCentralityVsNch", HistType::kTH2D, {axisCentEstimator, axisNch});
217-
registryData.add("hCentralityVsMultiplicity", "hCentralityVsMultiplicity", HistType::kTH2D, {axisCentEstimator, axisMult});
215+
registryData.add("hCentEstimator", "hCentEstimator", HistType::kTH1D, {101, 0.0f, 101.0f});
216+
registryData.add("hCentralityVsNch", "hCentralityVsNch", HistType::kTH2D, {101, 0.0f, 101.0f, axisNch});
217+
registryData.add("hCentralityVsMultiplicity", "hCentralityVsMultiplicity", HistType::kTH2D, {{101, 0.0f, 101.0f}, axisMult});
218218

219219
// Histograms for xi (data)
220220
registryData.add("hMassXipos", "hMassXipos", HistType::kTH3D, {centAxis, ptAxis, invMassXiAxis});
@@ -243,7 +243,7 @@ struct CascadeAnalysisLightIonsDerivedData {
243243
registryMC.get<TH1>(HIST("number_of_events_mc_rec"))->GetXaxis()->SetBinLabel(9, "kIsVertexTRDmatched");
244244
registryMC.get<TH1>(HIST("number_of_events_mc_rec"))->GetXaxis()->SetBinLabel(10, "kNoSameBunchPileup");
245245

246-
registryMC.add("number_of_events_mc_rec_vs_centrality", "number of events in mc_rec vs centrality", HistType::kTH2D, {{20, -0.5f, +19.5f}, axisCentEstimator});
246+
registryMC.add("number_of_events_mc_rec_vs_centrality", "number of events in mc_rec vs centrality", HistType::kTH2D, {{20, -0.5f, +19.5f}, {101, 0.0f, 101.0f}});
247247
registryMC.get<TH2>(HIST("number_of_events_mc_rec_vs_centrality"))->GetXaxis()->SetBinLabel(1, "All collisions");
248248
registryMC.get<TH2>(HIST("number_of_events_mc_rec_vs_centrality"))->GetXaxis()->SetBinLabel(2, "sel8 cut");
249249
registryMC.get<TH2>(HIST("number_of_events_mc_rec_vs_centrality"))->GetXaxis()->SetBinLabel(3, "posZ cut");
@@ -270,9 +270,9 @@ struct CascadeAnalysisLightIonsDerivedData {
270270
registryQC.add("hdcav0topvRec", "hdcav0topvRec", HistType::kTH1D, {{nBins, 0.0f, 2.0f}});
271271

272272
// Multiplicity Histograms
273-
registryMC.add("hCentEstimator_truerec", "hCentEstimator_truerec", HistType::kTH1D, {axisCentEstimator});
274-
registryMC.add("hCentralityVsNch_truerec", "hCentralityVsNch_truerec", HistType::kTH2D, {axisCentEstimator, axisNch});
275-
registryMC.add("hCentralityVsMultiplicity_truerec", "hCentralityVsMultiplicity_truerec", HistType::kTH2D, {axisCentEstimator, axisMult});
273+
registryMC.add("hCentEstimator_truerec", "hCentEstimator_truerec", HistType::kTH1D, {101, 0.0f, 101.0f});
274+
registryMC.add("hCentralityVsNch_truerec", "hCentralityVsNch_truerec", HistType::kTH2D, {{101, 0.0f, 101.0f}, axisNch});
275+
registryMC.add("hCentralityVsMultiplicity_truerec", "hCentralityVsMultiplicity_truerec", HistType::kTH2D, {{101, 0.0f, 101.0f}, axisMult});
276276

277277
// Histograms for xi (mc)
278278
registryMC.add("hMassXipos_truerec", "hMassXipos_truerec", HistType::kTH3D, {centAxis, ptAxis, invMassXiAxis});
@@ -296,12 +296,16 @@ struct CascadeAnalysisLightIonsDerivedData {
296296
registryMC.add("h2dGenXiPlusVsMultMC_RecoedEvt", "h2dGenXiPlusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
297297
registryMC.add("h2dGenXiMinusVsMultMC", "h2dGenXiMinusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
298298
registryMC.add("h2dGenXiPlusVsMultMC", "h2dGenXiPlusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
299+
registryMC.add("h2dGenXiMinus", "h2dGenXiMinus", HistType::kTH2D, {{101, 0.0f, 101.0f}, ptAxis});
300+
registryMC.add("h2dGenXiPlus", "h2dGenXiPlus", HistType::kTH2D, {{101, 0.0f, 101.0f}, ptAxis});
299301

300302
// Histograms for omega (mc)
301303
registryMC.add("h2dGenOmegaMinusVsMultMC_RecoedEvt", "h2dGenOmegaMinusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
302304
registryMC.add("h2dGenOmegaPlusVsMultMC_RecoedEvt", "h2dGenOmegaPlusVsMultMC_RecoedEvt", HistType::kTH2D, {axisNch, ptAxis});
303305
registryMC.add("h2dGenOmegaMinusVsMultMC", "h2dGenOmegaMinusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
304306
registryMC.add("h2dGenOmegaPlusVsMultMC", "h2dGenOmegaPlusVsMultMC", HistType::kTH2D, {axisNch, ptAxis});
307+
registryMC.add("h2dGenOmegaMinus", "h2dGenOmegaMinus", HistType::kTH2D, {{101, 0.0f, 101.0f}, ptAxis});
308+
registryMC.add("h2dGenOmegaPlus", "h2dGenOmegaPlus", HistType::kTH2D, {{101, 0.0f, 101.0f}, ptAxis});
305309

306310
// Histograms for event loss/splitting
307311
registryMC.add("hGenEvents", "hGenEvents", HistType::kTH2D, {{axisNch}, {2, -0.5f, +1.5f}});
@@ -751,6 +755,8 @@ struct CascadeAnalysisLightIonsDerivedData {
751755
multiplicity = collision.multNTracksGlobal();
752756
}
753757

758+
registryData.fill(HIST("number_of_events_data_vs_centrality"), 0, centrality);
759+
754760
// event selections
755761
if (applySel8 && !collision.sel8())
756762
return;
@@ -893,6 +899,8 @@ struct CascadeAnalysisLightIonsDerivedData {
893899
multiplicityMcRec = RecCol.multNTracksGlobal();
894900
}
895901

902+
registryMC.fill(HIST("number_of_events_mc_rec_vs_centrality"), 0, centralityMcRec);
903+
896904
// event selections
897905
if (applySel8 && !RecCol.sel8())
898906
continue;
@@ -1053,14 +1061,14 @@ struct CascadeAnalysisLightIonsDerivedData {
10531061
// Store the Zvtx
10541062
registryQC.fill(HIST("hVertexZGen"), mcCollision.posZ());
10551063

1056-
// float centrality = 100.5f;
1064+
float centrality = 100.5f;
10571065

10581066
if (listBestCollisionIdx[mcCollision.globalIndex()] > -1) {
1059-
// auto collision = RecCols.iteratorAt(listBestCollisionIdx[mcCollision.globalIndex()]);
1060-
// if (centralityEstimator == Option::kFT0C) centrality = collision.centFT0C();
1061-
// if (centralityEstimator == Option::kFT0M) centrality = collision.centFT0M();
1062-
// if (centralityEstimator == Option::kFV0A) centrality = collision.centFV0A();
1063-
// if (centralityEstimator == Option::kNGlobal) centrality = collision.centNGlobal();
1067+
auto collision = RecCols.iteratorAt(listBestCollisionIdx[mcCollision.globalIndex()]);
1068+
if (centralityEstimator == Option::kFT0C) centrality = collision.centFT0C();
1069+
if (centralityEstimator == Option::kFT0M) centrality = collision.centFT0M();
1070+
if (centralityEstimator == Option::kFV0A) centrality = collision.centFV0A();
1071+
if (centralityEstimator == Option::kNGlobal) centrality = collision.centNGlobal();
10641072

10651073
if (cascMC.pdgCode() == kXiMinus && std::abs(cascMC.rapidityMC(0)) < rapcut) {
10661074
registryMC.fill(HIST("h2dGenXiMinusVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), ptmc);
@@ -1077,15 +1085,19 @@ struct CascadeAnalysisLightIonsDerivedData {
10771085
}
10781086

10791087
if (cascMC.pdgCode() == kXiMinus && std::abs(cascMC.rapidityMC(0)) < rapcut) {
1088+
registryMC.fill(HIST("h2dGenXiMinus"), centralityMC, ptmc);
10801089
registryMC.fill(HIST("h2dGenXiMinusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
10811090
}
10821091
if (cascMC.pdgCode() == kXiPlusBar && std::abs(cascMC.rapidityMC(0)) < rapcut) {
1092+
registryMC.fill(HIST("h2dGenXiPlus"), centralityMC, ptmc);
10831093
registryMC.fill(HIST("h2dGenXiPlusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
10841094
}
10851095
if (cascMC.pdgCode() == kOmegaMinus && std::abs(cascMC.rapidityMC(2)) < rapcut) {
1096+
registryMC.fill(HIST("h2dGenOmegaMinus"), centralityMC, ptmc);
10861097
registryMC.fill(HIST("h2dGenOmegaMinusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
10871098
}
10881099
if (cascMC.pdgCode() == kOmegaPlusBar && std::abs(cascMC.rapidityMC(2)) < rapcut) {
1100+
registryMC.fill(HIST("h2dGenOmegaPlus"), centralityMC, ptmc);
10891101
registryMC.fill(HIST("h2dGenOmegaPlusVsMultMC"), mcCollision.multMCNParticlesEta05(), ptmc);
10901102
}
10911103
} // cascMC loop

0 commit comments

Comments
 (0)