diff --git a/PWGLF/TableProducer/Nuspex/trHeAnalysis.cxx b/PWGLF/TableProducer/Nuspex/trHeAnalysis.cxx index c6cca3bb498..a33295b6021 100644 --- a/PWGLF/TableProducer/Nuspex/trHeAnalysis.cxx +++ b/PWGLF/TableProducer/Nuspex/trHeAnalysis.cxx @@ -647,13 +647,14 @@ struct TrHeAnalysis { template float betheBlochAleph(Particle const& particle, T const& rigidity) { - double bg = particle.charge*rigidity / particle.mass; - double beta = bg / std::sqrt(1. + bg*bg); + double bg = particle.charge * rigidity / particle.mass; + double beta = bg / std::sqrt(1. + bg * bg); double aa = std::pow(beta, particle.betheParams[3]); double bb = std::pow(1. / bg, particle.betheParams[4]); - if ((particle.betheParams[2] + bb) <= 0) return 0; + if ((particle.betheParams[2] + bb) <= 0) + return 0; bb = std::log(particle.betheParams[2] + bb); - return std::pow(particle.charge, particle.chargeFactor) * 50 * (particle.betheParams[1] - aa - bb)*particle.betheParams[0] / aa; + return std::pow(particle.charge, particle.chargeFactor) * 50 * (particle.betheParams[1] - aa - bb) * particle.betheParams[0] / aa; } template @@ -684,7 +685,7 @@ struct TrHeAnalysis { { const float beta = track.beta(); const float rigidity = getRigidity(track); - float gamma = 1 / std::sqrt(1-beta*beta); + float gamma = 1 / std::sqrt(1 - beta * beta); float mass = (rigidity / std::sqrt(gamma * gamma - 1)); return mass; }