diff --git a/PWGLF/TableProducer/Nuspex/trHeAnalysis.cxx b/PWGLF/TableProducer/Nuspex/trHeAnalysis.cxx index d2b74f2af3d..0b1664d2b7f 100644 --- a/PWGLF/TableProducer/Nuspex/trHeAnalysis.cxx +++ b/PWGLF/TableProducer/Nuspex/trHeAnalysis.cxx @@ -711,36 +711,36 @@ struct TrHeAnalysis { return hePID ? track.tpcInnerParam() / 2 : track.tpcInnerParam(); } template -float getMassSquared(const T& track) -{ - if (cfgMassMethod == 0) { - float m = track.mass(); - return m * m; - } - if (cfgMassMethod == 1) { - const float beta = track.beta(); - const float rigidity = getRigidity(track); - float gamma = 1.f / std::sqrt(1.f - beta * beta); - float mass = rigidity / std::sqrt(gamma * gamma - 1.f); - return mass * mass; - } - if (cfgMassMethod == 2) { - const float rigidity = getRigidity(track); - float tofStartTime = track.evTimeForTrack(); - float tofTime = track.tofSignal(); - constexpr float CInCmPs = 2.99792458e-2f; - float length = track.length(); - float time = tofTime - tofStartTime; - if (time > 0.f && length > 0.f) { - float beta = length / (CInCmPs * time); + float getMassSquared(const T& track) + { + if (cfgMassMethod == 0) { + float m = track.mass(); + return m * m; + } + if (cfgMassMethod == 1) { + const float beta = track.beta(); + const float rigidity = getRigidity(track); float gamma = 1.f / std::sqrt(1.f - beta * beta); float mass = rigidity / std::sqrt(gamma * gamma - 1.f); return mass * mass; } + if (cfgMassMethod == 2) { + const float rigidity = getRigidity(track); + float tofStartTime = track.evTimeForTrack(); + float tofTime = track.tofSignal(); + constexpr float CInCmPs = 2.99792458e-2f; + float length = track.length(); + float time = tofTime - tofStartTime; + if (time > 0.f && length > 0.f) { + float beta = length / (CInCmPs * time); + float gamma = 1.f / std::sqrt(1.f - beta * beta); + float mass = rigidity / std::sqrt(gamma * gamma - 1.f); + return mass * mass; + } + return -1.f; + } return -1.f; } - return -1.f; -} }; WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)