]>
Commit | Line | Data |
---|---|---|
3551db50 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | /* $Id$ */ | |
17 | ||
4d18a639 | 18 | //////////////////////////////////////////////////////////////////////////// |
19 | // // | |
20 | // Class containing constant simulation parameters // | |
21 | // // | |
22 | // Request an instance with AliTRDSimParam::Instance() // | |
23 | // Then request the needed values // | |
24 | // // | |
25 | //////////////////////////////////////////////////////////////////////////// | |
3551db50 | 26 | |
a9151110 | 27 | #include <TMath.h> |
28 | ||
3551db50 | 29 | #include "AliRun.h" |
30 | ||
31 | #include "AliTRDSimParam.h" | |
a076fc2f | 32 | #include "AliTRDCommonParam.h" |
c93255fe | 33 | #include "AliLog.h" |
3551db50 | 34 | |
35 | ClassImp(AliTRDSimParam) | |
36 | ||
4d18a639 | 37 | AliTRDSimParam *AliTRDSimParam::fgInstance = 0; |
38 | Bool_t AliTRDSimParam::fgTerminated = kFALSE; | |
3551db50 | 39 | |
40 | //_ singleton implementation __________________________________________________ | |
41 | AliTRDSimParam* AliTRDSimParam::Instance() | |
42 | { | |
43 | // | |
44 | // Singleton implementation | |
45 | // Returns an instance of this class, it is created if neccessary | |
46 | // | |
47 | ||
4d18a639 | 48 | if (fgTerminated != kFALSE) { |
3551db50 | 49 | return 0; |
4d18a639 | 50 | } |
3551db50 | 51 | |
4d18a639 | 52 | if (fgInstance == 0) { |
3551db50 | 53 | fgInstance = new AliTRDSimParam(); |
4d18a639 | 54 | } |
55 | ||
3551db50 | 56 | return fgInstance; |
4d18a639 | 57 | |
3551db50 | 58 | } |
59 | ||
4d18a639 | 60 | //_ singleton implementation __________________________________________________ |
3551db50 | 61 | void AliTRDSimParam::Terminate() |
62 | { | |
63 | // | |
64 | // Singleton implementation | |
4d18a639 | 65 | // Deletes the instance of this class and sets the terminated flag, |
66 | // instances cannot be requested anymore | |
3551db50 | 67 | // This function can be called several times. |
68 | // | |
69 | ||
70 | fgTerminated = kTRUE; | |
71 | ||
4d18a639 | 72 | if (fgInstance != 0) { |
3551db50 | 73 | delete fgInstance; |
74 | fgInstance = 0; | |
75 | } | |
4d18a639 | 76 | |
3551db50 | 77 | } |
78 | ||
79 | //_____________________________________________________________________________ | |
80 | AliTRDSimParam::AliTRDSimParam() | |
4d18a639 | 81 | :TObject() |
82 | ,fGasGain(0.0) | |
83 | ,fNoise(0.0) | |
84 | ,fChipGain(0.0) | |
85 | ,fADCoutRange(0.0) | |
86 | ,fADCinRange(0.0) | |
4d18a639 | 87 | ,fADCbaseline(0) |
88 | ,fDiffusionOn(kFALSE) | |
89 | ,fElAttachOn(kFALSE) | |
90 | ,fElAttachProp(0.0) | |
91 | ,fTRFOn(kFALSE) | |
92 | ,fTRFsmp(0) | |
93 | ,fTRFbin(0) | |
94 | ,fTRFlo(0.0) | |
95 | ,fTRFhi(0.0) | |
96 | ,fTRFwid(0.0) | |
97 | ,fCTOn(kFALSE) | |
98 | ,fCTsmp(0) | |
4d18a639 | 99 | ,fPadCoupling(0.0) |
100 | ,fTimeCoupling(0.0) | |
101 | ,fTimeStructOn(kFALSE) | |
102 | ,fPRFOn(kFALSE) | |
966f6939 | 103 | ,fNTimeBins(0) |
876a928e | 104 | ,fNTBoverwriteOCDB(kFALSE) |
3551db50 | 105 | { |
106 | // | |
4d18a639 | 107 | // Default constructor |
3551db50 | 108 | // |
109 | ||
3551db50 | 110 | Init(); |
4d18a639 | 111 | |
3551db50 | 112 | } |
113 | ||
114 | //_____________________________________________________________________________ | |
115 | void AliTRDSimParam::Init() | |
116 | { | |
117 | // | |
4d18a639 | 118 | // Default initializiation |
3551db50 | 119 | // |
120 | ||
121 | // The default parameter for the digitization | |
b43a3e17 | 122 | fGasGain = 4000.0; |
123 | fChipGain = 12.4; | |
60243ea4 | 124 | fNoise = 1250.0; |
b43a3e17 | 125 | fADCoutRange = 1023.0; // 10-bit ADC |
126 | fADCinRange = 2000.0; // 2V input range | |
af3880b4 | 127 | fADCbaseline = 10; |
3551db50 | 128 | |
129 | // Diffusion on | |
b43a3e17 | 130 | fDiffusionOn = kTRUE; |
3551db50 | 131 | |
132 | // Propability for electron attachment | |
b43a3e17 | 133 | fElAttachOn = kFALSE; |
134 | fElAttachProp = 0.0; | |
3551db50 | 135 | |
136 | // The time response function | |
b43a3e17 | 137 | fTRFOn = kTRUE; |
3551db50 | 138 | |
139 | // The cross talk | |
b43a3e17 | 140 | fCTOn = kTRUE; |
3551db50 | 141 | |
3551db50 | 142 | // The pad coupling factor |
4d18a639 | 143 | // Use 0.46, instead of the theroetical value 0.3, since it reproduces better |
144 | // the test beam data, even tough it is not understood why. | |
b43a3e17 | 145 | fPadCoupling = 0.46; |
3551db50 | 146 | |
147 | // The time coupling factor (same number as for the TPC) | |
b43a3e17 | 148 | fTimeCoupling = 0.4; |
3551db50 | 149 | |
3551db50 | 150 | // Use drift time maps |
b43a3e17 | 151 | fTimeStructOn = kTRUE; |
3551db50 | 152 | |
cc7cef99 | 153 | // The pad response function |
b43a3e17 | 154 | fPRFOn = kTRUE; |
cc7cef99 | 155 | |
966f6939 | 156 | // The number of time bins |
29391df8 | 157 | fNTimeBins = 22; |
876a928e | 158 | fNTBoverwriteOCDB = kFALSE; |
966f6939 | 159 | |
3551db50 | 160 | ReInit(); |
4d18a639 | 161 | |
ab0a4106 | 162 | } |
3551db50 | 163 | |
164 | //_____________________________________________________________________________ | |
165 | AliTRDSimParam::~AliTRDSimParam() | |
166 | { | |
167 | // | |
4d18a639 | 168 | // Destructor |
3551db50 | 169 | // |
170 | ||
171 | if (fTRFsmp) { | |
172 | delete [] fTRFsmp; | |
173 | fTRFsmp = 0; | |
174 | } | |
175 | ||
176 | if (fCTsmp) { | |
177 | delete [] fCTsmp; | |
178 | fCTsmp = 0; | |
179 | } | |
4d18a639 | 180 | |
ab0a4106 | 181 | } |
3551db50 | 182 | |
183 | //_____________________________________________________________________________ | |
4d18a639 | 184 | AliTRDSimParam::AliTRDSimParam(const AliTRDSimParam &p) |
185 | :TObject(p) | |
186 | ,fGasGain(p.fGasGain) | |
187 | ,fNoise(p.fNoise) | |
188 | ,fChipGain(p.fChipGain) | |
189 | ,fADCoutRange(p.fADCoutRange) | |
190 | ,fADCinRange(p.fADCinRange) | |
4d18a639 | 191 | ,fADCbaseline(p.fADCbaseline) |
192 | ,fDiffusionOn(p.fDiffusionOn) | |
193 | ,fElAttachOn(p.fElAttachOn) | |
194 | ,fElAttachProp(p.fElAttachProp) | |
195 | ,fTRFOn(p.fTRFOn) | |
196 | ,fTRFsmp(0) | |
197 | ,fTRFbin(p.fTRFbin) | |
198 | ,fTRFlo(p.fTRFlo) | |
199 | ,fTRFhi(p.fTRFhi) | |
200 | ,fTRFwid(p.fTRFwid) | |
201 | ,fCTOn(p.fCTOn) | |
202 | ,fCTsmp(0) | |
4d18a639 | 203 | ,fPadCoupling(p.fPadCoupling) |
204 | ,fTimeCoupling(p.fTimeCoupling) | |
205 | ,fTimeStructOn(p.fTimeStructOn) | |
206 | ,fPRFOn(p.fPRFOn) | |
966f6939 | 207 | ,fNTimeBins(p.fNTimeBins) |
876a928e | 208 | ,fNTBoverwriteOCDB(p.fNTBoverwriteOCDB) |
3551db50 | 209 | { |
210 | // | |
4d18a639 | 211 | // Copy constructor |
3551db50 | 212 | // |
213 | ||
4d18a639 | 214 | Int_t iBin = 0; |
215 | ||
024c0422 | 216 | fTRFsmp = new Float_t[fTRFbin]; |
4d18a639 | 217 | for (iBin = 0; iBin < fTRFbin; iBin++) { |
024c0422 | 218 | fTRFsmp[iBin] = ((AliTRDSimParam &) p).fTRFsmp[iBin]; |
4d18a639 | 219 | } |
220 | ||
024c0422 | 221 | fCTsmp = new Float_t[fTRFbin]; |
4d18a639 | 222 | for (iBin = 0; iBin < fTRFbin; iBin++) { |
024c0422 | 223 | fCTsmp[iBin] = ((AliTRDSimParam &) p).fCTsmp[iBin]; |
4d18a639 | 224 | } |
3551db50 | 225 | |
4d18a639 | 226 | } |
3551db50 | 227 | |
228 | //_____________________________________________________________________________ | |
229 | AliTRDSimParam &AliTRDSimParam::operator=(const AliTRDSimParam &p) | |
230 | { | |
231 | // | |
232 | // Assignment operator | |
233 | // | |
234 | ||
089086d3 | 235 | if (this == &p) { |
236 | return *this; | |
b43a3e17 | 237 | } |
4d18a639 | 238 | |
1a0c151e | 239 | Init(); |
240 | ||
089086d3 | 241 | fGasGain = p.fGasGain; |
242 | fNoise = p.fNoise; | |
243 | fChipGain = p.fChipGain; | |
244 | fADCoutRange = p.fADCoutRange; | |
245 | fADCinRange = p.fADCinRange; | |
246 | fADCbaseline = p.fADCbaseline; | |
247 | fDiffusionOn = p.fDiffusionOn; | |
248 | fElAttachOn = p.fElAttachOn; | |
249 | fElAttachProp = p.fElAttachProp; | |
250 | fTRFOn = p.fTRFOn; | |
251 | fTRFsmp = 0; | |
252 | fTRFbin = p.fTRFbin; | |
253 | fTRFlo = p.fTRFlo; | |
254 | fTRFhi = p.fTRFhi; | |
255 | fTRFwid = p.fTRFwid; | |
256 | fCTOn = p.fCTOn; | |
257 | fCTsmp = 0; | |
258 | fPadCoupling = p.fPadCoupling; | |
259 | fTimeCoupling = p.fTimeCoupling; | |
260 | fTimeStructOn = p.fTimeStructOn; | |
261 | fPRFOn = p.fPRFOn; | |
262 | fNTimeBins = p.fNTimeBins; | |
263 | fNTBoverwriteOCDB = p.fNTBoverwriteOCDB; | |
264 | ||
265 | Int_t iBin = 0; | |
266 | ||
267 | if (fTRFsmp) { | |
268 | delete[] fTRFsmp; | |
269 | } | |
270 | fTRFsmp = new Float_t[fTRFbin]; | |
271 | for (iBin = 0; iBin < fTRFbin; iBin++) { | |
272 | fTRFsmp[iBin] = ((AliTRDSimParam &) p).fTRFsmp[iBin]; | |
273 | } | |
274 | ||
275 | if (fCTsmp) { | |
276 | delete[] fCTsmp; | |
277 | } | |
278 | fCTsmp = new Float_t[fTRFbin]; | |
279 | for (iBin = 0; iBin < fTRFbin; iBin++) { | |
280 | fCTsmp[iBin] = ((AliTRDSimParam &) p).fCTsmp[iBin]; | |
281 | } | |
282 | ||
3551db50 | 283 | return *this; |
4d18a639 | 284 | |
3551db50 | 285 | } |
286 | ||
287 | //_____________________________________________________________________________ | |
288 | void AliTRDSimParam::Copy(TObject &p) const | |
289 | { | |
290 | // | |
291 | // Copy function | |
292 | // | |
293 | ||
b43a3e17 | 294 | AliTRDSimParam *target = dynamic_cast<AliTRDSimParam *> (&p); |
4d18a639 | 295 | if (!target) { |
3551db50 | 296 | return; |
4d18a639 | 297 | } |
3551db50 | 298 | |
299 | target->fGasGain = fGasGain; | |
3551db50 | 300 | target->fNoise = fNoise; |
4d18a639 | 301 | target->fChipGain = fChipGain; |
3551db50 | 302 | target->fADCoutRange = fADCoutRange; |
303 | target->fADCinRange = fADCinRange; | |
3551db50 | 304 | target->fADCbaseline = fADCbaseline; |
3551db50 | 305 | target->fDiffusionOn = fDiffusionOn; |
3551db50 | 306 | target->fElAttachOn = fElAttachOn; |
307 | target->fElAttachProp = fElAttachProp; | |
3551db50 | 308 | target->fTRFOn = fTRFOn; |
3551db50 | 309 | target->fTRFbin = fTRFbin; |
310 | target->fTRFlo = fTRFlo; | |
311 | target->fTRFhi = fTRFhi; | |
312 | target->fTRFwid = fTRFwid; | |
3551db50 | 313 | target->fCTOn = fCTOn; |
4d18a639 | 314 | target->fPadCoupling = fPadCoupling; |
315 | target->fTimeCoupling = fTimeCoupling; | |
316 | target->fPRFOn = fPRFOn; | |
966f6939 | 317 | target->fNTimeBins = fNTimeBins; |
876a928e | 318 | target->fNTBoverwriteOCDB = fNTBoverwriteOCDB; |
4d18a639 | 319 | |
320 | if (target->fTRFsmp) { | |
321 | delete[] target->fTRFsmp; | |
322 | } | |
323 | target->fTRFsmp = new Float_t[fTRFbin]; | |
324 | for (Int_t iBin = 0; iBin < fTRFbin; iBin++) { | |
325 | target->fTRFsmp[iBin] = fTRFsmp[iBin]; | |
326 | } | |
327 | ||
328 | if (target->fCTsmp) { | |
3551db50 | 329 | delete[] target->fCTsmp; |
4d18a639 | 330 | } |
3551db50 | 331 | target->fCTsmp = new Float_t[fTRFbin]; |
332 | for (Int_t iBin = 0; iBin < fTRFbin; iBin++) { | |
333 | target->fCTsmp[iBin] = fCTsmp[iBin]; | |
334 | } | |
335 | ||
3551db50 | 336 | } |
337 | ||
338 | //_____________________________________________________________________________ | |
339 | void AliTRDSimParam::ReInit() | |
340 | { | |
341 | // | |
342 | // Reinitializes the parameter class after a change | |
343 | // | |
344 | ||
a076fc2f | 345 | if (AliTRDCommonParam::Instance()->IsXenon()) { |
f2979d08 | 346 | // The range and the binwidth for the sampled TRF |
347 | fTRFbin = 200; | |
348 | // Start 0.2 mus before the signal | |
349 | fTRFlo = -0.4; | |
350 | // End the maximum drift time after the signal | |
351 | fTRFhi = 3.58; | |
0a17cc30 | 352 | // Standard gas gain |
353 | fGasGain = 4000.0; | |
f2979d08 | 354 | } |
a076fc2f | 355 | else if (AliTRDCommonParam::Instance()->IsArgon()) { |
f2979d08 | 356 | // The range and the binwidth for the sampled TRF |
0a17cc30 | 357 | fTRFbin = 50; |
f2979d08 | 358 | // Start 0.2 mus before the signal |
0a17cc30 | 359 | fTRFlo = 0.02; |
f2979d08 | 360 | // End the maximum drift time after the signal |
0a17cc30 | 361 | fTRFhi = 1.98; |
362 | // Higher gas gain | |
363 | fGasGain = 8000.0; | |
f2979d08 | 364 | } |
365 | else { | |
366 | AliFatal("Not a valid gas mixture!"); | |
367 | exit(1); | |
368 | } | |
3551db50 | 369 | fTRFwid = (fTRFhi - fTRFlo) / ((Float_t) fTRFbin); |
370 | ||
371 | // Create the sampled TRF | |
372 | SampleTRF(); | |
3551db50 | 373 | |
4d18a639 | 374 | } |
915f999b | 375 | |
3551db50 | 376 | //_____________________________________________________________________________ |
377 | void AliTRDSimParam::SampleTRF() | |
378 | { | |
379 | // | |
915f999b | 380 | // Samples the new time response function. |
3551db50 | 381 | // |
382 | ||
f2979d08 | 383 | Int_t ipasa = 0; |
915f999b | 384 | |
f2979d08 | 385 | // Xenon |
386 | // From Antons measurements with Fe55 source, adjusted by C. Lippmann. | |
387 | // time bins are -0.4, -0.38, -0.36, ...., 3.54, 3.56, 3.58 microseconds | |
388 | const Int_t kNpasa = 200; // kNpasa should be equal to fTRFbin! | |
389 | Float_t xtalk[kNpasa]; | |
390 | Float_t signal[kNpasa] = { 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 | |
391 | , 0.0002, 0.0007, 0.0026, 0.0089, 0.0253, 0.0612, 0.1319 | |
392 | , 0.2416, 0.3913, 0.5609, 0.7295, 0.8662, 0.9581, 1.0000 | |
393 | , 0.9990, 0.9611, 0.8995, 0.8269, 0.7495, 0.6714, 0.5987 | |
394 | , 0.5334, 0.4756, 0.4249, 0.3811, 0.3433, 0.3110, 0.2837 | |
395 | , 0.2607, 0.2409, 0.2243, 0.2099, 0.1974, 0.1868, 0.1776 | |
396 | , 0.1695, 0.1627, 0.1566, 0.1509, 0.1457, 0.1407, 0.1362 | |
397 | , 0.1317, 0.1274, 0.1233, 0.1196, 0.1162, 0.1131, 0.1102 | |
398 | , 0.1075, 0.1051, 0.1026, 0.1004, 0.0979, 0.0956, 0.0934 | |
399 | , 0.0912, 0.0892, 0.0875, 0.0858, 0.0843, 0.0829, 0.0815 | |
400 | , 0.0799, 0.0786, 0.0772, 0.0757, 0.0741, 0.0729, 0.0718 | |
401 | , 0.0706, 0.0692, 0.0680, 0.0669, 0.0655, 0.0643, 0.0630 | |
402 | , 0.0618, 0.0607, 0.0596, 0.0587, 0.0576, 0.0568, 0.0558 | |
403 | , 0.0550, 0.0541, 0.0531, 0.0522, 0.0513, 0.0505, 0.0497 | |
404 | , 0.0490, 0.0484, 0.0474, 0.0465, 0.0457, 0.0449, 0.0441 | |
405 | , 0.0433, 0.0425, 0.0417, 0.0410, 0.0402, 0.0395, 0.0388 | |
406 | , 0.0381, 0.0374, 0.0368, 0.0361, 0.0354, 0.0348, 0.0342 | |
407 | , 0.0336, 0.0330, 0.0324, 0.0318, 0.0312, 0.0306, 0.0301 | |
408 | , 0.0296, 0.0290, 0.0285, 0.0280, 0.0275, 0.0270, 0.0265 | |
409 | , 0.0260, 0.0256, 0.0251, 0.0246, 0.0242, 0.0238, 0.0233 | |
410 | , 0.0229, 0.0225, 0.0221, 0.0217, 0.0213, 0.0209, 0.0206 | |
411 | , 0.0202, 0.0198, 0.0195, 0.0191, 0.0188, 0.0184, 0.0181 | |
412 | , 0.0178, 0.0175, 0.0171, 0.0168, 0.0165, 0.0162, 0.0159 | |
413 | , 0.0157, 0.0154, 0.0151, 0.0148, 0.0146, 0.0143, 0.0140 | |
414 | , 0.0138, 0.0135, 0.0133, 0.0131, 0.0128, 0.0126, 0.0124 | |
415 | , 0.0121, 0.0119, 0.0120, 0.0115, 0.0113, 0.0111, 0.0109 | |
416 | , 0.0107, 0.0105, 0.0103, 0.0101, 0.0100, 0.0098, 0.0096 | |
417 | , 0.0094, 0.0092, 0.0091, 0.0089, 0.0088, 0.0086, 0.0084 | |
418 | , 0.0083, 0.0081, 0.0080, 0.0078 }; | |
419 | signal[0] = 0.0; | |
420 | signal[1] = 0.0; | |
421 | signal[2] = 0.0; | |
4d18a639 | 422 | // With undershoot, positive peak corresponds to ~3% of the main signal: |
f2979d08 | 423 | for (ipasa = 3; ipasa < kNpasa; ipasa++) { |
4d18a639 | 424 | xtalk[ipasa] = 0.2 * (signal[ipasa-2] - signal[ipasa-3]); |
3551db50 | 425 | } |
4d18a639 | 426 | xtalk[0] = 0.0; |
427 | xtalk[1] = 0.0; | |
428 | xtalk[2] = 0.0; | |
915f999b | 429 | |
f2979d08 | 430 | // Argon |
431 | // Ar measurement with Fe55 source by Anton | |
432 | // time bins are 0.02, 0.06, 0.10, ...., 1.90, 1.94, 1.98 microseconds | |
433 | const Int_t kNpasaAr = 50; | |
434 | Float_t xtalkAr[kNpasaAr]; | |
435 | Float_t signalAr[kNpasaAr] = { -0.01, 0.01, 0.00, 0.00, 0.01 | |
436 | , -0.01, 0.01, 2.15, 22.28, 55.53 | |
437 | , 68.52, 58.21, 40.92, 27.12, 18.49 | |
438 | , 13.42, 10.48, 8.67, 7.49, 6.55 | |
439 | , 5.71, 5.12, 4.63, 4.22, 3.81 | |
440 | , 3.48, 3.20, 2.94, 2.77, 2.63 | |
441 | , 2.50, 2.37, 2.23, 2.13, 2.03 | |
442 | , 1.91, 1.83, 1.75, 1.68, 1.63 | |
443 | , 1.56, 1.49, 1.50, 1.49, 1.29 | |
444 | , 1.19, 1.21, 1.21, 1.20, 1.10 }; | |
445 | // Normalization to maximum | |
93ff78d6 | 446 | for (ipasa = 0; ipasa < kNpasaAr; ipasa++) { |
f2979d08 | 447 | signalAr[ipasa] /= 68.52; |
448 | } | |
449 | signalAr[0] = 0.0; | |
450 | signalAr[1] = 0.0; | |
451 | signalAr[2] = 0.0; | |
452 | // With undershoot, positive peak corresponds to ~3% of the main signal: | |
93ff78d6 | 453 | for (ipasa = 3; ipasa < kNpasaAr; ipasa++) { |
f2979d08 | 454 | xtalkAr[ipasa] = 0.2 * (signalAr[ipasa-2] - signalAr[ipasa-3]); |
455 | } | |
456 | xtalkAr[0] = 0.0; | |
457 | xtalkAr[1] = 0.0; | |
458 | xtalkAr[2] = 0.0; | |
4d18a639 | 459 | |
460 | if (fTRFsmp) { | |
461 | delete [] fTRFsmp; | |
462 | } | |
3551db50 | 463 | fTRFsmp = new Float_t[fTRFbin]; |
4d18a639 | 464 | |
465 | if (fCTsmp) { | |
466 | delete [] fCTsmp; | |
467 | } | |
3551db50 | 468 | fCTsmp = new Float_t[fTRFbin]; |
469 | ||
a076fc2f | 470 | if (AliTRDCommonParam::Instance()->IsXenon()) { |
0a17cc30 | 471 | if (fTRFbin != kNpasa) { |
472 | AliError("Array mismatch (xenon)\n\n"); | |
473 | } | |
474 | } | |
a076fc2f | 475 | else if (AliTRDCommonParam::Instance()->IsArgon()) { |
0a17cc30 | 476 | if (fTRFbin != kNpasaAr) { |
477 | AliError("Array mismatch (argon)\n\n"); | |
478 | } | |
479 | } | |
480 | ||
3551db50 | 481 | for (Int_t iBin = 0; iBin < fTRFbin; iBin++) { |
a076fc2f | 482 | if (AliTRDCommonParam::Instance()->IsXenon()) { |
f2979d08 | 483 | fTRFsmp[iBin] = signal[iBin]; |
484 | fCTsmp[iBin] = xtalk[iBin]; | |
485 | } | |
a076fc2f | 486 | else if (AliTRDCommonParam::Instance()->IsArgon()) { |
f2979d08 | 487 | fTRFsmp[iBin] = signalAr[iBin]; |
488 | fCTsmp[iBin] = xtalkAr[iBin]; | |
489 | } | |
3551db50 | 490 | } |
491 | ||
492 | } | |
493 | ||
494 | //_____________________________________________________________________________ | |
495 | Double_t AliTRDSimParam::TimeResponse(Double_t time) const | |
496 | { | |
497 | // | |
498 | // Applies the preamp shaper time response | |
915f999b | 499 | // (We assume a signal rise time of 0.2us = fTRFlo/2. |
3551db50 | 500 | // |
501 | ||
a076fc2f | 502 | Double_t rt = (time - .5*fTRFlo) / fTRFwid; |
503 | Int_t iBin = (Int_t) rt; | |
504 | Double_t dt = rt - iBin; | |
7093fb4f | 505 | if ((iBin >= 0) && (iBin+1 < fTRFbin)) { |
506 | return fTRFsmp[iBin] + (fTRFsmp[iBin+1] - fTRFsmp[iBin])*dt; | |
a076fc2f | 507 | } |
508 | else { | |
509 | return 0.0; | |
510 | } | |
511 | ||
3551db50 | 512 | } |
513 | ||
514 | //_____________________________________________________________________________ | |
515 | Double_t AliTRDSimParam::CrossTalk(Double_t time) const | |
516 | { | |
517 | // | |
518 | // Applies the pad-pad capacitive cross talk | |
519 | // | |
520 | ||
a076fc2f | 521 | Double_t rt = (time - fTRFlo) / fTRFwid; |
522 | Int_t iBin = (Int_t) rt; | |
523 | Double_t dt = rt - iBin; | |
7093fb4f | 524 | if ((iBin >= 0) && (iBin+1 < fTRFbin)) { |
525 | return fCTsmp[iBin] + (fCTsmp[iBin+1] - fCTsmp[iBin])*dt; | |
a076fc2f | 526 | } |
527 | else { | |
528 | return 0.0; | |
529 | } | |
530 | ||
3551db50 | 531 | } |