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"
14 class AliITSresponseSDD;
15 ///////////////////////////////////////////////////////
16 // Response for SDD //
17 ///////////////////////////////////////////////////////
19 class AliITSCalibrationSDD : public AliITSCalibration {
22 // Configuration methods
24 AliITSCalibrationSDD();
25 AliITSCalibrationSDD(const char *dataType);
26 virtual ~AliITSCalibrationSDD();
27 virtual void SetNoiseParam(Double_t /*n*/, Double_t /*b*/){
28 NotImplemented("SetNoiseParam");}
30 virtual void GetNoiseParam(Double_t &/*n*/, Double_t &/*b*/) const {
31 NotImplemented("GetNoiseParam");}
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;}
38 virtual void SetThresholds(Double_t mv, Double_t /* b */){
39 // Min value used in 2D - could be used as a threshold setting
41 virtual void Thresholds(Double_t & mv, Double_t & /* b */) const
43 virtual void GiveCompressParam(Int_t *x,Int_t ian) const;
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;}
61 void PrintGains() const;
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");}
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
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);
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);
93 virtual void SetMapA(Int_t wing,AliITSMapSDD* mapA) {
94 if(wing==0) fMapAW0=mapA;
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);
101 virtual void SetMapT(Int_t wing,AliITSMapSDD* mapT) {
102 if(wing==0) fMapTW0=mapT;
105 static Int_t GetMapTimeNBin() {return fgkMapTimeNBin;}
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 {
122 return fDriftVelParW0[0]+fDriftVelParW0[1]*nAnode+fDriftVelParW0[2]*nAnode*nAnode+fDriftVelParW0[3]*nAnode*nAnode*nAnode;
125 return fDriftVelParW1[0]+fDriftVelParW1[1]*nAnode+fDriftVelParW1[2]*nAnode*nAnode+fDriftVelParW1[3]*nAnode*nAnode*nAnode;
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];}
145 virtual void SetDo10to8(Bool_t bitcomp=kTRUE) {((AliITSresponseSDD*)fResponse)->SetDo10to8(bitcomp);}
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
167 Bool_t fIsDead; // module is dead or alive ?
168 TArrayI fBadChannels; //Array with bad anodes number (0-512)
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)
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
180 AliITSCalibrationSDD(const AliITSCalibrationSDD &ob); // copy constructor
181 AliITSCalibrationSDD& operator=(const AliITSCalibrationSDD & /* source */); // ass. op.
184 ClassDef(AliITSCalibrationSDD,7) // SDD response