@@ -115,8 +115,8 @@ static const std::vector<std::string> betheBlochParNames{
115115 " p0" , " p1" , " p2" , " p3" , " p4" , " resolution" };
116116constexpr float betheBlochDefault[nParticles][nBetheParams]{
117117 {0.248753 , 3.58634 , 0.0167065 , 2.29194 , 0.774344 ,
118- 0.07 }, // triton
119- {0.0274556 ,18.3054 , 3.99987e-05 , 3.17219 , 11.1775 ,
118+ 0.07 }, // triton
119+ {0.0274556 , 18.3054 , 3.99987e-05 , 3.17219 , 11.1775 ,
120120 0.07 }}; // Helion
121121} // namespace
122122using namespace o2 ;
@@ -647,13 +647,14 @@ struct TrHeAnalysis {
647647 template <class T >
648648 float betheBlochAleph (Particle const & particle, T const & rigidity)
649649 {
650- double bg = particle.charge * rigidity / particle.mass ;
651- double beta = bg / std::Sqrt (1 . + bg* bg);
650+ double bg = particle.charge * rigidity / particle.mass ;
651+ double beta = bg / std::Sqrt (1 . + bg * bg);
652652 double aa = std::Power (beta, particle.betheParams [3 ]);
653653 double bb = std::Power (1 . / bg, particle.betheParams [4 ]);
654- if ((particle.betheParams [2 ] + bb) <= 0 ) return 0 ;
654+ if ((particle.betheParams [2 ] + bb) <= 0 )
655+ return 0 ;
655656 bb = std::Log (particle.betheParams [2 ] + bb);
656- return std::pow (particle.charge , particle.chargeFactor ) * 50 * (particle.betheParams [1 ] - aa - bb)* particle.betheParams [0 ] / aa;
657+ return std::pow (particle.charge , particle.chargeFactor ) * 50 * (particle.betheParams [1 ] - aa - bb) * particle.betheParams [0 ] / aa;
657658 }
658659
659660 template <class T >
@@ -684,7 +685,7 @@ struct TrHeAnalysis {
684685 {
685686 const float beta = track.beta ();
686687 const float rigidity = getRigidity (track);
687- float gamma = 1 / std::Sqrt (1 - beta* beta);
688+ float gamma = 1 / std::Sqrt (1 - beta * beta);
688689 float mass = (rigidity / std::Sqrt (gamma * gamma - 1 ));
689690 return mass;
690691 }
0 commit comments