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