use std::hypot for vector norm calculation

Author: lubynets

PWGHF/TableProducer/treeCreatorLcToPKPi.cxx

@@ -875,11 +875,11 @@ struct HfTreeCreatorLcToPKPi {
     const float chi2Topo = candidate.kfChi2Topo();
     const float decayLength = candidate.kfDecayLength();
     const float dl = candidate.kfDecayLengthError();
-    const float pt = std::sqrt(candidate.kfPx() * candidate.kfPx() + candidate.kfPy() * candidate.kfPy());
+    const float pt = std::hypot(candidate.kfPx(), candidate.kfPy());
     const float deltaPt = std::sqrt(candidate.kfPx() * candidate.kfPx() * candidate.kfErrorPx() * candidate.kfErrorPx() +
                                     candidate.kfPy() * candidate.kfPy() * candidate.kfErrorPy() * candidate.kfErrorPy()) /
                           pt;
-    const float p = std::sqrt(pt * pt + candidate.kfPz() * candidate.kfPz());
+    const float p = std::hypot(pt, candidate.kfPz());
     const float deltaP = std::sqrt(pt * pt * deltaPt * deltaPt +
                                    candidate.kfPz() * candidate.kfPz() * candidate.kfErrorPz() * candidate.kfErrorPz()) /
                          p;

Tools/KFparticle/KFUtilities.h

@@ -298,7 +298,7 @@ float ldlFromKF(KFParticle kfpParticle, KFParticle PV)
   const float dxParticle = PV.GetX() - kfpParticle.GetX();
   const float dyParticle = PV.GetY() - kfpParticle.GetY();
   const float dzParticle = PV.GetZ() - kfpParticle.GetZ();
-  float lParticle = std::sqrt(dxParticle * dxParticle + dyParticle * dyParticle + dzParticle * dzParticle);
+  float lParticle = std::hypot(dxParticle, dyParticle, dzParticle);
   float dlParticle = (PV.GetCovariance(0) + kfpParticle.GetCovariance(0)) * dxParticle * dxParticle + (PV.GetCovariance(2) + kfpParticle.GetCovariance(2)) * dyParticle * dyParticle + (PV.GetCovariance(5) + kfpParticle.GetCovariance(5)) * dzParticle * dzParticle + 2 * ((PV.GetCovariance(1) + kfpParticle.GetCovariance(1)) * dxParticle * dyParticle + (PV.GetCovariance(3) + kfpParticle.GetCovariance(3)) * dxParticle * dzParticle + (PV.GetCovariance(4) + kfpParticle.GetCovariance(4)) * dyParticle * dzParticle);
   if (std::fabs(lParticle) < ArbitrarySmallNumber)
     lParticle = ArbitrarySmallNumber;
@@ -316,7 +316,7 @@ float ldlXYFromKF(KFParticle kfpParticle, KFParticle PV)
 {
   const float dxParticle = PV.GetX() - kfpParticle.GetX();
   const float dyParticle = PV.GetY() - kfpParticle.GetY();
-  float lParticle = std::sqrt(dxParticle * dxParticle + dyParticle * dyParticle);
+  float lParticle = std::hypot(dxParticle, dyParticle);
   float dlParticle = (PV.GetCovariance(0) + kfpParticle.GetCovariance(0)) * dxParticle * dxParticle + (PV.GetCovariance(2) + kfpParticle.GetCovariance(2)) * dyParticle * dyParticle + 2 * ((PV.GetCovariance(1) + kfpParticle.GetCovariance(1)) * dxParticle * dyParticle);
   if (std::fabs(lParticle) < ArbitrarySmallNumber)
     lParticle = ArbitrarySmallNumber;
@@ -366,7 +366,7 @@ float kfCalculateDistanceBetweenParticles(KFParticle track1, KFParticle track2)
   const float dx = params1[0] - params2[0];
   const float dy = params1[1] - params2[1];
   const float dz = params1[2] - params2[2];
-  return std::sqrt(dx * dx + dy * dy + dz * dz);
+  return std::hypot(dx, dy, dz);
 }
 
 /// @brief squared distance between two tracks normalised by its uncertainty evaluated in matrix form