1 #ifndef ALIITSCALIBRATIONSDD_H
2 #define ALIITSCALIBRATIONSDD_H
4 /* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
9 #include "AliITSCalibration.h"
10 #include "AliITSresponseSDD.h"
11 #include "AliITSsegmentationSDD.h"
15 class AliITSDriftSpeedArraySDD;
16 class AliITSresponseSDD;
17 ///////////////////////////////////////////////////////
18 // Response for SDD //
19 ///////////////////////////////////////////////////////
21 class AliITSCalibrationSDD : public AliITSCalibration {
24 // Configuration methods
26 AliITSCalibrationSDD();
27 AliITSCalibrationSDD(const char *dataType);
28 virtual ~AliITSCalibrationSDD();
29 virtual void SetNoiseParam(Double_t /*n*/, Double_t /*b*/){
30 NotImplemented("SetNoiseParam");}
32 virtual void GetNoiseParam(Double_t &/*n*/, Double_t &/*b*/) const {
33 NotImplemented("GetNoiseParam");}
35 virtual Float_t GetBaseline(Int_t anode) const {return fBaseline[anode];}
36 virtual void SetBaseline(Int_t anode,Double_t bas) {fBaseline[anode]=bas;}
37 virtual Float_t GetNoise(Int_t anode) const {return fNoise[anode];}
38 virtual void SetNoise(Int_t anode, Double_t noise) {fNoise[anode]=noise;}
40 virtual void SetThresholds(Double_t mv, Double_t /* b */){
41 // Min value used in 2D - could be used as a threshold setting
43 virtual void Thresholds(Double_t & mv, Double_t & /* b */) const
45 virtual void GiveCompressParam(Int_t *x,Int_t ian) const;
47 void SetNoiseAfterElectronics(Int_t anode,Double_t n=2.38){
48 // Noise after electronics (ADC units)
49 // 2.36 for ALICE from beam test measurements 2001
50 fNoiseAfterEl[anode]=n;}
51 Float_t GetNoiseAfterElectronics(Int_t anode) const {
52 // Noise after electronics (ADC units)
53 return fNoiseAfterEl[anode];}
54 //void SetDeadChannels(Int_t nchips=0, Int_t nchannels=0);
55 void SetDeadChannels(Int_t ndead=0){fDeadChannels=ndead; fBadChannels.Set(ndead);}
56 Int_t GetDeadChips() const { return fDeadChips; }
57 Int_t GetDeadChannels() const { return fDeadChannels; }
58 Float_t GetChannelGain(Int_t anode) const {return fGain[anode];}
59 virtual void SetGain(Int_t anode,Double_t g){fGain[anode]=g;}
62 Int_t GetWing(Int_t anode) const{
63 if(anode>=fgkChips*fgkChannels) return 1;
66 Int_t GetChipChannel(Int_t anode) const {return anode%fgkChannels;}
67 Int_t GetChip(Int_t anode) const {return anode/fgkChannels;}
68 Int_t GetAnodeNumber(Int_t iwing, Int_t ichip03, Int_t ichan) const {
69 if(iwing>=2 || ichip03>=4 || ichan>=64) return -1;
70 else return iwing*fgkChips*fgkChannels+ichip03*fgkChannels+ichan;
72 Int_t GetAnodeNumber(Int_t ichip07, Int_t ichan) const {
73 if(ichip07>=8 || ichan>=64) return -1;
74 else return ichip07*fgkChannels+ichan;
77 void PrintGains() const;
79 virtual void Print(ostream *os) const {AliITSCalibrationSDD::Print(os);}
80 virtual void Print(Option_t *option="") const {AliITSCalibrationSDD::Print(option);}
81 // not implemented virtual methods (devlared in the mother class
82 virtual void SetDetParam(Double_t *)
83 {NotImplemented("SetDetParam");}
84 virtual void GetDetParam(Double_t *) const
85 {NotImplemented("GetDetParam");}
86 virtual void SetNDetParam(Int_t /* n */)
87 {NotImplemented("SetNDetParam");}
88 virtual Int_t NDetParam() const
89 {NotImplemented("NDetParam"); return 0;}
90 virtual void SetSigmaSpread(Double_t, Double_t)
91 {NotImplemented("SetSigmaSpread");}
92 virtual void SigmaSpread(Double_t & /* p1 */,Double_t & /* p2 */) const
93 {NotImplemented("SigmaSpread");}
97 for(Int_t i=0;i<fgkChips*fgkWings;i++) fIsChipBad[i]=kTRUE;
99 virtual Bool_t IsBad() const { return fIsBad; }
100 void SetChipBad(Int_t nChip) {
101 fIsChipBad[nChip] = kTRUE;
103 virtual Bool_t IsChipBad(Int_t nChip) const {
104 return fIsChipBad[nChip];
106 Int_t Wings()const{return fgkWings;}//Total number of SDD wings
107 Int_t Chips() const{return fgkChips;} // Number of chips/module
108 Int_t Channels() const{ return fgkChannels;}//Number of channels/chip
110 virtual void SetBadChannel(Int_t i,Int_t anode);
111 Int_t GetBadChannel(Int_t i) const {return fBadChannels[i];}
112 Bool_t IsBadChannel(Int_t anode) const{
113 if(GetChannelGain(anode)==0) return kTRUE;
116 void SetUseCorrectionMaps(Bool_t useAnodeMap, Bool_t useDriftMap){
117 fUseACorrMap=useAnodeMap;
118 fUseTCorrMap=useDriftMap;
120 Float_t GetMapACell(Int_t i,Int_t j) const {
121 if(i<256) return fMapAW0->GetCellContent(i,j);
122 else return fMapAW1->GetCellContent(i-256,j);
124 virtual void SetMapA(Int_t wing,AliITSMapSDD* mapA) {
125 if(wing==0) fMapAW0=mapA;
128 Float_t GetMapTCell(Int_t i,Int_t j) const {
129 if(i<256) return fMapTW0->GetCellContent(i,j);
130 else return fMapTW1->GetCellContent(i-256,j);
132 virtual void SetMapT(Int_t wing,AliITSMapSDD* mapT) {
133 if(wing==0) fMapTW0=mapT;
137 virtual void SetDriftSpeed(Int_t wing, AliITSDriftSpeedArraySDD* arr){
138 if(wing==0) fDrSpeed0=arr;
142 virtual void SetElectronics(Int_t p1=1) {((AliITSresponseSDD*)fResponse)->SetElectronics(p1);}
143 virtual Int_t GetElectronics() const {return ((AliITSresponseSDD*)fResponse)->Electronics();}
144 virtual void SetMaxAdc(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetMaxAdc(p1);}
145 virtual Float_t GetMaxAdc() const {return ((AliITSresponseSDD*)fResponse)->MaxAdc();}
146 virtual void SetChargeLoss(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetChargeLoss(p1);}
147 virtual Float_t GetChargeLoss() const {return ((AliITSresponseSDD*)fResponse)->ChargeLoss();}
148 virtual void SetDynamicRange(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetDynamicRange(p1);}
149 virtual Float_t GetDynamicRange() const {return ((AliITSresponseSDD*)fResponse)->DynamicRange();}
151 virtual Float_t GetTimeOffset() const {return ((AliITSresponseSDD*)fResponse)->TimeOffset();}
152 virtual Float_t GetADC2keV() const {return ((AliITSresponseSDD*)fResponse)->ADC2keV();}
153 virtual Float_t GetDriftSpeedAtAnode(Float_t nAnode) const{
154 if(fDrSpeed0==0 || fDrSpeed1==0) AliFatal("Drift speed not set\n");
155 if(nAnode<256) return fDrSpeed0->GetDriftSpeed(0,nAnode);
156 else return fDrSpeed1->GetDriftSpeed(0,nAnode-256);
158 virtual void SetParamOptions(const char *opt1,const char *opt2) {((AliITSresponseSDD*)fResponse)->SetParamOptions(opt1,opt2);}
159 virtual void GetParamOptions(char *opt1,char *opt2) const {((AliITSresponseSDD*)fResponse)->ParamOptions(opt1,opt2);}
160 virtual Bool_t Do10to8() const {return ((AliITSresponseSDD*)fResponse)->Do10to8();}
161 virtual void SetZeroSupp (const char *opt) {((AliITSresponseSDD*)fResponse)->SetZeroSupp(opt);}
162 virtual const char *GetZeroSuppOption() const {return ((AliITSresponseSDD*)fResponse)->ZeroSuppOption();}
163 virtual void SetNSigmaIntegration(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetNSigmaIntegration(p1);}
164 virtual Float_t GetNSigmaIntegration() const {return ((AliITSresponseSDD*)fResponse)->NSigmaIntegration();}
165 virtual void SetNLookUp(Int_t p1) {((AliITSresponseSDD*)fResponse)->SetNLookUp(p1);}
166 virtual Int_t GetGausNLookUp() const {return ((AliITSresponseSDD*)fResponse)->GausNLookUp();}
167 virtual Float_t GetGausLookUp(Int_t i) const {return ((AliITSresponseSDD*)fResponse)->GausLookUp(i);}
168 virtual Int_t Convert8to10(Int_t signal) const {return ((AliITSresponseSDD*)fResponse)->Convert8to10(signal);}
169 virtual void SetJitterError(Double_t jitter=20) {((AliITSresponseSDD*)fResponse)->SetJitterError(jitter);}
170 virtual Float_t GetJitterError() const {return ((AliITSresponseSDD*)fResponse)->JitterError();}
171 virtual Float_t GetDriftPath(Float_t time, Float_t xAnode) const {return time*GetDriftSpeedAtAnode(xAnode);}
172 void GetCorrections(Float_t z, Float_t x, Float_t &devz, Float_t &devx, AliITSsegmentationSDD* seg);
173 virtual Float_t GetThresholdAnode(Int_t anode,Int_t nsigma=3) const {
174 return nsigma*fNoiseAfterEl[anode];}
176 virtual void SetDo10to8(Bool_t bitcomp=kTRUE) {((AliITSresponseSDD*)fResponse)->SetDo10to8(bitcomp);}
180 // these statis const should be move to AliITSsegmentationSDD
181 static const Int_t fgkWings = 2; // Number of wings per module
182 static const Int_t fgkChips = 4; // Number of chips/module
183 static const Int_t fgkChannels = 64; // Number of channels/chip
184 static const Float_t fgkTemperatureDefault; // default for fT (Kelvin)
185 static const Float_t fgkNoiseDefault; // default for fNoise
186 static const Float_t fgkBaselineDefault; // default for fBaseline
187 static const Float_t fgkMinValDefault; // default for fMinVal
188 static const Float_t fgkGainDefault; //default for gain
190 Int_t fDeadChips; // Number of dead chips
191 Int_t fDeadChannels; // Number of dead channels
192 Float_t fGain[fgkWings*fgkChips*fgkChannels]; //Array for channel gains
193 Float_t fNoise[fgkWings*fgkChips*fgkChannels]; // Noise array
194 Float_t fBaseline[fgkWings*fgkChips*fgkChannels]; // Baseline array
195 Float_t fNoiseAfterEl[fgkWings*fgkChips*fgkChannels]; // Noise after electronics
196 Float_t fMinVal; // Min value used in 2D zero-suppression algo
198 Bool_t fIsBad; // module is dead or alive ?
199 Bool_t fIsChipBad[fgkWings*fgkChips]; // chip is dead or alive ?
200 TArrayI fBadChannels; //Array with bad anodes number (0-512)
203 Bool_t fUseACorrMap; // flag for the use of correction maps (anode)
204 Bool_t fUseTCorrMap; // flag for the use of correction maps (drift)
206 AliITSMapSDD* fMapAW0; //! map of residuals on anode coord. wing 0
207 AliITSMapSDD* fMapAW1; //! map of residuals on anode coord. wing 1
208 AliITSMapSDD* fMapTW0; //! map of residuals on time coord. wing 0
209 AliITSMapSDD* fMapTW1; //! map of residuals on time coord. wing 1
210 AliITSDriftSpeedArraySDD* fDrSpeed0; //! drift speed for wing 0
211 AliITSDriftSpeedArraySDD* fDrSpeed1; //! drift speed for wing 1
214 AliITSCalibrationSDD(const AliITSCalibrationSDD &ob); // copy constructor
215 AliITSCalibrationSDD& operator=(const AliITSCalibrationSDD & /* source */); // ass. op.
218 ClassDef(AliITSCalibrationSDD,10) // SDD response