Test LorentzVector solution for MegaLinter

Author: ariffero

PWGUD/Tasks/FwdMuonsUPC.cxx

@@ -28,7 +28,7 @@
 #include "DataFormatsParameters/GRPECSObject.h"
 #include "PWGUD/DataModel/UDTables.h"
 
-#include "TLorentzVector.h"
+#include "LorentzVector.h"
 #include "TSystem.h"
 #include "TMath.h"
 #include "TRandom3.h"
@@ -461,7 +461,7 @@ struct FwdMuonsUPC {
   {
     float rAbs = fwdTrack.rAtAbsorberEnd();
     float pDca = fwdTrack.pDca();
-    TLorentzVector p;
+    LorentzVector p;
     auto mMu = particleMass(13);
     p.SetXYZM(fwdTrack.px(), fwdTrack.py(), fwdTrack.pz(), mMu);
     float eta = p.Eta();
@@ -480,10 +480,10 @@ struct FwdMuonsUPC {
   }
 
   // function to compute phi for azimuth anisotropy
-  void computePhiAnis(TLorentzVector p1, TLorentzVector p2, int sign1, float& phiAverage, float& phiCharge)
+  void computePhiAnis(LorentzVector p1, LorentzVector p2, int sign1, float& phiAverage, float& phiCharge)
   {
 
-    TLorentzVector tSum, tDiffAv, tDiffCh;
+    LorentzVector tSum, tDiffAv, tDiffCh;
     tSum = p1 + p2;
     if (sign1 > 0) {
       tDiffCh = p1 - p2;
@@ -560,11 +560,11 @@ struct FwdMuonsUPC {
       return;
 
     // form Lorentz vectors
-    TLorentzVector p1, p2;
+    LorentzVector p1, p2;
     auto mMu = particleMass(13);
     p1.SetXYZM(tr1.px(), tr1.py(), tr1.pz(), mMu);
     p2.SetXYZM(tr2.px(), tr2.py(), tr2.pz(), mMu);
-    TLorentzVector p = p1 + p2;
+    LorentzVector p = p1 + p2;
 
     // cut on pair kinematics
     // select mass
@@ -688,11 +688,11 @@ struct FwdMuonsUPC {
       LOGF(debug, "PDG codes: %d | %d", McPart1.pdgCode(), McPart2.pdgCode());
 
     // create Lorentz vectors
-    TLorentzVector p1, p2;
+    LorentzVector p1, p2;
     auto mMu = particleMass(13);
     p1.SetXYZM(McPart1.px(), McPart1.py(), McPart1.pz(), mMu);
     p2.SetXYZM(McPart2.px(), McPart2.py(), McPart2.pz(), mMu);
-    TLorentzVector p = p1 + p2;
+    LorentzVector p = p1 + p2;
 
     // cut on pair kinematics
     // select mass
@@ -803,11 +803,11 @@ struct FwdMuonsUPC {
       return;
 
     // form Lorentz vectors
-    TLorentzVector p1, p2;
+    LorentzVector p1, p2;
     auto mMu = particleMass(13);
     p1.SetXYZM(tr1.px(), tr1.py(), tr1.pz(), mMu);
     p2.SetXYZM(tr2.px(), tr2.py(), tr2.pz(), mMu);
-    TLorentzVector p = p1 + p2;
+    LorentzVector p = p1 + p2;
 
     // cut on pair kinematics (reco candidates)
     // select mass
@@ -832,10 +832,10 @@ struct FwdMuonsUPC {
     computePhiAnis(p1, p2, tr1.sign(), phiAverage, phiCharge);
 
     // gen particle
-    TLorentzVector p1Mc, p2Mc;
+    LorentzVector p1Mc, p2Mc;
     p1Mc.SetXYZM(McPart1.px(), McPart1.py(), McPart1.pz(), mMu);
     p2Mc.SetXYZM(McPart2.px(), McPart2.py(), McPart2.pz(), mMu);
-    TLorentzVector pMc = p1Mc + p2Mc;
+    LorentzVector pMc = p1Mc + p2Mc;
 
     // compute gen phi for azimuth anisotropy
     float phiGenAverage = 0;