|
49 | 49 | #include "ReconstructionDataFormats/Track.h" |
50 | 50 |
|
51 | 51 | #include "TGrid.h" |
| 52 | +#include <Math/Vector4D.h> |
52 | 53 | #include <TList.h> |
53 | 54 | #include <TPDGCode.h> |
54 | 55 | #include <TRandom3.h> |
55 | 56 | #include <TVector2.h> |
56 | 57 | #include <TVector3.h> |
57 | 58 |
|
| 59 | +#include <ROOT/Math/Boost.h> |
| 60 | +#include <ROOT/Math/PxPyPzMVector.h> |
58 | 61 | #include <fastjet/AreaDefinition.hh> |
59 | 62 | #include <fastjet/ClusterSequence.hh> |
60 | 63 | #include <fastjet/ClusterSequenceArea.hh> |
|
64 | 67 | #include <fastjet/tools/JetMedianBackgroundEstimator.hh> |
65 | 68 | #include <fastjet/tools/Subtractor.hh> |
66 | 69 |
|
67 | | -#include <Math/Vector4D.h> |
68 | | -#include <ROOT/Math/PxPyPzMVector.h> |
69 | | -#include <ROOT/Math/Boost.h> |
70 | | - |
71 | 70 | #include <chrono> |
72 | 71 | #include <cmath> |
73 | 72 | #include <memory> |
@@ -99,28 +98,31 @@ struct ReducedParticle { |
99 | 98 | double px; |
100 | 99 | double py; |
101 | 100 | double pz; |
102 | | - int pdgCode; |
103 | | - int mcIndex; |
104 | | - bool used; |
| 101 | + int pdgCode; |
| 102 | + int mcIndex; |
| 103 | + bool used; |
105 | 104 |
|
106 | 105 | // Pseudorapidity |
107 | | - double eta() const { |
| 106 | + double eta() const |
| 107 | + { |
108 | 108 | double p = std::sqrt(px * px + py * py + pz * pz); |
109 | 109 | if (p == std::abs(pz)) { |
110 | 110 | return (pz >= 0) ? 1e10 : -1e10; |
111 | 111 | } |
112 | 112 | return 0.5 * std::log((p + pz) / (p - pz)); |
113 | 113 | } |
114 | | - |
| 114 | + |
115 | 115 | // Azimuthal Angle |
116 | | - double phi() const { |
| 116 | + double phi() const |
| 117 | + { |
117 | 118 | double angle = PI + std::atan2(-py, -px); |
118 | 119 | return angle; |
119 | 120 | } |
120 | 121 |
|
121 | 122 | // Transverse Momentum |
122 | | - double pt() const { |
123 | | - return std::sqrt(px*px + py*py); |
| 123 | + double pt() const |
| 124 | + { |
| 125 | + return std::sqrt(px * px + py * py); |
124 | 126 | } |
125 | 127 | }; |
126 | 128 |
|
@@ -3170,17 +3172,17 @@ struct AntinucleiInJets { |
3170 | 3172 | registryMC.fill(HIST("genEventsCoalescence"), 2.5); |
3171 | 3173 |
|
3172 | 3174 | // Build deuterons |
3173 | | - for (int ip=0 ; ip<protons.size() ; ip++) { |
| 3175 | + for (int ip = 0; ip < protons.size(); ip++) { |
3174 | 3176 | auto& proton = protons[ip]; |
3175 | 3177 | if (proton.used) |
3176 | 3178 | continue; |
3177 | 3179 |
|
3178 | | - for (int in=0 ; in<neutrons.size() ; in++) { |
| 3180 | + for (int in = 0; in < neutrons.size(); in++) { |
3179 | 3181 | auto& neutron = neutrons[in]; |
3180 | 3182 | if (neutron.used) |
3181 | 3183 | continue; |
3182 | 3184 |
|
3183 | | - if (passDeuteronCoalescence(proton,neutron,coalescenceMomentum,mRand)) { |
| 3185 | + if (passDeuteronCoalescence(proton, neutron, coalescenceMomentum, mRand)) { |
3184 | 3186 |
|
3185 | 3187 | int sign = (proton.pdgCode > 0) ? +1 : -1; |
3186 | 3188 | int deuteronPdg = sign * o2::constants::physics::Pdg::kDeuteron; |
|
0 commit comments