3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 ////////////////////////////////////////////////
9 // C++ interface to Geant3 basic routines //
10 ////////////////////////////////////////////////
14 //______________________________________________________________
16 // Geant3 prototypes for commons
18 //______________________________________________________________
22 // COMMON/QUEST/IQUEST(100)
28 // COMMON/GCBANK/NZEBRA,GVERSN,ZVERSN,IXSTOR,IXDIV,IXCONS,FENDQ(16)
29 // + ,LMAIN,LR1,WS(KWBANK)
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
71 // COMMON/GCFLAG/IDEBUG,IDEMIN,IDEMAX,ITEST,IDRUN,IDEVT,IEORUN
72 // + ,IEOTRI,IEVENT,ISWIT(10),IFINIT(20),NEVENT,NRNDM(2)
90 // COMMON/GCKINE/IKINE,PKINE(10),ITRA,ISTAK,IVERT,IPART,ITRTYP
91 // + ,NAPART(5),AMASS,CHARGE,TLIFE,VERT(3),PVERT(4),IPAOLD
110 // COMMON/GCKING/KCASE,NGKINE,GKIN(5,MXGKIN),
111 // + TOFD(MXGKIN),IFLGK(MXGKIN)
116 Float_t gkin[MXGKIN][5];
117 Float_t tofd[MXGKIN];
122 // COMMON/GCKIN2/NGPHOT,XPHOT(11,MXPHOT)
126 Float_t xphot[MXPHOT][11];
130 // COMMON/GCKIN3/GPOS(3,MXGKIN)
132 Float_t gpos[MXGKIN][3];
136 // COMMON/GCMATE/NMAT,NAMATE(5),A,Z,DENS,RADL,ABSL
148 // COMMON/GCTMED/NUMED,NATMED(5),ISVOL,IFIELD,FIELDM,TMAXFD,STEMAX
149 // + ,DEEMAX,EPSIL,STMIN,CFIELD,PREC,IUPD,ISTPAR,NUMOLD
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
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)
224 Float_t gtran[15][3];
225 Float_t grmat[15][10];
231 // COMMON/GCSETS/IHSET,IHDET,ISET,IDET,IDTYPE,NVNAME,NUMBV(20)
243 // COMMON/GCNUM/NMATE ,NVOLUM,NROTM,NTMED,NTMULT,NTRACK,NPART
244 // + ,NSTMAX,NVERTX,NHEAD,NBIT
260 // COMMON/GCCUTS/CUTGAM,CUTELE,CUTNEU,CUTHAD,CUTMUO,BCUTE,BCUTM
261 // + ,DCUTE ,DCUTM ,PPCUTM,TOFMAX,GCUTS(5)
278 // COMMON/GCMULO/SINMUL(101),COSMUL(101),SQRMUL(101),OMCMOL,CHCMOL
279 // + ,EKMIN,EKMAX,NEKBIN,NEK1,EKINV,GEKA,GEKB,EKBIN(200),ELOW(200)
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
380 // COMMON/GCPHLT/ILABS,SLABS,SLLABS,ZINTLA,STEPLA
394 // COMMON/GCOPTI/IOPTIM
400 // COMMON/GCTLIT/THRIND,PMIN,DP,DNDL,JMIN,ITCKOV,IMCKOV,NPCKOV
413 // COMMON/GCVDMA/NVMANY,MANYLE(20),MANYNA(20,15),
414 // +MANYNU(20,15),NFMANY,MYCOUN,IMYSE,RAYTRA,VECCOS(3)
423 // COMMON/GCTPOL/POLAR(3), NAMEC1(MAXME1)
426 Int_t namec1[MAXME1];
429 /************************************************************************
431 * Commons for GEANE *
433 ************************************************************************/
437 // PARAMETER (MXPRED = 10)
438 // DOUBLE PRECISION ERDTRP
439 // REAL ERRIN, ERROUT, ERTRSP, ERXIN, ERXOUT, ERPIN,
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,
450 Double_t erdtrp[MXPRED*5*5];
452 Float_t errout[MXPRED*15];
453 Float_t ertrsp[MXPRED*5*5];
455 Float_t erxout[MXPRED*3];
457 Float_t erpout[MXPRED*3];
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
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];
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
518 /************************************************************************
520 * Commons for GEANE *
522 ************************************************************************/
525 class TGeant3 : public AliMC {
529 TGeant3(const char *title, Int_t nwgeant=0);
530 virtual ~TGeant3() {if(fVolNames) {
535 virtual void LoadAddress();
537 ///////////////////////////////////////////////////////////////////////
540 // Here are the service routines from the geometry //
541 // which could be implemented also in other geometries //
544 ///////////////////////////////////////////////////////////////////////
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 ©);
549 Int_t CurrentVolID(Int_t ©) const;
550 Int_t CurrentVolOffID(Int_t off, Int_t ©) 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,
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;}
597 virtual Int_t GetMedium() const;
598 virtual Float_t Edep() const;
599 virtual Float_t Etot() const;
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);
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);
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);
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);
616 /////////////////////////////////////////////////////////////////////////////////////////////
619 // Here are the interface functions with GEANT3.21 //
622 /////////////////////////////////////////////////////////////////////////////////////////////
624 // access functions to commons
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;}
653 // Access to GEANE commons
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;}
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);
669 virtual void Gtrig();
670 virtual void Gtrigc();
671 virtual void Gtrigi();
672 virtual void Gwork(Int_t nwork);
673 virtual void Gzinit();
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,
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,
704 virtual void Gstpar(Int_t itmed, const char *param, Float_t parval);
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);
715 // functions from GPHYS
716 virtual void Gphysi();
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);
732 // functions from GGEOM
733 virtual void Gdxyz(Int_t it);
734 virtual void Gdcxyz();
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);
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);
785 virtual void WriteEuclid(const char* filnam, const char* topvol,
786 Int_t number, Int_t nlevel);
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,
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,
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);
821 virtual void Vname(const char *name, char *vname);
823 virtual void InitLego();
825 // Routines from GEANE
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);
834 virtual void FinishGeometry();
835 virtual void BuildPhysics();
838 Int_t fNextVol; // Iterator for GeomIter
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
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
870 Ertrio_t *fErtrio; // ERTRIO common structure
871 Eropts_t *fEropts; // EROPTS common structure
872 Eroptc_t *fEroptc; // EROPTC common structure
873 Erwork_t *fErwork; // ERWORK common structure
875 //Put here all volume names
877 char (*fVolNames)[5]; //! Names of geant volumes as C++ chars
879 enum {kMaxParticles = 100};
881 Int_t fNPDGCodes; // Number of PDG codes known by G3
883 Int_t fPDGCode[kMaxParticles]; // Translation table of PDG codes
885 AliDecayer* fDecayer; // Pointer to decayer
887 AliMCProcess G3toVMC(Int_t iproc) const;
890 TGeant3(const TGeant3 &) {}
891 TGeant3 & operator=(const TGeant3&) {return *this;}
894 ClassDef(TGeant3,1) //C++ interface to Geant basic routines