4949#include " ReconstructionDataFormats/Track.h"
5050
5151#include " TGrid.h"
52+ #include < Math/Vector4D.h>
5253#include < TList.h>
5354#include < TPDGCode.h>
5455#include < TRandom3.h>
5556#include < TVector2.h>
5657#include < TVector3.h>
5758
59+ #include < ROOT/Math/Boost.h>
60+ #include < ROOT/Math/PxPyPzMVector.h>
5861#include < fastjet/AreaDefinition.hh>
5962#include < fastjet/ClusterSequence.hh>
6063#include < fastjet/ClusterSequenceArea.hh>
6467#include < fastjet/tools/JetMedianBackgroundEstimator.hh>
6568#include < fastjet/tools/Subtractor.hh>
6669
67- #include < Math/Vector4D.h>
68- #include < ROOT/Math/PxPyPzMVector.h>
69- #include < ROOT/Math/Boost.h>
70-
7170#include < chrono>
7271#include < cmath>
7372#include < memory>
@@ -99,12 +98,13 @@ struct ReducedParticle {
9998 double px;
10099 double py;
101100 double pz;
102- int pdgCode;
103- int mcIndex;
104- bool used;
101+ int pdgCode;
102+ int mcIndex;
103+ bool used;
105104
106105 // Pseudorapidity
107- double eta () const {
106+ double eta () const
107+ {
108108 double p = std::sqrt (px * px + py * py + pz * pz);
109109 if (p == std::abs (pz)) {
110110 return (pz >= 0 ) ? 1e10 : -1e10 ;
@@ -113,14 +113,16 @@ struct ReducedParticle {
113113 }
114114
115115 // Azimuthal Angle
116- double phi () const {
116+ double phi () const
117+ {
117118 double angle = PI + std::atan2 (-py, -px);
118119 return angle;
119120 }
120121
121122 // 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);
124126 }
125127};
126128
@@ -3170,17 +3172,17 @@ struct AntinucleiInJets {
31703172 registryMC.fill (HIST (" genEventsCoalescence" ), 2.5 );
31713173
31723174 // Build deuterons
3173- for (int ip= 0 ; ip< protons.size () ; ip++) {
3175+ for (int ip = 0 ; ip < protons.size (); ip++) {
31743176 auto & proton = protons[ip];
31753177 if (proton.used )
31763178 continue ;
31773179
3178- for (int in= 0 ; in< neutrons.size () ; in++) {
3180+ for (int in = 0 ; in < neutrons.size (); in++) {
31793181 auto & neutron = neutrons[in];
31803182 if (neutron.used )
31813183 continue ;
31823184
3183- if (passDeuteronCoalescence (proton,neutron,coalescenceMomentum,mRand )) {
3185+ if (passDeuteronCoalescence (proton, neutron, coalescenceMomentum, mRand )) {
31843186
31853187 int sign = (proton.pdgCode > 0 ) ? +1 : -1 ;
31863188 int deuteronPdg = sign * o2::constants::physics::Pdg::kDeuteron ;
0 commit comments