Bad modules and chip treatment added (E. Fragiacomo)
[u/mrichter/AliRoot.git] / ITS / AliITSCalibrationSSD.h
1 #ifndef ALIITSCALIBRATIONSSD_H
2 #define ALIITSCALIBRATIONSSD_H
3 /* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice                               */ 
5 #include "AliITSCalibration.h"
6 #include "AliITSNoiseSSD.h"
7 #include "AliITSPedestalSSD.h"
8 #include "AliITSGainSSD.h"
9 #include "AliITSBadChannelsSSD.h"
10 #include "AliITSresponseSSD.h"
11 #include "TArrayF.h"
12 #include "TArrayI.h"
13
14 /* $Id$ */
15 //////////////////////////////////////////////
16 // Response class for SSD                   //
17 //                                          //
18 //////////////////////////////////////////////
19 class AliITSCalibrationSSD : public AliITSCalibration {
20
21  public:
22     AliITSCalibrationSSD();
23     AliITSCalibrationSSD(const char *dataType);
24     virtual ~AliITSCalibrationSSD();
25
26     virtual  void   SetNoiseParam(Double_t np, Double_t nn) {
27         // set noise par
28         fNoiseP=np; fNoiseN=nn;
29     }
30
31     virtual void    GetNoiseParam(Double_t &np, Double_t &nn) const {
32         // get noise par
33         np=fNoiseP; nn=fNoiseN;
34     }
35
36     void AddNoiseP(Int_t c, Float_t n) { fNoise->AddNoiseP(c,n);}       
37     TArrayF GetNoiseP() {return fNoise->GetNoiseP(); }
38     Float_t GetNoiseP(Int_t n) {return fNoise->GetNoiseP(n); }
39     void AddNoiseN(Int_t c, Float_t n) { fNoise->AddNoiseN(c,n);}
40     TArrayF GetNoiseN() {return fNoise->GetNoiseN(); }
41     Float_t GetNoiseN(Int_t n) {return fNoise->GetNoiseN(n); }
42     void SetNoise( AliITSNoiseSSD* noise) {fNoise=noise;}
43
44     void AddPedestalP(Int_t c, Float_t n) { fPedestal->AddPedestalP(c,n);}       
45     TArrayF GetPedestalP() {return fPedestal->GetPedestalP(); }
46     Float_t GetPedestalP(Int_t n) {return fPedestal->GetPedestalP(n); }
47     void AddPedestalN(Int_t c, Float_t n) { fPedestal->AddPedestalN(c,n);}
48     TArrayF GetPedestalN() {return fPedestal->GetPedestalN(); }
49     Float_t GetPedestalN(Int_t n) {return fPedestal->GetPedestalN(n); }
50     void SetPedestal( AliITSPedestalSSD* pedestal) {fPedestal=pedestal;}
51
52     void AddGainP(Int_t c, Float_t n) { fGain->AddGainP(c,n);}       
53     TArrayF GetGainP() {return fGain->GetGainP(); }
54     Float_t GetGainP(Int_t n) {return fGain->GetGainP(n); }
55     void AddGainN(Int_t c, Float_t n) { fGain->AddGainN(c,n);}
56     TArrayF GetGainN() {return fGain->GetGainN(); }
57     Float_t GetGainN(Int_t n) {return fGain->GetGainN(n); }
58     void SetGain( AliITSGainSSD* gain) {fGain=gain;}
59
60     void   SetBad() {
61       fIsBad = kTRUE;
62       for(Int_t i=0;i<fgkChipsPerModule;i++) fIsChipBad[i]=kTRUE;
63     }
64     virtual Bool_t IsBad() const { return fIsBad; }
65     void   SetChipBad(Int_t nChip) {
66       fIsChipBad[nChip] = kTRUE;
67     }
68     virtual Bool_t IsChipBad(Int_t nChip) const {
69       return fIsChipBad[nChip];
70     }
71     Int_t ChipsPerModule() const{return fgkChipsPerModule;} // Number of chips/module
72     Int_t ChannelsPerChip() const{ return fgkChannelsPerChip;}//Number of channels/chip
73
74     TArrayI GetBadPChannelsList() { return fBadChannels->GetBadPChannelsList(); } const
75     TArrayI GetBadNChannelsList() { return fBadChannels->GetBadNChannelsList(); } const
76     void SetBadChannels( AliITSBadChannelsSSD* badchannels) {fBadChannels=badchannels;}
77
78     void SetNoisePThreshold(Int_t threshold) { fNoisePThreshold = threshold;}
79     void AddNoisyPChannel(Int_t c, Int_t n) { fNoisyPChannelsList.AddAt(n,c);}
80     TArrayI GetNoisyPChannelsList() const {return fNoisyPChannelsList; }
81     void SetNoiseNThreshold(Int_t threshold) { fNoiseNThreshold = threshold;}
82     void AddNoisyNChannel(Int_t c, Int_t n) { fNoisyNChannelsList.AddAt(n,c);}
83     TArrayI GetNoisyNChannelsList() const {return fNoisyNChannelsList; }
84
85     void SetNDeadPChannelsList(Int_t n) { fDeadPChannelsList.Set(n); }
86     void AddDeadPChannel(Int_t c, Int_t n) { fDeadPChannelsList.AddAt(n,c);}
87     TArrayI GetDeadPChannelsList() const {return fDeadPChannelsList; }
88     void SetNDeadNChannelsList(Int_t n) { fDeadNChannelsList.Set(n); }
89     void AddDeadNChannel(Int_t c, Int_t n) { fDeadNChannelsList.AddAt(n,c);}
90     TArrayI GetDeadNChannelsList() const {return fDeadNChannelsList; }
91     //
92
93
94     virtual  void   SetNDetParam(Int_t npar) {
95         // set number of param
96         fNPar=npar;
97     }
98     virtual  void   SetDetParam(Double_t *par);
99
100     // Parameters options
101     virtual Int_t   NDetParam() const {
102         // number of param
103         return fNPar;
104     }
105     virtual void    GetDetParam(Double_t *dpar) const; 
106
107     virtual void    SetSigmaSpread(Double_t p1, Double_t p2) {
108         // Set sigmas of the charge spread function: Pside-Nside
109         // square of (microns)
110         fSigmaP=p1; fSigmaN=p2;
111     }
112
113     virtual void    SigmaSpread(Double_t &sP, Double_t &sN) const {
114       // Get sigmas for the charge spread 
115       sP=fSigmaP; sN=fSigmaN;
116     }
117
118     virtual void   SetThresholds(Double_t /* a */, Double_t /* b */)
119       {NotImplemented("SetThresholds");}
120     virtual void   Thresholds(Double_t & /* a */, Double_t & /* b */) const 
121       {NotImplemented("Thresholds");}
122   
123     virtual void SetParamOptions(const char *opt1, const char *opt2) {((AliITSresponseSSD*)fResponse)->SetParamOptions(opt1,opt2);}
124     virtual void GetParamOptions(char *opt1,char *opt2) const {((AliITSresponseSSD*)fResponse)->ParamOptions(opt1,opt2);}
125     virtual void SetADCpereV(Double_t a=120./24888.9) {((AliITSresponseSSD*)fResponse)->SetADCpereV(a);}
126     virtual Double_t GetDEvToADC(Double_t eV) const {return ((AliITSresponseSSD*)fResponse)->DEvToADC(eV);}
127     virtual Int_t IEvToADC(Double_t eV) const {return ((AliITSresponseSSD*)fResponse)->IEvToADC(eV);} 
128
129     virtual void SetKeVperADC(Double_t a=86.4/120.0) {((AliITSresponseSSD*)fResponse)->SetKeVperADC(a);}
130     virtual Double_t ADCToKeV(Double_t adc) const {return ((AliITSresponseSSD*)fResponse)->ADCToKeV(adc);}
131
132   virtual Double_t GetCouplingPR() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingPR();}
133   virtual Double_t GetCouplingPL() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingPL();}
134   virtual Double_t GetCouplingNR() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingNR();}
135   virtual Double_t GetCouplingNL() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingNL();}
136   virtual void SetCouplings(Double_t pr, Double_t pl, Double_t nr, Double_t nl) 
137     { ((AliITSresponseSSD*)fResponse)->SetCouplings(pr,pl,nr,nl);}
138
139   virtual Int_t GetZSThreshold() const {return ((AliITSresponseSSD*)fResponse)->GetZSThreshold();}
140   virtual void SetZSThreshold(Int_t zsth) 
141     { ((AliITSresponseSSD*)fResponse)->SetZSThreshold(zsth);}
142
143 protected:
144
145     static const Int_t fgkChipsPerModule  = 12;    // Number of chips/module
146     static const Int_t fgkChannelsPerChip = 128;   // Number of channels/chip
147
148     static const Double_t fgkNoiseNDefault; // default for fNoiseN
149     static const Double_t fgkNoisePDefault; // default for fNoiseP
150     static const Int_t fgkNParDefault; // default for fNPar
151     static const Double_t fgkSigmaPDefault; //default for fSigmaP
152     static const Double_t fgkSigmaNDefault; //default for fSigmaP
153
154     Int_t   fNPar;            // Number of detector param 
155     Double_t *fDetPar;         //[fNPar] Array of parameters
156
157     Double_t fNoiseP;          // Noise on Pside
158     Double_t fNoiseN;          // Noise on Nside
159     Double_t fSigmaP;          // Sigma charge spread on Pside
160     Double_t fSigmaN;          // Sigma charge spread on Nside
161     
162     AliITSNoiseSSD *fNoise;
163     AliITSPedestalSSD *fPedestal;
164     AliITSGainSSD *fGain;
165     AliITSBadChannelsSSD *fBadChannels;
166
167     Bool_t   fIsBad;                         // module is dead or alive ?
168     Bool_t   fIsChipBad[fgkChipsPerModule];  // chip is dead or alive ?
169
170     TArrayF fGainP;           // Gain for P side channels
171     TArrayF fGainN;           // Gain for N side channels
172
173     TArrayF fNoisP;           // Noise for P side channels
174     TArrayF fNoisN;           // Noise for N side channels
175
176     TArrayF fPedP;
177     TArrayF fPedN;
178
179     Float_t fNoisePThreshold;     // need to decide if channel is noisy  
180     TArrayI  fNoisyPChannelsList; // list of P side noisy channels
181     Float_t fNoiseNThreshold;     // need to decide if channel is noisy  
182     TArrayI  fNoisyNChannelsList; // list of N side noisy channels
183
184     TArrayI  fDeadNChannelsList;  // list of P side dead channels
185     TArrayI  fDeadPChannelsList;  // list of N side dead channels
186
187  private:
188     AliITSCalibrationSSD(const AliITSCalibrationSSD &source); // copy constructor
189     AliITSCalibrationSSD& operator=(const AliITSCalibrationSSD &source); // ass. op.
190
191     ClassDef(AliITSCalibrationSSD,3) //Response class for SSD
192 };
193 #endif