8de976fa66332b0eb7b8d8f8bb2fa20d07aa8f74
[u/mrichter/AliRoot.git] / ITS / AliITSsimulationSPD.h
1 #ifndef ALIITSSIMULATIONSPD_H
2 #define ALIITSSIMULATIONSPD_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice                               */
5
6 /*
7 $Id$
8 */
9 #include "AliITSresponseSPD.h"
10 #include "AliITSsegmentationSPD.h"
11 #include "AliITSsimulation.h"
12
13 class AliITSMapA2;
14 class AliITSpList;
15 class AliITSmodule;
16
17 //-------------------------------------------------------------------
18
19 class AliITSsimulationSPD : public AliITSsimulation {
20
21  public:        
22     AliITSsimulationSPD(); // Default constructor
23     // Standard constructor
24     AliITSsimulationSPD(AliITSsegmentation *seg, AliITSresponse *res);
25     ~AliITSsimulationSPD();// destructor
26     AliITSsimulationSPD(const AliITSsimulationSPD &source); // copy constructo
27     // assignment operator
28     AliITSsimulationSPD& operator=(const AliITSsimulationSPD &source);
29     // Initilizes the variables
30     void Init(AliITSsegmentationSPD *seg, AliITSresponseSPD *resp);
31
32     // Sum digitize module
33     // Create maps to build the lists of tracks for each summable digit
34     void InitSimulationModule(Int_t module,Int_t events);
35 //    // Add summable digits to module maps.
36 //    void AddSDigitsToModule(TClonesArray *sdig,Int_t mask);
37     // Digitize module from the sum of summable digits.
38     void FinishSDigitiseModule();
39     void SDigitiseModule(AliITSmodule *mod, Int_t dummy0,Int_t dummy1);
40     // digitize module. Also need to digitize modules with only noise.
41     void DigitiseModule(AliITSmodule *mod,Int_t dummy0, Int_t dummy1);
42     // sum digits to Digits.
43     void SDigitsToDigits(Int_t module,AliITSpList *pList);
44     // updates the Map of signal, adding the energy  (ene) released by
45     // the current track
46     void UpdateMapSignal(Int_t row,Int_t col,Int_t trk,Int_t hit,Int_t mod,
47                          Double_t ene,AliITSpList *pList);
48     // updates the Map of noise, adding the energy  (ene) give my noise
49     void UpdateMapNoise(Int_t row,Int_t col,Int_t mod,Double_t ene,
50                         AliITSpList *pList);
51     // Loops over all hits to produce Analog/floting point digits. This
52     // is also the first task in producing standard digits.
53     void HitsToAnalogDigits(AliITSmodule *mod,Int_t *frowpixel,
54                             Int_t *fcolpixel,Double_t *fenepixel,
55                             AliITSpList *pList);
56     //  Steering function to determine the digits associated to a given
57     // hit (hitpos)
58     // The digits are created by charge sharing (ChargeSharing) and by
59     // capacitive coupling (SetCoupling). At all the created digits is
60     // associated the track number of the hit (ntrack)
61     void HitToDigit(AliITSmodule *mod, Int_t hitpos,Int_t *frowpixel,
62                     Int_t *fcolpixel, Double_t *fenepixel,AliITSpList *pList);
63     //  Take into account the geometrical charge sharing when the track
64     //  crosses more than one pixel.
65     void ChargeSharing(Float_t x1l,Float_t z1l,Float_t x2l,Float_t z2l,
66                        Int_t c1,Int_t r1,Int_t c2,Int_t r2,Float_t etot,
67                        Int_t &npixel,Int_t *frowpixel,Int_t *fcolpixel,
68                        Double_t *fenepixel);
69     //  Take into account the coupling between adiacent pixels.
70     //  The parameters probcol and probrow are the fractions of the
71     //  signal in one pixel shared in the two adjacent pixels along
72     //  the column and row direction, respectively.
73     void SetCoupling(Int_t row,Int_t col,Int_t ntrack,Int_t idhit,Int_t module,
74                      AliITSpList *pList);
75     // The pixels are fired if the energy deposited inside them is above
76     // the threshold parameter ethr. Fired pixed are interpreted as digits
77     // and stored in the file digitfilename. One also needs to write out
78     // cases when there is only noise (nhits==0).
79     void CreateDigit(Int_t module,AliITSpList *pList);
80     //  Set the electronic noise and threshold non-uniformities to all the
81     //  pixels in a detector.
82     //  The parameter fSigma is the squared sum of the sigma due to noise
83     //  and the sigma of the threshold distribution among pixels.
84     void SetFluctuations(AliITSpList *pList,Int_t module);
85     //  Apply a mask to the SPD module. 1% of the pixel channels are
86     //  masked. When the database will be ready, the masked pixels
87     //  should be read from it.
88     void SetMask();
89     // Create Histograms
90     void CreateHistograms();
91     // Reset histograms for this detector
92     void ResetHistograms();
93     // Fills the Summable digits Tree
94     void WriteSDigits(AliITSpList *pList);
95     // Fills fMap2A from the pList of Summable digits
96     void FillMapFrompList(AliITSpList *pList);
97     // get hist array
98     TObjArray*  GetHistArray() {return fHis;}
99
100  private:
101     // Getters for data kept in fSegmentation and fResponse.
102     // Returns the Threshold in electrons
103     Double_t GetThreshold(){Float_t a=0.0,b=0.0;
104     ((AliITSresponseSPD*)fResponse)->Thresholds(a,b); return a;}
105     // Returns the threshold and rms noise.
106     void GetThresholds(Float_t &t,Float_t &s){
107     ((AliITSresponseSPD*)fResponse)->Thresholds(t,s);}
108     // Returns the couplings Columb and Row.
109     void GetCouplings(Float_t &cc,Float_t cr){
110         ((AliITSresponseSPD*)fResponse)->GetNoiseParam(cc,cr);}
111     // Returns the number of pixels in x
112     Int_t GetNPixelsX(){return ((AliITSsegmentationSPD*)fSegmentation)->Npx();}
113     // Returns the number of pixels in z
114     Int_t GetNPixelsZ(){return ((AliITSsegmentationSPD*)fSegmentation)->Npz();}
115
116  private:
117 /*    Float_t      fThresh;   //! Threshold from fResponse
118     Float_t      fSigma;    //! Noise from fResponse
119     Float_t      fCouplCol; //! Coupling along columns from fResponse
120     Float_t      fCouplRow; //! Coupling along rows from fResponse
121     Int_t        fNPixelsX; //! NPixelsX from fSegmentation
122     Int_t        fNPixelsZ; //! NPixelsZ from fSegmentation
123 */
124 //    AliITSpList  *fList;    //! Internal map and sdigits
125 //    Int_t         fModule;  //! Module number being simulated
126 //    Int_t         fEvent;   //! Event number being simulated
127     AliITSMapA2  *fMapA2;   //! MapA2 for Local internal use only
128     TObjArray    *fHis;     //! just in case for histogramming for Local
129                             // internal use only
130
131     ClassDef(AliITSsimulationSPD,1)  // Simulation of SPD clusters
132
133 };
134
135 #endif