Time added (A. Kisiel)

diff --git a/STEER/AliAODMCParticle.cxx b/STEER/AliAODMCParticle.cxx
index bbe75d763db4b88a1f06e57a5628767088aaf33b..ee25c4de4b039ba03947e055715343a696b03ace 100644 (file)
--- a/STEER/AliAODMCParticle.cxx
+++ b/STEER/AliAODMCParticle.cxx
@@ -46,7 +46,8 @@ AliVParticle(),
   fE(0),
   fVx(0),
   fVy(0),
-  fVz(0)
+  fVz(0),
+  fVt(0)
 {
   // Default Constructor
   fDaughter[0] =   fDaughter[1] = 0;
@@ -65,7 +66,8 @@ AliAODMCParticle::AliAODMCParticle(AliMCParticle* mcpart, Int_t label,Int_t flag
     fE(mcpart->Particle()->Energy()),
     fVx(mcpart->Particle()->Vx()),
     fVy(mcpart->Particle()->Vy()),
-    fVz(mcpart->Particle()->Vz())
+    fVz(mcpart->Particle()->Vz()),
+    fVt(mcpart->Particle()->T())
 {
   fDaughter[0] =  mcpart->GetFirstDaughter(); 
   fDaughter[1] =  mcpart->GetLastDaughter();
@@ -86,7 +88,8 @@ AliAODMCParticle::AliAODMCParticle(const AliAODMCParticle& mcPart) :
     fE(mcPart.fE),
     fVx(mcPart.fVx),
     fVy(mcPart.fVy),
-    fVz(mcPart.fVz)
+    fVz(mcPart.fVz),
+    fVt(mcPart.fVt)
 {
   // Copy constructor
   fDaughter[0] = mcPart.fDaughter[0]; 
@@ -110,6 +113,7 @@ AliAODMCParticle& AliAODMCParticle::operator=(const AliAODMCParticle& mcPart)
     fVx         = mcPart.fVx;
     fVy         = mcPart.fVy;
     fVz         = mcPart.fVz;
+    fVt         = mcPart.fVt;
     fDaughter[0] = mcPart.fDaughter[0]; 
     fDaughter[1] = mcPart.fDaughter[1]; 
   }  
@@ -146,6 +150,6 @@ void AliAODMCParticle::Print(const Option_t */*opt*/) const {
   else{
     Printf(">>> PDG (%d) : %s",fPdgCode,"Unknown");
   }
-  Printf(">>  P(%3.3f,%3.3f,%3.3f) V((%3.3f,%3.3f,%3.3f)",fPx,fPy,fPz,fVx,fVy,fVz);  
+  Printf(">>  P(%3.3f,%3.3f,%3.3f) V((%3.3f,%3.3f,%3.3f,%3.3f)",fPx,fPy,fPz,fVx,fVy,fVz,fVt);  
   Printf(">   Mother %d, First Daughter %d Last Daughter %d Process %d",fMother,fDaughter[0],fDaughter[1],TObject::GetUniqueID());
 }

diff --git a/STEER/AliAODMCParticle.h b/STEER/AliAODMCParticle.h
index 9e38e6d7c84adb65d31f4b5b6cf163d8c2f27183..fbf30db500b38df6459eb4ada51943acd976c5fe 100644 (file)
--- a/STEER/AliAODMCParticle.h
+++ b/STEER/AliAODMCParticle.h
@@ -48,6 +48,7 @@ class AliAODMCParticle: public AliVParticle {
     virtual Double_t Yv() const;
     virtual Double_t Zv() const;
     virtual Bool_t   XvYvZv(Double_t x[3]) const;  
+    virtual Double_t T() const;
 
     virtual Double_t E()          const;
     virtual Double_t M()          const;
@@ -57,8 +58,8 @@ class AliAODMCParticle: public AliVParticle {
     
     virtual Short_t Charge()      const;
 
-    virtual Int_t      Label()       const;
-    virtual Int_t      GetLabel() const {return Label();} 
+    virtual Int_t   Label()       const;
+    virtual Int_t   GetLabel()    const {return Label();} 
            
     // PID
     virtual const Double_t *PID() const {return 0;} // return PID object (to be defined, still)
@@ -108,6 +109,7 @@ class AliAODMCParticle: public AliVParticle {
   Double32_t       fVx;                   // [0.,0.,12] x of production vertex
   Double32_t       fVy;                   // [0.,0.,12] y of production vertex
   Double32_t       fVz;                   // [0.,0.,12] z of production vertex
+  Double32_t       fVt;                   // [0.,0.,12] t of production vertex
 
   // Copy the uniquID to another data member? unique ID is correctly handled 
   // via TOBject Copy construct but not by AliVParticle ctor (no passing of 
@@ -140,12 +142,13 @@ inline Double_t AliAODMCParticle::Pt()        const {return TMath::Sqrt(fPx*fPx+
 inline Double_t AliAODMCParticle::P()         const {return TMath::Sqrt(fPx*fPx+fPy*fPy+fPz*fPz); }
 inline Double_t AliAODMCParticle::OneOverPt() const {return 1. / Pt();}
 inline Bool_t   AliAODMCParticle::PxPyPz(Double_t p[3]) const { p[0] = fPx; p[1] = fPy; p[2] = fPz; return kTRUE; }
-inline Double_t AliAODMCParticle::Phi()       const { return TMath::Pi()+TMath::ATan2(-fPy,-fPx); }  // note that Phi() returns an angle between 0 and 2pi
-inline Double_t AliAODMCParticle::Theta()     const { return (fPz==0)?TMath::PiOver2():TMath::ACos(fPz/P()); }
+inline Double_t AliAODMCParticle::Phi()       const {return TMath::Pi()+TMath::ATan2(-fPy,-fPx); }  // note that Phi() returns an angle between 0 and 2pi
+inline Double_t AliAODMCParticle::Theta()     const {return (fPz==0)?TMath::PiOver2():TMath::ACos(fPz/P()); }
 inline Double_t AliAODMCParticle::Xv()        const {return fVx;}
 inline Double_t AliAODMCParticle::Yv()        const {return fVy;}
 inline Double_t AliAODMCParticle::Zv()        const {return fVz;}
 inline Bool_t   AliAODMCParticle::XvYvZv(Double_t x[3]) const { x[0] = fVx; x[1] = fVy; x[2] = fVz; return kTRUE; }
+inline Double_t AliAODMCParticle::T()         const {return fVt;}
 inline Double_t AliAODMCParticle::E()         const {return fE;}
 inline Double_t AliAODMCParticle::Eta()       const {  
   Double_t pmom = P();