|
1 | | -//Code modified to obtain histograms for subsample division |
2 | | - |
3 | 1 | // Copyright 2019-2020 CERN and copyright holders of ALICE O2. |
4 | 2 | // See https://alice-o2.web.cern.ch/copyright for details of the copyright holders. |
5 | 3 | // All rights not expressly granted are reserved. |
@@ -98,7 +96,7 @@ struct AntinucleiInJets { |
98 | 96 | HistogramRegistry registryQC{"registryQC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; |
99 | 97 | HistogramRegistry registryMult{"registryMult", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; |
100 | 98 | HistogramRegistry registryCorr{"registryCorr", {}, OutputObjHandlingPolicy::AnalysisObject, true, true}; |
101 | | - |
| 99 | + |
102 | 100 | // Random generator for subsample assignment |
103 | 101 | TRandom3 mRand; |
104 | 102 |
|
@@ -210,7 +208,7 @@ struct AntinucleiInJets { |
210 | 208 | if (setMCDefaultItsParams) { |
211 | 209 | itsResponse.setMCDefaultParameters(); |
212 | 210 | } |
213 | | - |
| 211 | + |
214 | 212 | // Initialize random seed using high-resolution clock to ensure unique sequences across parallel Grid jobs |
215 | 213 | auto time_seed = std::chrono::high_resolution_clock::now().time_since_epoch().count(); |
216 | 214 | mRand.SetSeed(time_seed); |
@@ -469,39 +467,37 @@ struct AntinucleiInJets { |
469 | 467 | const AxisSpec nBarP2Axis{100, 0.0, 100.0, "N_{#bar{p}}^{i} #times N_{#bar{p}}^{j}"}; |
470 | 468 | const AxisSpec nBarDnBarPAxis{100, 0.0, 100.0, "N_{#bar{d}}^{i} #times N_{#bar{p}}^{j}"}; |
471 | 469 | const AxisSpec subsampleAxis{20, 0, 20, "Subsample Index"}; |
472 | | - |
473 | 470 |
|
474 | 471 | // Event counter |
475 | 472 | registryCorr.add("eventCounter", "number of events", HistType::kTH1F, {{20, 0, 20, "counter"}}); |
476 | | - |
477 | 473 | registryCorr.add("eventCounter_centrality_fullEvent", "Number of events per centrality (Full Event)", HistType::kTH2F, {multiplicityAxis, subsampleAxis}); |
478 | 474 | //registryCorr.add("eventCounter_centrality_jet", "Number of events per centrality (Jet)", HistType::kTH1F, {multiplicityAxis}); |
479 | 475 | //registryCorr.add("eventCounter_centrality_ue", "Number of events per centrality (Underlying Event)", HistType::kTH1F, {multiplicityAxis}); |
480 | | - |
| 476 | + |
481 | 477 | // Correlation histograms: antiproton vs. antideuteron number vs. event multiplicity |
482 | 478 | //registryCorr.add("rho_jet", "rho_jet", HistType::kTH3F, {nAntideuteronsAxis, nAntiprotonsAxis, multiplicityAxis}); |
483 | 479 | //registryCorr.add("rho_ue", "rho_ue", HistType::kTH3F, {nAntideuteronsAxis, nAntiprotonsAxis, multiplicityAxis}); |
484 | 480 | registryCorr.add("rho_fullEvent", "rho_fullEvent", HistType::kTHnSparseD, {nAntideuteronsAxis, nAntiprotonsAxis, multiplicityAxis, subsampleAxis}); |
485 | | - |
| 481 | + |
486 | 482 | // Correlation histograms: net antiproton vs. net antideuteron numbers |
487 | 483 | //registryCorr.add("rho_netP_netD_jet", "rho_netP_netD_jet", HistType::kTH2F, {nAntideuteronsAxis, nAntiprotonsAxis}); |
488 | 484 | //registryCorr.add("rho_netP_netD_ue", "rho_netP_netD_ue", HistType::kTH2F, {nAntideuteronsAxis, nAntiprotonsAxis}); |
489 | 485 | registryCorr.add("rho_netP_netD_fullEvent", "rho_netP_netD_fullEvent", HistType::kTH3F, {nAntideuteronsAxis, nAntiprotonsAxis, subsampleAxis}); |
490 | | - |
| 486 | + |
491 | 487 | // Efficiency histograms jet |
492 | 488 | //registryCorr.add("q1d_jet", "q1d_jet", HistType::kTH3F, {nAntideuteronsAxis, ptPerNucleonAxis, multiplicityAxis}); |
493 | 489 | //registryCorr.add("q1p_jet", "q1p_jet", HistType::kTH3F, {nAntiprotonsAxis, ptPerNucleonAxis, multiplicityAxis}); |
494 | 490 | //registryCorr.add("q1d_square_jet", "q1d_square_jet", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarD2Axis, multiplicityAxis}); |
495 | 491 | //registryCorr.add("q1p_square_jet", "q1p_square_jet", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarP2Axis, multiplicityAxis}); |
496 | 492 | //registryCorr.add("q1d_q1p_jet", "q1d_q1p_jet", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarDnBarPAxis, multiplicityAxis}); |
497 | 493 |
|
498 | | - // Efficiency histograms UE |
499 | | - //registryCorr.add("q1d_ue", "q1d_ue", HistType::kTH3F, {nAntideuteronsAxis, ptPerNucleonAxis, multiplicityAxis}); |
500 | | - //registryCorr.add("q1p_ue", "q1p_ue", HistType::kTH3F, {nAntiprotonsAxis, ptPerNucleonAxis, multiplicityAxis}); |
501 | | - //registryCorr.add("q1d_square_ue", "q1d_square_ue", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarD2Axis, multiplicityAxis}); |
502 | | - //registryCorr.add("q1p_square_ue", "q1p_square_ue", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarP2Axis, multiplicityAxis}); |
503 | | - //registryCorr.add("q1d_q1p_ue", "q1d_q1p_ue", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarDnBarPAxis, multiplicityAxis}); |
504 | | - |
| 494 | + // Efficiency histograms UE |
| 495 | + //registryCorr.add("q1d_ue", "q1d_ue", HistType::kTH3F, {nAntideuteronsAxis, ptPerNucleonAxis, multiplicityAxis}); |
| 496 | + //registryCorr.add("q1p_ue", "q1p_ue", HistType::kTH3F, {nAntiprotonsAxis, ptPerNucleonAxis, multiplicityAxis}); |
| 497 | + //registryCorr.add("q1d_square_ue", "q1d_square_ue", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarD2Axis, multiplicityAxis}); |
| 498 | + //registryCorr.add("q1p_square_ue", "q1p_square_ue", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarP2Axis, multiplicityAxis}); |
| 499 | + //registryCorr.add("q1d_q1p_ue", "q1d_q1p_ue", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarDnBarPAxis, multiplicityAxis}); |
| 500 | + |
505 | 501 | // Efficiency histograms full event |
506 | 502 | registryCorr.add("q1d_fullEvent", "q1d_fullEvent", HistType::kTHnSparseD, {nAntideuteronsAxis, ptPerNucleonAxis, multiplicityAxis, subsampleAxis}); |
507 | 503 | registryCorr.add("q1p_fullEvent", "q1p_fullEvent", HistType::kTHnSparseD, {nAntiprotonsAxis, ptPerNucleonAxis, multiplicityAxis, subsampleAxis}); |
|
0 commit comments