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 "AliITSsegmentationSDD.h"
13 class AliITSCorrMapSDD;
14 class AliITSDriftSpeedArraySDD;
16 ///////////////////////////////////////////////////////
17 // Response for SDD //
18 ///////////////////////////////////////////////////////
20 class AliITSCalibrationSDD : public AliITSCalibration {
23 // Configuration methods
25 AliITSCalibrationSDD();
26 AliITSCalibrationSDD(const char *dataType);
27 virtual ~AliITSCalibrationSDD();
29 virtual Float_t GetBaseline(Int_t anode) const {return fBaseline[anode];}
30 virtual void SetBaseline(Int_t anode,Double_t bas) {fBaseline[anode]=bas;}
31 virtual Float_t GetNoise(Int_t anode) const {return fNoise[anode];}
32 virtual void SetNoise(Int_t anode, Double_t noise) {fNoise[anode]=noise;}
34 virtual void GiveCompressParam(Int_t *x) const;
36 void SetZSLowThreshold(Int_t iWing, Int_t thr=25){fZSTL[iWing]=thr;}
37 void SetZSHighThreshold(Int_t iWing, Int_t thr=29){fZSTH[iWing]=thr;}
38 Int_t GetZSLowThreshold(Int_t iWing) const {return fZSTL[iWing];}
39 Int_t GetZSHighThreshold(Int_t iWing) const {return fZSTH[iWing];}
41 void SetNoiseAfterElectronics(Int_t anode,Double_t n=2.38){
42 // Noise after electronics (ADC units)
43 // 2.36 for ALICE from beam test measurements 2001
44 fNoiseAfterEl[anode]=n;}
45 Float_t GetNoiseAfterElectronics(Int_t anode) const {
46 // Noise after electronics (ADC units)
47 return fNoiseAfterEl[anode];}
48 //void SetDeadChannels(Int_t nchips=0, Int_t nchannels=0);
49 void SetDeadChannels(Int_t ndead=0){fDeadChannels=ndead; fBadChannels.Set(ndead);}
50 Int_t GetDeadChips() const { return fDeadChips; }
51 Int_t GetDeadChannels() const { return fDeadChannels; }
52 Float_t GetChannelGain(Int_t anode) const {return fGain[anode];}
53 virtual void SetGain(Int_t anode,Double_t g){fGain[anode]=g;}
56 Int_t GetWing(Int_t anode) const{
57 if(anode>=fgkChips*fgkChannels) return 1;
60 Int_t GetChipChannel(Int_t anode) const {return anode%fgkChannels;}
61 Int_t GetChip(Int_t anode) const {return anode/fgkChannels;}
62 Int_t GetAnodeNumber(Int_t iwing, Int_t ichip03, Int_t ichan) const {
63 if(iwing>=2 || ichip03>=4 || ichan>=64) return -1;
64 else return iwing*fgkChips*fgkChannels+ichip03*fgkChannels+ichan;
66 Int_t GetAnodeNumber(Int_t ichip07, Int_t ichan) const {
67 if(ichip07>=8 || ichan>=64) return -1;
68 else return ichip07*fgkChannels+ichan;
71 void PrintGains() const;
73 virtual void Print(ostream *os) const {AliITSCalibrationSDD::Print(os);}
74 virtual void Print(Option_t *option="") const {AliITSCalibrationSDD::Print(option);}
75 // not implemented virtual methods (devlared in the mother class
76 virtual void SetDetParam(Double_t *)
77 {NotImplemented("SetDetParam");}
78 virtual void GetDetParam(Double_t *) const
79 {NotImplemented("GetDetParam");}
80 virtual void SetNDetParam(Int_t /* n */)
81 {NotImplemented("SetNDetParam");}
82 virtual Int_t NDetParam() const
83 {NotImplemented("NDetParam"); return 0;}
84 virtual void SetSigmaSpread(Double_t, Double_t)
85 {NotImplemented("SetSigmaSpread");}
86 virtual void SigmaSpread(Double_t & /* p1 */,Double_t & /* p2 */) const
87 {NotImplemented("SigmaSpread");}
91 for(Int_t i=0;i<fgkChips*fgkWings;i++) fIsChipBad[i]=kTRUE;
93 virtual Bool_t IsBad() const { return fIsBad; }
94 void SetChipBad(Int_t nChip) {
95 fIsChipBad[nChip] = kTRUE;
97 virtual Bool_t IsChipBad(Int_t nChip) const {
98 return fIsChipBad[nChip];
100 Int_t Wings()const{return fgkWings;}//Total number of SDD wings
101 Int_t Chips() const{return fgkChips;} // Number of chips/module
102 Int_t Channels() const{ return fgkChannels;}//Number of channels/chip
103 Int_t NOfAnodes() const {return fgkChannels*fgkChips*fgkWings;}
105 virtual void SetBadChannel(Int_t i,Int_t anode);
106 Int_t GetBadChannel(Int_t i) const {return fBadChannels[i];}
107 Bool_t IsBadChannel(Int_t anode) const{
108 if(GetChannelGain(anode)==0) return kTRUE;
111 Float_t GetMapACell(Int_t i,Int_t j) const {
112 if(i<256) return fMapAW0->GetCellContent(i,j);
113 else return fMapAW1->GetCellContent(i-256,j);
115 virtual void SetMapA(Int_t wing,AliITSCorrMapSDD* mapA) {
116 if(wing==0) fMapAW0=mapA;
119 Float_t GetMapTCell(Int_t i,Int_t j) const {
120 if(i<256) return fMapTW0->GetCellContent(i,j);
121 else return fMapTW1->GetCellContent(i-256,j);
123 virtual void SetMapT(Int_t wing,AliITSCorrMapSDD* mapT) {
124 if(wing==0) fMapTW0=mapT;
128 virtual void SetDriftSpeed(Int_t wing, AliITSDriftSpeedArraySDD* arr){
129 if(wing==0) fDrSpeed0=arr;
134 virtual Float_t GetDriftSpeedAtAnode(Float_t nAnode) const{
135 if(fDrSpeed0==0 || fDrSpeed1==0) AliFatal("Drift speed not set\n");
136 if(nAnode<256) return fDrSpeed0->GetDriftSpeed(0,nAnode);
137 else return fDrSpeed1->GetDriftSpeed(0,nAnode-256);
140 virtual void SetZeroSupp(Bool_t opt=kTRUE) {fZeroSupp=opt;}
141 virtual Bool_t GetZeroSupp() const {return fZeroSupp;}
143 virtual void SetAMAt40MHz() {fAMAt20MHz=kFALSE;}
144 virtual void SetAMAt20MHz() {fAMAt20MHz=kTRUE;}
145 virtual Bool_t IsAMAt20MHz() const {return fAMAt20MHz;}
147 void GetCorrections(Float_t z, Float_t x, Float_t &devz, Float_t &devx, AliITSsegmentationSDD* seg);
148 virtual Float_t GetThresholdAnode(Int_t anode, Double_t nsigma=2.2) const {
149 return nsigma*fNoiseAfterEl[anode];}
155 // these statis const should be move to AliITSsegmentationSDD
156 static const Int_t fgkWings = 2; // Number of wings per module
157 static const Int_t fgkChips = 4; // Number of chips/module
158 static const Int_t fgkChannels = 64; // Number of channels/chip
159 static const Float_t fgkTemperatureDefault; // default for fT (Kelvin)
160 static const Float_t fgkNoiseDefault; // default for fNoise
161 static const Float_t fgkBaselineDefault; // default for fBaseline
162 static const Float_t fgkGainDefault; //default for gain
164 Bool_t fZeroSupp; // zero suppression
165 Bool_t fAMAt20MHz; // flag for Analog memory of Pascal at 20 MHz
166 Int_t fDeadChips; // Number of dead chips
167 Int_t fDeadChannels; // Number of dead channels
168 Float_t fGain[fgkWings*fgkChips*fgkChannels]; //Array for channel gains
169 Float_t fNoise[fgkWings*fgkChips*fgkChannels]; // Noise array
170 Float_t fBaseline[fgkWings*fgkChips*fgkChannels]; // Baseline array
171 Float_t fNoiseAfterEl[fgkWings*fgkChips*fgkChannels]; // Noise after electronics
173 Int_t fZSTL[2]; // Low threshold in 2D zero-suppression (2 hybrids)
174 Int_t fZSTH[2]; // High threshold in 2D zero-suppression (2 hybrids)
176 Bool_t fIsBad; // module is dead or alive ?
177 Bool_t fIsChipBad[fgkWings*fgkChips]; // chip is dead or alive ?
178 TArrayI fBadChannels; //Array with bad anodes number (0-512)
181 AliITSCorrMapSDD* fMapAW0; //! map of residuals on anode coord. wing 0
182 AliITSCorrMapSDD* fMapAW1; //! map of residuals on anode coord. wing 1
183 AliITSCorrMapSDD* fMapTW0; //! map of residuals on time coord. wing 0
184 AliITSCorrMapSDD* fMapTW1; //! map of residuals on time coord. wing 1
185 AliITSDriftSpeedArraySDD* fDrSpeed0; //! drift speed for wing 0
186 AliITSDriftSpeedArraySDD* fDrSpeed1; //! drift speed for wing 1
189 AliITSCalibrationSDD(const AliITSCalibrationSDD &ob); // copy constructor
190 AliITSCalibrationSDD& operator=(const AliITSCalibrationSDD & /* source */); // ass. op.
193 ClassDef(AliITSCalibrationSDD,17)