1 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
2 * See cxx source for full Copyright notice */
6 // Author: I. Hrivnacova
10 // Fluka implementation of the MonteCarlo interface.
15 #include "TMCProcess.h"
16 #include "TVirtualMC.h"
19 class TG4GeometryManager;
23 class TFluka: public TVirtualMC
26 TFluka(const char* name, const char* title);
29 // TFluka(const TFluka& right);
33 // methods for building/management of geometry
34 // ------------------------------------------------
37 // functions from GBASE
38 virtual void FinishGeometry();
41 // functions from GCONS
42 virtual void Gfmate(Int_t imat, char *name, Float_t &a, Float_t &z,
43 Float_t &dens, Float_t &radl, Float_t &absl,
44 Float_t* ubuf, Int_t& nbuf);
46 // detector composition
47 virtual void Material(Int_t& kmat, const char* name, Float_t a,
48 Float_t z, Float_t dens, Float_t radl, Float_t absl,
49 Float_t* buf, Int_t nwbuf);
50 virtual void Mixture(Int_t& kmat, const char *name, Float_t *a,
51 Float_t *z, Float_t dens, Int_t nlmat, Float_t *wmat);
52 virtual void Medium(Int_t& kmed, const char *name, Int_t nmat,
53 Int_t isvol, Int_t ifield, Float_t fieldm, Float_t tmaxfd,
54 Float_t stemax, Float_t deemax, Float_t epsil,
55 Float_t stmin, Float_t* ubuf, Int_t nbuf);
56 virtual void Matrix(Int_t& krot, Float_t thetaX, Float_t phiX,
57 Float_t thetaY, Float_t phiY, Float_t thetaZ,
59 virtual void Gstpar(Int_t itmed, const char *param, Float_t parval);
60 virtual void SetCerenkov(Int_t itmed, Int_t npckov, Float_t *ppckov,
61 Float_t *absco, Float_t *effic, Float_t *rindex);
63 // functions from GGEOM
64 virtual Int_t Gsvolu(const char *name, const char *shape, Int_t nmed,
65 Float_t *upar, Int_t np);
66 virtual void Gsdvn(const char *name, const char *mother, Int_t ndiv,
68 virtual void Gsdvn2(const char *name, const char *mother, Int_t ndiv,
69 Int_t iaxis, Float_t c0i, Int_t numed);
70 virtual void Gsdvt(const char *name, const char *mother, Float_t step,
71 Int_t iaxis, Int_t numed, Int_t ndvmx);
72 virtual void Gsdvt2(const char *name, const char *mother, Float_t step,
73 Int_t iaxis, Float_t c0, Int_t numed, Int_t ndvmx);
74 virtual void Gsord(const char *name, Int_t iax);
75 virtual void Gspos(const char *name, Int_t nr, const char *mother,
76 Float_t x, Float_t y, Float_t z, Int_t irot,
78 virtual void Gsposp(const char *name, Int_t nr, const char *mother,
79 Float_t x, Float_t y, Float_t z, Int_t irot,
80 const char *konly, Float_t *upar, Int_t np);
81 virtual void Gsbool(const char* onlyVolName, const char* manyVolName) {}
84 virtual void WriteEuclid(const char* fileName, const char* topVol,
85 Int_t number, Int_t nlevel);
88 virtual Int_t VolId(const Text_t* volName) const;
89 virtual const char* VolName(Int_t id) const;
90 virtual Int_t NofVolumes() const;
91 virtual Int_t VolId2Mate(Int_t id) const;
94 // methods for physics management
95 // ------------------------------------------------
98 virtual void BuildPhysics();
101 virtual void SetCut(const char* cutName, Float_t cutValue);
102 virtual void SetProcess(const char* flagName, Int_t flagValue);
103 virtual Float_t Xsec(char* reac, Float_t energy, Int_t part, Int_t mate);
104 virtual void SetExternalDecayer(AliDecayer* decayer); //NEW
107 virtual AliDecayer* Decayer() const; //NEW
109 // particle table usage
110 virtual Int_t IdFromPDG(Int_t pdgID) const;
111 virtual Int_t PDGFromId(Int_t mcID) const;
112 virtual void DefineParticles();
115 // methods for step management
116 // ------------------------------------------------
120 virtual void StopTrack();
121 virtual void StopEvent();
124 virtual void SetMaxStep(Float_t);
125 virtual void SetMaxNStep(Int_t);
126 virtual void SetUserDecay(Int_t); //NEW
129 // tracking volume(s)
130 virtual Int_t CurrentVolID(Int_t& copyNo) const;
131 virtual Int_t CurrentVolOffID(Int_t off, Int_t& copyNo) const;
132 virtual const char* CurrentVolName() const;
133 virtual const char* CurrentVolOffName(Int_t off) const;
134 virtual Int_t CurrentMaterial(Float_t &a, Float_t &z,
135 Float_t &dens, Float_t &radl, Float_t &absl) const;
136 virtual void Gmtod(Float_t* xm, Float_t* xd, Int_t iflag);
137 virtual void Gdtom(Float_t* xd, Float_t* xm, Int_t iflag);
138 virtual Float_t MaxStep() const;
139 virtual Int_t GetMaxNStep() const;
140 virtual Int_t GetMedium() const;
143 // dynamic properties
144 virtual void TrackPosition(TLorentzVector& position) const;
145 virtual void TrackPosition(Double_t &x, Double_t &y, Double_t &z) const;
146 virtual void TrackMomentum(TLorentzVector& momentum) const;
147 virtual void TrackMomentum(Double_t &px, Double_t &py, Double_t &pz, Double_t &etot) const;
148 virtual void TrackVertexPosition(TLorentzVector& position) const;
149 virtual void TrackVertexMomentum(TLorentzVector& momentum) const;
150 virtual Float_t TrackStep() const;
151 virtual Float_t TrackLength() const;
152 virtual Float_t TrackTime() const;
153 virtual Float_t Edep() const;
155 virtual Int_t TrackPid() const;
156 virtual Float_t TrackCharge() const;
157 virtual Float_t TrackMass() const;
158 virtual Float_t Etot() const;
161 virtual Bool_t IsTrackInside() const;
162 virtual Bool_t IsTrackEntering() const;
163 virtual Bool_t IsTrackExiting() const;
164 virtual Bool_t IsTrackOut() const;
165 virtual Bool_t IsTrackDisappeared() const;
166 virtual Bool_t IsTrackStop() const;
167 virtual Bool_t IsTrackAlive() const;
168 virtual Bool_t IsNewTrack() const;
171 virtual Int_t NSecondaries() const;
172 virtual void GetSecondary(Int_t isec, Int_t& particleId,
173 TLorentzVector& position, TLorentzVector& momentum);
174 virtual TMCProcess ProdProcess(Int_t isec) const;
175 virtual Int_t StepProcesses(TArrayI &proc) const;
178 // methods for visualization
179 // ------------------------------------------------
181 // functions for drawing
182 virtual void DrawOneSpec(const char* name);
183 virtual void Gsatt(const char* name, const char* att, Int_t val);
184 virtual void Gdraw(const char* name, Float_t theta, Float_t phi,
185 Float_t psi, Float_t u0, Float_t v0,
186 Float_t ul, Float_t vl);
190 // Geant3 specific methods
191 // !!! need to be transformed to common interface
193 virtual void Gdopt(const char* name , const char* value);
194 virtual void SetClipBox(const char *name, Float_t xmin, Float_t xmax,
195 Float_t ymin, Float_t ymax, Float_t zmin, Float_t zmax);
196 virtual void DefaultRange();
197 virtual void Gdhead(Int_t isel, const char* name, Float_t chrsiz);
198 virtual void Gdman(Float_t u, Float_t v, const char* type);
199 virtual void SetColors();
200 virtual void Gtreve();
201 virtual void GtreveRoot();
202 virtual void Gckmat(Int_t itmed, char* natmed);
203 virtual void InitLego();
204 virtual void Gfpart(Int_t ipart, char *name, Int_t& itrtyp,
205 Float_t& amass, Float_t& charge, Float_t& tlife);
206 virtual void Gspart(Int_t ipart, const char *name, Int_t itrtyp,
207 Float_t amass, Float_t charge, Float_t tlife);
210 // methods for run control
211 // ------------------------------------------------
215 virtual void ProcessEvent();
216 virtual void ProcessRun(Int_t nofEvents);
218 // UI control methods
221 void ProcessGeantMacro(const char* macroName);
222 void ProcessGeantCommand(const char* commandPath);
225 virtual Int_t CurrentEvent() const;
226 virtual AliMCGeomType GetMCGeomType() const { return kGeant4; }
230 TFluka(const TFluka& right);
233 TFluka& operator=(const TFluka& right);
237 TG4GeometryManager* fGeometryManager; //geometry manager