3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
9 // This is the ALICE implementation of TVirtualMCApplication
10 // for simulation with different implementations
11 // of the Virtual MonteCarlo
17 #include <TMCProcess.h>
18 #include <TVirtualMCApplication.h>
24 class AliMC : public TVirtualMCApplication {
27 AliMC(const char *name, const char *title);
28 AliMC(const AliMC &mc);
31 AliMC& operator= (const AliMC &mc) {
32 // Assignment operator
40 virtual void ConstructGeometry();
41 virtual Bool_t MisalignGeometry();
42 virtual void ConstructOpGeometry();
43 virtual void InitGeometry();
44 virtual void SetAllAlignableVolumes();
45 virtual void GeneratePrimaries();
46 virtual void BeginEvent();
47 virtual void BeginPrimary();
48 virtual void PreTrack();
49 virtual void Stepping();
50 virtual void PostTrack();
51 virtual void FinishPrimary();
52 virtual void FinishEvent();
53 virtual Double_t TrackingZmax() const {return fTrZmax;}
54 virtual Double_t TrackingRmax() const {return fTrRmax;}
55 virtual void Field(const Double_t* x, Double_t* b) const;
56 virtual Int_t DetFromMate(Int_t i) const { return (*fImedia)[i];}
59 virtual AliGenerator* Generator() const {return fGenerator;}
60 virtual void SetGenerator(AliGenerator *generator);
61 virtual void ResetGenerator(AliGenerator *generator);
64 virtual void ReadTransPar();
65 virtual void MediaTable();
66 virtual void EnergySummary();
67 virtual void FinishRun();
68 void AddEnergyDeposit(Int_t id, Float_t edep)
69 {fEventEnergy[id]+=edep;}
70 virtual void ResetHits();
71 virtual void TrackingLimits( Float_t rmax=1.e10, Float_t zmax=1.e10)
72 {fTrRmax=rmax; fTrZmax=zmax;}
73 virtual void DecayLimits( Float_t rmin = -1., Float_t rmax = -1., Int_t pdg = 0)
74 {fRDecayMin = rmin; fRDecayMax = rmax; fDecayPdg = pdg;}
77 virtual void SetTransPar(const char *filename="$(ALICE_ROOT)/data/galice.cuts");
78 virtual void Browse(TBrowser *b);
79 AliMCQA *GetMCQA() const {return fMCQA;}
81 virtual void AddHit(Int_t id, Int_t track, Int_t *vol, Float_t *hits) const;
82 virtual void AddDigit(Int_t id, Int_t *tracks, Int_t *digits) const;
83 virtual void AddHitList(TCollection *hitList) {if (hitList) fHitLists->Add(hitList);}
84 Int_t GetCurrentTrackNumber() const;
85 virtual void DumpPart (Int_t i) const;
86 virtual void DumpPStack () const;
87 TList* GetHitLists() const {return fHitLists ;}
88 Int_t GetNtrack() const;
89 virtual Int_t GetPrimary(Int_t track) const;
90 TObjArray *Particles() const;
91 TParticle *Particle(Int_t i) const;
92 virtual void PushTrack(Int_t done, Int_t parent, Int_t pdg,
93 Float_t *pmom, Float_t *vpos, Float_t *polar,
94 Float_t tof, TMCProcess mech, Int_t &ntr,
95 Float_t weight = 1, Int_t is = 0) const;
96 virtual void PushTrack(Int_t done, Int_t parent, Int_t pdg,
97 Double_t px, Double_t py, Double_t pz, Double_t e,
98 Double_t vx, Double_t vy, Double_t vz, Double_t tof,
99 Double_t polx, Double_t poly, Double_t polz,
100 TMCProcess mech, Int_t &ntr, Float_t weight=1,
102 virtual void SetHighWaterMark(Int_t nt) const;
104 virtual void KeepTrack(Int_t itra) const;
105 virtual void FlagTrack(Int_t track) const;
106 virtual void SetCurrentTrack(Int_t track) const;
107 // Track reference related
108 virtual void AddTrackReference(Int_t label);
109 TClonesArray *TrackReferences() const {return fTrackReferences;}
110 virtual void RemapTrackReferencesIDs(Int_t *map); //remaping track references MI
111 virtual void ResetTrackReferences();
112 virtual void FixParticleDecaytime();
115 void Copy (TObject &mc) const;
117 AliGenerator *fGenerator; // Generator used in the MC
118 TArrayF fEventEnergy; //! Energy deposit for current event
119 TArrayF fSummEnergy; //! Energy per event in each volume
120 TArrayF fSum2Energy; //! Energy squared per event in each volume
121 Float_t fTrRmax; // Maximum radius for tracking
122 Float_t fTrZmax; // Maximu z for tracking
123 Float_t fRDecayMax; // Maximum radius for decay
124 Float_t fRDecayMin; // Minimum radius for decay
125 Int_t fDecayPdg; // PDG code of particle with forced decay length
126 TArrayI *fImedia; //! Array of correspondence between media and detectors
127 TString fTransParName; // Name of the transport parameters file
128 AliMCQA *fMCQA; // Pointer to MC Quality assurance class
130 TList *fHitLists; //! Lists of hits to be remapped by PurifyKine
131 TClonesArray *fTrackReferences; //!list of track references - for one primary track only -MI