@@ -3304,35 +3304,36 @@ void VarManager::FillTrackCollisionMC(T1 const& track, T2 const& MotherTrack, C
33043304}
33053305
33063306template <int candidateType, typename T1>
3307- void VarManager::FillTrackCollisionMC (T1 const & track, const std::array<double , 3 >& collPos, float massHyp, float * values){
3307+ void VarManager::FillTrackCollisionMC (T1 const & track, const std::array<double , 3 >& collPos, float massHyp, float * values)
3308+ {
33083309
3309- if (!values) {
3310+ if (!values) {
33103311 values = fgValues;
33113312 }
33123313
33133314 float m = o2::constants::physics::MassJPsi;
3314- if (massHyp) m = massHyp;
3315+ if (massHyp)
3316+ m = massHyp;
33153317
33163318 // displaced vertex is computed with decay product (track) and vertex collPos
3317- values[kMCVertexingLxy ] = (collPos[0 ]- track.vx ()) * (collPos[0 ]- track.vx ()) +
3318- (collPos[1 ]- track.vy ()) * (collPos[1 ]- track.vy ());
3319- values[kMCVertexingLz ] = (collPos[2 ]- track.vz ()) * (collPos[2 ]- track.vz ());
3319+ values[kMCVertexingLxy ] = (collPos[0 ] - track.vx ()) * (collPos[0 ] - track.vx ()) +
3320+ (collPos[1 ] - track.vy ()) * (collPos[1 ] - track.vy ());
3321+ values[kMCVertexingLz ] = (collPos[2 ] - track.vz ()) * (collPos[2 ] - track.vz ());
33203322 values[kMCVertexingLxyz ] = values[kMCVertexingLxy ] + values[kMCVertexingLz ];
33213323 values[kMCVertexingLxy ] = std::sqrt (values[kMCVertexingLxy ]);
33223324 values[kMCVertexingLz ] = std::sqrt (values[kMCVertexingLz ]);
33233325 values[kMCVertexingLxyz ] = std::sqrt (values[kMCVertexingLxyz ]);
3324- values[kMCVertexingTauz ] = (collPos[2 ]- track.vz ()) * m / (TMath::Abs (track.pz ()) * o2::constants::physics::LightSpeedCm2NS);
3326+ values[kMCVertexingTauz ] = (collPos[2 ] - track.vz ()) * m / (TMath::Abs (track.pz ()) * o2::constants::physics::LightSpeedCm2NS);
33253327 values[kMCVertexingTauxy ] = values[kMCVertexingLxy ] * m / (track.pt () * o2::constants::physics::LightSpeedCm2NS);
33263328
3327- values[kMCVertexingLzProjected ] = ((track.vz ()- collPos[2 ]) * track.pz ()) / TMath::Abs (track.pz ());
3328- values[kMCVertexingLxyProjected ] = (((track.vx ()- collPos[0 ]) * track.px ()) + ((track.vy ()- collPos[1 ]) * track.py ())) / TMath::Abs (track.pt ());
3329- values[kMCVertexingLxyzProjected ] = (((track.vx ()- collPos[0 ]) * track.px ()) + ((track.vy ()- collPos[1 ]) * track.py ()) + ((track.vz ()- collPos[2 ]) * track.pz ())) / track.p ();
3329+ values[kMCVertexingLzProjected ] = ((track.vz () - collPos[2 ]) * track.pz ()) / TMath::Abs (track.pz ());
3330+ values[kMCVertexingLxyProjected ] = (((track.vx () - collPos[0 ]) * track.px ()) + ((track.vy () - collPos[1 ]) * track.py ())) / TMath::Abs (track.pt ());
3331+ values[kMCVertexingLxyzProjected ] = (((track.vx () - collPos[0 ]) * track.px ()) + ((track.vy () - collPos[1 ]) * track.py ()) + ((track.vz () - collPos[2 ]) * track.pz ())) / track.p ();
33303332 values[kMCVertexingTauxyProjected ] = values[kMCVertexingLxyProjected ] * m / (track.pt ());
33313333 values[kMCVertexingTauzProjected ] = values[kMCVertexingLzProjected ] * m / TMath::Abs (track.pz ());
33323334 values[kMCVertexingTauxyzProjected ] = values[kMCVertexingLxyzProjected ] * m / (track.p ());
33333335}
33343336
3335-
33363337template <int pairType, typename T, typename T1>
33373338void VarManager::FillEnergyCorrelatorsMC (T const & track, T1 const & t1, float * values, float Translow, float Transhigh)
33383339{
0 commit comments