Some calculations simplified.
authorvkucera <vit.kucera@cern.ch>
Wed, 29 Oct 2014 22:49:33 +0000 (23:49 +0100)
committermverweij <marta.verweij@cern.ch>
Mon, 3 Nov 2014 14:10:20 +0000 (15:10 +0100)
PWGJE/EMCALJetTasks/AliEmcalJet.cxx
PWGJE/EMCALJetTasks/AliEmcalJet.h

index b84160b..e68a91e 100644 (file)
@@ -396,8 +396,7 @@ void AliEmcalJet::GetMom(TLorentzVector& vec) const
 {
   // Return momentum as four vector.
 
-  Double_t p = fPt * TMath::CosH(fEta);
-  vec.SetPtEtaPhiE(fPt, fEta, fPhi, TMath::Sqrt(p * p + fM * fM));
+  vec.SetPtEtaPhiE(fPt, fEta, fPhi, E());
 }
 
 //__________________________________________________________________________________________________
@@ -478,10 +477,11 @@ Double_t AliEmcalJet::GetZ(const Double_t trkPx, const Double_t trkPy, const Dou
 {
   // Get the z of a constituent inside of a jet
 
-  Double_t pJetSq = this->Px() * this->Px() + this->Py() * this->Py() + this->Pz() * this->Pz();
+  Double_t pJetSq = P();
+  pJetSq *= pJetSq;
 
   if(pJetSq > 1e-6)
-  { return (trkPx * this->Px() + trkPy * this->Py() + trkPz * this->Pz()) / pJetSq ; }
+  { return (trkPx * Px() + trkPy * Py() + trkPz * Pz()) / pJetSq ; }
   else
   { AliWarning(Form("%s: strange, pjet*pjet seems to be zero pJetSq: %f", GetName(), pJetSq)); return -1; }
 }
index c2d7d5e..d913636 100644 (file)
@@ -47,10 +47,10 @@ class AliEmcalJet : public AliVParticle
   Double_t          OneOverPt()                  const { return 1./fPt;  }
   Double_t          Phi()                        const { return fPhi;    }
   Double_t          Theta()                      const { return 2*TMath::ATan(TMath::Exp(-fEta));         }
-  Double_t          E()                          const { Double_t p=P(); return TMath::Sqrt(M()*M()+p*p); }
+  Double_t          E()                          const { Double_t p=P(); return TMath::Sqrt(fM*fM+p*p); }
   Double_t          M()                          const { return fM; }
   Double_t          Eta()                        const { return fEta;    }
-  Double_t          Y()                          const { return 0.5*TMath::Log((E()+Pz())/(E()-Pz()));    }
+  Double_t          Y()                          const { Double_t e = E(); Double_t pz = Pz(); return 0.5*TMath::Log((e+pz)/(e-pz));    }
   Short_t           Charge()                     const { return 0;       }
   Int_t             GetLabel()                   const { return fLabel;  }
   Int_t             PdgCode()                    const { return 0;       }