PDG definition codes in AliPDG.h -- to go to Root later??
[u/mrichter/AliRoot.git] / STEER / GParticle.h
CommitLineData
fe4da5cc 1#ifndef ROOT_GParticle
2#define ROOT_GParticle
3
4#include "TObject.h"
5#include "TMath.h"
6
7#define Keep_Bit 1
8#define Children_Bit 2
9#define Done_Bit 4
10
11class GParticle : public TObject {
12
13private:
14
15 Int_t fKS; // status of particle ( LUJETS K[1] )
16 Int_t fKF; // Geant3 particle type
17 Int_t fParent; // parent id
18 Int_t fFirstChild; // id of first child
19 Int_t fLastChild; // id of last child
20
21 Float_t fPx; // X momentum [GeV/c]
22 Float_t fPy; // Y momentum [GeV/c]
23 Float_t fPz; // Z momentum [GeV/c]
24 Float_t fEnergy; // Energy [GeV]
25 Float_t fMass; // Mass [Gev/c^2]
26
27 Float_t fVx; // X vertex [cm]
28 Float_t fVy; // Y vertex [cm]
29 Float_t fVz; // Z vertex [cm]
30
31 Float_t fPolx; // X component of polarisation
32 Float_t fPoly; // Y component of polarisation
33 Float_t fPolz; // Z component of polarisation
34
35 Float_t fTime; // time of procuction
36 Float_t fLifeTime; // proper lifetime
37 Float_t fProcessTime; // Cpu Time to process this track
38
39 char fOrigin[11]; // generation mechanism
40 Float_t fWgt;
41 Float_t fPT;
42 Float_t fTheta;
43 Float_t fPhi;
44 Float_t fEta;
45
46public:
47 GParticle() { }
48
49 GParticle(const GParticle &p);
50
51 GParticle(const Int_t kS, const Int_t kF,
52 const Int_t parent, const Int_t firstchild, const Int_t lastchild,
53 const Float_t px, const Float_t py, const Float_t pz,
54 const Float_t energy, const Float_t mass,
55 const Float_t vx, const Float_t vy, const Float_t vz,
56 const Float_t polx, const Float_t poly, const Float_t polz,
57 const Float_t time, const Float_t lifetime,
58 const char *Origin="Unknown",Float_t wgt=1.) :
59
60 fKS(kS),
61 fKF(kF),
62 fParent(parent),
63 fFirstChild(firstchild),
64 fLastChild(lastchild),
65 fPx(px),
66 fPy(py),
67 fPz(pz),
68 fEnergy(energy),
69 fMass(mass),
70 fVx(vx),
71 fVy(vy),
72 fVz(vz),
73 fPolx(polx),
74 fPoly(poly),
75 fPolz(polz),
76 fTime(time),
77 fLifeTime(lifetime),
78 fProcessTime(-0.1),
79 fWgt(wgt)
80 {Int_t p=0;
81 while((fOrigin[p]=Origin[p]) && p++<10);fOrigin[10]=0;
82 fPT =TMath::Sqrt(fPx*fPx+fPy*fPy);
83 fTheta =Float_t(TMath::ATan2(Double_t(fPT),Double_t(fPz)));
84 fPhi =Float_t(TMath::ATan2(Double_t(fPy),Double_t(fPx)))+TMath::Pi();
85 fEta =-TMath::Log(TMath::Tan(fTheta/2));
86 }
87
88 virtual ~GParticle() { }
89
90 inline Int_t GetKS() const {return fKS;}
91 inline Int_t GetKF() const {return fKF;}
92 inline Int_t GetParent() const {return fParent;}
93 inline Int_t GetFirstChild() const {return fFirstChild;}
94 inline Int_t GetLastChild() const {return fLastChild;}
95
96 inline Float_t GetPx() const {return fPx;}
97 inline Float_t GetPy() const {return fPy;}
98 inline Float_t GetPz() const {return fPz;}
99
100 inline Float_t GetPolx() const {return fPolx;}
101 inline Float_t GetPoly() const {return fPoly;}
102 inline Float_t GetPolz() const {return fPolz;}
103
104 inline Float_t GetEnergy() const {return fEnergy;}
105 inline Float_t GetMass() const {return fMass;}
106 inline Float_t GetMomentum() const {return TMath::Sqrt(fPx*fPx+fPy*fPy+fPz*fPz);}
107 inline void GetOrigin(char *Origin) {strcpy(Origin,fOrigin);}
108 inline Float_t GetPT() const {return fPT;}
109 inline Float_t GetTheta() const {return fTheta;}
110 inline Float_t GetEta() const {return fEta;}
111 inline Float_t GetPhi() const {return fPhi;}
112
113 inline Float_t GetVx() const {return fVx;}
114 inline Float_t GetVy() const {return fVy;}
115 inline Float_t GetVz() const {return fVz;}
116 inline Float_t GetTime() const {return fTime;}
117 inline Float_t GetLifeTime() const {return fLifeTime;}
118 inline Float_t GetWgt() const {return fWgt;}
119 inline void SetWgt(Float_t wgt) {fWgt=wgt;}
120 inline Float_t GetProcessTime() const {return fProcessTime;}
121 virtual const Text_t *GetName() const;
122 virtual const Text_t *GetTitle() const;
123
124 inline void SetKS(Int_t KS) {fKS=KS;}
125 inline void SetFirstChild(const Int_t firstchild) {fFirstChild=firstchild;}
126 inline void SetLastChild(const Int_t lastchild) {fLastChild=lastchild;}
127 inline void SetParent(const Int_t parent) {fParent=parent;}
128 inline void SetProcessTime(Float_t cputime) {fProcessTime=cputime;}
129
130 GParticle &operator=(const GParticle &p);
131
132 ClassDef(GParticle,1) // MonteCarlo Particle Class
133};
134
135inline GParticle &GParticle::operator = (const GParticle &p)
136{
137 fKS = p.fKS;
138 fKF = p.fKF;
139 fParent = p.fParent;
140 fFirstChild = p.fFirstChild;
141 fLastChild = p.fLastChild;
142 fPx = p.fPx;
143 fPy = p.fPy;
144 fPz = p.fPz;
145 fEnergy = p.fEnergy;
146 fMass = p.fMass;
147 fVx = p.fVx;
148 fVy = p.fVy;
149 fVz = p.fVz;
150 fPolx = p.fPolx;
151 fPoly = p.fPoly;
152 fPolz = p.fPolz;
153 fTime = p.fTime;
154 fLifeTime = p.fLifeTime;
155 fProcessTime = p.fProcessTime;
156 fWgt = p.fWgt;
157 fPT = p.fPT;
158 fTheta = p.fTheta;
159 fPhi = p.fPhi;
160 fEta = p.fEta;
161 strcpy(fOrigin,p.fOrigin);
162 return *this;
163}
164
165
166inline GParticle::GParticle(const GParticle &p) { *this = p; }
167
168#endif
169
170