]>
Commit | Line | Data |
---|---|---|
1 | #ifndef ALIITSRESPONSESDD_H | |
2 | #define ALIITSRESPONSESDD_H | |
3 | ||
4 | #include "TArrayF.h" | |
5 | #include <TString.h> | |
6 | #include <iostream.h> | |
7 | #include "AliITSresponse.h" | |
8 | ||
9 | // response for SDD | |
10 | ||
11 | class AliITSresponseSDD : public AliITSresponse { | |
12 | public: | |
13 | // | |
14 | // Configuration methods | |
15 | // | |
16 | ||
17 | AliITSresponseSDD(); | |
18 | AliITSresponseSDD(const char *dataType); | |
19 | virtual ~AliITSresponseSDD(); | |
20 | ||
21 | void SetElectronics(Int_t p1=1) { | |
22 | // Electronics: Pascal (1) or OLA (2) | |
23 | fElectronics=p1; | |
24 | } | |
25 | ||
26 | Int_t Electronics() { | |
27 | // Electronics: 1 = Pascal; 2 = OLA | |
28 | return fElectronics; | |
29 | } | |
30 | ||
31 | void SetMaxAdc(Float_t p1=1024.) { | |
32 | // Adc-count saturation value | |
33 | fMaxAdc=p1; | |
34 | } | |
35 | Float_t MaxAdc() { | |
36 | // Get maximum Adc-count value | |
37 | return fMaxAdc; | |
38 | } | |
39 | ||
40 | void SetChargeLoss(Float_t p1=0.0) { | |
41 | // Set Linear Charge Loss Steepness // 0.01 for 20% | |
42 | fChargeLoss=p1; | |
43 | } | |
44 | Float_t ChargeLoss() { | |
45 | // Get Charge Loss Coefficient | |
46 | return fChargeLoss; | |
47 | } | |
48 | ||
49 | void SetDynamicRange(Float_t p1=132.) { | |
50 | // Set Dynamic Range | |
51 | fDynamicRange=p1; | |
52 | } | |
53 | Float_t DynamicRange() { | |
54 | // Get Dynamic Range | |
55 | return fDynamicRange; | |
56 | } | |
57 | ||
58 | void SetDiffCoeff(Float_t p1=3.23,Float_t p2=30.) { | |
59 | // Diffusion coefficients | |
60 | fDiffCoeff=p1; | |
61 | fDiffCoeff1=p2; | |
62 | } | |
63 | void DiffCoeff(Float_t&diff,Float_t&diff1) { | |
64 | // Get diffusion coefficients | |
65 | diff = fDiffCoeff; | |
66 | diff1 = fDiffCoeff1; | |
67 | } | |
68 | ||
69 | void SetDriftSpeed(Float_t p1=7.3) { | |
70 | // Drift velocity | |
71 | fDriftSpeed=p1; | |
72 | } | |
73 | Float_t DriftSpeed() { | |
74 | // drift speed | |
75 | return fDriftSpeed; | |
76 | } | |
77 | ||
78 | void SetTemperature(Float_t p1=23.) { | |
79 | // Temperature | |
80 | fTemperature=p1; | |
81 | } | |
82 | Float_t Temperature() { | |
83 | // Get temperature | |
84 | return fTemperature; | |
85 | } | |
86 | ||
87 | void SetDataType(const char *data="simulated") { | |
88 | // Type of data - real or simulated | |
89 | fDataType=data; | |
90 | } | |
91 | const char *DataType() const { | |
92 | // Get data type | |
93 | return fDataType.Data(); | |
94 | } | |
95 | ||
96 | void SetParamOptions(const char *opt1="same",const char *opt2="same"){ | |
97 | // Parameters: "same" or read from "file" | |
98 | fParam1=opt1; fParam2=opt2; | |
99 | } | |
100 | void ParamOptions(char *opt1,char *opt2) { | |
101 | // options | |
102 | strcpy(opt1,fParam1.Data()); strcpy(opt2,fParam2.Data()); | |
103 | } | |
104 | ||
105 | void SetNoiseParam(Float_t n=0., Float_t b=20.){ | |
106 | // Noise and baseline // 8.3 for ALICE with beam test measurements | |
107 | fNoise=n; fBaseline=b; | |
108 | } | |
109 | void SetNoiseAfterElectronics(Float_t n=0.){ | |
110 | // Noise after electronics (ADC units) // 1.6 for ALICE from beam test measurements | |
111 | fNoiseAfterEl=n; | |
112 | } | |
113 | void GetNoiseParam(Float_t &n, Float_t &b) { | |
114 | // get noise param | |
115 | n=fNoise; b=fBaseline; | |
116 | } | |
117 | Float_t GetNoiseAfterElectronics(){ | |
118 | // Noise after electronics (ADC units) | |
119 | return fNoiseAfterEl; | |
120 | } | |
121 | ||
122 | void SetDo10to8(Bool_t bitcomp=kTRUE) { | |
123 | // set the option for 10 to 8 bit compression | |
124 | fBitComp = bitcomp; | |
125 | } | |
126 | ||
127 | Bool_t Do10to8() { | |
128 | // get 10 to 8 compression option | |
129 | return fBitComp; | |
130 | } | |
131 | ||
132 | void SetZeroSupp (const char *opt="1D") { | |
133 | // Zero-suppression option - could be 1D, 2D or non-ZS | |
134 | fOption=opt; | |
135 | } | |
136 | const char *ZeroSuppOption() const { | |
137 | // Get zero-suppression option | |
138 | return fOption.Data(); | |
139 | } | |
140 | void SetMinVal(Int_t mv=4) { | |
141 | // Min value used in 2D - could be used as a threshold setting | |
142 | fMinVal = mv; | |
143 | } | |
144 | Int_t MinVal() { | |
145 | // min val | |
146 | return fMinVal; | |
147 | } | |
148 | ||
149 | void SetFilenames(const char *f1="",const char *f2="",const char *f3="") { | |
150 | // Set filenames - input, output, parameters .... | |
151 | fFileName1=f1; fFileName2=f2; fFileName3=f3; | |
152 | } | |
153 | void Filenames(char *input,char *baseline,char *param) { | |
154 | // Filenames | |
155 | strcpy(input,fFileName1.Data()); strcpy(baseline,fFileName2.Data()); | |
156 | strcpy(param,fFileName3.Data()); | |
157 | } | |
158 | ||
159 | ||
160 | void SetOutputOption(Bool_t write=kFALSE) { | |
161 | // set output option | |
162 | fWrite = write; | |
163 | } | |
164 | Bool_t OutputOption() { | |
165 | // output option | |
166 | return fWrite; | |
167 | } | |
168 | // | |
169 | // Compression parameters | |
170 | void SetCompressParam(Int_t cp[8]); | |
171 | void GiveCompressParam(Int_t *x); | |
172 | ||
173 | // | |
174 | // Detector type response methods | |
175 | void SetNSigmaIntegration(Float_t p1=3.) { | |
176 | // Set number of sigmas over which cluster disintegration is performed | |
177 | fNsigmas=p1; | |
178 | } | |
179 | Float_t NSigmaIntegration() { | |
180 | // Get number of sigmas over which cluster disintegration is performed | |
181 | return fNsigmas; | |
182 | } | |
183 | void SetNLookUp(Int_t p1=121) { | |
184 | // Set number of sigmas over which cluster disintegration is performed | |
185 | fNcomps=p1; | |
186 | fGaus = new TArrayF(fNcomps+1); | |
187 | for(Int_t i=0; i<=fNcomps; i++) { | |
188 | Float_t x = -fNsigmas + (2.*i*fNsigmas)/(fNcomps-1); | |
189 | (*fGaus)[i] = exp(-((x*x)/2)); | |
190 | // cout << "fGaus[" << i << "]: " << fGaus->At(i) << endl; | |
191 | } | |
192 | } | |
193 | // Get number of intervals in which the gaussian lookup table is divided | |
194 | Int_t GausNLookUp() {return fNcomps;} | |
195 | ||
196 | Float_t IntPH(Float_t eloss) { | |
197 | // Pulse height from scored quantity (eloss) | |
198 | return 0.; | |
199 | } | |
200 | Float_t IntXZ(AliITSsegmentation *) { | |
201 | // Charge disintegration | |
202 | return 0.; | |
203 | } | |
204 | Float_t GausLookUp(Int_t i) { | |
205 | if(i<0 || i>=fNcomps) return 0.; | |
206 | return fGaus->At(i); | |
207 | } | |
208 | void SetDeadChannels(Int_t nmodules=0, Int_t nchips=0, Int_t nchannels=0); | |
209 | void PrintGains(); | |
210 | void Print(); | |
211 | ||
212 | ||
213 | private: | |
214 | ||
215 | AliITSresponseSDD(const AliITSresponseSDD &source); // copy constructor | |
216 | AliITSresponseSDD& operator=(const AliITSresponseSDD &source); // ass. op. | |
217 | ||
218 | protected: | |
219 | ||
220 | static const Int_t fModules = 520; // Total number of SDD modules | |
221 | static const Int_t fChips = 4; // Number of chips/module | |
222 | static const Int_t fChannels = 64; // Number of channels/chip | |
223 | Float_t fGain[fModules][fChips][fChannels]; // Array for channel gains | |
224 | ||
225 | Int_t fCPar[8]; // Hardware compression parameters | |
226 | Float_t fNoise; // Noise | |
227 | Float_t fBaseline; // Baseline | |
228 | Float_t fNoiseAfterEl; // Noise after electronics | |
229 | Float_t fDynamicRange; // Set Dynamic Range | |
230 | Float_t fChargeLoss; // Set Linear Coefficient for Charge Loss | |
231 | Float_t fTemperature; // Temperature | |
232 | Float_t fDriftSpeed; // Drift velocity | |
233 | Int_t fElectronics; // Electronics | |
234 | ||
235 | Float_t fMaxAdc; // Adc saturation value | |
236 | Float_t fDiffCoeff; // Diffusion Coefficient (scaling the time) | |
237 | Float_t fDiffCoeff1; // Diffusion Coefficient (constant term) | |
238 | Float_t fNsigmas; // Number of sigmas over which charge disintegration | |
239 | // is performed | |
240 | TArrayF *fGaus; // Gaussian lookup table for signal generation | |
241 | Int_t fNcomps; // Number of samplings along the gaussian | |
242 | ||
243 | Int_t fMinVal; // Min value used in 2D zero-suppression algo | |
244 | ||
245 | Bool_t fWrite; // Write option for the compression algorithms | |
246 | Bool_t fBitComp; // 10 to 8 bit compression option | |
247 | ||
248 | TString fOption; // Zero-suppresion option (1D, 2D or none) | |
249 | TString fParam1; // Read baselines from file option | |
250 | TString fParam2; // Read compression algo thresholds from file | |
251 | ||
252 | TString fDataType; // data type - real or simulated | |
253 | TString fFileName1; // input keys : run, module # | |
254 | TString fFileName2; // baseline & noise val or output coded // signal or monitored bgr. | |
255 | TString fFileName3; // param values or output coded signal | |
256 | ||
257 | ClassDef(AliITSresponseSDD,2) // SDD response | |
258 | }; | |
259 | #endif | |
260 |