]>
Commit | Line | Data |
---|---|---|
ef42d733 | 1 | #ifndef TGEANT3_H |
2 | #define TGEANT3_H | |
3da30618 | 3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | /* $Id$ */ | |
7 | ||
fe4da5cc | 8 | //////////////////////////////////////////////// |
9 | // C++ interface to Geant3 basic routines // | |
10 | //////////////////////////////////////////////// | |
11 | ||
12 | #include <AliMC.h> | |
13 | ||
14 | //______________________________________________________________ | |
15 | // | |
16 | // Geant3 prototypes for commons | |
17 | // | |
18 | //______________________________________________________________ | |
19 | // | |
20 | ||
21 | //----------QUEST | |
22 | // COMMON/QUEST/IQUEST(100) | |
23 | typedef struct { | |
24 | Int_t iquest[100]; | |
25 | } Quest_t; | |
26 | ||
7ac3f11b | 27 | //----------GCBANK |
28 | // COMMON/GCBANK/NZEBRA,GVERSN,ZVERSN,IXSTOR,IXDIV,IXCONS,FENDQ(16) | |
29 | // + ,LMAIN,LR1,WS(KWBANK) | |
30 | typedef struct { | |
31 | Int_t nzebra; | |
32 | Float_t gversn; | |
33 | Float_t zversn; | |
34 | Int_t ixstor; | |
35 | Int_t ixdiv; | |
36 | Int_t ixcons; | |
37 | Float_t fendq[16]; | |
38 | Int_t lmain; | |
39 | Int_t lr1; | |
40 | } Gcbank_t; | |
41 | ||
fe4da5cc | 42 | //----------GCLINK |
43 | // COMMON/GCLINK/JDIGI ,JDRAW ,JHEAD ,JHITS ,JKINE ,JMATE ,JPART | |
44 | // + ,JROTM ,JRUNG ,JSET ,JSTAK ,JGSTAT,JTMED ,JTRACK,JVERTX | |
45 | // + ,JVOLUM,JXYZ ,JGPAR ,JGPAR2,JSKLT | |
46 | typedef struct { | |
47 | Int_t jdigi; | |
48 | Int_t jdraw; | |
49 | Int_t jhead; | |
50 | Int_t jhits; | |
51 | Int_t jkine; | |
52 | Int_t jmate; | |
53 | Int_t jpart; | |
54 | Int_t jrotm; | |
55 | Int_t jrung; | |
56 | Int_t jset; | |
57 | Int_t jstak; | |
58 | Int_t jgstat; | |
59 | Int_t jtmed; | |
60 | Int_t jtrack; | |
61 | Int_t jvertx; | |
62 | Int_t jvolum; | |
63 | Int_t jxyz; | |
64 | Int_t jgpar; | |
65 | Int_t jgpar2; | |
66 | Int_t jsklt; | |
67 | } Gclink_t; | |
68 | ||
69 | ||
70 | //----------GCFLAG | |
71 | // COMMON/GCFLAG/IDEBUG,IDEMIN,IDEMAX,ITEST,IDRUN,IDEVT,IEORUN | |
72 | // + ,IEOTRI,IEVENT,ISWIT(10),IFINIT(20),NEVENT,NRNDM(2) | |
73 | typedef struct { | |
74 | Int_t idebug; | |
75 | Int_t idemin; | |
76 | Int_t idemax; | |
77 | Int_t itest; | |
78 | Int_t idrun; | |
79 | Int_t idevt; | |
80 | Int_t ieorun; | |
81 | Int_t ieotri; | |
82 | Int_t ievent; | |
83 | Int_t iswit[10]; | |
84 | Int_t ifinit[20]; | |
85 | Int_t nevent; | |
86 | Int_t nrndm[2]; | |
87 | } Gcflag_t; | |
88 | ||
89 | //----------GCKINE | |
90 | // COMMON/GCKINE/IKINE,PKINE(10),ITRA,ISTAK,IVERT,IPART,ITRTYP | |
91 | // + ,NAPART(5),AMASS,CHARGE,TLIFE,VERT(3),PVERT(4),IPAOLD | |
92 | typedef struct { | |
93 | Int_t ikine; | |
94 | Float_t pkine[10]; | |
95 | Int_t itra; | |
96 | Int_t istak; | |
97 | Int_t ivert; | |
98 | Int_t ipart; | |
99 | Int_t itrtyp; | |
100 | Int_t napart[5]; | |
101 | Float_t amass; | |
102 | Float_t charge; | |
103 | Float_t tlife; | |
104 | Float_t vert[3]; | |
105 | Float_t pvert[4]; | |
106 | Int_t ipaold; | |
107 | } Gckine_t; | |
108 | ||
109 | //----------GCKING | |
110 | // COMMON/GCKING/KCASE,NGKINE,GKIN(5,MXGKIN), | |
111 | // + TOFD(MXGKIN),IFLGK(MXGKIN) | |
112 | #define MXGKIN 100 | |
113 | typedef struct { | |
114 | Int_t kcase; | |
115 | Int_t ngkine; | |
116 | Float_t gkin[MXGKIN][5]; | |
aee8290b | 117 | Float_t tofd[MXGKIN]; |
fe4da5cc | 118 | Int_t iflgk[MXGKIN]; |
119 | } Gcking_t; | |
120 | ||
121 | //----------GCKIN2 | |
122 | // COMMON/GCKIN2/NGPHOT,XPHOT(11,MXPHOT) | |
123 | #define MXPHOT 800 | |
124 | typedef struct { | |
125 | Int_t ngphot; | |
126 | Float_t xphot[MXPHOT][11]; | |
127 | } Gckin2_t; | |
128 | ||
129 | //----------GCKIN3 | |
130 | // COMMON/GCKIN3/GPOS(3,MXGKIN) | |
131 | typedef struct { | |
132 | Float_t gpos[MXGKIN][3]; | |
133 | } Gckin3_t; | |
134 | ||
135 | //----------GCMATE | |
136 | // COMMON/GCMATE/NMAT,NAMATE(5),A,Z,DENS,RADL,ABSL | |
137 | typedef struct { | |
138 | Int_t nmat; | |
139 | Int_t namate[5]; | |
140 | Float_t a; | |
141 | Float_t z; | |
142 | Float_t dens; | |
143 | Float_t radl; | |
144 | Float_t absl; | |
145 | } Gcmate_t; | |
146 | ||
147 | //----------GCTMED | |
148 | // COMMON/GCTMED/NUMED,NATMED(5),ISVOL,IFIELD,FIELDM,TMAXFD,STEMAX | |
149 | // + ,DEEMAX,EPSIL,STMIN,CFIELD,PREC,IUPD,ISTPAR,NUMOLD | |
150 | typedef struct { | |
151 | Int_t numed; | |
152 | Int_t natmed[5]; | |
153 | Int_t isvol; | |
154 | Int_t ifield; | |
155 | Float_t fieldm; | |
156 | Float_t tmaxfd; | |
157 | Float_t stemax; | |
158 | Float_t deemax; | |
159 | Float_t epsil; | |
160 | Float_t stmin; | |
161 | Float_t cfield; | |
162 | Float_t prec; | |
163 | Int_t iupd; | |
164 | Int_t istpar; | |
165 | Int_t numold; | |
166 | } Gctmed_t; | |
167 | ||
168 | //----------GCTRAK | |
169 | #define MAXMEC 30 | |
170 | // PARAMETER (MAXMEC=30) | |
171 | // COMMON/GCTRAK/VECT(7),GETOT,GEKIN,VOUT(7),NMEC,LMEC(MAXMEC) | |
172 | // + ,NAMEC(MAXMEC),NSTEP ,MAXNST,DESTEP,DESTEL,SAFETY,SLENG | |
173 | // + ,STEP ,SNEXT ,SFIELD,TOFG ,GEKRAT,UPWGHT,IGNEXT,INWVOL | |
174 | // + ,ISTOP ,IGAUTO,IEKBIN, ILOSL, IMULL,INGOTO,NLDOWN,NLEVIN | |
175 | // + ,NLVSAV,ISTORY | |
176 | typedef struct { | |
177 | Float_t vect[7]; | |
178 | Float_t getot; | |
179 | Float_t gekin; | |
180 | Int_t vout[7]; | |
181 | Int_t nmec; | |
182 | Int_t lmec[MAXMEC]; | |
183 | Int_t namec[MAXMEC]; | |
184 | Int_t nstep; | |
185 | Int_t maxnst; | |
186 | Float_t destep; | |
187 | Float_t destel; | |
188 | Float_t safety; | |
189 | Float_t sleng; | |
190 | Float_t step; | |
191 | Float_t snext; | |
192 | Float_t sfield; | |
193 | Float_t tofg; | |
194 | Float_t gekrat; | |
195 | Float_t upwght; | |
196 | Int_t ignext; | |
197 | Int_t inwvol; | |
198 | Int_t istop; | |
199 | Int_t igauto; | |
200 | Int_t iekbin; | |
201 | Int_t ilosl; | |
202 | Int_t imull; | |
203 | Int_t ingoto; | |
204 | Int_t nldown; | |
205 | Int_t nlevin; | |
206 | Int_t nlsav; | |
207 | Int_t istory; | |
208 | } Gctrak_t; | |
209 | ||
210 | //----------GCVOLU | |
211 | // COMMON/GCVOLU/NLEVEL,NAMES(15),NUMBER(15), | |
212 | // +LVOLUM(15),LINDEX(15),INFROM,NLEVMX,NLDEV(15),LINMX(15), | |
213 | // +GTRAN(3,15),GRMAT(10,15),GONLY(15),GLX(3) | |
214 | typedef struct { | |
215 | Int_t nlevel; | |
216 | Int_t names[15]; | |
217 | Int_t number[15]; | |
218 | Int_t lvolum[15]; | |
219 | Int_t lindex[15]; | |
220 | Int_t infrom; | |
221 | Int_t nlevmx; | |
222 | Int_t nldev[15]; | |
223 | Int_t linmx[15]; | |
224 | Float_t gtran[15][3]; | |
225 | Float_t grmat[15][10]; | |
226 | Float_t gonly[15]; | |
227 | Float_t glx[3]; | |
228 | } Gcvolu_t; | |
229 | ||
230 | //----------GCSETS | |
231 | // COMMON/GCSETS/IHSET,IHDET,ISET,IDET,IDTYPE,NVNAME,NUMBV(20) | |
232 | typedef struct { | |
233 | Int_t ihset; | |
234 | Int_t ihdet; | |
235 | Int_t iset; | |
236 | Int_t idet; | |
237 | Int_t idtype; | |
238 | Int_t nvname; | |
239 | Int_t numbv[20]; | |
240 | } Gcsets_t; | |
241 | ||
242 | //----------GCNUM | |
243 | // COMMON/GCNUM/NMATE ,NVOLUM,NROTM,NTMED,NTMULT,NTRACK,NPART | |
244 | // + ,NSTMAX,NVERTX,NHEAD,NBIT | |
245 | typedef struct { | |
246 | Int_t nmate; | |
247 | Int_t nvolum; | |
248 | Int_t nrotm; | |
249 | Int_t ntmed; | |
250 | Int_t ntmult; | |
251 | Int_t ntrack; | |
252 | Int_t npart; | |
253 | Int_t nstmax; | |
254 | Int_t nvertx; | |
255 | Int_t nhead; | |
256 | Int_t nbit; | |
257 | } Gcnum_t; | |
258 | ||
259 | //----------GCCUTS | |
260 | // COMMON/GCCUTS/CUTGAM,CUTELE,CUTNEU,CUTHAD,CUTMUO,BCUTE,BCUTM | |
261 | // + ,DCUTE ,DCUTM ,PPCUTM,TOFMAX,GCUTS(5) | |
262 | typedef struct { | |
263 | Float_t cutgam; | |
264 | Float_t cutele; | |
265 | Float_t cutneu; | |
266 | Float_t cuthad; | |
267 | Float_t cutmuo; | |
268 | Float_t bcute; | |
269 | Float_t bcutm; | |
270 | Float_t dcute; | |
271 | Float_t dcutm; | |
272 | Float_t ppcutm; | |
273 | Float_t tofmax; | |
274 | Float_t gcuts[5]; | |
275 | } Gccuts_t; | |
276 | ||
9e7a32d1 | 277 | //----------GCMULO |
278 | // COMMON/GCMULO/SINMUL(101),COSMUL(101),SQRMUL(101),OMCMOL,CHCMOL | |
279 | // + ,EKMIN,EKMAX,NEKBIN,NEK1,EKINV,GEKA,GEKB,EKBIN(200),ELOW(200) | |
280 | typedef struct { | |
281 | Float_t sinmul[101]; | |
282 | Float_t cosmul[101]; | |
283 | Float_t sqrmul[101]; | |
284 | Float_t omcmol; | |
285 | Float_t chcmol; | |
286 | Float_t ekmin; | |
287 | Float_t ekmax; | |
288 | Int_t nekbin; | |
289 | Int_t nek1; | |
290 | Float_t ekinv; | |
291 | Float_t geka; | |
292 | Float_t gekb; | |
293 | Float_t ekbin[200]; | |
294 | Float_t elow[200]; | |
295 | } Gcmulo_t; | |
296 | ||
fe4da5cc | 297 | //----------GCPHYS |
298 | // COMMON/GCPHYS/IPAIR,SPAIR,SLPAIR,ZINTPA,STEPPA | |
299 | // + ,ICOMP,SCOMP,SLCOMP,ZINTCO,STEPCO | |
300 | // + ,IPHOT,SPHOT,SLPHOT,ZINTPH,STEPPH | |
301 | // + ,IPFIS,SPFIS,SLPFIS,ZINTPF,STEPPF | |
302 | // + ,IDRAY,SDRAY,SLDRAY,ZINTDR,STEPDR | |
303 | // + ,IANNI,SANNI,SLANNI,ZINTAN,STEPAN | |
304 | // + ,IBREM,SBREM,SLBREM,ZINTBR,STEPBR | |
305 | // + ,IHADR,SHADR,SLHADR,ZINTHA,STEPHA | |
306 | // + ,IMUNU,SMUNU,SLMUNU,ZINTMU,STEPMU | |
307 | // + ,IDCAY,SDCAY,SLIFE ,SUMLIF,DPHYS1 | |
308 | // + ,ILOSS,SLOSS,SOLOSS,STLOSS,DPHYS2 | |
309 | // + ,IMULS,SMULS,SOMULS,STMULS,DPHYS3 | |
310 | // + ,IRAYL,SRAYL,SLRAYL,ZINTRA,STEPRA | |
fe4da5cc | 311 | typedef struct { |
312 | Int_t ipair; | |
313 | Float_t spair; | |
314 | Float_t slpair; | |
315 | Float_t zintpa; | |
316 | Float_t steppa; | |
317 | Int_t icomp; | |
318 | Float_t scomp; | |
319 | Float_t slcomp; | |
320 | Float_t zintco; | |
321 | Float_t stepco; | |
322 | Int_t iphot; | |
323 | Float_t sphot; | |
324 | Float_t slphot; | |
325 | Float_t zintph; | |
326 | Float_t stepph; | |
327 | Int_t ipfis; | |
328 | Float_t spfis; | |
329 | Float_t slpfis; | |
330 | Float_t zintpf; | |
331 | Float_t steppf; | |
332 | Int_t idray; | |
333 | Float_t sdray; | |
334 | Float_t sldray; | |
335 | Float_t zintdr; | |
336 | Float_t stepdr; | |
337 | Int_t ianni; | |
338 | Float_t sanni; | |
339 | Float_t slanni; | |
340 | Float_t zintan; | |
341 | Float_t stepan; | |
342 | Int_t ibrem; | |
343 | Float_t sbrem; | |
344 | Float_t slbrem; | |
345 | Float_t zintbr; | |
346 | Float_t stepbr; | |
347 | Int_t ihadr; | |
348 | Float_t shadr; | |
349 | Float_t slhadr; | |
350 | Float_t zintha; | |
351 | Float_t stepha; | |
352 | Int_t imunu; | |
353 | Float_t smunu; | |
354 | Float_t slmunu; | |
355 | Float_t zintmu; | |
356 | Float_t stepmu; | |
357 | Int_t idcay; | |
358 | Float_t sdcay; | |
359 | Float_t slife; | |
360 | Float_t sumlif; | |
361 | Float_t dphys1; | |
362 | Int_t iloss; | |
363 | Float_t sloss; | |
364 | Float_t soloss; | |
365 | Float_t stloss; | |
366 | Float_t dphys2; | |
367 | Int_t imuls; | |
368 | Float_t smuls; | |
369 | Float_t somuls; | |
370 | Float_t stmuls; | |
371 | Float_t dphys3; | |
372 | Int_t irayl; | |
373 | Float_t srayl; | |
374 | Float_t slrayl; | |
375 | Float_t zintra; | |
376 | Float_t stepra; | |
377 | } Gcphys_t; | |
fd91b664 | 378 | |
379 | //----------GCPHLT | |
380 | // COMMON/GCPHLT/ILABS,SLABS,SLLABS,ZINTLA,STEPLA | |
381 | // + ,ISYNC | |
382 | // + ,ISTRA | |
383 | typedef struct { | |
384 | Int_t ilabs; | |
385 | Float_t slabs; | |
386 | Float_t sllabs; | |
387 | Float_t zintla; | |
388 | Float_t stepla; | |
389 | Int_t isync; | |
390 | Int_t istra; | |
391 | } Gcphlt_t; | |
fe4da5cc | 392 | |
393 | //----------GCOPTI | |
394 | // COMMON/GCOPTI/IOPTIM | |
395 | typedef struct { | |
396 | Int_t ioptim; | |
397 | } Gcopti_t; | |
398 | ||
399 | //----------GCTLIT | |
400 | // COMMON/GCTLIT/THRIND,PMIN,DP,DNDL,JMIN,ITCKOV,IMCKOV,NPCKOV | |
401 | typedef struct { | |
402 | Float_t thrind; | |
403 | Float_t pmin; | |
404 | Float_t dp; | |
405 | Float_t dndl; | |
406 | Int_t jmin; | |
407 | Int_t itckov; | |
408 | Int_t imckov; | |
409 | Int_t npckov; | |
410 | } Gctlit_t; | |
411 | ||
412 | //----------GCVDMA | |
413 | // COMMON/GCVDMA/NVMANY,MANYLE(20),MANYNA(20,15), | |
414 | // +MANYNU(20,15),NFMANY,MYCOUN,IMYSE,RAYTRA,VECCOS(3) | |
415 | typedef struct { | |
416 | Int_t vdma[624]; | |
417 | Float_t raytra; | |
418 | Float_t veccos[3]; | |
419 | } Gcvdma_t; | |
420 | ||
421 | //----------GCTPOL | |
422 | #define MAXME1 30 | |
423 | // COMMON/GCTPOL/POLAR(3), NAMEC1(MAXME1) | |
424 | typedef struct { | |
425 | Float_t polar[3]; | |
426 | Int_t namec1[MAXME1]; | |
427 | } Gctpol_t; | |
428 | ||
6991054d | 429 | /************************************************************************ |
430 | * * | |
431 | * Commons for GEANE * | |
432 | * * | |
433 | ************************************************************************/ | |
434 | ||
435 | //------------ERTRIO | |
436 | // INTEGER MXPRED | |
437 | // PARAMETER (MXPRED = 10) | |
438 | // DOUBLE PRECISION ERDTRP | |
439 | // REAL ERRIN, ERROUT, ERTRSP, ERXIN, ERXOUT, ERPIN, | |
440 | // + ERPOUT | |
441 | // INTEGER NEPRED, INLIST, ILPRED, IEPRED | |
442 | // COMMON /ERTRIO/ ERDTRP(5,5,MXPRED), ERRIN(15), ERROUT(15,MXPRED), | |
443 | // + ERTRSP(5,5,MXPRED), ERXIN( 3), ERXOUT( 3,MXPRED), | |
444 | // + ERPIN(3), ERPOUT(3,MXPRED), NEPRED,INLIST,ILPRED, | |
445 | // + IEPRED(MXPRED) | |
446 | // | |
447 | ||
448 | #define MXPRED 10 | |
449 | typedef struct { | |
450 | Double_t erdtrp[MXPRED*5*5]; | |
451 | Float_t errin[5]; | |
452 | Float_t errout[MXPRED*15]; | |
453 | Float_t ertrsp[MXPRED*5*5]; | |
454 | Float_t erxin[3]; | |
455 | Float_t erxout[MXPRED*3]; | |
456 | Float_t erpin[3]; | |
457 | Float_t erpout[MXPRED*3]; | |
458 | Int_t nepred; | |
459 | Int_t inlist; | |
460 | Int_t ilpred; | |
461 | Int_t iepred; | |
462 | } Ertrio_t; | |
463 | ||
464 | //-----------EROTPS | |
465 | // CHARACTER*8 CHOPTI | |
466 | // LOGICAL LEEXAC, LELENG, LEONLY, LEPLAN, LEPOIN, LEVOLU | |
467 | // REAL ERPLI, ERPLO, ERLENG | |
468 | // INTEGER NAMEER, NUMVER, IOVLER | |
469 | // COMMON /EROPTS/ ERPLI(3,2), ERPLO(3,4,MXPRED), ERLENG(MXPRED), | |
470 | // + NAMEER(MXPRED), NUMVER(MXPRED), IOVLER(MXPRED), | |
471 | // + LEEXAC, LELENG, LEONLY, LEPLAN, LEPOIN, LEVOLU | |
472 | // COMMON /EROPTC/CHOPTI | |
473 | ||
474 | typedef struct { | |
475 | Float_t erpli[3*2]; | |
476 | Float_t erplo[MXPRED*3*4]; | |
477 | Float_t erleng[MXPRED]; | |
478 | Int_t nameer[MXPRED]; | |
479 | Int_t numver[MXPRED]; | |
480 | Int_t iovler[MXPRED]; | |
481 | Bool_t leexac; | |
482 | Bool_t leleng; | |
483 | Bool_t leonly; | |
484 | Bool_t leplan; | |
485 | Bool_t lepoin; | |
486 | Bool_t levolu; | |
487 | } Eropts_t; | |
488 | ||
489 | typedef struct { | |
490 | char chopti[8]; | |
491 | } Eroptc_t; | |
492 | ||
493 | //-------ERWORK | |
494 | // DOUBLE PRECISION EI, EF, ASDSC | |
495 | // COMMON /ERWORK/ EI(15), EF(15), ASDSC(5,5), | |
496 | // + XI(3), PPI(3), HI(9), | |
497 | // + XF(3), PF(3), HF(9), | |
498 | // + CHTR, DEDX2, BACKTR, CUTEK, TLGCM2, TLRAD | |
499 | ||
500 | typedef struct { | |
501 | Double_t ei[15]; | |
502 | Double_t ef[15]; | |
503 | Double_t asdsc[5*5]; | |
504 | Float_t xi[3]; | |
505 | Float_t ppi[3]; | |
506 | Float_t hi[9]; | |
507 | Float_t xf[3]; | |
508 | Float_t pf[3]; | |
509 | Float_t hf[9]; | |
510 | Float_t chtr; | |
511 | Float_t dedx2; | |
512 | Float_t backtr; | |
513 | Float_t cutek; | |
514 | Float_t tlgcm2; | |
515 | Float_t tlrad; | |
516 | } Erwork_t; | |
517 | ||
518 | /************************************************************************ | |
519 | * * | |
520 | * Commons for GEANE * | |
521 | * * | |
522 | ************************************************************************/ | |
fe4da5cc | 523 | |
524 | class TGeant3 : public AliMC { | |
525 | ||
fe4da5cc | 526 | public: |
527 | TGeant3(); | |
528 | TGeant3(const char *title, Int_t nwgeant=0); | |
9006703b | 529 | virtual ~TGeant3() {if(fVolNames) { |
530 | delete [] fVolNames; | |
531 | fVolNames=0;} | |
532 | } | |
ef42d733 | 533 | |
fe4da5cc | 534 | virtual void LoadAddress(); |
535 | ||
536 | /////////////////////////////////////////////////////////////////////// | |
537 | // // | |
538 | // // | |
539 | // Here are the service routines from the geometry // | |
540 | // which could be implemented also in other geometries // | |
541 | // // | |
542 | // // | |
543 | /////////////////////////////////////////////////////////////////////// | |
544 | ||
545 | void GeomIter(); | |
546 | Int_t CurrentMaterial(Float_t &a, Float_t &z, Float_t &dens, Float_t &radl, Float_t &absl) const; | |
547 | Int_t NextVolUp(Text_t *name, Int_t ©); | |
0a6d8768 | 548 | Int_t CurrentVolID(Int_t ©) const; |
549 | Int_t CurrentVolOffID(Int_t off, Int_t ©) const; | |
550 | const char* CurrentVolName() const; | |
551 | const char *CurrentVolOffName(Int_t off) const; | |
229004c4 | 552 | Int_t VolId(const Text_t *name) const; |
1578254f | 553 | Int_t IdFromPDG(Int_t pdg) const; |
554 | Int_t PDGFromId(Int_t pdg) const; | |
555 | void DefineParticles(); | |
099385a4 | 556 | const char* VolName(Int_t id) const; |
ef42d733 | 557 | Float_t Xsec(char* reac, Float_t energy, Int_t part, Int_t mate); |
0a6d8768 | 558 | void TrackPosition(TLorentzVector &xyz) const; |
559 | void TrackMomentum(TLorentzVector &xyz) const; | |
1f97a957 | 560 | Int_t NofVolumes() const; |
fe4da5cc | 561 | Float_t TrackTime() const; |
562 | Float_t TrackCharge() const; | |
563 | Float_t TrackMass() const; | |
564 | Float_t TrackStep() const; | |
565 | Float_t TrackLength() const; | |
566 | Int_t TrackPid() const; | |
a01a8b12 | 567 | Bool_t IsNewTrack() const; |
0a6d8768 | 568 | Bool_t IsTrackInside() const; |
569 | Bool_t IsTrackEntering() const; | |
570 | Bool_t IsTrackExiting() const; | |
571 | Bool_t IsTrackOut() const; | |
572 | Bool_t IsTrackDisappeared() const; | |
573 | Bool_t IsTrackStop() const; | |
574 | Bool_t IsTrackAlive() const; | |
6a935c13 | 575 | Int_t NSecondaries() const; |
576 | Int_t CurrentEvent() const; | |
577 | const char* ProdProcess() const; | |
ef42d733 | 578 | void GetSecondary(Int_t isec, Int_t& ipart, TLorentzVector &x, |
579 | TLorentzVector &p); | |
fe4da5cc | 580 | void StopTrack(); |
581 | void StopEvent(); | |
582 | Float_t MaxStep() const; | |
fe4da5cc | 583 | void SetMaxStep(Float_t maxstep); |
584 | void SetMaxNStep(Int_t maxnstp); | |
585 | Int_t GetMaxNStep() const; | |
229004c4 | 586 | void SetCut(const char* cutName, Float_t cutValue); |
587 | void SetProcess(const char* flagName, Int_t flagValue); | |
42dcee26 | 588 | // void GetParticle(const Int_t pdg, char *name, Float_t &mass) const; |
fe4da5cc | 589 | virtual Int_t GetMedium() const; |
590 | virtual Float_t Edep() const; | |
591 | virtual Float_t Etot() const; | |
592 | virtual void Rndm(Float_t* r, const Int_t n) const; | |
ef42d733 | 593 | |
594 | virtual void Material(Int_t& kmat, const char* name, Float_t a, Float_t z, | |
595 | Float_t dens, Float_t radl, Float_t absl, | |
596 | Float_t* buf=0, Int_t nwbuf=0); | |
597 | ||
598 | virtual void Mixture(Int_t& kmat, const char* name, Float_t* a,Float_t* z, | |
599 | Float_t dens, Int_t nlmat, Float_t* wmat); | |
600 | ||
601 | virtual void Medium(Int_t& kmed, const char* name, Int_t nmat, Int_t isvol, | |
602 | Int_t ifield, Float_t fieldm, Float_t tmaxfd, | |
603 | Float_t stemax, Float_t deemax, Float_t epsil, | |
604 | Float_t stmin, Float_t* ubuf=0, Int_t nbuf=0); | |
605 | ||
606 | virtual void Matrix(Int_t& krot, Float_t thex, Float_t phix, Float_t they, | |
607 | Float_t phiy, Float_t thez, Float_t phiz); | |
fe4da5cc | 608 | |
609 | ///////////////////////////////////////////////////////////////////////////////////////////// | |
610 | // // | |
611 | // // | |
612 | // Here are the interface functions with GEANT3.21 // | |
613 | // // | |
614 | // // | |
615 | ///////////////////////////////////////////////////////////////////////////////////////////// | |
616 | ||
617 | // access functions to commons | |
618 | ||
619 | virtual Quest_t* Quest() const {return fQuest;} | |
7ac3f11b | 620 | virtual Gcbank_t* Gcbank() const {return fGcbank;} |
fe4da5cc | 621 | virtual Gclink_t* Gclink() const {return fGclink;} |
622 | virtual Gccuts_t* Gccuts() const {return fGccuts;} | |
9e7a32d1 | 623 | virtual Gcmulo_t* Gcmulo() const {return fGcmulo;} |
fe4da5cc | 624 | virtual Gcmate_t* Gcmate() const {return fGcmate;} |
625 | virtual Gctpol_t* Gctpol() const {return fGctpol;} | |
626 | virtual Gcnum_t* Gcnum() const {return fGcnum;} | |
627 | virtual Gcsets_t* Gcsets() const {return fGcsets;} | |
628 | virtual Gcopti_t* Gcopti() const {return fGcopti;} | |
629 | virtual Gctlit_t* Gctlit() const {return fGctlit;} | |
630 | virtual Gcvdma_t* Gcvdma() const {return fGcvdma;} | |
631 | virtual Gcvolu_t* Gcvolu() const {return fGcvolu;} | |
632 | virtual Gckine_t* Gckine() const {return fGckine;} | |
633 | virtual Gcflag_t* Gcflag() const {return fGcflag;} | |
634 | virtual Gctmed_t* Gctmed() const {return fGctmed;} | |
635 | virtual Gcphys_t* Gcphys() const {return fGcphys;} | |
fd91b664 | 636 | virtual Gcphlt_t* Gcphlt() const {return fGcphlt;} |
fe4da5cc | 637 | virtual Gcking_t* Gcking() const {return fGcking;} |
638 | virtual Gckin2_t* Gckin2() const {return fGckin2;} | |
639 | virtual Gckin3_t* Gckin3() const {return fGckin3;} | |
640 | virtual Gctrak_t* Gctrak() const {return fGctrak;} | |
24380615 | 641 | virtual Int_t* Iq() const {return fZiq;} |
642 | virtual Int_t* Lq() const {return fZlq;} | |
643 | virtual Float_t* Q() const {return fZq;} | |
fe4da5cc | 644 | |
645 | ||
6991054d | 646 | // Access to GEANE commons |
647 | ||
648 | virtual Ertrio_t* Ertrio() const {return fErtrio;} | |
649 | virtual Eropts_t* Eropts() const {return fEropts;} | |
650 | virtual Eroptc_t* Eroptc() const {return fEroptc;} | |
651 | virtual Erwork_t* Erwork() const {return fErwork;} | |
652 | ||
653 | ||
654 | ||
fe4da5cc | 655 | // functions from GBASE |
656 | virtual void Gpcxyz(); | |
657 | virtual void Ggclos(); | |
658 | virtual void Gfile(const char *filename, const char *option="I"); | |
659 | virtual void Glast(); | |
660 | virtual void Gprint(const char *name); | |
661 | virtual void Grun(); | |
662 | virtual void Gtrig(); | |
663 | virtual void Gtrigc(); | |
664 | virtual void Gtrigi(); | |
665 | virtual void Gwork(Int_t nwork); | |
666 | virtual void Gzinit(); | |
667 | ||
668 | // functions from GCONS | |
669 | virtual void Gfmate(Int_t imat, char *name, Float_t &a, Float_t &z, Float_t &dens, | |
670 | Float_t &radl, Float_t &absl, Float_t* ubuf, Int_t& nbuf); | |
671 | virtual void Gfpart(Int_t ipart, char *name, Int_t &itrtyp, | |
672 | Float_t &amass, Float_t &charge, Float_t &tlife); | |
673 | virtual void Gftmed(Int_t numed, char *name, Int_t &nmat, Int_t &isvol, | |
674 | Int_t &ifield, Float_t &fieldm, Float_t &tmaxfd, | |
675 | Float_t &stemax, Float_t &deemax, Float_t &epsil, | |
98be8990 | 676 | Float_t &stmin, Float_t *buf=0, Int_t *nbuf=0); |
677 | virtual void Gftmat(Int_t imate, Int_t ipart, char *chmeca, Int_t kdim, | |
678 | Float_t* tkin, Float_t* value, Float_t* pcut, | |
679 | Int_t &ixst); | |
680 | virtual Float_t Gbrelm(Float_t z, Float_t t, Float_t cut); | |
681 | virtual Float_t Gprelm(Float_t z, Float_t t, Float_t cut); | |
fe4da5cc | 682 | virtual void Gmate(); |
683 | virtual void Gpart(); | |
684 | virtual void Gsckov(Int_t itmed, Int_t npckov, Float_t *ppckov, | |
685 | Float_t *absco, Float_t *effic, Float_t *rindex); | |
686 | virtual void Gsdk(Int_t ipart, Float_t *bratio, Int_t *mode); | |
687 | virtual void Gsmate(Int_t imat, const char *name, Float_t a, Float_t z, | |
688 | Float_t dens, Float_t radl, Float_t absl); | |
689 | virtual void Gsmixt(Int_t imat, const char *name, Float_t *a, Float_t *z, | |
690 | Float_t dens, Int_t nlmat, Float_t *wmat); | |
691 | virtual void Gspart(Int_t ipart, const char *name, Int_t itrtyp, | |
692 | Float_t amass, Float_t charge, Float_t tlife); | |
693 | virtual void Gstmed(Int_t numed, const char *name, Int_t nmat, Int_t isvol, | |
694 | Int_t ifield, Float_t fieldm, Float_t tmaxfd, | |
695 | Float_t stemax, Float_t deemax, Float_t epsil, | |
696 | Float_t stmin); | |
697 | virtual void Gstpar(Int_t itmed, const char *param, Float_t parval); | |
698 | ||
699 | // functions from GKINE | |
700 | virtual void Gfkine(Int_t itra, Float_t *vert, Float_t *pvert, | |
701 | Int_t &ipart, Int_t &nvert); | |
702 | virtual void Gfvert(Int_t nvtx, Float_t *v, Int_t &ntbeam, Int_t &nttarg, Float_t &tofg); | |
703 | virtual Int_t Gskine(Float_t *plab, Int_t ipart, Int_t nv, Float_t *ubuf=0, Int_t nwbuf=0); | |
704 | virtual Int_t Gsvert(Float_t *v, Int_t ntbeam, Int_t nttarg, Float_t *ubuf=0, Int_t nwbuf=0); | |
705 | ||
706 | // functions from GPHYS | |
707 | virtual void Gphysi(); | |
708 | ||
709 | // functions from GTRAK | |
710 | virtual void Gdebug(); | |
711 | virtual void Gekbin(); | |
712 | virtual void Gfinds(); | |
713 | virtual void Gsking(Int_t igk); | |
714 | virtual void Gskpho(Int_t igk); | |
715 | virtual void Gsstak(Int_t iflag); | |
716 | virtual void Gsxyz(); | |
717 | virtual void Gtrack(); | |
718 | virtual void Gtreve(); | |
aee8290b | 719 | virtual void GtreveRoot(); |
fe4da5cc | 720 | virtual void Grndm(Float_t *rvec, const Int_t len) const; |
721 | virtual void Grndmq(Int_t &is1, Int_t &is2, const Int_t iseq, const Text_t *chopt); | |
722 | ||
723 | // functions from GGEOM | |
ef42d733 | 724 | virtual void Gdxyz(Int_t it); |
fe4da5cc | 725 | virtual void Gdcxyz(); |
726 | ||
727 | // functions from GGEOM | |
728 | virtual void Gdtom(Float_t *xd, Float_t *xm, Int_t iflag); | |
729 | virtual void Glmoth(const char* iudet, Int_t iunum, Int_t &nlev, | |
730 | Int_t *lvols, Int_t *lindx); | |
731 | virtual void Gmedia(Float_t *x, Int_t &numed); | |
732 | virtual void Gmtod(Float_t *xm, Float_t *xd, Int_t iflag); | |
733 | virtual void Gsdvn(const char *name, const char *mother, Int_t ndiv, Int_t iaxis); | |
734 | virtual void Gsdvn2(const char *name, const char *mother, Int_t ndiv, Int_t iaxis, Float_t c0i, Int_t numed); | |
735 | virtual void Gsdvs(const char *name, const char *mother, Float_t step, Int_t iaxis, Int_t numed); | |
736 | virtual void Gsdvs2(const char *name, const char *mother, Float_t step, Int_t iaxis, Float_t c0, Int_t numed); | |
737 | virtual void Gsdvt(const char *name, const char *mother, Float_t step, Int_t iaxis, Int_t numed, Int_t ndvmx); | |
738 | virtual void Gsdvt2(const char *name, const char *mother, Float_t step, Int_t iaxis, | |
739 | Float_t c0, Int_t numed, Int_t ndvmx); | |
740 | virtual void Gsord(const char *name, Int_t iax); | |
741 | virtual void Gspos(const char *name, Int_t nr, const char *mother, | |
742 | Float_t x, Float_t y, Float_t z, Int_t irot, const char *konly="ONLY"); | |
743 | virtual void Gsposp(const char *name, Int_t nr, const char *mother, | |
744 | Float_t x, Float_t y, Float_t z, Int_t irot, const char *konly, Float_t *upar, Int_t np); | |
745 | virtual void Gsrotm(Int_t nmat, Float_t theta1, Float_t phi1, Float_t theta2, Float_t phi2, | |
746 | Float_t theta3, Float_t phi3); | |
747 | virtual void Gprotm(Int_t nmat=0); | |
748 | virtual Int_t Gsvolu(const char *name, const char *shape, Int_t nmed, | |
749 | Float_t *upar, Int_t np); | |
750 | virtual void Gsatt(const char *name, const char *att, Int_t val); | |
751 | virtual void Gfpara(const char *name, Int_t number, Int_t intext, Int_t& npar, | |
752 | Int_t& natt, Float_t* par, Float_t* att); | |
ef42d733 | 753 | virtual void Gckpar(Int_t ish, Int_t npar, Float_t *par); |
754 | virtual void Gckmat(Int_t itmed, char *natmed); | |
fe4da5cc | 755 | |
756 | // functions from GDRAW | |
757 | virtual void DefaultRange(); | |
758 | virtual void InitHIGZ(); | |
759 | virtual void Gdopen(Int_t view); | |
760 | virtual void Gdclose(); | |
761 | virtual void Gdelete(Int_t view); | |
762 | virtual void Gdshow(Int_t view); | |
763 | virtual void Gdopt(const char *name,const char *value); | |
764 | virtual void Gdraw(const char *name,Float_t theta=30, Float_t phi=30, Float_t psi=0,Float_t u0=10,Float_t v0=10,Float_t ul=0.01,Float_t vl=0.01); | |
765 | virtual void Gdrawc(const char *name,Int_t axis=1, Float_t cut=0,Float_t u0=10,Float_t v0=10,Float_t ul=0.01,Float_t vl=0.01); | |
766 | virtual void Gdrawx(const char *name,Float_t cutthe, Float_t cutphi, Float_t cutval, | |
767 | Float_t theta=30, Float_t phi=30,Float_t u0=10,Float_t v0=10,Float_t ul=0.01,Float_t vl=0.01); | |
768 | virtual void Gdhead(Int_t isel, const char *name, Float_t chrsiz=0.6); | |
769 | virtual void Gdman(Float_t u0, Float_t v0, const char *type="MAN"); | |
770 | virtual void Gdspec(const char *name); | |
771 | virtual void DrawOneSpec(const char *name); | |
772 | virtual void Gdtree(const char *name,Int_t levmax=15,Int_t ispec=0); | |
773 | virtual void GdtreeParent(const char *name,Int_t levmax=15,Int_t ispec=0); | |
774 | ||
ef42d733 | 775 | virtual void WriteEuclid(const char* filnam, const char* topvol, |
776 | Int_t number, Int_t nlevel); | |
fe4da5cc | 777 | |
778 | virtual void SetABAN(Int_t par=1); | |
779 | virtual void SetANNI(Int_t par=1); | |
780 | virtual void SetAUTO(Int_t par=1); | |
781 | virtual void SetBOMB(Float_t bomb=1); | |
782 | virtual void SetBREM(Int_t par=1); | |
783 | virtual void SetCKOV(Int_t par=1); | |
784 | virtual void SetClipBox(const char *name,Float_t xmin=-9999,Float_t xmax=0, Float_t ymin=-9999,Float_t ymax=0,Float_t zmin=-9999,Float_t zmax=0); | |
785 | virtual void SetCOMP(Int_t par=1); | |
786 | virtual void SetCUTS(Float_t cutgam,Float_t cutele,Float_t cutneu,Float_t cuthad, | |
787 | Float_t cutmuo ,Float_t bcute ,Float_t bcutm ,Float_t dcute , | |
788 | Float_t dcutm ,Float_t ppcutm, Float_t tofmax); | |
789 | virtual void SetDCAY(Int_t par=1); | |
790 | virtual void SetDEBU(Int_t emin=1, Int_t emax=999, Int_t emod=1); | |
791 | virtual void SetDRAY(Int_t par=1); | |
9e7a32d1 | 792 | virtual void SetERAN(Float_t ekmin=1.e-5, Float_t ekmax=1.e4, |
793 | Int_t nekbin=90); | |
fe4da5cc | 794 | virtual void SetHADR(Int_t par=1); |
795 | virtual void SetKINE(Int_t kine, Float_t xk1=0, Float_t xk2=0, Float_t xk3=0, Float_t xk4=0, | |
796 | Float_t xk5=0, Float_t xk6=0, Float_t xk7=0, Float_t xk8=0, Float_t xk9=0, | |
797 | Float_t xk10=0); | |
798 | virtual void SetLOSS(Int_t par=2); | |
799 | virtual void SetMULS(Int_t par=1); | |
800 | virtual void SetMUNU(Int_t par=1); | |
801 | virtual void SetOPTI(Int_t par=2); | |
802 | virtual void SetPAIR(Int_t par=1); | |
803 | virtual void SetPFIS(Int_t par=1); | |
804 | virtual void SetPHOT(Int_t par=1); | |
805 | virtual void SetRAYL(Int_t par=1); | |
fd91b664 | 806 | virtual void SetSTRA(Int_t par=0); |
fe4da5cc | 807 | virtual void SetSWIT(Int_t sw, Int_t val=1); |
808 | virtual void SetTRIG(Int_t nevents=1); | |
7ac3f11b | 809 | virtual void SetUserDecay(Int_t ipart); |
fe4da5cc | 810 | |
811 | virtual void Vname(const char *name, char *vname); | |
812 | ||
813 | virtual void InitLego(); | |
6991054d | 814 | |
815 | // Routines from GEANE | |
816 | ||
817 | virtual void Ertrgo(); | |
818 | virtual void Ertrak(const Float_t *const x1, const Float_t *const p1, | |
819 | const Float_t *x2, const Float_t *p2, | |
820 | Int_t ipa, Option_t *chopt); | |
fe4da5cc | 821 | |
875c717b | 822 | // Control Methods |
823 | ||
824 | virtual void FinishGeometry(); | |
825 | virtual void BuildPhysics(); | |
826 | ||
ef42d733 | 827 | protected: |
828 | Int_t fNextVol; // Iterator for GeomIter | |
829 | ||
830 | //--------------Declarations for ZEBRA--------------------- | |
831 | Int_t *fZiq; // Good Old IQ of Zebra | |
832 | Int_t *fZlq; // Good Old LQ of Zebra | |
833 | Float_t *fZq; // Good Old Q of Zebra | |
834 | ||
835 | Quest_t *fQuest; // QUEST common structure | |
836 | Gcbank_t *fGcbank; // GCBANK common structure | |
837 | Gclink_t *fGclink; // GCLINK common structure | |
838 | Gccuts_t *fGccuts; // GCCUTS common structure | |
839 | Gcmulo_t *fGcmulo; // GCMULO common structure | |
840 | Gcmate_t *fGcmate; // GCMATE common structure | |
841 | Gctpol_t *fGctpol; // GCTPOL common structure | |
842 | Gcnum_t *fGcnum; // GCNUM common structure | |
843 | Gcsets_t *fGcsets; // GCSETS common structure | |
844 | Gcopti_t *fGcopti; // GCOPTI common structure | |
845 | Gctlit_t *fGctlit; // GCTLIT common structure | |
846 | Gcvdma_t *fGcvdma; // GCVDMA common structure | |
847 | Gcvolu_t *fGcvolu; // GCVOLU common structure | |
848 | Gckine_t *fGckine; // GCKINE common structure | |
849 | Gcflag_t *fGcflag; // GCFLAG common structure | |
850 | Gctmed_t *fGctmed; // GCTMED common structure | |
851 | Gcphys_t *fGcphys; // GCPHYS common structure | |
852 | Gcphlt_t *fGcphlt; // GCPHLT common structure | |
853 | Gcking_t *fGcking; // GCKING common structure | |
854 | Gckin2_t *fGckin2; // GCKIN2 common structure | |
855 | Gckin3_t *fGckin3; // GCKIN3 common structure | |
856 | Gctrak_t *fGctrak; // GCTRAK common structure | |
857 | ||
858 | ||
859 | // commons for GEANE | |
860 | Ertrio_t *fErtrio; // ERTRIO common structure | |
861 | Eropts_t *fEropts; // EROPTS common structure | |
862 | Eroptc_t *fEroptc; // EROPTC common structure | |
863 | Erwork_t *fErwork; // ERWORK common structure | |
864 | ||
865 | //Put here all volume names | |
866 | ||
867 | char (*fVolNames)[5]; //! Names of geant volumes as C++ chars | |
868 | ||
869 | enum {kMaxParticles = 100}; | |
870 | ||
871 | Int_t fNPDGCodes; // Number of PDG codes known by G3 | |
872 | ||
873 | Int_t fPDGCode[kMaxParticles]; // Translation table of PDG codes | |
874 | ||
875 | private: | |
876 | TGeant3(const TGeant3 &) {} | |
877 | TGeant3 & operator=(const TGeant3&) {return *this;} | |
878 | ||
879 | ||
fe4da5cc | 880 | ClassDef(TGeant3,1) //C++ interface to Geant basic routines |
881 | }; | |
882 | ||
883 | #endif |