time offset and charge conversion to energy unit moved from Cluster finder V2 to...
[u/mrichter/AliRoot.git] / ITS / AliITSCalibrationSDD.h
1 #ifndef ALIITSCALIBRATIONSDD_H
2 #define ALIITSCALIBRATIONSDD_H
3  
4 /* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6  
7 /* $Id$ */
8
9 #include "AliITSCalibration.h"
10 #include "AliITSresponseSDD.h"
11 #include "TArrayI.h"
12
13 class AliITSMapSDD;
14 class AliITSresponseSDD;
15 ///////////////////////////////////////////////////////
16 //  Response for SDD                                 //
17 ///////////////////////////////////////////////////////
18
19 class AliITSCalibrationSDD : public AliITSCalibration {
20   public:
21     //
22     // Configuration methods
23     //
24     AliITSCalibrationSDD();
25     AliITSCalibrationSDD(const char *dataType);
26     virtual ~AliITSCalibrationSDD();
27     virtual void  SetNoiseParam(Double_t /*n*/, Double_t /*b*/){
28       NotImplemented("SetNoiseParam");}
29  
30     virtual void  GetNoiseParam(Double_t &/*n*/, Double_t &/*b*/) const {
31       NotImplemented("GetNoiseParam");}
32
33     virtual Float_t GetBaseline(Int_t anode) const {return fBaseline[anode];}
34     virtual void SetBaseline(Int_t anode,Double_t bas) {fBaseline[anode]=bas;}
35     virtual Float_t GetNoise(Int_t anode) const {return fNoise[anode];}
36     virtual void SetNoise(Int_t anode, Double_t noise) {fNoise[anode]=noise;}
37
38     virtual void   SetThresholds(Double_t  mv, Double_t /* b */){
39        // Min value used in 2D - could be used as a threshold setting
40         fMinVal = mv;}
41     virtual void   Thresholds(Double_t &  mv, Double_t & /* b */) const 
42       {mv = fMinVal;}
43     virtual void  GiveCompressParam(Int_t *x,Int_t ian) const;
44
45     void  SetNoiseAfterElectronics(Int_t anode,Double_t n=2.38){
46         // Noise after electronics (ADC units)
47         // 2.36 for ALICE from beam test measurements 2001
48         fNoiseAfterEl[anode]=n;}
49     Float_t  GetNoiseAfterElectronics(Int_t anode) const {
50         // Noise after electronics (ADC units)
51         return fNoiseAfterEl[anode];} 
52     //void SetDeadChannels(Int_t nchips=0, Int_t nchannels=0);
53     void SetDeadChannels(Int_t ndead=0){fDeadChannels=ndead; fBadChannels.Set(ndead);}
54     Int_t GetDeadChips() const { return fDeadChips; }
55     Int_t GetDeadChannels() const { return fDeadChannels; }
56     Float_t Gain(Int_t wing,Int_t chip,Int_t ch)const 
57         {return fGain[wing][chip][ch]; }
58     virtual void SetGain(Double_t g,Int_t wing,Int_t chip, Int_t ch) 
59       {fGain[wing][chip][ch]=g;}
60     
61     void    PrintGains() const;
62     void    Print();
63     virtual void Print(ostream *os) const {AliITSCalibrationSDD::Print(os);}
64     virtual void Print(Option_t *option="") const {AliITSCalibrationSDD::Print(option);}
65     // not implemented virtual methods (devlared in the mother class
66     virtual  void   SetDetParam(Double_t *) 
67       {NotImplemented("SetDetParam");}
68     virtual void   GetDetParam(Double_t *) const 
69       {NotImplemented("GetDetParam");}
70     virtual  void   SetNDetParam(Int_t /* n */)
71       {NotImplemented("SetNDetParam");}
72     virtual Int_t  NDetParam() const
73       {NotImplemented("NDetParam"); return 0;}
74     virtual void    SetSigmaSpread(Double_t, Double_t) 
75       {NotImplemented("SetSigmaSpread");}
76     virtual void    SigmaSpread(Double_t & /* p1 */,Double_t & /* p2 */) const 
77       {NotImplemented("SigmaSpread");}
78
79     void   SetDead() { fIsDead = kTRUE; };
80     Bool_t IsDead() const { return fIsDead; };
81     Int_t Wings()const{return fgkWings;}//Total number of SDD wings
82     Int_t Chips() const{return fgkChips;} // Number of chips/module
83     Int_t Channels() const{ return fgkChannels;}//Number of channels/chip
84     
85     virtual void SetBadChannel(Int_t i,Int_t anode);
86     Int_t GetBadChannel(Int_t i) const {return fBadChannels[i];}
87     Bool_t IsBadChannel(Int_t anode); 
88
89     Float_t GetMapACell(Int_t i,Int_t j) const {
90       if(i<256) return fMapAW0->GetCellContent(i,j);
91       else return fMapAW1->GetCellContent(i,j);
92     }
93     virtual void SetMapA(Int_t wing,AliITSMapSDD* mapA) {
94       if(wing==0) fMapAW0=mapA;
95       else fMapAW1=mapA;
96     } 
97     Float_t GetMapTCell(Int_t i,Int_t j) const {
98       if(i<256) return fMapTW0->GetCellContent(i,j);
99       else return fMapTW1->GetCellContent(i,j);
100     }
101     virtual void SetMapT(Int_t wing,AliITSMapSDD* mapT) {
102       if(wing==0) fMapTW0=mapT;
103       else fMapTW1=mapT;
104     } 
105     static Int_t GetMapTimeNBin() {return fgkMapTimeNBin;} 
106
107     virtual void SetElectronics(Int_t p1=1) {((AliITSresponseSDD*)fResponse)->SetElectronics(p1);}
108     virtual Int_t GetElectronics() const {return ((AliITSresponseSDD*)fResponse)->Electronics();}
109     virtual void SetMaxAdc(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetMaxAdc(p1);}
110     virtual Float_t GetMaxAdc() const {return ((AliITSresponseSDD*)fResponse)->MaxAdc();} 
111     virtual void SetChargeLoss(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetChargeLoss(p1);}
112     virtual Float_t GetChargeLoss() const {return ((AliITSresponseSDD*)fResponse)->ChargeLoss();}
113     virtual void SetDynamicRange(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetDynamicRange(p1);}
114     virtual Float_t GetDynamicRange() const {return ((AliITSresponseSDD*)fResponse)->DynamicRange();} 
115     virtual void SetDriftSpeed(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetDriftSpeed(p1);}
116     virtual void SetDriftSpeedParam(Int_t iWing, Float_t* p);
117     virtual Float_t GetTimeOffset() const {return ((AliITSresponseSDD*)fResponse)->TimeOffset();}
118     virtual Float_t GetADC2keV() const {return ((AliITSresponseSDD*)fResponse)->ADC2keV();}
119     virtual Float_t GetDriftSpeed() const {return ((AliITSresponseSDD*)fResponse)->DriftSpeed();}
120     virtual Float_t GetDriftSpeedAtAnode(Float_t nAnode) const {
121       if(nAnode<256){
122         return fDriftVelParW0[0]+fDriftVelParW0[1]*nAnode+fDriftVelParW0[2]*nAnode*nAnode+fDriftVelParW0[3]*nAnode*nAnode*nAnode;
123       }else{
124         nAnode-=256;
125         return fDriftVelParW1[0]+fDriftVelParW1[1]*nAnode+fDriftVelParW1[2]*nAnode*nAnode+fDriftVelParW1[3]*nAnode*nAnode*nAnode;
126       }
127     }
128     virtual void SetParamOptions(const char *opt1,const char *opt2) {((AliITSresponseSDD*)fResponse)->SetParamOptions(opt1,opt2);}
129     virtual void GetParamOptions(char *opt1,char *opt2) const {((AliITSresponseSDD*)fResponse)->ParamOptions(opt1,opt2);}
130     virtual Bool_t Do10to8() const {return ((AliITSresponseSDD*)fResponse)->Do10to8();}
131     virtual void SetZeroSupp (const char *opt) {((AliITSresponseSDD*)fResponse)->SetZeroSupp(opt);} 
132     virtual const char *GetZeroSuppOption() const {return ((AliITSresponseSDD*)fResponse)->ZeroSuppOption();}
133     virtual void SetNSigmaIntegration(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetNSigmaIntegration(p1);}
134     virtual Float_t GetNSigmaIntegration() const {return ((AliITSresponseSDD*)fResponse)->NSigmaIntegration();}
135     virtual void SetNLookUp(Int_t p1) {((AliITSresponseSDD*)fResponse)->SetNLookUp(p1);}
136     virtual Int_t GetGausNLookUp() const {return ((AliITSresponseSDD*)fResponse)->GausNLookUp();}
137     virtual Float_t GetGausLookUp(Int_t i) const {return ((AliITSresponseSDD*)fResponse)->GausLookUp(i);}
138     virtual Int_t Convert8to10(Int_t signal) const {return ((AliITSresponseSDD*)fResponse)->Convert8to10(signal);}
139     virtual void  SetJitterError(Double_t jitter=20) {((AliITSresponseSDD*)fResponse)->SetJitterError(jitter);}
140     virtual Float_t GetJitterError() const {return ((AliITSresponseSDD*)fResponse)->JitterError();}
141     virtual Float_t GetDriftPath(Float_t time,Float_t /*anodecoord*/) const {return time*GetDriftSpeed();}
142     virtual Float_t GetThresholdAnode(Int_t anode,Int_t nsigma=3) const {
143       return nsigma*fNoiseAfterEl[anode];}
144
145     virtual void  SetDo10to8(Bool_t bitcomp=kTRUE) {((AliITSresponseSDD*)fResponse)->SetDo10to8(bitcomp);}
146  protected:
147
148
149     // these statis const should be move to AliITSsegmentationSDD
150     static const Int_t fgkWings = 2;     // Number of wings per module
151     static const Int_t fgkChips = 4;        // Number of chips/module
152     static const Int_t fgkChannels = 64;    // Number of channels/chip
153     static const Float_t fgkTemperatureDefault; // default for fT (Kelvin)
154     static const Float_t fgkNoiseDefault; // default for fNoise
155     static const Float_t fgkBaselineDefault; // default for fBaseline
156     static const Float_t fgkMinValDefault; // default for fMinVal
157     static const Float_t fgkGainDefault; //default for gain
158     static const Int_t fgkMapTimeNBin = 72; //map granularity along drift direction
159     Int_t fDeadChips;                     // Number of dead chips
160     Int_t fDeadChannels;                  // Number of dead channels
161     Float_t fGain[fgkWings][fgkChips][fgkChannels];//Array for channel gains
162     Float_t fNoise[fgkWings*fgkChips*fgkChannels];          // Noise array
163     Float_t fBaseline[fgkWings*fgkChips*fgkChannels];       // Baseline array
164     Float_t fNoiseAfterEl[fgkWings*fgkChips*fgkChannels];   // Noise after electronics
165     Float_t fMinVal;        // Min value used in 2D zero-suppression algo
166
167     Bool_t   fIsDead;  // module is dead or alive ?
168     TArrayI  fBadChannels; //Array with bad anodes number (0-512) 
169
170     Float_t fDriftVelParW0[4];  // Coeff. of pol3 fit to drift speed vs. anode (wing0)
171     Float_t fDriftVelParW1[4];  // Coeff. of pol3 fit to drift speed vs. anode (wing1)
172     
173     AliITSMapSDD* fMapAW0;     //! map of residuals on anode coord. wing 0
174     AliITSMapSDD* fMapAW1;     //! map of residuals on anode coord. wing 1
175     AliITSMapSDD* fMapTW0;     //! map of residuals on time coord. wing 0
176     AliITSMapSDD* fMapTW1;     //! map of residuals on time coord. wing 1
177
178
179  private:
180     AliITSCalibrationSDD(const AliITSCalibrationSDD &ob); // copy constructor
181     AliITSCalibrationSDD& operator=(const AliITSCalibrationSDD & /* source */); // ass. op.
182
183
184     ClassDef(AliITSCalibrationSDD,7) // SDD response 
185     
186     };
187 #endif