]>
Commit | Line | Data |
---|---|---|
5d12ce38 | 1 | #ifndef ALIMC_H |
2 | #define ALIMC_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 | ||
2057aecb | 8 | // |
9 | // This is the ALICE implementation of TVirtualMCApplication | |
10 | // for simulation with different implementations | |
11 | // of the Virtual MonteCarlo | |
12 | // | |
13 | ||
5d12ce38 | 14 | #include <TArrayF.h> |
15 | #include <TArrayI.h> | |
024a7e64 | 16 | #include <TMCProcess.h> |
5d12ce38 | 17 | #include <TVirtualMCApplication.h> |
18 | ||
d6276823 | 19 | class TParticle; |
5d12ce38 | 20 | class AliGenerator; |
21 | class AliMCQA; | |
22 | ||
23 | class AliMC : public TVirtualMCApplication { | |
24 | public: | |
25 | AliMC(); | |
26 | AliMC(const char *name, const char *title); | |
27 | AliMC(const AliMC &mc); | |
28 | virtual ~AliMC(); | |
29 | ||
d0d4a6b3 | 30 | AliMC& operator= (const AliMC &mc) { |
31 | // Assignment operator | |
32 | mc.Copy(*this); | |
33 | return *this; | |
34 | } | |
5d12ce38 | 35 | |
36 | // | |
37 | // MC Application | |
38 | // | |
39 | virtual void ConstructGeometry(); | |
40 | virtual void InitGeometry(); | |
41 | virtual void GeneratePrimaries(); | |
42 | virtual void BeginEvent(); | |
43 | virtual void BeginPrimary(); | |
44 | virtual void PreTrack(); | |
45 | virtual void Stepping(); | |
46 | virtual void PostTrack(); | |
47 | virtual void FinishPrimary(); | |
48 | virtual void FinishEvent(); | |
49 | virtual Double_t TrackingZmax() const {return fTrZmax;} | |
50 | virtual Double_t TrackingRmax() const {return fTrRmax;} | |
51 | virtual void Field(const Double_t* x, Double_t* b) const; | |
52 | virtual Int_t DetFromMate(Int_t i) const { return (*fImedia)[i];} | |
53 | // | |
54 | ||
55 | virtual AliGenerator* Generator() const {return fGenerator;} | |
56 | virtual void SetGenerator(AliGenerator *generator); | |
57 | virtual void ResetGenerator(AliGenerator *generator); | |
58 | ||
59 | // | |
60 | virtual void ReadTransPar(); | |
61 | virtual void MediaTable(); | |
62 | virtual void EnergySummary(); | |
63 | virtual void FinishRun(); | |
64 | void AddEnergyDeposit(Int_t id, Float_t edep) | |
65 | {fEventEnergy[id]+=edep;} | |
66 | virtual void ResetHits(); | |
67 | virtual void TrackingLimits( Float_t rmax=1.e10, Float_t zmax=1.e10) | |
0561efeb | 68 | {fTrRmax=rmax; fTrZmax=zmax;} |
69 | virtual void DecayLimits( Float_t rmin = -1., Float_t rmax = -1., Int_t pdg = 0) | |
70 | {fRDecayMin = rmin; fRDecayMax = rmax; fDecayPdg = pdg;} | |
71 | ||
5d12ce38 | 72 | virtual void Init(); |
73 | virtual void SetTransPar(const char *filename="$(ALICE_ROOT)/data/galice.cuts"); | |
0054628d | 74 | virtual void Browse(TBrowser *b); |
5d12ce38 | 75 | AliMCQA *GetMCQA() const {return fMCQA;} |
76 | //PH | |
77 | virtual void AddHit(Int_t id, Int_t track, Int_t *vol, Float_t *hits) const; | |
78 | virtual void AddDigit(Int_t id, Int_t *tracks, Int_t *digits) const; | |
4fc0efd3 | 79 | virtual void AddHitList(TCollection *hitList) {if (hitList) fHitLists->Add(hitList);} |
5d12ce38 | 80 | Int_t GetCurrentTrackNumber() const; |
81 | virtual void DumpPart (Int_t i) const; | |
82 | virtual void DumpPStack () const; | |
83 | TList* GetHitLists() const {return fHitLists ;} | |
84 | Int_t GetNtrack() const; | |
85 | virtual Int_t GetPrimary(Int_t track) const; | |
86 | TObjArray *Particles() const; | |
87 | TParticle *Particle(Int_t i) const; | |
88 | virtual void PushTrack(Int_t done, Int_t parent, Int_t pdg, | |
89 | Float_t *pmom, Float_t *vpos, Float_t *polar, | |
90 | Float_t tof, TMCProcess mech, Int_t &ntr, | |
2057aecb | 91 | Float_t weight = 1, Int_t is = 0) const; |
5d12ce38 | 92 | virtual void PushTrack(Int_t done, Int_t parent, Int_t pdg, |
93 | Double_t px, Double_t py, Double_t pz, Double_t e, | |
94 | Double_t vx, Double_t vy, Double_t vz, Double_t tof, | |
95 | Double_t polx, Double_t poly, Double_t polz, | |
96 | TMCProcess mech, Int_t &ntr, Float_t weight=1, | |
2057aecb | 97 | Int_t is = 0) const; |
98 | virtual void SetHighWaterMark(Int_t nt) const; | |
5d12ce38 | 99 | |
2057aecb | 100 | virtual void KeepTrack(Int_t itra) const; |
101 | virtual void FlagTrack(Int_t track) const; | |
102 | virtual void SetCurrentTrack(Int_t track) const; | |
5d12ce38 | 103 | // Track reference related |
104 | virtual void AddTrackReference(Int_t label); | |
105 | TClonesArray *TrackReferences() const {return fTrackReferences;} | |
106 | virtual void RemapTrackReferencesIDs(Int_t *map); //remaping track references MI | |
107 | virtual void ResetTrackReferences(); | |
0561efeb | 108 | virtual void FixParticleDecaytime(); |
5d12ce38 | 109 | |
110 | private: | |
6c4904c2 | 111 | void Copy (TObject &mc) const; |
5d12ce38 | 112 | AliGenerator *fGenerator; // Generator used in the MC |
113 | TArrayF fEventEnergy; //! Energy deposit for current event | |
114 | TArrayF fSummEnergy; //! Energy per event in each volume | |
115 | TArrayF fSum2Energy; //! Energy squared per event in each volume | |
116 | Float_t fTrRmax; // Maximum radius for tracking | |
117 | Float_t fTrZmax; // Maximu z for tracking | |
0561efeb | 118 | Float_t fRDecayMax; // Maximum radius for decay |
119 | Float_t fRDecayMin; // Minimum radius for decay | |
120 | Int_t fDecayPdg; // PDG code of particle with forced decay length | |
5d12ce38 | 121 | TArrayI *fImedia; //! Array of correspondence between media and detectors |
122 | TString fTransParName; // Name of the transport parameters file | |
123 | AliMCQA *fMCQA; // Pointer to MC Quality assurance class | |
124 | //PH | |
125 | TList *fHitLists; //! Lists of hits to be remapped by PurifyKine | |
126 | TClonesArray *fTrackReferences; //!list of track references - for one primary track only -MI | |
127 | ||
21bf7095 | 128 | ClassDef(AliMC,2) |
5d12ce38 | 129 | }; |
130 | ||
131 | ||
132 | #endif |