|
13 | 13 | /// \file multiplicityPt.cxx |
14 | 14 | /// \brief Analysis to do PID with MC - Full correction factors for pions, kaons, protons |
15 | 15 |
|
16 | | -#include "PWGLF/DataModel/LFParticleIdentification.h" |
17 | | -#include "PWGLF/DataModel/mcCentrality.h" |
18 | | -#include "PWGLF/DataModel/spectraTOF.h" |
19 | | -#include "PWGLF/Utils/inelGt.h" |
20 | | - |
21 | | -#include "Common/Core/RecoDecay.h" |
22 | 16 | #include "Common/Core/TrackSelection.h" |
23 | 17 | #include "Common/Core/TrackSelectionDefaults.h" |
24 | 18 | #include "Common/DataModel/Centrality.h" |
25 | 19 | #include "Common/DataModel/EventSelection.h" |
26 | | -#include "Common/DataModel/McCollisionExtra.h" |
27 | 20 | #include "Common/DataModel/Multiplicity.h" |
28 | | -#include "Common/DataModel/PIDResponseTOF.h" |
29 | 21 | #include "Common/DataModel/PIDResponseTPC.h" |
30 | 22 | #include "Common/DataModel/TrackSelectionTables.h" |
31 | 23 |
|
|
34 | 26 | #include <Framework/ASoAHelpers.h> |
35 | 27 | #include <Framework/AnalysisDataModel.h> |
36 | 28 | #include <Framework/AnalysisTask.h> |
| 29 | +#include <Framework/Configurable.h> |
37 | 30 | #include <Framework/HistogramRegistry.h> |
38 | 31 | #include <Framework/Logger.h> |
39 | 32 | #include <Framework/O2DatabasePDGPlugin.h> |
40 | 33 | #include <Framework/runDataProcessing.h> |
41 | 34 | #include <ReconstructionDataFormats/Track.h> |
42 | 35 |
|
43 | 36 | #include <TF1.h> |
| 37 | +#include <TH2.h> |
| 38 | +#include <TH3.h> |
44 | 39 | #include <TMCProcess.h> |
45 | 40 | #include <TPDGCode.h> |
46 | 41 |
|
47 | | -#include <algorithm> |
48 | 42 | #include <array> |
49 | 43 | #include <cmath> |
50 | | -#include <map> |
51 | 44 | #include <memory> |
52 | | -#include <numeric> |
53 | | -#include <set> |
54 | 45 | #include <string> |
55 | 46 | #include <vector> |
56 | 47 |
|
@@ -81,31 +72,31 @@ using TracksMC = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, |
81 | 72 | aod::TrackSelection, aod::McTrackLabels, |
82 | 73 | aod::pidTPCPi, aod::pidTPCKa, aod::pidTPCPr>; |
83 | 74 |
|
84 | | -static constexpr int NCentHists{10}; |
85 | | -static constexpr int NPartHists{5}; |
86 | | -std::array<std::shared_ptr<TH3>, NCentHists> hDedxVsMomentumVsCentPos{}; |
87 | | -std::array<std::shared_ptr<TH3>, NCentHists> hDedxVsMomentumVsCentNeg{}; |
88 | | -std::array<std::shared_ptr<TH3>, NCentHists + 1> hDedxVspTMomentumVsCent{}; |
89 | | -std::array<std::shared_ptr<TH2>, NCentHists + 1> hMomentumVsEtaPos{}; |
90 | | -std::array<std::shared_ptr<TH2>, NCentHists + 1> hMomentumVsEtaNeg{}; |
91 | | -std::array<std::shared_ptr<TH2>, NCentHists + 1> hpTVsEtaPos{}; |
92 | | -std::array<std::shared_ptr<TH2>, NCentHists + 1> hpTVsEtaNeg{}; |
93 | | -// Total counts |
94 | | -std::array<std::shared_ptr<TH2>, NCentHists + 1> hTotalMomPosCent{}; |
95 | | -std::array<std::shared_ptr<TH2>, NCentHists + 1> hTotalMomNegCent{}; |
96 | | -std::array<std::shared_ptr<TH2>, NCentHists + 1> hTotalPtPosCent{}; |
97 | | -std::array<std::shared_ptr<TH2>, NCentHists + 1> hTotalPtNegCent{}; |
98 | | -// Counts for particles |
99 | | -std::array<std::array<std::shared_ptr<TH2>, NCentHists + 1>, NPartHists> hFracMomPosCent{}; |
100 | | -std::array<std::array<std::shared_ptr<TH2>, NCentHists + 1>, NPartHists> hFracMomNegCent{}; |
101 | | -std::array<std::array<std::shared_ptr<TH2>, NCentHists + 1>, NPartHists> hFracPtPosCent{}; |
102 | | -std::array<std::array<std::shared_ptr<TH2>, NCentHists + 1>, NPartHists> hFracPtNegCent{}; |
103 | | - |
104 | 75 | struct MultiplicityPt { |
105 | 76 |
|
106 | 77 | // ── Services ────────────────────────────────────────────── |
107 | | - Service<o2::framework::O2DatabasePDG> pdg; |
108 | | - Service<ccdb::BasicCCDBManager> ccdb; |
| 78 | + Service<o2::framework::O2DatabasePDG> pdg{}; |
| 79 | + Service<ccdb::BasicCCDBManager> ccdb{}; |
| 80 | + |
| 81 | + static constexpr int NCentHists{10}; |
| 82 | + static constexpr int NPartHists{5}; |
| 83 | + std::array<std::shared_ptr<TH3>, NCentHists> hDedxVsMomentumVsCentPos{}; |
| 84 | + std::array<std::shared_ptr<TH3>, NCentHists> hDedxVsMomentumVsCentNeg{}; |
| 85 | + std::array<std::shared_ptr<TH3>, NCentHists + 1> hDedxVspTMomentumVsCent{}; |
| 86 | + std::array<std::shared_ptr<TH2>, NCentHists + 1> hMomentumVsEtaPos{}; |
| 87 | + std::array<std::shared_ptr<TH2>, NCentHists + 1> hMomentumVsEtaNeg{}; |
| 88 | + std::array<std::shared_ptr<TH2>, NCentHists + 1> hpTVsEtaPos{}; |
| 89 | + std::array<std::shared_ptr<TH2>, NCentHists + 1> hpTVsEtaNeg{}; |
| 90 | + // Total counts |
| 91 | + std::array<std::shared_ptr<TH2>, NCentHists + 1> hTotalMomPosCent{}; |
| 92 | + std::array<std::shared_ptr<TH2>, NCentHists + 1> hTotalMomNegCent{}; |
| 93 | + std::array<std::shared_ptr<TH2>, NCentHists + 1> hTotalPtPosCent{}; |
| 94 | + std::array<std::shared_ptr<TH2>, NCentHists + 1> hTotalPtNegCent{}; |
| 95 | + // Counts for particles |
| 96 | + std::array<std::array<std::shared_ptr<TH2>, NCentHists + 1>, NPartHists> hFracMomPosCent{}; |
| 97 | + std::array<std::array<std::shared_ptr<TH2>, NCentHists + 1>, NPartHists> hFracMomNegCent{}; |
| 98 | + std::array<std::array<std::shared_ptr<TH2>, NCentHists + 1>, NPartHists> hFracPtPosCent{}; |
| 99 | + std::array<std::array<std::shared_ptr<TH2>, NCentHists + 1>, NPartHists> hFracPtNegCent{}; |
109 | 100 |
|
110 | 101 | // ── Constant values ────────────────────────────────── |
111 | 102 | static constexpr int CentBinMax = 100; |
|
0 commit comments