- Double_t E(UInt_t pdg) const;
- Double_t Y(UInt_t pdg) const {return 0.5*TMath::Log((E(pdg)+Pz())/(E(pdg)-Pz()+1.e-13));}
- Double_t DecayLength(Double_t point[3]) const;
+ Double_t E2(UInt_t pdg) const;
+ Double_t E(UInt_t pdg) const {return TMath::Sqrt(E2(pdg));}
+ Double_t Y(UInt_t pdg) const {Double_t e=E(pdg); return 0.5*TMath::Log((e+Pz())/(e-Pz()+1.e-13));}
+ Double_t DecayLength2(Double_t point[3]) const;
+ Double_t DecayLength(Double_t point[3]) const {return TMath::Sqrt(DecayLength2(point));}
+ Double_t DecayLength2(AliAODVertex *vtx1) const
+ {return GetSecondaryVtx()->Distance2ToVertex(vtx1);}