]>
Commit | Line | Data |
---|---|---|
1 | #ifndef ALIGENERATOR_H | |
2 | #define ALIGENERATOR_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | /* $Id$ */ | |
7 | ||
8 | /////////////////////////////////////////////////////////// | |
9 | // // | |
10 | // Class to generate the particles for the MC // | |
11 | // The base class is empty // | |
12 | // // | |
13 | /////////////////////////////////////////////////////////// | |
14 | ||
15 | #include "TLorentzVector.h" | |
16 | #include "TArrayF.h" | |
17 | #include "AliRndm.h" | |
18 | #include "AliMCProcess.h" | |
19 | class TGenerator; | |
20 | class AliStack; | |
21 | ||
22 | ||
23 | typedef enum { kNoSmear, kPerEvent, kPerTrack } VertexSmear_t; | |
24 | typedef enum { kExternal, kInternal} VertexSource_t; | |
25 | ||
26 | ||
27 | class AliGenerator : public TNamed, public AliRndm | |
28 | { | |
29 | ||
30 | public: | |
31 | AliGenerator(); | |
32 | AliGenerator(Int_t npart); | |
33 | AliGenerator(const AliGenerator &gen); | |
34 | virtual ~AliGenerator(); | |
35 | virtual void Init(); | |
36 | void Copy(AliGenerator &gen) const; | |
37 | virtual void SetOrigin(Float_t ox, Float_t oy, Float_t oz); | |
38 | virtual void SetOrigin(const TLorentzVector &o); | |
39 | virtual void SetSigma(Float_t sx, Float_t sy, Float_t sz); | |
40 | virtual void SetMomentumRange(Float_t pmin=0, Float_t pmax=1.e10); | |
41 | virtual void SetPtRange(Float_t ptmin=0, Float_t ptmax=1.e10); | |
42 | virtual void SetPhiRange(Float_t phimin=-180., Float_t phimax=180); | |
43 | virtual void SetYRange(Float_t ymin=-100, Float_t ymax=100); | |
44 | virtual void SetVRange(Float_t vxmin, Float_t vxmax, | |
45 | Float_t vymin, Float_t vymax, | |
46 | Float_t vzmin, Float_t vzmax); | |
47 | virtual void SetNumberParticles(Int_t npart=100) {fNpart=npart;} | |
48 | virtual Int_t NumberParticles() const {return fNpart;} | |
49 | virtual void SetThetaRange(Float_t thetamin=0, Float_t thetamax=180); | |
50 | virtual void Generate()=0; | |
51 | virtual void SetParentWeight(Float_t wgt) {fParentWeight=wgt;} | |
52 | virtual void SetChildWeight(Float_t wgt) {fChildWeight=wgt;} | |
53 | virtual void SetAnalog(Int_t flag=1) {fAnalog=flag;} | |
54 | virtual void SetVertexSmear(VertexSmear_t smear) {fVertexSmear = smear;} | |
55 | virtual void SetCutVertexZ(Float_t cut=999999.) {fCutVertexZ = cut;} | |
56 | virtual void SetVertexSource(VertexSource_t smear) {fVertexSource = kInternal;} | |
57 | virtual void SetTrackingFlag(Int_t flag=1) {fTrackIt=flag;} | |
58 | void Vertex(); | |
59 | void VertexExternal(); | |
60 | virtual void VertexInternal(); | |
61 | virtual void FinishRun(){;} | |
62 | ||
63 | virtual void SetMC(TGenerator *theMC) | |
64 | {if (!fgMCEvGen) fgMCEvGen =theMC;} | |
65 | ||
66 | AliGenerator & operator=(const AliGenerator &gen); | |
67 | ||
68 | // Getters | |
69 | ||
70 | virtual void GetOrigin(Float_t &ox, Float_t &oy, Float_t &oz) const | |
71 | {ox=fOrigin.At(0);oy=fOrigin.At(1);oz=fOrigin.At(2);} | |
72 | virtual void GetOrigin(TLorentzVector &o) const | |
73 | {o[0]=fOrigin.At(0);o[1]=fOrigin.At(1);o[2]=fOrigin.At(2);o[3]=0;} | |
74 | ||
75 | void SetStack (AliStack *stack) {fStack = stack;} | |
76 | AliStack* GetStack(){return fStack;} | |
77 | // Comminication with stack | |
78 | protected: | |
79 | virtual void SetTrack(Int_t done, Int_t parent, Int_t pdg, | |
80 | Float_t *pmom, Float_t *vpos, Float_t *polar, | |
81 | Float_t tof, AliMCProcess mech, Int_t &ntr, | |
82 | Float_t weight = 1, Int_t is = 0); | |
83 | virtual void SetTrack(Int_t done, Int_t parent, Int_t pdg, | |
84 | Double_t px, Double_t py, Double_t pz, Double_t e, | |
85 | Double_t vx, Double_t vy, Double_t vz, Double_t tof, | |
86 | Double_t polx, Double_t poly, Double_t polz, | |
87 | AliMCProcess mech, Int_t &ntr, Float_t weight = 1, Int_t is = 0); | |
88 | virtual void KeepTrack(Int_t itrack); | |
89 | virtual void SetHighWaterMark(Int_t nt); | |
90 | ||
91 | protected: | |
92 | static TGenerator* fgMCEvGen; //Pointer to the generator | |
93 | Float_t fThetaMin; //Minimum theta of generation in radians | |
94 | Float_t fThetaMax; //Maximum theta of generation in radians | |
95 | Float_t fPhiMin; //Minimum phi of generation in radians | |
96 | Float_t fPhiMax; //Maximum phi of generation in radians | |
97 | Float_t fPMin; //Minimum momentum of generation in GeV/c | |
98 | Float_t fPMax; //Minimum momentum of generation in GeV/c | |
99 | Float_t fPtMin; //Minimum transverse momentum | |
100 | Float_t fPtMax; //Maximum transverse momentum | |
101 | Float_t fYMin; //Minimum rapidity | |
102 | Float_t fYMax; //Maximum rapidity | |
103 | TArrayF fVMin; //Minimum Decaylength | |
104 | TArrayF fVMax; //Minimum Decaylength | |
105 | Int_t fNpart; //Maximum number of particles per event | |
106 | Float_t fParentWeight; //Parent Weight | |
107 | Float_t fChildWeight; //ChildWeight | |
108 | Int_t fAnalog; //Flag for anolog or pt-weighted generation | |
109 | // | |
110 | VertexSmear_t fVertexSmear; //Vertex Smearing mode | |
111 | VertexSource_t fVertexSource; //Vertex source (internal/external) | |
112 | Float_t fCutVertexZ; // Vertex cut in units of sigma_z | |
113 | Int_t fTrackIt; // if 1, Track final state particles | |
114 | TArrayF fOrigin; // Origin of event | |
115 | TArrayF fOsigma; // Sigma of the Origin of event | |
116 | TArrayF fVertex; //! Vertex of current event | |
117 | TArrayF fEventVertex; //!The current event vertex | |
118 | ||
119 | AliStack* fStack; //! Local pointer to stack | |
120 | /*************************************************************************/ | |
121 | enum {kThetaRange = BIT(14), | |
122 | kVertexRange = BIT(15), | |
123 | kPhiRange = BIT(16), | |
124 | kPtRange = BIT(17), | |
125 | kYRange = BIT(18), | |
126 | kMomentumRange = BIT(19) | |
127 | }; | |
128 | ClassDef(AliGenerator,2) // Base class for event generators | |
129 | }; | |
130 | ||
131 | #endif | |
132 | ||
133 | ||
134 | ||
135 | ||
136 | ||
137 | ||
138 | ||
139 | ||
140 | ||
141 | ||
142 | ||
143 | ||
144 | ||
145 |