1- // Code modified to obtain histograms for subsample division
1+ // Code modified to obtain histograms for subsample division
22
33// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
44// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
6666#include < fastjet/tools/JetMedianBackgroundEstimator.hh>
6767#include < fastjet/tools/Subtractor.hh>
6868
69+ #include < chrono>
6970#include < cmath>
7071#include < memory>
7172#include < random>
7273#include < string>
7374#include < vector>
74- #include < chrono>
7575
7676using namespace std ;
7777using namespace o2 ;
@@ -98,7 +98,7 @@ struct AntinucleiInJets {
9898 HistogramRegistry registryQC{" registryQC" , {}, OutputObjHandlingPolicy::AnalysisObject, true , true };
9999 HistogramRegistry registryMult{" registryMult" , {}, OutputObjHandlingPolicy::AnalysisObject, true , true };
100100 HistogramRegistry registryCorr{" registryCorr" , {}, OutputObjHandlingPolicy::AnalysisObject, true , true };
101-
101+
102102 // Random generator for subsample assignment
103103 TRandom3 mRand ;
104104
@@ -210,7 +210,7 @@ struct AntinucleiInJets {
210210 if (setMCDefaultItsParams) {
211211 itsResponse.setMCDefaultParameters ();
212212 }
213-
213+
214214 // Initialize random seed using high-resolution clock to ensure unique sequences across parallel Grid jobs
215215 auto time_seed = std::chrono::high_resolution_clock::now ().time_since_epoch ().count ();
216216 mRand .SetSeed (time_seed);
@@ -469,39 +469,38 @@ struct AntinucleiInJets {
469469 const AxisSpec nBarP2Axis{100 , 0.0 , 100.0 , " N_{#bar{p}}^{i} #times N_{#bar{p}}^{j}" };
470470 const AxisSpec nBarDnBarPAxis{100 , 0.0 , 100.0 , " N_{#bar{d}}^{i} #times N_{#bar{p}}^{j}" };
471471 const AxisSpec subsampleAxis{20 , 0 , 20 , " Subsample Index" };
472-
473472
474473 // Event counter
475474 registryCorr.add (" eventCounter" , " number of events" , HistType::kTH1F , {{20 , 0 , 20 , " counter" }});
476-
475+
477476 registryCorr.add (" eventCounter_centrality_fullEvent" , " Number of events per centrality (Full Event)" , HistType::kTH2F , {multiplicityAxis, subsampleAxis});
478- // registryCorr.add("eventCounter_centrality_jet", "Number of events per centrality (Jet)", HistType::kTH1F, {multiplicityAxis});
479- // registryCorr.add("eventCounter_centrality_ue", "Number of events per centrality (Underlying Event)", HistType::kTH1F, {multiplicityAxis});
480-
477+ // registryCorr.add("eventCounter_centrality_jet", "Number of events per centrality (Jet)", HistType::kTH1F, {multiplicityAxis});
478+ // registryCorr.add("eventCounter_centrality_ue", "Number of events per centrality (Underlying Event)", HistType::kTH1F, {multiplicityAxis});
479+
481480 // Correlation histograms: antiproton vs. antideuteron number vs. event multiplicity
482- // registryCorr.add("rho_jet", "rho_jet", HistType::kTH3F, {nAntideuteronsAxis, nAntiprotonsAxis, multiplicityAxis});
483- // registryCorr.add("rho_ue", "rho_ue", HistType::kTH3F, {nAntideuteronsAxis, nAntiprotonsAxis, multiplicityAxis});
481+ // registryCorr.add("rho_jet", "rho_jet", HistType::kTH3F, {nAntideuteronsAxis, nAntiprotonsAxis, multiplicityAxis});
482+ // registryCorr.add("rho_ue", "rho_ue", HistType::kTH3F, {nAntideuteronsAxis, nAntiprotonsAxis, multiplicityAxis});
484483 registryCorr.add (" rho_fullEvent" , " rho_fullEvent" , HistType::kTHnSparseD , {nAntideuteronsAxis, nAntiprotonsAxis, multiplicityAxis, subsampleAxis});
485-
484+
486485 // Correlation histograms: net antiproton vs. net antideuteron numbers
487- // registryCorr.add("rho_netP_netD_jet", "rho_netP_netD_jet", HistType::kTH2F, {nAntideuteronsAxis, nAntiprotonsAxis});
488- // registryCorr.add("rho_netP_netD_ue", "rho_netP_netD_ue", HistType::kTH2F, {nAntideuteronsAxis, nAntiprotonsAxis});
486+ // registryCorr.add("rho_netP_netD_jet", "rho_netP_netD_jet", HistType::kTH2F, {nAntideuteronsAxis, nAntiprotonsAxis});
487+ // registryCorr.add("rho_netP_netD_ue", "rho_netP_netD_ue", HistType::kTH2F, {nAntideuteronsAxis, nAntiprotonsAxis});
489488 registryCorr.add (" rho_netP_netD_fullEvent" , " rho_netP_netD_fullEvent" , HistType::kTH3F , {nAntideuteronsAxis, nAntiprotonsAxis, subsampleAxis});
490-
489+
491490 // Efficiency histograms jet
492- // registryCorr.add("q1d_jet", "q1d_jet", HistType::kTH3F, {nAntideuteronsAxis, ptPerNucleonAxis, multiplicityAxis});
493- // registryCorr.add("q1p_jet", "q1p_jet", HistType::kTH3F, {nAntiprotonsAxis, ptPerNucleonAxis, multiplicityAxis});
494- // registryCorr.add("q1d_square_jet", "q1d_square_jet", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarD2Axis, multiplicityAxis});
495- // registryCorr.add("q1p_square_jet", "q1p_square_jet", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarP2Axis, multiplicityAxis});
496- // registryCorr.add("q1d_q1p_jet", "q1d_q1p_jet", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarDnBarPAxis, multiplicityAxis});
497-
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-
491+ // registryCorr.add("q1d_jet", "q1d_jet", HistType::kTH3F, {nAntideuteronsAxis, ptPerNucleonAxis, multiplicityAxis});
492+ // registryCorr.add("q1p_jet", "q1p_jet", HistType::kTH3F, {nAntiprotonsAxis, ptPerNucleonAxis, multiplicityAxis});
493+ // registryCorr.add("q1d_square_jet", "q1d_square_jet", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarD2Axis, multiplicityAxis});
494+ // registryCorr.add("q1p_square_jet", "q1p_square_jet", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarP2Axis, multiplicityAxis});
495+ // registryCorr.add("q1d_q1p_jet", "q1d_q1p_jet", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarDnBarPAxis, multiplicityAxis});
496+
497+ // Efficiency histograms UE
498+ // registryCorr.add("q1d_ue", "q1d_ue", HistType::kTH3F, {nAntideuteronsAxis, ptPerNucleonAxis, multiplicityAxis});
499+ // registryCorr.add("q1p_ue", "q1p_ue", HistType::kTH3F, {nAntiprotonsAxis, ptPerNucleonAxis, multiplicityAxis});
500+ // registryCorr.add("q1d_square_ue", "q1d_square_ue", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarD2Axis, multiplicityAxis});
501+ // registryCorr.add("q1p_square_ue", "q1p_square_ue", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarP2Axis, multiplicityAxis});
502+ // registryCorr.add("q1d_q1p_ue", "q1d_q1p_ue", HistType::kTHnSparseD, {ptPerNucleonAxis, ptPerNucleonAxis, nBarDnBarPAxis, multiplicityAxis});
503+
505504 // Efficiency histograms full event
506505 registryCorr.add (" q1d_fullEvent" , " q1d_fullEvent" , HistType::kTHnSparseD , {nAntideuteronsAxis, ptPerNucleonAxis, multiplicityAxis, subsampleAxis});
507506 registryCorr.add (" q1p_fullEvent" , " q1p_fullEvent" , HistType::kTHnSparseD , {nAntiprotonsAxis, ptPerNucleonAxis, multiplicityAxis, subsampleAxis});
0 commit comments