1 #ifndef AliMCParticle_H
2 #define AliMCParticle_H
3 /* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 //-------------------------------------------------------------------------
9 // AliVParticle realisation for MC Particles
10 // Author: Andreas Morsch, CERN
11 //-------------------------------------------------------------------------
14 #include <TParticle.h>
15 #include <TParticlePDG.h>
16 #include <TRefArray.h>
18 #include "AliTrackReference.h"
19 #include "AliVParticle.h"
21 class AliMCParticle: public AliVParticle {
24 AliMCParticle(TParticle* part, TRefArray* rarray = 0, Int_t label=-1);
25 virtual ~AliMCParticle();
26 AliMCParticle(const AliMCParticle& mcPart);
27 AliMCParticle& operator=(const AliMCParticle& mcPart);
30 virtual Double_t Px() const;
31 virtual Double_t Py() const;
32 virtual Double_t Pz() const;
33 virtual Double_t Pt() const;
34 virtual Double_t P() const;
35 virtual Bool_t PxPyPz(Double_t p[3]) const;
37 virtual Double_t OneOverPt() const;
38 virtual Double_t Phi() const;
39 virtual Double_t Theta() const;
41 virtual Double_t Xv() const;
42 virtual Double_t Yv() const;
43 virtual Double_t Zv() const;
44 virtual Bool_t XvYvZv(Double_t x[3]) const;
46 virtual Double_t E() const;
47 virtual Double_t M() const;
49 virtual Double_t Eta() const;
50 virtual Double_t Y() const;
52 virtual Short_t Charge() const;
54 virtual Int_t Label() const;
55 virtual TParticle* Particle() const {return fParticle;}
58 virtual const Double_t *PID() const {return 0;} // return PID object (to be defined, still)
61 Int_t GetNumberOfTrackReferences() {return fNTrackRef;}
62 AliTrackReference* GetTrackReference(Int_t i)
63 {return dynamic_cast<AliTrackReference*>((*fTrackReferences)[i]);}
65 TParticle *fParticle; // The wrapped TParticle
66 TRefArray *fTrackReferences; // Reference array to track references
67 Int_t fNTrackRef; // Number of track references
68 Int_t fLabel; // fParticle Label in the Stack
70 ClassDef(AliMCParticle,0) // AliVParticle realisation for MCParticles
73 inline Double_t AliMCParticle::Px() const {return fParticle->Px();}
74 inline Double_t AliMCParticle::Py() const {return fParticle->Py();}
75 inline Double_t AliMCParticle::Pz() const {return fParticle->Pz();}
76 inline Double_t AliMCParticle::Pt() const {return fParticle->Pt();}
77 inline Double_t AliMCParticle::P() const {return fParticle->P(); }
78 inline Double_t AliMCParticle::OneOverPt() const {return 1. / fParticle->Pt();}
79 inline Bool_t AliMCParticle::PxPyPz(Double_t p[3]) const { p[0] = Px(); p[1] = Py(); p[2] = Pz(); return kTRUE; }
80 inline Double_t AliMCParticle::Phi() const {return fParticle->Phi();}
81 inline Double_t AliMCParticle::Theta() const {return fParticle->Theta();}
82 inline Double_t AliMCParticle::Xv() const {return fParticle->Vx();}
83 inline Double_t AliMCParticle::Yv() const {return fParticle->Vy();}
84 inline Double_t AliMCParticle::Zv() const {return fParticle->Vz();}
85 inline Bool_t AliMCParticle::XvYvZv(Double_t x[3]) const { x[0] = Xv(); x[1] = Yv(); x[2] = Zv(); return kTRUE; }
86 inline Double_t AliMCParticle::E() const {return fParticle->Energy();}
87 inline Double_t AliMCParticle::Eta() const {return fParticle->Eta();}
89 inline Double_t AliMCParticle::M() const
91 TParticlePDG* pdg = fParticle->GetPDG();
95 return (fParticle->GetCalcMass());
100 inline Double_t AliMCParticle::Y() const
103 Double_t pz = TMath::Abs(Pz());
105 return 0.5*TMath::Log((e+pz)/(e-pz));
111 inline Short_t AliMCParticle::Charge() const
113 TParticlePDG* pdg = fParticle->GetPDG();
115 return (Short_t (pdg->Charge()));
121 inline Int_t AliMCParticle::Label() const {return fLabel;}