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;
45 virtual Double_t T() const;
47 virtual Double_t E() const;
48 virtual Double_t M() const;
50 virtual Double_t Eta() const;
51 virtual Double_t Y() const;
53 virtual Short_t Charge() const;
55 virtual Int_t Label() const;
56 virtual Int_t GetLabel() const {return Label();}
57 virtual Int_t PdgCode() const {return fParticle->GetPdgCode();}
58 virtual TParticle* Particle() const {return fParticle;}
61 virtual const Double_t *PID() const {return 0;} // return PID object (to be defined, still)
64 Int_t GetNumberOfTrackReferences() const {return fNTrackRef;}
65 AliTrackReference* GetTrackReference(Int_t i)
66 {return dynamic_cast<AliTrackReference*>((*fTrackReferences)[i]);}
68 // "Trackable" criteria
69 Float_t GetTPCTrackLength(Float_t bz, Float_t ptmin, Int_t &counter, Float_t deadWidth);
71 Int_t GetMother() const {return fMother;}
72 Int_t GetFirstDaughter() const {return fFirstDaughter;}
73 Int_t GetLastDaughter() const {return fLastDaughter;}
74 void SetMother(Int_t idx) {fMother = idx;}
75 void SetFirstDaughter(Int_t idx) {fFirstDaughter = idx;}
76 void SetLastDaughter(Int_t idx) {fLastDaughter = idx;}
77 void SetLabel(Int_t label) {fLabel = label;}
80 TParticle *fParticle; // The wrapped TParticle
81 TRefArray *fTrackReferences; // Reference array to track references
82 Int_t fNTrackRef; // Number of track references
83 Int_t fLabel; // fParticle Label in the Stack
84 Int_t fMother; // Mother particles
85 Int_t fFirstDaughter; // First daughter
86 Int_t fLastDaughter; // LastDaughter
88 ClassDef(AliMCParticle,0) // AliVParticle realisation for MCParticles
91 inline Double_t AliMCParticle::Px() const {return fParticle->Px();}
92 inline Double_t AliMCParticle::Py() const {return fParticle->Py();}
93 inline Double_t AliMCParticle::Pz() const {return fParticle->Pz();}
94 inline Double_t AliMCParticle::Pt() const {return fParticle->Pt();}
95 inline Double_t AliMCParticle::P() const {return fParticle->P(); }
96 inline Double_t AliMCParticle::OneOverPt() const {return 1. / fParticle->Pt();}
97 inline Bool_t AliMCParticle::PxPyPz(Double_t p[3]) const { p[0] = Px(); p[1] = Py(); p[2] = Pz(); return kTRUE; }
98 inline Double_t AliMCParticle::Phi() const {return fParticle->Phi();}
99 inline Double_t AliMCParticle::Theta() const {return fParticle->Theta();}
100 inline Double_t AliMCParticle::Xv() const {return fParticle->Vx();}
101 inline Double_t AliMCParticle::Yv() const {return fParticle->Vy();}
102 inline Double_t AliMCParticle::Zv() const {return fParticle->Vz();}
103 inline Bool_t AliMCParticle::XvYvZv(Double_t x[3]) const { x[0] = Xv(); x[1] = Yv(); x[2] = Zv(); return kTRUE; }
104 inline Double_t AliMCParticle::T() const {return fParticle->T();}
105 inline Double_t AliMCParticle::E() const {return fParticle->Energy();}
106 inline Double_t AliMCParticle::Eta() const {return fParticle->Eta();}
109 inline Double_t AliMCParticle::M() const
111 TParticlePDG* pdg = fParticle->GetPDG();
113 return (pdg->Mass());
115 return (fParticle->GetCalcMass());
120 inline Double_t AliMCParticle::Y() const
125 if (e != TMath::Abs(pz)) {
126 return 0.5*TMath::Log((e+pz)/(e-pz));
132 inline Short_t AliMCParticle::Charge() const
134 TParticlePDG* pdg = fParticle->GetPDG();
136 return (Short_t (pdg->Charge()));
142 inline Int_t AliMCParticle::Label() const {return fLabel;}