]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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]; | |
117 | Int_t tofd[MXGKIN]; | |
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 | ||
277 | //----------GCPHYS | |
278 | // COMMON/GCPHYS/IPAIR,SPAIR,SLPAIR,ZINTPA,STEPPA | |
279 | // + ,ICOMP,SCOMP,SLCOMP,ZINTCO,STEPCO | |
280 | // + ,IPHOT,SPHOT,SLPHOT,ZINTPH,STEPPH | |
281 | // + ,IPFIS,SPFIS,SLPFIS,ZINTPF,STEPPF | |
282 | // + ,IDRAY,SDRAY,SLDRAY,ZINTDR,STEPDR | |
283 | // + ,IANNI,SANNI,SLANNI,ZINTAN,STEPAN | |
284 | // + ,IBREM,SBREM,SLBREM,ZINTBR,STEPBR | |
285 | // + ,IHADR,SHADR,SLHADR,ZINTHA,STEPHA | |
286 | // + ,IMUNU,SMUNU,SLMUNU,ZINTMU,STEPMU | |
287 | // + ,IDCAY,SDCAY,SLIFE ,SUMLIF,DPHYS1 | |
288 | // + ,ILOSS,SLOSS,SOLOSS,STLOSS,DPHYS2 | |
289 | // + ,IMULS,SMULS,SOMULS,STMULS,DPHYS3 | |
290 | // + ,IRAYL,SRAYL,SLRAYL,ZINTRA,STEPRA | |
291 | // COMMON/GCPHLT/ILABS,SLABS,SLLABS,ZINTLA,STEPLA | |
292 | // + ,ISYNC | |
293 | // + ,ISTRA | |
294 | typedef struct { | |
295 | Int_t ipair; | |
296 | Float_t spair; | |
297 | Float_t slpair; | |
298 | Float_t zintpa; | |
299 | Float_t steppa; | |
300 | Int_t icomp; | |
301 | Float_t scomp; | |
302 | Float_t slcomp; | |
303 | Float_t zintco; | |
304 | Float_t stepco; | |
305 | Int_t iphot; | |
306 | Float_t sphot; | |
307 | Float_t slphot; | |
308 | Float_t zintph; | |
309 | Float_t stepph; | |
310 | Int_t ipfis; | |
311 | Float_t spfis; | |
312 | Float_t slpfis; | |
313 | Float_t zintpf; | |
314 | Float_t steppf; | |
315 | Int_t idray; | |
316 | Float_t sdray; | |
317 | Float_t sldray; | |
318 | Float_t zintdr; | |
319 | Float_t stepdr; | |
320 | Int_t ianni; | |
321 | Float_t sanni; | |
322 | Float_t slanni; | |
323 | Float_t zintan; | |
324 | Float_t stepan; | |
325 | Int_t ibrem; | |
326 | Float_t sbrem; | |
327 | Float_t slbrem; | |
328 | Float_t zintbr; | |
329 | Float_t stepbr; | |
330 | Int_t ihadr; | |
331 | Float_t shadr; | |
332 | Float_t slhadr; | |
333 | Float_t zintha; | |
334 | Float_t stepha; | |
335 | Int_t imunu; | |
336 | Float_t smunu; | |
337 | Float_t slmunu; | |
338 | Float_t zintmu; | |
339 | Float_t stepmu; | |
340 | Int_t idcay; | |
341 | Float_t sdcay; | |
342 | Float_t slife; | |
343 | Float_t sumlif; | |
344 | Float_t dphys1; | |
345 | Int_t iloss; | |
346 | Float_t sloss; | |
347 | Float_t soloss; | |
348 | Float_t stloss; | |
349 | Float_t dphys2; | |
350 | Int_t imuls; | |
351 | Float_t smuls; | |
352 | Float_t somuls; | |
353 | Float_t stmuls; | |
354 | Float_t dphys3; | |
355 | Int_t irayl; | |
356 | Float_t srayl; | |
357 | Float_t slrayl; | |
358 | Float_t zintra; | |
359 | Float_t stepra; | |
360 | } Gcphys_t; | |
361 | ||
362 | //----------GCOPTI | |
363 | // COMMON/GCOPTI/IOPTIM | |
364 | typedef struct { | |
365 | Int_t ioptim; | |
366 | } Gcopti_t; | |
367 | ||
368 | //----------GCTLIT | |
369 | // COMMON/GCTLIT/THRIND,PMIN,DP,DNDL,JMIN,ITCKOV,IMCKOV,NPCKOV | |
370 | typedef struct { | |
371 | Float_t thrind; | |
372 | Float_t pmin; | |
373 | Float_t dp; | |
374 | Float_t dndl; | |
375 | Int_t jmin; | |
376 | Int_t itckov; | |
377 | Int_t imckov; | |
378 | Int_t npckov; | |
379 | } Gctlit_t; | |
380 | ||
381 | //----------GCVDMA | |
382 | // COMMON/GCVDMA/NVMANY,MANYLE(20),MANYNA(20,15), | |
383 | // +MANYNU(20,15),NFMANY,MYCOUN,IMYSE,RAYTRA,VECCOS(3) | |
384 | typedef struct { | |
385 | Int_t vdma[624]; | |
386 | Float_t raytra; | |
387 | Float_t veccos[3]; | |
388 | } Gcvdma_t; | |
389 | ||
390 | //----------GCTPOL | |
391 | #define MAXME1 30 | |
392 | // COMMON/GCTPOL/POLAR(3), NAMEC1(MAXME1) | |
393 | typedef struct { | |
394 | Float_t polar[3]; | |
395 | Int_t namec1[MAXME1]; | |
396 | } Gctpol_t; | |
397 | ||
6991054d | 398 | /************************************************************************ |
399 | * * | |
400 | * Commons for GEANE * | |
401 | * * | |
402 | ************************************************************************/ | |
403 | ||
404 | //------------ERTRIO | |
405 | // INTEGER MXPRED | |
406 | // PARAMETER (MXPRED = 10) | |
407 | // DOUBLE PRECISION ERDTRP | |
408 | // REAL ERRIN, ERROUT, ERTRSP, ERXIN, ERXOUT, ERPIN, | |
409 | // + ERPOUT | |
410 | // INTEGER NEPRED, INLIST, ILPRED, IEPRED | |
411 | // COMMON /ERTRIO/ ERDTRP(5,5,MXPRED), ERRIN(15), ERROUT(15,MXPRED), | |
412 | // + ERTRSP(5,5,MXPRED), ERXIN( 3), ERXOUT( 3,MXPRED), | |
413 | // + ERPIN(3), ERPOUT(3,MXPRED), NEPRED,INLIST,ILPRED, | |
414 | // + IEPRED(MXPRED) | |
415 | // | |
416 | ||
417 | #define MXPRED 10 | |
418 | typedef struct { | |
419 | Double_t erdtrp[MXPRED*5*5]; | |
420 | Float_t errin[5]; | |
421 | Float_t errout[MXPRED*15]; | |
422 | Float_t ertrsp[MXPRED*5*5]; | |
423 | Float_t erxin[3]; | |
424 | Float_t erxout[MXPRED*3]; | |
425 | Float_t erpin[3]; | |
426 | Float_t erpout[MXPRED*3]; | |
427 | Int_t nepred; | |
428 | Int_t inlist; | |
429 | Int_t ilpred; | |
430 | Int_t iepred; | |
431 | } Ertrio_t; | |
432 | ||
433 | //-----------EROTPS | |
434 | // CHARACTER*8 CHOPTI | |
435 | // LOGICAL LEEXAC, LELENG, LEONLY, LEPLAN, LEPOIN, LEVOLU | |
436 | // REAL ERPLI, ERPLO, ERLENG | |
437 | // INTEGER NAMEER, NUMVER, IOVLER | |
438 | // COMMON /EROPTS/ ERPLI(3,2), ERPLO(3,4,MXPRED), ERLENG(MXPRED), | |
439 | // + NAMEER(MXPRED), NUMVER(MXPRED), IOVLER(MXPRED), | |
440 | // + LEEXAC, LELENG, LEONLY, LEPLAN, LEPOIN, LEVOLU | |
441 | // COMMON /EROPTC/CHOPTI | |
442 | ||
443 | typedef struct { | |
444 | Float_t erpli[3*2]; | |
445 | Float_t erplo[MXPRED*3*4]; | |
446 | Float_t erleng[MXPRED]; | |
447 | Int_t nameer[MXPRED]; | |
448 | Int_t numver[MXPRED]; | |
449 | Int_t iovler[MXPRED]; | |
450 | Bool_t leexac; | |
451 | Bool_t leleng; | |
452 | Bool_t leonly; | |
453 | Bool_t leplan; | |
454 | Bool_t lepoin; | |
455 | Bool_t levolu; | |
456 | } Eropts_t; | |
457 | ||
458 | typedef struct { | |
459 | char chopti[8]; | |
460 | } Eroptc_t; | |
461 | ||
462 | //-------ERWORK | |
463 | // DOUBLE PRECISION EI, EF, ASDSC | |
464 | // COMMON /ERWORK/ EI(15), EF(15), ASDSC(5,5), | |
465 | // + XI(3), PPI(3), HI(9), | |
466 | // + XF(3), PF(3), HF(9), | |
467 | // + CHTR, DEDX2, BACKTR, CUTEK, TLGCM2, TLRAD | |
468 | ||
469 | typedef struct { | |
470 | Double_t ei[15]; | |
471 | Double_t ef[15]; | |
472 | Double_t asdsc[5*5]; | |
473 | Float_t xi[3]; | |
474 | Float_t ppi[3]; | |
475 | Float_t hi[9]; | |
476 | Float_t xf[3]; | |
477 | Float_t pf[3]; | |
478 | Float_t hf[9]; | |
479 | Float_t chtr; | |
480 | Float_t dedx2; | |
481 | Float_t backtr; | |
482 | Float_t cutek; | |
483 | Float_t tlgcm2; | |
484 | Float_t tlrad; | |
485 | } Erwork_t; | |
486 | ||
487 | /************************************************************************ | |
488 | * * | |
489 | * Commons for GEANE * | |
490 | * * | |
491 | ************************************************************************/ | |
fe4da5cc | 492 | |
493 | class TGeant3 : public AliMC { | |
494 | ||
495 | private: | |
496 | Int_t fNextVol; // Iterator for GeomIter | |
497 | ||
498 | //--------------Declarations for ZEBRA--------------------- | |
499 | Int_t *fZiq, *fZlq; | |
500 | Float_t *fZq; | |
501 | ||
502 | Quest_t *fQuest; | |
7ac3f11b | 503 | Gcbank_t *fGcbank; |
fe4da5cc | 504 | Gclink_t *fGclink; |
505 | Gccuts_t *fGccuts; | |
506 | Gcmate_t *fGcmate; | |
507 | Gctpol_t *fGctpol; | |
508 | Gcnum_t *fGcnum; | |
509 | Gcsets_t *fGcsets; | |
510 | Gcopti_t *fGcopti; | |
511 | Gctlit_t *fGctlit; | |
512 | Gcvdma_t *fGcvdma; | |
513 | Gcvolu_t *fGcvolu; | |
514 | Gckine_t *fGckine; | |
515 | Gcflag_t *fGcflag; | |
516 | Gctmed_t *fGctmed; | |
517 | Gcphys_t *fGcphys; | |
518 | Gcking_t *fGcking; | |
519 | Gckin2_t *fGckin2; | |
520 | Gckin3_t *fGckin3; | |
521 | Gctrak_t *fGctrak; | |
522 | ||
6991054d | 523 | |
524 | // commons for GEANE | |
525 | Ertrio_t *fErtrio; | |
526 | Eropts_t *fEropts; | |
527 | Eroptc_t *fEroptc; | |
528 | Erwork_t *fErwork; | |
529 | ||
1578254f | 530 | enum {kMaxParticles = 100}; |
fe4da5cc | 531 | |
1578254f | 532 | Int_t fNPDGCodes; |
533 | ||
534 | Int_t fPDGCode[kMaxParticles]; | |
fe4da5cc | 535 | |
536 | public: | |
537 | TGeant3(); | |
538 | TGeant3(const char *title, Int_t nwgeant=0); | |
539 | virtual ~TGeant3() {} | |
540 | virtual void LoadAddress(); | |
541 | ||
542 | /////////////////////////////////////////////////////////////////////// | |
543 | // // | |
544 | // // | |
545 | // Here are the service routines from the geometry // | |
546 | // which could be implemented also in other geometries // | |
547 | // // | |
548 | // // | |
549 | /////////////////////////////////////////////////////////////////////// | |
550 | ||
551 | void GeomIter(); | |
552 | Int_t CurrentMaterial(Float_t &a, Float_t &z, Float_t &dens, Float_t &radl, Float_t &absl) const; | |
553 | Int_t NextVolUp(Text_t *name, Int_t ©); | |
0a6d8768 | 554 | Int_t CurrentVolID(Int_t ©) const; |
555 | Int_t CurrentVolOffID(Int_t off, Int_t ©) const; | |
556 | const char* CurrentVolName() const; | |
557 | const char *CurrentVolOffName(Int_t off) const; | |
fe4da5cc | 558 | Int_t VolId(Text_t *name) const; |
1578254f | 559 | Int_t IdFromPDG(Int_t pdg) const; |
560 | Int_t PDGFromId(Int_t pdg) const; | |
561 | void DefineParticles(); | |
099385a4 | 562 | const char* VolName(Int_t id) const; |
d5a6b28b | 563 | Float_t Xsec(char*, Float_t, Int_t, Int_t); |
0a6d8768 | 564 | void TrackPosition(TLorentzVector &xyz) const; |
565 | void TrackMomentum(TLorentzVector &xyz) const; | |
1f97a957 | 566 | Int_t NofVolumes() const; |
fe4da5cc | 567 | Float_t TrackTime() const; |
568 | Float_t TrackCharge() const; | |
569 | Float_t TrackMass() const; | |
570 | Float_t TrackStep() const; | |
571 | Float_t TrackLength() const; | |
572 | Int_t TrackPid() const; | |
0a6d8768 | 573 | Bool_t IsTrackInside() const; |
574 | Bool_t IsTrackEntering() const; | |
575 | Bool_t IsTrackExiting() const; | |
576 | Bool_t IsTrackOut() const; | |
577 | Bool_t IsTrackDisappeared() const; | |
578 | Bool_t IsTrackStop() const; | |
579 | Bool_t IsTrackAlive() const; | |
fe4da5cc | 580 | Int_t NSecondaries() const; |
581 | Int_t CurrentEvent() const; | |
582 | void ProdProcess(char*) const; | |
583 | void GetSecondary(Int_t, Int_t&, Float_t*, Float_t*); | |
584 | void StopTrack(); | |
585 | void StopEvent(); | |
586 | Float_t MaxStep() const; | |
587 | void SetColors(); | |
588 | void SetMaxStep(Float_t maxstep); | |
589 | void SetMaxNStep(Int_t maxnstp); | |
590 | Int_t GetMaxNStep() const; | |
42dcee26 | 591 | // void GetParticle(const Int_t pdg, char *name, Float_t &mass) const; |
fe4da5cc | 592 | virtual Int_t GetMedium() const; |
593 | virtual Float_t Edep() const; | |
594 | virtual Float_t Etot() const; | |
595 | virtual void Rndm(Float_t* r, const Int_t n) const; | |
596 | virtual void Material(Int_t&, const char*, Float_t, Float_t, Float_t, Float_t, | |
597 | Float_t, Float_t* buf=0, Int_t nwbuf=0); | |
598 | virtual void Mixture(Int_t&, const char*, Float_t*, Float_t*, Float_t, Int_t, Float_t*); | |
599 | virtual void Medium(Int_t&, const char*, Int_t, Int_t, Int_t, Float_t, Float_t, | |
600 | Float_t, Float_t, Float_t, Float_t, Float_t* ubuf=0, Int_t nbuf=0); | |
601 | virtual void Matrix(Int_t&, Float_t, Float_t, Float_t, Float_t, Float_t, Float_t); | |
602 | ||
603 | ///////////////////////////////////////////////////////////////////////////////////////////// | |
604 | // // | |
605 | // // | |
606 | // Here are the interface functions with GEANT3.21 // | |
607 | // // | |
608 | // // | |
609 | ///////////////////////////////////////////////////////////////////////////////////////////// | |
610 | ||
611 | // access functions to commons | |
612 | ||
613 | virtual Quest_t* Quest() const {return fQuest;} | |
7ac3f11b | 614 | virtual Gcbank_t* Gcbank() const {return fGcbank;} |
fe4da5cc | 615 | virtual Gclink_t* Gclink() const {return fGclink;} |
616 | virtual Gccuts_t* Gccuts() const {return fGccuts;} | |
617 | virtual Gcmate_t* Gcmate() const {return fGcmate;} | |
618 | virtual Gctpol_t* Gctpol() const {return fGctpol;} | |
619 | virtual Gcnum_t* Gcnum() const {return fGcnum;} | |
620 | virtual Gcsets_t* Gcsets() const {return fGcsets;} | |
621 | virtual Gcopti_t* Gcopti() const {return fGcopti;} | |
622 | virtual Gctlit_t* Gctlit() const {return fGctlit;} | |
623 | virtual Gcvdma_t* Gcvdma() const {return fGcvdma;} | |
624 | virtual Gcvolu_t* Gcvolu() const {return fGcvolu;} | |
625 | virtual Gckine_t* Gckine() const {return fGckine;} | |
626 | virtual Gcflag_t* Gcflag() const {return fGcflag;} | |
627 | virtual Gctmed_t* Gctmed() const {return fGctmed;} | |
628 | virtual Gcphys_t* Gcphys() const {return fGcphys;} | |
629 | virtual Gcking_t* Gcking() const {return fGcking;} | |
630 | virtual Gckin2_t* Gckin2() const {return fGckin2;} | |
631 | virtual Gckin3_t* Gckin3() const {return fGckin3;} | |
632 | virtual Gctrak_t* Gctrak() const {return fGctrak;} | |
24380615 | 633 | virtual Int_t* Iq() const {return fZiq;} |
634 | virtual Int_t* Lq() const {return fZlq;} | |
635 | virtual Float_t* Q() const {return fZq;} | |
fe4da5cc | 636 | |
637 | ||
6991054d | 638 | // Access to GEANE commons |
639 | ||
640 | virtual Ertrio_t* Ertrio() const {return fErtrio;} | |
641 | virtual Eropts_t* Eropts() const {return fEropts;} | |
642 | virtual Eroptc_t* Eroptc() const {return fEroptc;} | |
643 | virtual Erwork_t* Erwork() const {return fErwork;} | |
644 | ||
645 | ||
646 | ||
fe4da5cc | 647 | // functions from GBASE |
648 | virtual void Gpcxyz(); | |
649 | virtual void Ggclos(); | |
650 | virtual void Gfile(const char *filename, const char *option="I"); | |
651 | virtual void Glast(); | |
652 | virtual void Gprint(const char *name); | |
653 | virtual void Grun(); | |
654 | virtual void Gtrig(); | |
655 | virtual void Gtrigc(); | |
656 | virtual void Gtrigi(); | |
657 | virtual void Gwork(Int_t nwork); | |
658 | virtual void Gzinit(); | |
659 | ||
660 | // functions from GCONS | |
661 | virtual void Gfmate(Int_t imat, char *name, Float_t &a, Float_t &z, Float_t &dens, | |
662 | Float_t &radl, Float_t &absl, Float_t* ubuf, Int_t& nbuf); | |
663 | virtual void Gfpart(Int_t ipart, char *name, Int_t &itrtyp, | |
664 | Float_t &amass, Float_t &charge, Float_t &tlife); | |
665 | virtual void Gftmed(Int_t numed, char *name, Int_t &nmat, Int_t &isvol, | |
666 | Int_t &ifield, Float_t &fieldm, Float_t &tmaxfd, | |
667 | Float_t &stemax, Float_t &deemax, Float_t &epsil, | |
668 | Float_t &stmin, Float_t *buf=0, Int_t *nbuf=0); | |
669 | virtual void Gmate(); | |
670 | virtual void Gpart(); | |
671 | virtual void Gsckov(Int_t itmed, Int_t npckov, Float_t *ppckov, | |
672 | Float_t *absco, Float_t *effic, Float_t *rindex); | |
673 | virtual void Gsdk(Int_t ipart, Float_t *bratio, Int_t *mode); | |
674 | virtual void Gsmate(Int_t imat, const char *name, Float_t a, Float_t z, | |
675 | Float_t dens, Float_t radl, Float_t absl); | |
676 | virtual void Gsmixt(Int_t imat, const char *name, Float_t *a, Float_t *z, | |
677 | Float_t dens, Int_t nlmat, Float_t *wmat); | |
678 | virtual void Gspart(Int_t ipart, const char *name, Int_t itrtyp, | |
679 | Float_t amass, Float_t charge, Float_t tlife); | |
680 | virtual void Gstmed(Int_t numed, const char *name, Int_t nmat, Int_t isvol, | |
681 | Int_t ifield, Float_t fieldm, Float_t tmaxfd, | |
682 | Float_t stemax, Float_t deemax, Float_t epsil, | |
683 | Float_t stmin); | |
684 | virtual void Gstpar(Int_t itmed, const char *param, Float_t parval); | |
685 | ||
686 | // functions from GKINE | |
687 | virtual void Gfkine(Int_t itra, Float_t *vert, Float_t *pvert, | |
688 | Int_t &ipart, Int_t &nvert); | |
689 | virtual void Gfvert(Int_t nvtx, Float_t *v, Int_t &ntbeam, Int_t &nttarg, Float_t &tofg); | |
690 | virtual Int_t Gskine(Float_t *plab, Int_t ipart, Int_t nv, Float_t *ubuf=0, Int_t nwbuf=0); | |
691 | virtual Int_t Gsvert(Float_t *v, Int_t ntbeam, Int_t nttarg, Float_t *ubuf=0, Int_t nwbuf=0); | |
692 | ||
693 | // functions from GPHYS | |
694 | virtual void Gphysi(); | |
695 | ||
696 | // functions from GTRAK | |
697 | virtual void Gdebug(); | |
698 | virtual void Gekbin(); | |
699 | virtual void Gfinds(); | |
700 | virtual void Gsking(Int_t igk); | |
701 | virtual void Gskpho(Int_t igk); | |
702 | virtual void Gsstak(Int_t iflag); | |
703 | virtual void Gsxyz(); | |
704 | virtual void Gtrack(); | |
705 | virtual void Gtreve(); | |
1578254f | 706 | virtual void Gtreve_root(); |
fe4da5cc | 707 | virtual void Grndm(Float_t *rvec, const Int_t len) const; |
708 | virtual void Grndmq(Int_t &is1, Int_t &is2, const Int_t iseq, const Text_t *chopt); | |
709 | ||
710 | // functions from GGEOM | |
711 | virtual void Gdxyz(Int_t ); | |
712 | virtual void Gdcxyz(); | |
713 | ||
714 | // functions from GGEOM | |
715 | virtual void Gdtom(Float_t *xd, Float_t *xm, Int_t iflag); | |
716 | virtual void Glmoth(const char* iudet, Int_t iunum, Int_t &nlev, | |
717 | Int_t *lvols, Int_t *lindx); | |
718 | virtual void Gmedia(Float_t *x, Int_t &numed); | |
719 | virtual void Gmtod(Float_t *xm, Float_t *xd, Int_t iflag); | |
720 | virtual void Gsdvn(const char *name, const char *mother, Int_t ndiv, Int_t iaxis); | |
721 | virtual void Gsdvn2(const char *name, const char *mother, Int_t ndiv, Int_t iaxis, Float_t c0i, Int_t numed); | |
722 | virtual void Gsdvs(const char *name, const char *mother, Float_t step, Int_t iaxis, Int_t numed); | |
723 | virtual void Gsdvs2(const char *name, const char *mother, Float_t step, Int_t iaxis, Float_t c0, Int_t numed); | |
724 | virtual void Gsdvt(const char *name, const char *mother, Float_t step, Int_t iaxis, Int_t numed, Int_t ndvmx); | |
725 | virtual void Gsdvt2(const char *name, const char *mother, Float_t step, Int_t iaxis, | |
726 | Float_t c0, Int_t numed, Int_t ndvmx); | |
727 | virtual void Gsord(const char *name, Int_t iax); | |
728 | virtual void Gspos(const char *name, Int_t nr, const char *mother, | |
729 | Float_t x, Float_t y, Float_t z, Int_t irot, const char *konly="ONLY"); | |
730 | virtual void Gsposp(const char *name, Int_t nr, const char *mother, | |
731 | Float_t x, Float_t y, Float_t z, Int_t irot, const char *konly, Float_t *upar, Int_t np); | |
732 | virtual void Gsrotm(Int_t nmat, Float_t theta1, Float_t phi1, Float_t theta2, Float_t phi2, | |
733 | Float_t theta3, Float_t phi3); | |
734 | virtual void Gprotm(Int_t nmat=0); | |
735 | virtual Int_t Gsvolu(const char *name, const char *shape, Int_t nmed, | |
736 | Float_t *upar, Int_t np); | |
737 | virtual void Gsatt(const char *name, const char *att, Int_t val); | |
738 | virtual void Gfpara(const char *name, Int_t number, Int_t intext, Int_t& npar, | |
739 | Int_t& natt, Float_t* par, Float_t* att); | |
740 | virtual void Gckpar(Int_t, Int_t, Float_t*); | |
741 | virtual void Gckmat(Int_t, char*); | |
742 | ||
743 | // functions from GDRAW | |
744 | virtual void DefaultRange(); | |
745 | virtual void InitHIGZ(); | |
746 | virtual void Gdopen(Int_t view); | |
747 | virtual void Gdclose(); | |
748 | virtual void Gdelete(Int_t view); | |
749 | virtual void Gdshow(Int_t view); | |
750 | virtual void Gdopt(const char *name,const char *value); | |
751 | 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); | |
752 | 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); | |
753 | virtual void Gdrawx(const char *name,Float_t cutthe, Float_t cutphi, Float_t cutval, | |
754 | 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); | |
755 | virtual void Gdhead(Int_t isel, const char *name, Float_t chrsiz=0.6); | |
756 | virtual void Gdman(Float_t u0, Float_t v0, const char *type="MAN"); | |
757 | virtual void Gdspec(const char *name); | |
758 | virtual void DrawOneSpec(const char *name); | |
759 | virtual void Gdtree(const char *name,Int_t levmax=15,Int_t ispec=0); | |
760 | virtual void GdtreeParent(const char *name,Int_t levmax=15,Int_t ispec=0); | |
761 | ||
762 | virtual void WriteEuclid(const char*, const char*, Int_t, Int_t); | |
763 | ||
764 | virtual void SetABAN(Int_t par=1); | |
765 | virtual void SetANNI(Int_t par=1); | |
766 | virtual void SetAUTO(Int_t par=1); | |
767 | virtual void SetBOMB(Float_t bomb=1); | |
768 | virtual void SetBREM(Int_t par=1); | |
769 | virtual void SetCKOV(Int_t par=1); | |
770 | 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); | |
771 | virtual void SetCOMP(Int_t par=1); | |
772 | virtual void SetCUTS(Float_t cutgam,Float_t cutele,Float_t cutneu,Float_t cuthad, | |
773 | Float_t cutmuo ,Float_t bcute ,Float_t bcutm ,Float_t dcute , | |
774 | Float_t dcutm ,Float_t ppcutm, Float_t tofmax); | |
775 | virtual void SetDCAY(Int_t par=1); | |
776 | virtual void SetDEBU(Int_t emin=1, Int_t emax=999, Int_t emod=1); | |
777 | virtual void SetDRAY(Int_t par=1); | |
778 | virtual void SetHADR(Int_t par=1); | |
779 | virtual void SetKINE(Int_t kine, Float_t xk1=0, Float_t xk2=0, Float_t xk3=0, Float_t xk4=0, | |
780 | Float_t xk5=0, Float_t xk6=0, Float_t xk7=0, Float_t xk8=0, Float_t xk9=0, | |
781 | Float_t xk10=0); | |
782 | virtual void SetLOSS(Int_t par=2); | |
783 | virtual void SetMULS(Int_t par=1); | |
784 | virtual void SetMUNU(Int_t par=1); | |
785 | virtual void SetOPTI(Int_t par=2); | |
786 | virtual void SetPAIR(Int_t par=1); | |
787 | virtual void SetPFIS(Int_t par=1); | |
788 | virtual void SetPHOT(Int_t par=1); | |
789 | virtual void SetRAYL(Int_t par=1); | |
790 | virtual void SetSWIT(Int_t sw, Int_t val=1); | |
791 | virtual void SetTRIG(Int_t nevents=1); | |
7ac3f11b | 792 | virtual void SetUserDecay(Int_t ipart); |
fe4da5cc | 793 | |
794 | virtual void Vname(const char *name, char *vname); | |
795 | ||
796 | virtual void InitLego(); | |
6991054d | 797 | |
798 | // Routines from GEANE | |
799 | ||
800 | virtual void Ertrgo(); | |
801 | virtual void Ertrak(const Float_t *const x1, const Float_t *const p1, | |
802 | const Float_t *x2, const Float_t *p2, | |
803 | Int_t ipa, Option_t *chopt); | |
fe4da5cc | 804 | |
805 | ClassDef(TGeant3,1) //C++ interface to Geant basic routines | |
806 | }; | |
807 | ||
808 | #endif |