Adding possibility to set ISTRA, PAI model
[u/mrichter/AliRoot.git] / TGeant3 / TGeant3.h
1 #ifndef TGeant3_H 
2 #define TGeant3_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 //  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  
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
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 //----------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
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
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; 
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;
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
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  ************************************************************************/
523
524 class TGeant3 : public AliMC { 
525
526 protected:
527   Int_t fNextVol;    // Iterator for GeomIter
528
529 //--------------Declarations for ZEBRA--------------------- 
530   Int_t *fZiq, *fZlq; 
531   Float_t *fZq; 
532
533   Quest_t  *fQuest; 
534   Gcbank_t *fGcbank;
535   Gclink_t *fGclink; 
536   Gccuts_t *fGccuts; 
537   Gcmulo_t *fGcmulo; 
538   Gcmate_t *fGcmate; 
539   Gctpol_t *fGctpol; 
540   Gcnum_t  *fGcnum; 
541   Gcsets_t *fGcsets; 
542   Gcopti_t *fGcopti; 
543   Gctlit_t *fGctlit; 
544   Gcvdma_t *fGcvdma; 
545   Gcvolu_t *fGcvolu; 
546   Gckine_t *fGckine; 
547   Gcflag_t *fGcflag; 
548   Gctmed_t *fGctmed; 
549   Gcphys_t *fGcphys; 
550   Gcphlt_t *fGcphlt; 
551   Gcking_t *fGcking; 
552   Gckin2_t *fGckin2; 
553   Gckin3_t *fGckin3; 
554   Gctrak_t *fGctrak; 
555
556
557   // commons for GEANE
558   Ertrio_t *fErtrio;
559   Eropts_t *fEropts;
560   Eroptc_t *fEroptc;
561   Erwork_t *fErwork;
562
563   //Put here all volume names
564
565   char (*fVolNames)[5];           //! Names of geant volumes as C++ chars
566
567   enum {kMaxParticles = 100};
568
569   Int_t fNPDGCodes;
570
571   Int_t fPDGCode[kMaxParticles];
572
573 public: 
574   TGeant3(); 
575   TGeant3(const char *title, Int_t nwgeant=0); 
576   virtual ~TGeant3() {if(fVolNames) {
577     delete [] fVolNames;
578     fVolNames=0;}
579   } 
580   virtual void LoadAddress(); 
581  
582 ///////////////////////////////////////////////////////////////////////
583 //                                                                   //
584 //                                                                   //
585 //     Here are the service routines from the geometry               //
586 //     which could be implemented also in other geometries           //
587 //                                                                   //
588 //                                                                   //
589 ///////////////////////////////////////////////////////////////////////
590
591   void  GeomIter();
592   Int_t CurrentMaterial(Float_t &a, Float_t &z, Float_t &dens, Float_t &radl, Float_t &absl) const;
593   Int_t NextVolUp(Text_t *name, Int_t &copy);
594   Int_t CurrentVolID(Int_t &copy) const;
595   Int_t CurrentVolOffID(Int_t off, Int_t &copy) const;
596   const char* CurrentVolName() const;
597   const char *CurrentVolOffName(Int_t off) const;
598   Int_t VolId(const Text_t *name) const;
599   Int_t IdFromPDG(Int_t pdg) const;
600   Int_t PDGFromId(Int_t pdg) const;
601   void  DefineParticles();
602   const char* VolName(Int_t id) const;
603   Float_t Xsec(char*, Float_t, Int_t, Int_t);
604   void  TrackPosition(TLorentzVector &xyz) const;
605   void  TrackMomentum(TLorentzVector &xyz) const;  
606   Int_t NofVolumes() const;
607   Float_t TrackTime() const;  
608   Float_t TrackCharge() const;
609   Float_t TrackMass() const;
610   Float_t TrackStep() const;
611   Float_t TrackLength() const;
612   Int_t   TrackPid() const;
613   Bool_t IsNewTrack() const;
614   Bool_t IsTrackInside() const;
615   Bool_t IsTrackEntering() const;
616   Bool_t IsTrackExiting() const;
617   Bool_t IsTrackOut() const;
618   Bool_t IsTrackDisappeared() const;
619   Bool_t IsTrackStop() const;
620   Bool_t IsTrackAlive() const;
621   Int_t  NSecondaries() const;
622   Int_t  CurrentEvent() const;
623   const char*  ProdProcess() const;
624   void   GetSecondary(Int_t, Int_t&, TLorentzVector&, TLorentzVector&);
625   void   StopTrack();
626   void   StopEvent();
627   Float_t MaxStep() const;
628   void  SetMaxStep(Float_t maxstep);
629   void  SetMaxNStep(Int_t maxnstp);
630   Int_t GetMaxNStep() const;
631   void SetCut(const char* cutName, Float_t cutValue);
632   void SetProcess(const char* flagName, Int_t flagValue);
633   //  void GetParticle(const Int_t pdg, char *name, Float_t &mass) const;
634   virtual Int_t GetMedium() const;
635   virtual Float_t Edep() const;
636   virtual Float_t Etot() const;
637   virtual void    Rndm(Float_t* r, const Int_t n) const;
638   virtual void    Material(Int_t&, const char*, Float_t, Float_t, Float_t, Float_t,
639                             Float_t, Float_t* buf=0, Int_t nwbuf=0);
640   virtual void    Mixture(Int_t&, const char*, Float_t*, Float_t*, Float_t, Int_t, Float_t*);
641   virtual void    Medium(Int_t&, const char*, Int_t, Int_t, Int_t, Float_t, Float_t, 
642                    Float_t, Float_t, Float_t, Float_t, Float_t* ubuf=0, Int_t nbuf=0);
643   virtual void    Matrix(Int_t&, Float_t, Float_t, Float_t, Float_t, Float_t, Float_t);
644
645 /////////////////////////////////////////////////////////////////////////////////////////////
646 //                                                                                         //
647 //                                                                                         //
648 //     Here are the interface functions with GEANT3.21                                     //
649 //                                                                                         //
650 //                                                                                         //
651 /////////////////////////////////////////////////////////////////////////////////////////////
652
653   // access functions to commons
654  
655   virtual Quest_t* Quest() const {return fQuest;}
656   virtual Gcbank_t* Gcbank() const {return fGcbank;}
657   virtual Gclink_t* Gclink() const {return fGclink;}
658   virtual Gccuts_t* Gccuts() const {return fGccuts;}
659   virtual Gcmulo_t* Gcmulo() const {return fGcmulo;}
660   virtual Gcmate_t* Gcmate() const {return fGcmate;}
661   virtual Gctpol_t* Gctpol() const {return fGctpol;}
662   virtual Gcnum_t* Gcnum() const {return fGcnum;}
663   virtual Gcsets_t* Gcsets() const {return fGcsets;}
664   virtual Gcopti_t* Gcopti() const {return fGcopti;}
665   virtual Gctlit_t* Gctlit() const {return fGctlit;}
666   virtual Gcvdma_t* Gcvdma() const {return fGcvdma;}
667   virtual Gcvolu_t* Gcvolu() const {return fGcvolu;}
668   virtual Gckine_t* Gckine() const {return fGckine;}
669   virtual Gcflag_t* Gcflag() const {return fGcflag;}
670   virtual Gctmed_t* Gctmed() const {return fGctmed;}
671   virtual Gcphys_t* Gcphys() const {return fGcphys;}
672   virtual Gcphlt_t* Gcphlt() const {return fGcphlt;}
673   virtual Gcking_t* Gcking() const {return fGcking;}
674   virtual Gckin2_t* Gckin2() const {return fGckin2;}
675   virtual Gckin3_t* Gckin3() const {return fGckin3;}
676   virtual Gctrak_t* Gctrak() const {return fGctrak;}
677   virtual Int_t* Iq() const {return fZiq;}
678   virtual Int_t* Lq() const {return fZlq;}
679   virtual Float_t* Q() const {return fZq;}
680
681
682   // Access to GEANE commons
683
684   virtual Ertrio_t* Ertrio() const {return fErtrio;}
685   virtual Eropts_t* Eropts() const {return fEropts;}
686   virtual Eroptc_t* Eroptc() const {return fEroptc;}
687   virtual Erwork_t* Erwork() const {return fErwork;}
688
689
690
691       // functions from GBASE 
692    virtual  void  Gpcxyz(); 
693    virtual  void  Ggclos(); 
694    virtual  void  Gfile(const char *filename, const char *option="I"); 
695    virtual  void  Glast(); 
696    virtual  void  Gprint(const char *name); 
697    virtual  void  Grun(); 
698    virtual  void  Gtrig(); 
699    virtual  void  Gtrigc(); 
700    virtual  void  Gtrigi(); 
701    virtual  void  Gwork(Int_t nwork); 
702    virtual  void  Gzinit(); 
703  
704       // functions from GCONS 
705    virtual  void  Gfmate(Int_t imat, char *name, Float_t &a, Float_t &z, Float_t &dens, 
706                          Float_t &radl, Float_t &absl, Float_t* ubuf, Int_t& nbuf); 
707    virtual  void  Gfpart(Int_t ipart, char *name, Int_t &itrtyp,  
708                          Float_t &amass, Float_t &charge, Float_t &tlife); 
709    virtual  void  Gftmed(Int_t numed, char *name, Int_t &nmat, Int_t &isvol,  
710                          Int_t &ifield, Float_t &fieldm, Float_t &tmaxfd, 
711                          Float_t &stemax, Float_t &deemax, Float_t &epsil, 
712                          Float_t &stmin, Float_t *buf=0, Int_t *nbuf=0);
713    virtual  void  Gftmat(Int_t imate, Int_t ipart, char *chmeca, Int_t kdim,
714                          Float_t* tkin, Float_t* value, Float_t* pcut,
715                          Int_t &ixst);
716    virtual  Float_t Gbrelm(Float_t z, Float_t t, Float_t cut);
717    virtual  Float_t Gprelm(Float_t z, Float_t t, Float_t cut);
718    virtual  void  Gmate(); 
719    virtual  void  Gpart(); 
720    virtual  void  Gsckov(Int_t itmed, Int_t npckov, Float_t *ppckov,
721                          Float_t *absco, Float_t *effic, Float_t *rindex); 
722    virtual  void  Gsdk(Int_t ipart, Float_t *bratio, Int_t *mode); 
723    virtual  void  Gsmate(Int_t imat, const char *name, Float_t a, Float_t z,  
724                          Float_t dens, Float_t radl, Float_t absl); 
725    virtual  void  Gsmixt(Int_t imat, const char *name, Float_t *a, Float_t *z,  
726                          Float_t dens, Int_t nlmat, Float_t *wmat); 
727    virtual  void  Gspart(Int_t ipart, const char *name, Int_t itrtyp,  
728                          Float_t amass, Float_t charge, Float_t tlife); 
729    virtual  void  Gstmed(Int_t numed, const char *name, Int_t nmat, Int_t isvol,  
730                          Int_t ifield, Float_t fieldm, Float_t tmaxfd, 
731                          Float_t stemax, Float_t deemax, Float_t epsil, 
732                          Float_t stmin); 
733    virtual  void  Gstpar(Int_t itmed, const char *param, Float_t parval); 
734  
735       // functions from GKINE 
736    virtual  void  Gfkine(Int_t itra, Float_t *vert, Float_t *pvert, 
737                          Int_t &ipart, Int_t &nvert); 
738    virtual  void  Gfvert(Int_t nvtx, Float_t *v, Int_t &ntbeam, Int_t &nttarg, Float_t &tofg); 
739    virtual  Int_t Gskine(Float_t *plab, Int_t ipart, Int_t nv, Float_t *ubuf=0, Int_t nwbuf=0); 
740    virtual  Int_t Gsvert(Float_t *v, Int_t ntbeam, Int_t nttarg, Float_t *ubuf=0, Int_t nwbuf=0); 
741  
742       // functions from GPHYS 
743    virtual  void  Gphysi(); 
744  
745       // functions from GTRAK 
746    virtual  void  Gdebug(); 
747    virtual  void  Gekbin(); 
748    virtual  void  Gfinds(); 
749    virtual  void  Gsking(Int_t igk); 
750    virtual  void  Gskpho(Int_t igk); 
751    virtual  void  Gsstak(Int_t iflag); 
752    virtual  void  Gsxyz(); 
753    virtual  void  Gtrack(); 
754    virtual  void  Gtreve(); 
755    virtual  void  Gtreve_root(); 
756    virtual  void  Grndm(Float_t *rvec, const Int_t len) const; 
757    virtual  void  Grndmq(Int_t &is1, Int_t &is2, const Int_t iseq, const Text_t *chopt); 
758  
759       // functions from GGEOM 
760    virtual  void  Gdxyz(Int_t ); 
761    virtual  void  Gdcxyz(); 
762
763       // functions from GGEOM 
764    virtual  void  Gdtom(Float_t *xd, Float_t *xm, Int_t iflag); 
765    virtual  void  Glmoth(const char* iudet, Int_t iunum, Int_t &nlev, 
766                          Int_t *lvols, Int_t *lindx); 
767    virtual  void  Gmedia(Float_t *x, Int_t &numed); 
768    virtual  void  Gmtod(Float_t *xm, Float_t *xd, Int_t iflag); 
769    virtual  void  Gsdvn(const char *name, const char *mother, Int_t ndiv, Int_t iaxis); 
770    virtual  void  Gsdvn2(const char *name, const char *mother, Int_t ndiv, Int_t iaxis, Float_t c0i, Int_t numed); 
771    virtual  void  Gsdvs(const char *name, const char *mother, Float_t step, Int_t iaxis, Int_t numed); 
772    virtual  void  Gsdvs2(const char *name, const char *mother, Float_t step, Int_t iaxis, Float_t c0, Int_t numed); 
773    virtual  void  Gsdvt(const char *name, const char *mother, Float_t step, Int_t iaxis, Int_t numed, Int_t ndvmx); 
774    virtual  void  Gsdvt2(const char *name, const char *mother, Float_t step, Int_t iaxis,
775                          Float_t c0, Int_t numed, Int_t ndvmx); 
776    virtual  void  Gsord(const char *name, Int_t iax); 
777    virtual  void  Gspos(const char *name, Int_t nr, const char *mother,  
778                          Float_t x, Float_t y, Float_t z, Int_t irot, const char *konly="ONLY"); 
779    virtual  void  Gsposp(const char *name, Int_t nr, const char *mother,  
780                          Float_t x, Float_t y, Float_t z, Int_t irot, const char *konly, Float_t *upar, Int_t np); 
781    virtual  void  Gsrotm(Int_t nmat, Float_t theta1, Float_t phi1, Float_t theta2, Float_t phi2, 
782                          Float_t theta3, Float_t phi3); 
783    virtual  void  Gprotm(Int_t nmat=0); 
784    virtual  Int_t Gsvolu(const char *name, const char *shape, Int_t nmed,  
785                          Float_t *upar, Int_t np); 
786    virtual  void  Gsatt(const char *name, const char *att, Int_t val);
787    virtual  void  Gfpara(const char *name, Int_t number, Int_t intext, Int_t& npar,
788                          Int_t& natt, Float_t* par, Float_t* att);
789    virtual  void  Gckpar(Int_t, Int_t, Float_t*);
790    virtual  void  Gckmat(Int_t, char*);
791     
792       // functions from GDRAW 
793    virtual  void  DefaultRange();
794    virtual  void  InitHIGZ();
795    virtual  void  Gdopen(Int_t view);
796    virtual  void  Gdclose();
797    virtual  void  Gdelete(Int_t view);
798    virtual  void  Gdshow(Int_t view);
799    virtual  void  Gdopt(const char *name,const char *value);
800    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);
801    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);
802    virtual  void  Gdrawx(const char *name,Float_t cutthe, Float_t cutphi, Float_t cutval,
803                          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);
804    virtual  void  Gdhead(Int_t isel, const char *name, Float_t chrsiz=0.6);   
805    virtual  void  Gdman(Float_t u0, Float_t v0, const char *type="MAN");
806    virtual  void  Gdspec(const char *name);
807    virtual  void  DrawOneSpec(const char *name);
808    virtual  void  Gdtree(const char *name,Int_t levmax=15,Int_t ispec=0);
809    virtual  void  GdtreeParent(const char *name,Int_t levmax=15,Int_t ispec=0);
810
811    virtual  void  WriteEuclid(const char*, const char*, Int_t, Int_t);
812
813    virtual  void  SetABAN(Int_t par=1);
814    virtual  void  SetANNI(Int_t par=1);
815    virtual  void  SetAUTO(Int_t par=1);
816    virtual  void  SetBOMB(Float_t bomb=1);
817    virtual  void  SetBREM(Int_t par=1);
818    virtual  void  SetCKOV(Int_t par=1);
819    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);
820    virtual  void  SetCOMP(Int_t par=1);
821    virtual  void  SetCUTS(Float_t cutgam,Float_t cutele,Float_t cutneu,Float_t cuthad,
822                       Float_t cutmuo ,Float_t bcute ,Float_t bcutm ,Float_t dcute ,
823                       Float_t dcutm ,Float_t ppcutm, Float_t tofmax);
824    virtual  void  SetDCAY(Int_t par=1);
825    virtual  void  SetDEBU(Int_t emin=1, Int_t emax=999, Int_t emod=1);
826    virtual  void  SetDRAY(Int_t par=1);
827    virtual  void  SetERAN(Float_t ekmin=1.e-5, Float_t ekmax=1.e4,
828                           Int_t nekbin=90);
829    virtual  void  SetHADR(Int_t par=1);
830    virtual  void  SetKINE(Int_t kine, Float_t xk1=0, Float_t xk2=0, Float_t xk3=0, Float_t xk4=0,
831                          Float_t xk5=0, Float_t xk6=0, Float_t xk7=0, Float_t xk8=0, Float_t xk9=0,
832                          Float_t xk10=0);
833    virtual  void  SetLOSS(Int_t par=2);
834    virtual  void  SetMULS(Int_t par=1);
835    virtual  void  SetMUNU(Int_t par=1);
836    virtual  void  SetOPTI(Int_t par=2);
837    virtual  void  SetPAIR(Int_t par=1);
838    virtual  void  SetPFIS(Int_t par=1);
839    virtual  void  SetPHOT(Int_t par=1);
840    virtual  void  SetRAYL(Int_t par=1);
841    virtual  void  SetSTRA(Int_t par=0);
842    virtual  void  SetSWIT(Int_t sw, Int_t val=1);
843    virtual  void  SetTRIG(Int_t nevents=1);
844    virtual  void  SetUserDecay(Int_t ipart);
845
846    virtual  void  Vname(const char *name, char *vname);
847
848    virtual  void  InitLego();
849
850   // Routines from GEANE
851
852     virtual void Ertrgo();
853     virtual void Ertrak(const Float_t *const x1, const Float_t *const p1, 
854                         const Float_t *x2, const Float_t *p2,
855                         Int_t ipa,  Option_t *chopt);
856         
857   // Control Methods
858
859   virtual void FinishGeometry();
860   virtual void BuildPhysics();
861
862    ClassDef(TGeant3,1)  //C++ interface to Geant basic routines 
863 }; 
864
865 #endif