|
18 | 18 | /// Also contains checks and monte-carlo (efficiency, purity, mc-true correlation,...) |
19 | 19 | /// End goal of studying correlations between direct photons and jets |
20 | 20 |
|
21 | | -#define ETA_MAX_DEFAULT 0.8 |
22 | | -#define DPHI_SCALE constants::math::TwoPI - constants::math::PIHalf |
23 | | -#define DETA_SCALE 4 * ETA_MAX_DEFAULT - 2 * ETA_MAX_DEFAULT |
24 | | - |
25 | 21 | #include <cmath> |
26 | 22 | #include <deque> |
27 | 23 | #include <string> |
|
51 | 47 | #include "PWGEM/PhotonMeson/Utils/PCMUtilities.h" |
52 | 48 |
|
53 | 49 | #include "PWGJE/DataModel/Jet.h" |
| 50 | +#include "PWGJE/DataModel/PhotonChargedTriggerCorrelation.h" |
54 | 51 |
|
55 | | -// derived data for correlations (on-the-fly) |
56 | | -// should be more efficient due to usage in correlations and event mixing |
57 | | - |
58 | | -namespace o2::aod |
59 | | -{ |
60 | | - |
61 | | -// basic correlation particle columns |
62 | | -namespace corr_particle |
63 | | -{ |
64 | | -DECLARE_SOA_INDEX_COLUMN_FULL(JetCollision, jetCollision, int, JCollisions, ""); |
65 | | -DECLARE_SOA_INDEX_COLUMN_FULL(JetMcCollision, jetMcCollision, int, JMcCollisions, ""); |
66 | | -DECLARE_SOA_COLUMN(Pt, pt, float); |
67 | | -DECLARE_SOA_COLUMN(Phi, phi, float); |
68 | | -DECLARE_SOA_COLUMN(Eta, eta, float); |
69 | | -} // namespace corr_particle |
70 | | - |
71 | | -// reco |
72 | | - |
73 | | -// collision extension |
74 | | -namespace collision_extra_corr |
75 | | -{ |
76 | | -DECLARE_SOA_COLUMN(SelEv, selEv, bool); |
77 | | -DECLARE_SOA_COLUMN(TrigEv, trigEv, bool); |
78 | | -} // namespace collision_extra_corr |
79 | | -DECLARE_SOA_TABLE(CollisionsExtraCorr, "AOD", "COLLISIONSEXTRACORR", |
80 | | - collision_extra_corr::SelEv, collision_extra_corr::TrigEv); |
81 | | - |
82 | | -// trigger |
83 | | -namespace trigger |
84 | | -{ |
85 | | -DECLARE_SOA_INDEX_COLUMN_FULL(JetTrack, jetTrack, int, JetTracks, ""); |
86 | | -} // namespace trigger |
87 | | -DECLARE_SOA_TABLE(Triggers, "AOD", "TRIGGERS", |
88 | | - o2::soa::Index<>, corr_particle::JetCollisionId, trigger::JetTrackId, |
89 | | - corr_particle::Pt, corr_particle::Phi, corr_particle::Eta); |
90 | | -using Trigger = Triggers::iterator; |
91 | | - |
92 | | -// hadrons (global tracks) |
93 | | -namespace hadron |
94 | | -{ |
95 | | -DECLARE_SOA_INDEX_COLUMN_FULL(JetTrack, jetTrack, int, JetTracks, ""); |
96 | | -} // namespace hadron |
97 | | -DECLARE_SOA_TABLE(Hadrons, "AOD", "HADRONS", |
98 | | - o2::soa::Index<>, corr_particle::JetCollisionId, hadron::JetTrackId, |
99 | | - corr_particle::Pt, corr_particle::Phi, corr_particle::Eta); |
100 | | -using Hadron = Hadrons::iterator; |
101 | | - |
102 | | -// pipm |
103 | | -namespace pipm |
104 | | -{ |
105 | | -DECLARE_SOA_INDEX_COLUMN_FULL(JetTrack, jetTrack, int, JetTracks, ""); |
106 | | -} // namespace pipm |
107 | | -DECLARE_SOA_TABLE(Pipms, "AOD", "PIPMS", |
108 | | - o2::soa::Index<>, corr_particle::JetCollisionId, pipm::JetTrackId, |
109 | | - corr_particle::Pt, corr_particle::Phi, corr_particle::Eta); |
110 | | -using Pipm = Pipms::iterator; |
111 | | - |
112 | | -// photonPCM |
113 | | -namespace photon_pcm |
114 | | -{ |
115 | | -DECLARE_SOA_INDEX_COLUMN_FULL(V0PhotonKF, v0PhotonKF, int, V0PhotonsKF, ""); |
116 | | -DECLARE_SOA_COLUMN(PosTrackId, posTrackId, int); |
117 | | -DECLARE_SOA_COLUMN(NegTrackId, negTrackId, int); |
118 | | -} // namespace photon_pcm |
119 | | -DECLARE_SOA_TABLE(PhotonPCMs, "AOD", "PHOTONPCMS", |
120 | | - o2::soa::Index<>, corr_particle::JetCollisionId, photon_pcm::V0PhotonKFId, |
121 | | - photon_pcm::PosTrackId, photon_pcm::NegTrackId, |
122 | | - corr_particle::Pt, corr_particle::Phi, corr_particle::Eta); |
123 | | -using PhotonPCM = PhotonPCMs::iterator; |
124 | | - |
125 | | -// photonPCM pairs (pi0) |
126 | | -namespace photon_pcm_pair |
127 | | -{ |
128 | | -DECLARE_SOA_INDEX_COLUMN_FULL(V0PhotonKF1, v0PhotonKF1, int, V0PhotonsKF, "_1"); |
129 | | -DECLARE_SOA_INDEX_COLUMN_FULL(V0PhotonKF2, v0PhotonKF2, int, V0PhotonsKF, "_2"); |
130 | | -DECLARE_SOA_COLUMN(PosTrack1Id, posTrack1Id, int); |
131 | | -DECLARE_SOA_COLUMN(NegTrack1Id, negTrack1Id, int); |
132 | | -DECLARE_SOA_COLUMN(PosTrack2Id, posTrack2Id, int); |
133 | | -DECLARE_SOA_COLUMN(NegTrack2Id, negTrack2Id, int); |
134 | | -DECLARE_SOA_COLUMN(Mgg, mgg, float); |
135 | | -} // namespace photon_pcm_pair |
136 | | -DECLARE_SOA_TABLE(PhotonPCMPairs, "AOD", "PHOTONPCMPAIRS", |
137 | | - o2::soa::Index<>, corr_particle::JetCollisionId, photon_pcm_pair::V0PhotonKF1Id, photon_pcm_pair::V0PhotonKF2Id, |
138 | | - photon_pcm_pair::PosTrack1Id, photon_pcm_pair::NegTrack1Id, photon_pcm_pair::PosTrack2Id, photon_pcm_pair::NegTrack2Id, |
139 | | - corr_particle::Pt, corr_particle::Phi, corr_particle::Eta, photon_pcm_pair::Mgg); |
140 | | -using PhotonPCMPair = PhotonPCMPairs::iterator; |
141 | | - |
142 | | -// mc |
143 | | - |
144 | | -// mcCollision extension |
145 | | -namespace mc_collision_extra_corr |
146 | | -{ |
147 | | -DECLARE_SOA_COLUMN(TrigEv, trigEv, bool); |
148 | | -} // namespace mc_collision_extra_corr |
149 | | -DECLARE_SOA_TABLE(McCollisionsExtraCorr, "AOD", "MCCOLLISIONSEXTRACORR", |
150 | | - mc_collision_extra_corr::TrigEv); |
151 | | - |
152 | | -// trigger |
153 | | -namespace trigger_particle |
154 | | -{ |
155 | | -DECLARE_SOA_INDEX_COLUMN_FULL(JetMcParticle, jetMcParticle, int, JetParticles, ""); |
156 | | -} // namespace trigger_particle |
157 | | -DECLARE_SOA_TABLE(TriggerParticles, "AOD", "TRIGGERPARTICLES", |
158 | | - o2::soa::Index<>, corr_particle::JetMcCollisionId, trigger_particle::JetMcParticleId, |
159 | | - corr_particle::Pt, corr_particle::Phi, corr_particle::Eta); |
160 | | -using TriggerParticle = TriggerParticles::iterator; |
161 | | -} // namespace o2::aod |
| 52 | +const double absEtaMax = 0.8; |
| 53 | +#define DPHI_SCALE constants::math::TwoPI - constants::math::PIHalf |
| 54 | +#define DETA_SCALE 4 * absEtaMax - 2 * absEtaMax |
162 | 55 |
|
163 | 56 | using namespace o2; |
164 | 57 | using namespace o2::framework; |
@@ -189,7 +82,7 @@ struct CorrelationTableProducer { |
189 | 82 | Configurable<double> zPvMax{"zPvMax", 7, "maximum absZ primary-vertex cut"}; |
190 | 83 | Configurable<int> occupancyMin{"occupancyMin", 0, "minimum occupancy cut"}; |
191 | 84 | Configurable<int> occupancyMax{"occupancyMax", 2000, "maximum occupancy cut"}; |
192 | | - Configurable<double> etaMax{"etaMax", ETA_MAX_DEFAULT, "maximum absEta cut"}; |
| 85 | + Configurable<double> etaMax{"etaMax", 1 * absEtaMax, "maximum absEta cut"}; |
193 | 86 |
|
194 | 87 | Configurable<std::string> eventSelections{"eventSelections", "sel8", "JE framework - event selection"}; |
195 | 88 | Configurable<std::string> trackSelections{"trackSelections", "globalTracks", "JE framework - track selections"}; |
|
0 commit comments