@@ -711,36 +711,36 @@ struct TrHeAnalysis {
711711 return hePID ? track.tpcInnerParam () / 2 : track.tpcInnerParam ();
712712 }
713713 template <class T >
714- float getMassSquared (const T& track)
715- {
716- if (cfgMassMethod == 0 ) {
717- float m = track.mass ();
718- return m * m;
719- }
720- if (cfgMassMethod == 1 ) {
721- const float beta = track.beta ();
722- const float rigidity = getRigidity (track);
723- float gamma = 1 .f / std::sqrt (1 .f - beta * beta);
724- float mass = rigidity / std::sqrt (gamma * gamma - 1 .f );
725- return mass * mass;
726- }
727- if (cfgMassMethod == 2 ) {
728- const float rigidity = getRigidity (track);
729- float tofStartTime = track.evTimeForTrack ();
730- float tofTime = track.tofSignal ();
731- constexpr float CInCmPs = 2 .99792458e-2f ;
732- float length = track.length ();
733- float time = tofTime - tofStartTime;
734- if (time > 0 .f && length > 0 .f ) {
735- float beta = length / (CInCmPs * time);
714+ float getMassSquared (const T& track)
715+ {
716+ if (cfgMassMethod == 0 ) {
717+ float m = track.mass ();
718+ return m * m;
719+ }
720+ if (cfgMassMethod == 1 ) {
721+ const float beta = track.beta ();
722+ const float rigidity = getRigidity (track);
736723 float gamma = 1 .f / std::sqrt (1 .f - beta * beta);
737724 float mass = rigidity / std::sqrt (gamma * gamma - 1 .f );
738725 return mass * mass;
739726 }
727+ if (cfgMassMethod == 2 ) {
728+ const float rigidity = getRigidity (track);
729+ float tofStartTime = track.evTimeForTrack ();
730+ float tofTime = track.tofSignal ();
731+ constexpr float CInCmPs = 2 .99792458e-2f ;
732+ float length = track.length ();
733+ float time = tofTime - tofStartTime;
734+ if (time > 0 .f && length > 0 .f ) {
735+ float beta = length / (CInCmPs * time);
736+ float gamma = 1 .f / std::sqrt (1 .f - beta * beta);
737+ float mass = rigidity / std::sqrt (gamma * gamma - 1 .f );
738+ return mass * mass;
739+ }
740+ return -1 .f ;
741+ }
740742 return -1 .f ;
741743 }
742- return -1 .f ;
743- }
744744};
745745
746746WorkflowSpec defineDataProcessing (ConfigContext const & cfgc)
0 commit comments