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/PxPyPzMVector.h>
60+ #include < ROOT/Math/Boost.h>
61+
5862#include < fastjet/AreaDefinition.hh>
5963#include < fastjet/ClusterSequence.hh>
6064#include < fastjet/ClusterSequenceArea.hh>
6468#include < fastjet/tools/JetMedianBackgroundEstimator.hh>
6569#include < fastjet/tools/Subtractor.hh>
6670
67- #include < Math/Vector4D.h>
68- #include < ROOT/Math/PxPyPzMVector.h>
69- #include < ROOT/Math/Boost.h>
70-
7171#include < chrono>
7272#include < cmath>
7373#include < memory>
@@ -99,12 +99,13 @@ struct ReducedParticle {
9999 double px;
100100 double py;
101101 double pz;
102- int pdgCode;
103- int mcIndex;
104- bool used;
102+ int pdgCode;
103+ int mcIndex;
104+ bool used;
105105
106106 // Pseudorapidity
107- double eta () const {
107+ double eta () const
108+ {
108109 double p = std::sqrt (px * px + py * py + pz * pz);
109110 if (p == std::abs (pz)) {
110111 return (pz >= 0 ) ? 1e10 : -1e10 ;
@@ -113,13 +114,15 @@ struct ReducedParticle {
113114 }
114115
115116 // Azimuthal Angle
116- double phi () const {
117+ double phi () const
118+ {
117119 double angle = PI + std::atan2 (-py, -px);
118120 return angle;
119121 }
120122
121123 // Transverse Momentum
122- double pt () const {
124+ double pt () const
125+ {
123126 return std::sqrt (px*px + py*py);
124127 }
125128};
@@ -3170,17 +3173,17 @@ struct AntinucleiInJets {
31703173 registryMC.fill (HIST (" genEventsCoalescence" ), 2.5 );
31713174
31723175 // Build deuterons
3173- for (int ip= 0 ; ip< protons.size () ; ip++) {
3176+ for (int ip = 0 ; ip < protons.size () ; ip++) {
31743177 auto & proton = protons[ip];
31753178 if (proton.used )
31763179 continue ;
31773180
3178- for (int in= 0 ; in< neutrons.size () ; in++) {
3181+ for (int in = 0 ; in < neutrons.size () ; in++) {
31793182 auto & neutron = neutrons[in];
31803183 if (neutron.used )
31813184 continue ;
31823185
3183- if (passDeuteronCoalescence (proton,neutron,coalescenceMomentum,mRand )) {
3186+ if (passDeuteronCoalescence (proton, neutron, coalescenceMomentum, mRand )) {
31843187
31853188 int sign = (proton.pdgCode > 0 ) ? +1 : -1 ;
31863189 int deuteronPdg = sign * o2::constants::physics::Pdg::kDeuteron ;
0 commit comments