]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TRD/AliTRDparameter.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Functions renamed to get a prefix PHOS
[u/mrichter/AliRoot.git] / TRD / AliTRDparameter.cxx
CommitLineData
17b26de4 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/*
17$Log$
5443e65e 18Revision 1.5 2002/04/30 08:30:40 cblume
19gAlice now only read by AliRunDigitizer. Therefore it is just deleted in AliTRDmerge.C
20
c57e2264 21Revision 1.4 2002/04/12 12:13:23 cblume
22Add Jiris changes
23
8e64dd77 24Revision 1.3 2002/03/28 14:59:07 cblume
25Coding conventions
26
0a29d0f1 27Revision 1.2 2002/03/28 10:00:36 hristov
28Some additional initialisation
29
67689223 30Revision 1.1 2002/03/25 20:01:18 cblume
31Introduce parameter class
32
17b26de4 33*/
34
35///////////////////////////////////////////////////////////////////////////////
36// //
37// TRD parameter class //
38// //
39///////////////////////////////////////////////////////////////////////////////
40
41#include "AliRun.h"
42#include "AliMagF.h"
43
44#include "AliTRD.h"
45#include "AliTRDparameter.h"
46#include "AliTRDgeometryFull.h"
47
48ClassImp(AliTRDparameter)
49
50//_____________________________________________________________________________
51AliTRDparameter::AliTRDparameter():TNamed()
52{
53 //
54 // AliTRDparameter default constructor
55 //
56
57 fGeo = 0;
58 fPRFsmp = 0;
59 fTRFsmp = 0;
60 fCTsmp = 0;
61 fGasGain = 0.0;
62 fNoise = 0.0;
63 fChipGain = 0.0;
64 fADCoutRange = 0.0;
65 fADCinRange = 0.0;
66 fADCthreshold = 0;
67 fDiffusionOn = 0;
68 fDiffusionT = 0.0;
69 fDiffusionL = 0.0;
70 fElAttachOn = 0;
71 fElAttachProp = 0.0;
72 fExBOn = 0;
73 fOmegaTau = 0.0;
74 fPRFOn = 0;
75 fTRFOn = 0;
76 fCTOn = 0;
77 fTCOn = 0;
78 fDriftVelocity = 0.0;
79 fPadCoupling = 0.0;
80 fTimeCoupling = 0.0;
81 fTimeBinWidth = 0.0;
82 fField = 0.0;
83 fTiltingAngle = 0.0;
84 fPRFbin = 0;
85 fPRFlo = 0.0;
86 fPRFhi = 0.0;
87 fPRFwid = 0.0;
88 fPRFpad = 0;
89 fTRFbin = 0;
90 fTRFlo = 0.0;
91 fTRFhi = 0.0;
92 fTRFwid = 0.0;
93 fTCnexp = 0;
94
95 fLUTOn = 0;
96 fLUT = 0;
97 fClusMaxThresh = 0;
98 fClusSigThresh = 0;
99
100}
101
102//_____________________________________________________________________________
103AliTRDparameter::AliTRDparameter(const Text_t *name, const Text_t *title)
104 :TNamed(name,title)
105{
106 //
107 // AliTRDparameter constructor
108 //
109
110 fGeo = new AliTRDgeometryFull();
67689223 111 fPRFsmp = 0;
112 fTRFsmp = 0;
113 fCTsmp = 0;
114 fGasGain = 0.0;
115 fNoise = 0.0;
116 fChipGain = 0.0;
117 fADCoutRange = 0.0;
118 fADCinRange = 0.0;
119 fADCthreshold = 0;
120 fDiffusionOn = 0;
121 fDiffusionT = 0.0;
122 fDiffusionL = 0.0;
123 fElAttachOn = 0;
124 fElAttachProp = 0.0;
125 fExBOn = 0;
126 fOmegaTau = 0.0;
127 fPRFOn = 0;
128 fTRFOn = 0;
129 fCTOn = 0;
130 fTCOn = 0;
131 fDriftVelocity = 0.0;
132 fPadCoupling = 0.0;
133 fTimeCoupling = 0.0;
134 fTimeBinWidth = 0.0;
135 fField = 0.0;
136 fTiltingAngle = 0.0;
137 fPRFbin = 0;
138 fPRFlo = 0.0;
139 fPRFhi = 0.0;
140 fPRFwid = 0.0;
141 fPRFpad = 0;
142 fTRFbin = 0;
143 fTRFlo = 0.0;
144 fTRFhi = 0.0;
145 fTRFwid = 0.0;
146 fTCnexp = 0;
147
148 fLUTOn = 0;
149 fLUT = 0;
150 fClusMaxThresh = 0;
151 fClusSigThresh = 0;
17b26de4 152
153 Init();
154
155}
156
157
158//_____________________________________________________________________________
159AliTRDparameter::AliTRDparameter(const AliTRDparameter &p)
160{
161 //
162 // AliTRDparameter copy constructor
163 //
164
165 ((AliTRDparameter &) p).Copy(*this);
166
167}
168
169///_____________________________________________________________________________
170AliTRDparameter::~AliTRDparameter()
171{
172 //
173 // AliTRDparameter destructor
174 //
175
176 if (fTRFsmp) {
177 delete [] fTRFsmp;
178 fTRFsmp = 0;
179 }
180
181 if (fPRFsmp) {
182 delete [] fPRFsmp;
183 fPRFsmp = 0;
184 }
185
186 if (fCTsmp) {
187 delete [] fCTsmp;
188 fCTsmp = 0;
189 }
190
191 if (fLUT) {
192 delete [] fLUT;
193 fLUT = 0;
194 }
195
196 if (fGeo) {
197 delete fGeo;
198 fGeo = 0;
199 }
200
201}
202
203//_____________________________________________________________________________
204AliTRDparameter &AliTRDparameter::operator=(const AliTRDparameter &p)
205{
206 //
207 // Assignment operator
208 //
209
210 if (this != &p) ((AliTRDparameter &) p).Copy(*this);
211 return *this;
212
213}
214
215//_____________________________________________________________________________
216void AliTRDparameter::Copy(TObject &p)
217{
218 //
219 // Copy function
220 //
221
222 Int_t iBin;
223
224 ((AliTRDparameter &) p).fGasGain = fGasGain;
225 ((AliTRDparameter &) p).fNoise = fNoise;
226 ((AliTRDparameter &) p).fChipGain = fChipGain;
227 ((AliTRDparameter &) p).fADCoutRange = fADCoutRange;
228 ((AliTRDparameter &) p).fADCinRange = fADCinRange;
229 ((AliTRDparameter &) p).fADCthreshold = fADCthreshold;
230 ((AliTRDparameter &) p).fDiffusionOn = fDiffusionOn;
231 ((AliTRDparameter &) p).fDiffusionT = fDiffusionT;
232 ((AliTRDparameter &) p).fDiffusionL = fDiffusionL;
233 ((AliTRDparameter &) p).fElAttachOn = fElAttachOn;
234 ((AliTRDparameter &) p).fElAttachProp = fElAttachProp;
235 ((AliTRDparameter &) p).fExBOn = fExBOn;
236 ((AliTRDparameter &) p).fOmegaTau = fOmegaTau;
237 ((AliTRDparameter &) p).fLorentzFactor = fLorentzFactor;
238 ((AliTRDparameter &) p).fDriftVelocity = fDriftVelocity;
239 ((AliTRDparameter &) p).fPadCoupling = fPadCoupling;
240 ((AliTRDparameter &) p).fTimeCoupling = fTimeCoupling;
241 ((AliTRDparameter &) p).fTimeBinWidth = fTimeBinWidth;
242 ((AliTRDparameter &) p).fField = fField;
243 ((AliTRDparameter &) p).fPRFOn = fPRFOn;
244 ((AliTRDparameter &) p).fTRFOn = fTRFOn;
245 ((AliTRDparameter &) p).fCTOn = fCTOn;
246 ((AliTRDparameter &) p).fTCOn = fTCOn;
247 ((AliTRDparameter &) p).fTiltingAngle = fTiltingAngle;
248 ((AliTRDparameter &) p).fPRFbin = fPRFbin;
249 ((AliTRDparameter &) p).fPRFlo = fPRFlo;
250 ((AliTRDparameter &) p).fPRFhi = fPRFhi;
251 ((AliTRDparameter &) p).fPRFwid = fPRFwid;
252 ((AliTRDparameter &) p).fPRFpad = fPRFpad;
253 if (((AliTRDparameter &) p).fPRFsmp) delete [] ((AliTRDparameter &) p).fPRFsmp;
254 ((AliTRDparameter &) p).fPRFsmp = new Float_t[fPRFbin];
255 for (iBin = 0; iBin < fPRFbin; iBin++) {
256 ((AliTRDparameter &) p).fPRFsmp[iBin] = fPRFsmp[iBin];
257 }
258 ((AliTRDparameter &) p).fTRFbin = fTRFbin;
259 ((AliTRDparameter &) p).fTRFlo = fTRFlo;
260 ((AliTRDparameter &) p).fTRFhi = fTRFhi;
261 ((AliTRDparameter &) p).fTRFwid = fTRFwid;
262 if (((AliTRDparameter &) p).fTRFsmp) delete [] ((AliTRDparameter &) p).fTRFsmp;
263 ((AliTRDparameter &) p).fTRFsmp = new Float_t[fTRFbin];
264 for (iBin = 0; iBin < fTRFbin; iBin++) {
265 ((AliTRDparameter &) p).fTRFsmp[iBin] = fTRFsmp[iBin];
266 }
267 if (((AliTRDparameter &) p).fCTsmp) delete [] ((AliTRDparameter &) p).fCTsmp;
268 ((AliTRDparameter &) p).fCTsmp = new Float_t[fTRFbin];
269 for (iBin = 0; iBin < fTRFbin; iBin++) {
270 ((AliTRDparameter &) p).fCTsmp[iBin] = fCTsmp[iBin];
271 }
272 ((AliTRDparameter &) p).fTCnexp = fTCnexp;
273
274 ((AliTRDparameter &) p).fLUTOn = fLUTOn;
275 ((AliTRDparameter &) p).fLUTbin = fLUTbin;
276 if (((AliTRDparameter &) p).fLUT) delete [] ((AliTRDparameter &) p).fLUT;
277 ((AliTRDparameter &) p).fLUT = new Float_t[fLUTbin];
278 for (iBin = 0; iBin < fLUTbin; iBin++) {
279 ((AliTRDparameter &) p).fLUT[iBin] = fLUT[iBin];
280 }
281 ((AliTRDparameter &) p).fClusMaxThresh = fClusMaxThresh;
282 ((AliTRDparameter &) p).fClusSigThresh = fClusSigThresh;
283
284}
285
286//_____________________________________________________________________________
287void AliTRDparameter::Init()
288{
289 //
290 // Initializes the parameter
291 //
292 // The maximum number of pads
293 // and the position of pad 0,0,0
294 //
295 // chambers seen from the top:
296 // +----------------------------+
297 // | |
298 // | | ^
299 // | | rphi|
300 // | | |
301 // |0 | |
302 // +----------------------------+ +------>
303 // z
304 // chambers seen from the side: ^
305 // +----------------------------+ drift|
306 // |0 | |
307 // | | |
308 // +----------------------------+ +------>
309 // z
310 //
311 // IMPORTANT: time bin 0 is now the first one in the drift region
312 // closest to the readout !!!
313 //
314
315 //
316 // ----------------------------------------------------------------------------
317 // The pad definition
318 // ----------------------------------------------------------------------------
319 //
320
321 // The pad size in column direction (rphi-direction)
322 SetColPadSize(0,0.65);
323 SetColPadSize(1,0.68);
324 SetColPadSize(2,0.71);
325 SetColPadSize(3,0.74);
326 SetColPadSize(4,0.77);
327 SetColPadSize(5,0.80);
328
329 // The pad row (z-direction)
330 SetNRowPad();
331
332 // The number of time bins. Default is 100 ns timbin size
333 SetNTimeBin(15);
334
335 // Additional time bins before and after the drift region.
336 // Default is to only sample the drift region
337 SetExpandTimeBin(0,0);
338
339 //
340 // ----------------------------------------------------------------------------
341 // The digitization parameter
342 // ----------------------------------------------------------------------------
343 //
344
345 // The default parameter for the digitization
346 fGasGain = 4500.;
347 fChipGain = 12.4;
348 fNoise = 1000.;
349 fADCoutRange = 1023.; // 10-bit ADC
350 fADCinRange = 1000.; // 1V input range
351 fADCthreshold = 1;
352
353 // The drift velocity (cm / mus)
354 fDriftVelocity = 1.5;
355
356 // Diffusion on
357 fDiffusionOn = 1;
358
359 // E x B effects
360 fExBOn = 0;
361
362 // Propability for electron attachment
363 fElAttachOn = 0;
364 fElAttachProp = 0.0;
365
366 // The pad response function
367 fPRFOn = 1;
368
369 // The time response function
370 fTRFOn = 1;
371
372 // The cross talk
373 fCTOn = 0;
374
375 // The tail cancelation
376 fTCOn = 1;
377
378 // The number of exponentials
379 fTCnexp = 2;
380
381 // The pad coupling factor (same number as for the TPC)
382 fPadCoupling = 0.5;
383
384 // The time coupling factor (same number as for the TPC)
385 fTimeCoupling = 0.4;
386
387 // The tilting angle for the readout pads
5443e65e 388 //SetTiltingAngle(5.0);
389 SetTiltingAngle(0.0);
17b26de4 390
391 // The magnetic field strength in Tesla
c57e2264 392 //fField = 0.2 * gAlice->Field()->Factor();
393 fField = 0.4;
17b26de4 394
395 //
396 // ----------------------------------------------------------------------------
397 // The clusterization parameter
398 // ----------------------------------------------------------------------------
399 //
400
401 // The default parameter for the clustering
402 fClusMaxThresh = 3;
403 fClusSigThresh = 1;
404
405 // Use the LUT
406 fLUTOn = 1;
407
408 ReInit();
409
410}
411
412//_____________________________________________________________________________
413void AliTRDparameter::ReInit()
414{
415 //
416 // Reinitializes the parameter class after a change
417 //
418
419 // Calculate the time bin width in ns
420 fTimeBinWidth = fTimeBinSize / fDriftVelocity * 1000.0;
421
422 // The range and the binwidth for the sampled TRF
423 fTRFbin = 100;
424 // Start 0.2 mus before the signal
425 fTRFlo = -0.2 * fDriftVelocity;
426 // End the maximum driftlength after the signal
427 fTRFhi = AliTRDgeometry::DrThick()
428 + fTimeAfter * fTimeBinSize;
429 fTRFwid = (fTRFhi - fTRFlo) / ((Float_t) fTRFbin);
430
431 // Transverse and longitudinal diffusion coefficients (Xe/CO2)
432 fDiffusionT = GetDiffusionT(fDriftVelocity,fField);
433 fDiffusionL = GetDiffusionL(fDriftVelocity,fField);
434
435 // omega * tau.= tan(Lorentz-angle)
436 fOmegaTau = GetOmegaTau(fDriftVelocity,fField);
437
438 // The Lorentz factor
439 if (fExBOn) {
440 fLorentzFactor = 1.0 / (1.0 + fOmegaTau*fOmegaTau);
441 }
442 else {
443 fLorentzFactor = 1.0;
444 }
445
446 // Create the sampled PRF
447 SamplePRF();
448
449 // Create the sampled TRF
450 SampleTRF();
451
452 // Create the LUT
453 FillLUT();
454
455}
456
457//_____________________________________________________________________________
458void AliTRDparameter::SetNRowPad(const Int_t p, const Int_t c, const Int_t npad)
459{
460 //
461 // Redefines the number of pads in raw direction for
462 // a given plane and chamber number
463 //
464
465 for (Int_t isect = 0; isect < AliTRDgeometry::Nsect(); isect++) {
466
467 fRowMax[p][c][isect] = npad;
468
469 fRowPadSize[p][c][isect] = (fGeo->GetChamberLength(p,c)
470 - 2.* AliTRDgeometry::RpadW())
471 / ((Float_t) npad);
472
473 }
474
475}
476
477//_____________________________________________________________________________
478void AliTRDparameter::SetNRowPad()
479{
480 //
481 // Defines the number of pads in row direction
482 //
483
484 Int_t isect;
485 Int_t icham;
486 Int_t iplan;
487
488 Int_t rowMax[kNplan][kNcham] = { { 16, 16, 12, 16, 16 }
489 , { 16, 16, 12, 16, 16 }
490 , { 16, 16, 12, 16, 16 }
491 , { 16, 16, 12, 16, 16 }
492 , { 14, 16, 12, 16, 14 }
493 , { 13, 16, 12, 16, 13 } };
494
495 Float_t rpadW = AliTRDgeometry::RpadW();
496
497 for (isect = 0; isect < kNsect; isect++) {
498 for (icham = 0; icham < kNcham; icham++) {
499 for (iplan = 0; iplan < kNplan; iplan++) {
500
501 fRowMax[iplan][icham][isect] = rowMax[iplan][icham];
502
503 fRowPadSize[iplan][icham][isect] = (fGeo->GetChamberLength(iplan,icham)
504 - 2.*rpadW)
505 / ((Float_t) rowMax[iplan][icham]);
506
507 Float_t row0 = rpadW - fGeo->GetChamberLength(iplan,0)
508 - fGeo->GetChamberLength(iplan,1)
509 - fGeo->GetChamberLength(iplan,2) / 2.;
510 for (Int_t ic = 0; ic < icham; ic++) {
511 row0 += fGeo->GetChamberLength(iplan,icham);
512 }
513 fRow0[iplan][icham][isect] = row0;
514
515 }
516 }
517 }
518
519}
520
521//_____________________________________________________________________________
522void AliTRDparameter::SetColPadSize(const Int_t p, const Float_t s)
523{
524 //
525 // Redefines the pad size in column direction
526 //
527
528 Float_t cpadW = AliTRDgeometry::CpadW();
529
530 fColPadSize[p] = s;
531 fCol0[p] = - fGeo->GetChamberWidth(p)/2. + cpadW;
532 fColMax[p] = ((Int_t) ((fGeo->GetChamberWidth(p) - 2.*cpadW) / s));
533
534}
535
536//_____________________________________________________________________________
537void AliTRDparameter::SetNTimeBin(const Int_t nbin)
538{
539 //
540 // Redefines the number of time bins in the drift region.
541 // The time bin width is defined by the length of the
542 // drift region divided by <nbin>.
543 //
544
545 fTimeMax = nbin;
546 fTimeBinSize = AliTRDgeometry::DrThick() / ((Float_t) fTimeMax);
547 for (Int_t iplan = 0; iplan < AliTRDgeometry::Nplan(); iplan++) {
548 fTime0[iplan] = AliTRDgeometry::Rmin()
549 + AliTRDgeometry::CraHght()
550 + AliTRDgeometry::CdrHght()
551 + iplan * (AliTRDgeometry::Cheight()
552 + AliTRDgeometry::Cspace());
553 }
554
555}
556
557//_____________________________________________________________________________
0a29d0f1 558Float_t AliTRDparameter::CrossTalk(Float_t time) const
17b26de4 559{
560 //
561 // Applies the pad-pad capacitive cross talk
562 //
563
564 Int_t iBin = ((Int_t) ((time - fTRFlo) / fTRFwid));
565 if ((iBin >= 0) && (iBin < fTRFbin)) {
566 return fCTsmp[iBin];
567 }
568 else {
569 return 0.0;
570 }
571
572}
573
574//_____________________________________________________________________________
575Int_t AliTRDparameter::Diffusion(Float_t driftlength, Float_t *xyz)
576{
577 //
578 // Applies the diffusion smearing to the position of a single electron
579 //
580
581 Float_t driftSqrt = TMath::Sqrt(driftlength);
582 Float_t sigmaT = driftSqrt * fDiffusionT;
583 Float_t sigmaL = driftSqrt * fDiffusionL;
584 xyz[0] = gRandom->Gaus(xyz[0], sigmaL * fLorentzFactor);
585 xyz[1] = gRandom->Gaus(xyz[1], sigmaT * fLorentzFactor);
586 xyz[2] = gRandom->Gaus(xyz[2], sigmaT);
587
588 return 1;
589
590}
591
592//_____________________________________________________________________________
0a29d0f1 593Int_t AliTRDparameter::ExB(Float_t driftlength, Float_t *xyz) const
17b26de4 594{
595 //
596 // Applies E x B effects to the position of a single electron
597 //
598
599 xyz[0] = xyz[0];
600 xyz[1] = xyz[1] + fOmegaTau * driftlength;
601 xyz[2] = xyz[2];
602
603 return 1;
604
605}
606
607//_____________________________________________________________________________
608Int_t AliTRDparameter::PadResponse(Float_t signal, Float_t dist
0a29d0f1 609 , Int_t plane, Float_t *pad) const
17b26de4 610{
611 //
612 // Applies the pad response
613 //
614
615 const Int_t kNplan = AliTRDgeometry::kNplan;
616
617 Int_t iBin = ((Int_t) (( - dist - fPRFlo) / fPRFwid));
618 Int_t iOff = plane * fPRFbin;
619
620 Int_t iBin0 = iBin - fPRFpad + iOff;
621 Int_t iBin1 = iBin + iOff;
622 Int_t iBin2 = iBin + fPRFpad + iOff;
623
624 pad[0] = 0.0;
625 pad[1] = 0.0;
626 pad[2] = 0.0;
627 if ((iBin1 >= 0) && (iBin1 < (fPRFbin*kNplan))) {
628
629 if (iBin0 >= 0) {
630 pad[0] = signal * fPRFsmp[iBin0];
631 }
632 pad[1] = signal * fPRFsmp[iBin1];
633 if (iBin2 < (fPRFbin*kNplan)) {
634 pad[2] = signal * fPRFsmp[iBin2];
635 }
636
637 return 1;
638
639 }
640 else {
641
642 return 0;
643
644 }
645
646}
647
648//_____________________________________________________________________________
0a29d0f1 649Float_t AliTRDparameter::TimeResponse(Float_t time) const
17b26de4 650{
651 //
652 // Applies the preamp shaper time response
653 //
654
655 Int_t iBin = ((Int_t) ((time - fTRFlo) / fTRFwid));
656 if ((iBin >= 0) && (iBin < fTRFbin)) {
657 return fTRFsmp[iBin];
658 }
659 else {
660 return 0.0;
661 }
662
663}
664
665//_____________________________________________________________________________
666Float_t AliTRDparameter::Col0Tilted(Float_t col0, Float_t rowOffset
667 , Int_t plane)
668{
669 //
670 // Calculates col0 for tilted pads
671 //
672
673 Float_t diff = fTiltingAngle * rowOffset;
674 return (col0 + TMath::Power(-1.0,plane) * diff);
675
676}
677
678//_____________________________________________________________________________
679void AliTRDparameter::SampleTRF()
680{
681 //
682 // Samples the time response function
683 //
684 // New TRF from Venelin Angelov, simulated with CADENCE
685 // Pad-ground capacitance = 25 pF
686 // Pad-pad cross talk capacitance = 6 pF
687 //
688
689 Int_t ipos1;
690 Int_t ipos2;
691 Float_t diff;
692
693 const Int_t kNpasa = 252;
694
695 Float_t time[kNpasa] = { -0.220000, -0.210000, -0.200000, -0.190000
696 , -0.180000, -0.170000, -0.160000, -0.150000
697 , -0.140000, -0.130000, -0.120000, -0.110000
698 , -0.100000, -0.090000, -0.080000, -0.070000
699 , -0.060000, -0.050000, -0.040000, -0.030000
700 , -0.020000, -0.010000, -0.000000, 0.010000
701 , 0.020000, 0.030000, 0.040000, 0.050000
702 , 0.060000, 0.070000, 0.080000, 0.090000
703 , 0.100000, 0.110000, 0.120000, 0.130000
704 , 0.140000, 0.150000, 0.160000, 0.170000
705 , 0.180000, 0.190000, 0.200000, 0.210000
706 , 0.220000, 0.230000, 0.240000, 0.250000
707 , 0.260000, 0.270000, 0.280000, 0.290000
708 , 0.300000, 0.310000, 0.320000, 0.330000
709 , 0.340000, 0.350000, 0.360000, 0.370000
710 , 0.380000, 0.390000, 0.400000, 0.410000
711 , 0.420000, 0.430000, 0.440000, 0.450000
712 , 0.460000, 0.470000, 0.480000, 0.490000
713 , 0.500000, 0.510000, 0.520000, 0.530000
714 , 0.540000, 0.550000, 0.560000, 0.570000
715 , 0.580000, 0.590000, 0.600000, 0.610000
716 , 0.620000, 0.630000, 0.640000, 0.650000
717 , 0.660000, 0.670000, 0.680000, 0.690000
718 , 0.700000, 0.710000, 0.720000, 0.730000
719 , 0.740000, 0.750000, 0.760000, 0.770000
720 , 0.780000, 0.790000, 0.800000, 0.810000
721 , 0.820000, 0.830000, 0.840000, 0.850000
722 , 0.860000, 0.870000, 0.880000, 0.890000
723 , 0.900000, 0.910000, 0.920000, 0.930000
724 , 0.940000, 0.950000, 0.960000, 0.970000
725 , 0.980000, 0.990000, 1.000000, 1.010000
726 , 1.020000, 1.030000, 1.040000, 1.050000
727 , 1.060000, 1.070000, 1.080000, 1.090000
728 , 1.100000, 1.110000, 1.120000, 1.130000
729 , 1.140000, 1.150000, 1.160000, 1.170000
730 , 1.180000, 1.190000, 1.200000, 1.210000
731 , 1.220000, 1.230000, 1.240000, 1.250000
732 , 1.260000, 1.270000, 1.280000, 1.290000
733 , 1.300000, 1.310000, 1.320000, 1.330000
734 , 1.340000, 1.350000, 1.360000, 1.370000
735 , 1.380000, 1.390000, 1.400000, 1.410000
736 , 1.420000, 1.430000, 1.440000, 1.450000
737 , 1.460000, 1.470000, 1.480000, 1.490000
738 , 1.500000, 1.510000, 1.520000, 1.530000
739 , 1.540000, 1.550000, 1.560000, 1.570000
740 , 1.580000, 1.590000, 1.600000, 1.610000
741 , 1.620000, 1.630000, 1.640000, 1.650000
742 , 1.660000, 1.670000, 1.680000, 1.690000
743 , 1.700000, 1.710000, 1.720000, 1.730000
744 , 1.740000, 1.750000, 1.760000, 1.770000
745 , 1.780000, 1.790000, 1.800000, 1.810000
746 , 1.820000, 1.830000, 1.840000, 1.850000
747 , 1.860000, 1.870000, 1.880000, 1.890000
748 , 1.900000, 1.910000, 1.920000, 1.930000
749 , 1.940000, 1.950000, 1.960000, 1.970000
750 , 1.980000, 1.990000, 2.000000, 2.010000
751 , 2.020000, 2.030000, 2.040000, 2.050000
752 , 2.060000, 2.070000, 2.080000, 2.090000
753 , 2.100000, 2.110000, 2.120000, 2.130000
754 , 2.140000, 2.150000, 2.160000, 2.170000
755 , 2.180000, 2.190000, 2.200000, 2.210000
756 , 2.220000, 2.230000, 2.240000, 2.250000
757 , 2.260000, 2.270000, 2.280000, 2.290000 };
758
759 Float_t signal[kNpasa] = { 0.000000, 0.000000, 0.000000, 0.000000
760 , 0.000000, 0.000000, 0.000000, 0.000396
761 , 0.005096, 0.022877, 0.061891, 0.126614
762 , 0.215798, 0.324406, 0.444507, 0.566817
763 , 0.683465, 0.787089, 0.873159, 0.937146
764 , 0.979049, 0.999434, 1.000000, 0.983579
765 , 0.954134, 0.913364, 0.866365, 0.813703
766 , 0.759910, 0.706116, 0.653454, 0.603624
767 , 0.556625, 0.514156, 0.475085, 0.439977
768 , 0.408834, 0.380578, 0.355549, 0.333352
769 , 0.313647, 0.296093, 0.280351, 0.266195
770 , 0.253397, 0.241789, 0.231257, 0.221574
771 , 0.212627, 0.204417, 0.196772, 0.189581
772 , 0.182956, 0.176784, 0.171008, 0.165515
773 , 0.160419, 0.155606, 0.151076, 0.146716
774 , 0.142639, 0.138845, 0.135221, 0.131767
775 , 0.128482, 0.125368, 0.122424, 0.119592
776 , 0.116931, 0.114326, 0.111891, 0.109513
777 , 0.107248, 0.105096, 0.103058, 0.101019
778 , 0.099151, 0.097282, 0.095527, 0.093715
779 , 0.092129, 0.090544, 0.088958, 0.087429
780 , 0.086014, 0.084598, 0.083239, 0.081880
781 , 0.080634, 0.079388, 0.078143, 0.077010
782 , 0.075878, 0.074745, 0.073669, 0.072593
783 , 0.071574, 0.070612, 0.069649, 0.068686
784 , 0.067780, 0.066874, 0.066025, 0.065176
785 , 0.064326, 0.063533, 0.062684, 0.061948
786 , 0.061212, 0.060419, 0.059740, 0.059003
787 , 0.058324, 0.057644, 0.057022, 0.056342
788 , 0.055663, 0.055096, 0.054473, 0.053851
789 , 0.053284, 0.052718, 0.052152, 0.051585
790 , 0.051019, 0.050566, 0.050000, 0.049490
791 , 0.048981, 0.048528, 0.048018, 0.047508
792 , 0.047055, 0.046602, 0.046149, 0.045696
793 , 0.045300, 0.044904, 0.044451, 0.044054
794 , 0.043658, 0.043205, 0.042865, 0.042469
795 , 0.042072, 0.041733, 0.041336, 0.040997
796 , 0.040657, 0.040260, 0.039921, 0.039581
797 , 0.039241, 0.038958, 0.038618, 0.038335
798 , 0.037995, 0.037656, 0.037373, 0.037089
799 , 0.036806, 0.036467, 0.036183, 0.035900
800 , 0.035617, 0.035334, 0.035108, 0.034824
801 , 0.034541, 0.034315, 0.034032, 0.033805
802 , 0.033522, 0.033296, 0.033069, 0.032786
803 , 0.032559, 0.032333, 0.032106, 0.031880
804 , 0.031653, 0.031427, 0.031200, 0.030974
805 , 0.030804, 0.030578, 0.030351, 0.030125
806 , 0.029955, 0.029785, 0.029558, 0.029332
807 , 0.029162, 0.028992, 0.028766, 0.028596
808 , 0.028426, 0.028199, 0.028086, 0.027860
809 , 0.027746, 0.027633, 0.027463, 0.027293
810 , 0.027180, 0.027067, 0.026954, 0.026954
811 , 0.026840, 0.026727, 0.026727, 0.026614
812 , 0.026614, 0.026614, 0.026557, 0.026501
813 , 0.026501, 0.026501, 0.026501, 0.026501
814 , 0.026501, 0.026501, 0.026501, 0.026387
815 , 0.026387, 0.026387, 0.026387, 0.026387
816 , 0.026387, 0.026387, 0.026387, 0.026387
817 , 0.026387, 0.026387, 0.026387, 0.026387
818 , 0.026387, 0.026274, 0.026274, 0.026274
819 , 0.026274, 0.026274, 0.026274, 0.026274
820 , 0.026274, 0.026274, 0.026274, 0.026274
821 , 0.026274, 0.026274, 0.026274, 0.026161 };
822
823 Float_t xtalk[kNpasa] = { 0.000000, 0.000000, 0.000000, 0.000000
824 , 0.000000, 0.000000, 0.000000, 0.000113
825 , 0.000793, 0.003058, 0.007305, 0.013194
826 , 0.019706, 0.025821, 0.030634, 0.033465
827 , 0.034145, 0.032729, 0.029615, 0.025198
828 , 0.019989, 0.014496, 0.009003, 0.003964
829 , -0.000510, -0.004190, -0.007191, -0.009400
830 , -0.010872, -0.011835, -0.012288, -0.012288
831 , -0.012005, -0.011495, -0.010872, -0.010136
832 , -0.009343, -0.008607, -0.007871, -0.007191
833 , -0.006512, -0.005946, -0.005379, -0.004926
834 , -0.004473, -0.004077, -0.003737, -0.003398
835 , -0.003114, -0.002831, -0.002605, -0.002378
836 , -0.002208, -0.002039, -0.001869, -0.001699
837 , -0.001585, -0.001472, -0.001359, -0.001246
838 , -0.001132, -0.001019, -0.001019, -0.000906
839 , -0.000906, -0.000793, -0.000793, -0.000680
840 , -0.000680, -0.000680, -0.000566, -0.000566
841 , -0.000566, -0.000566, -0.000453, -0.000453
842 , -0.000453, -0.000453, -0.000453, -0.000453
843 , -0.000340, -0.000340, -0.000340, -0.000340
844 , -0.000340, -0.000340, -0.000340, -0.000340
845 , -0.000340, -0.000340, -0.000340, -0.000340
846 , -0.000340, -0.000227, -0.000227, -0.000227
847 , -0.000227, -0.000227, -0.000227, -0.000227
848 , -0.000227, -0.000227, -0.000227, -0.000227
849 , -0.000227, -0.000227, -0.000227, -0.000227
850 , -0.000227, -0.000227, -0.000227, -0.000227
851 , -0.000227, -0.000227, -0.000227, -0.000227
852 , -0.000227, -0.000227, -0.000227, -0.000227
853 , -0.000227, -0.000227, -0.000227, -0.000227
854 , -0.000227, -0.000227, -0.000227, -0.000227
855 , -0.000227, -0.000227, -0.000227, -0.000113
856 , -0.000113, -0.000113, -0.000113, -0.000113
857 , -0.000113, -0.000113, -0.000113, -0.000113
858 , -0.000113, -0.000113, -0.000113, -0.000113
859 , -0.000113, -0.000113, -0.000113, -0.000113
860 , -0.000113, -0.000113, -0.000113, -0.000113
861 , -0.000113, -0.000113, -0.000113, -0.000113
862 , -0.000113, -0.000113, -0.000113, -0.000113
863 , -0.000113, -0.000113, -0.000113, -0.000113
864 , -0.000113, -0.000113, -0.000113, -0.000113
865 , -0.000113, -0.000113, -0.000113, -0.000113
866 , -0.000113, -0.000113, -0.000113, -0.000113
867 , -0.000113, -0.000113, -0.000113, -0.000113
868 , -0.000113, -0.000113, -0.000113, -0.000113
869 , -0.000113, -0.000113, -0.000113, -0.000113
870 , -0.000113, -0.000113, -0.000113, -0.000113
871 , -0.000113, -0.000113, -0.000113, -0.000113
872 , -0.000113, -0.000113, -0.000113, -0.000113
873 , -0.000113, -0.000113, -0.000113, -0.000113
874 , -0.000113, -0.000113, -0.000113, -0.000113
875 , -0.000113, -0.000113, -0.000113, 0.000000
876 , 0.000000, 0.000000, 0.000000, 0.000000
877 , 0.000000, 0.000000, 0.000000, 0.000000
878 , 0.000000, 0.000000, 0.000000, 0.000000
879 , 0.000000, 0.000000, 0.000000, 0.000000
880 , 0.000000, 0.000000, 0.000000, 0.000000
881 , 0.000000, 0.000000, 0.000000, 0.000000
882 , 0.000000, 0.000000, 0.000000, 0.000000
883 , 0.000000, 0.000000, 0.000000, 0.000000
884 , 0.000000, 0.000000, 0.000000, 0.000000
885 , 0.000000, 0.000000, 0.000000, 0.000000 };
886
887 // increase CrossTalk to measurements
888 for (Int_t ipasa = 0; ipasa < kNpasa; ipasa++) {
889 xtalk[ipasa] *= 1.75;
890 }
891
892 if (fTRFsmp) delete [] fTRFsmp;
893 fTRFsmp = new Float_t[fTRFbin];
894 if (fCTsmp) delete [] fCTsmp;
895 fCTsmp = new Float_t[fTRFbin];
896
897 Float_t loTRF = TMath::Max(fTRFlo / fDriftVelocity,time[0]);
898 Float_t hiTRF = TMath::Min(fTRFhi / fDriftVelocity,time[kNpasa-1]);
899 Float_t binWidth = (hiTRF - loTRF) / ((Float_t) fTRFbin);
900
901 // Take the linear interpolation
902 for (Int_t iBin = 0; iBin < fTRFbin; iBin++) {
903
904 Float_t bin = (((Float_t) iBin) + 0.5) * binWidth + loTRF;
905 ipos1 = ipos2 = 0;
906 diff = 0;
907 do {
908 diff = bin - time[ipos2++];
909 } while (diff > 0);
910 ipos2--;
911 if (ipos2 >= kNpasa) ipos2 = kNpasa - 1;
912 ipos1 = ipos2 - 1;
913
914 fTRFsmp[iBin] = signal[ipos2]
915 + diff * (signal[ipos2] - signal[ipos1])
916 / ( time[ipos2] - time[ipos1]);
917
918 fCTsmp[iBin] = xtalk[ipos2]
919 + diff * (xtalk[ipos2] - xtalk[ipos1])
920 / ( time[ipos2] - time[ipos1]);
921
922 }
923
924}
925
926//_____________________________________________________________________________
927void AliTRDparameter::SamplePRF()
928{
929 //
930 // Samples the pad response function
931 //
932
933 const Int_t kNplan = AliTRDgeometry::kNplan;
934 const Int_t kPRFbin = 61;
935
936 Float_t prf[kNplan][kPRFbin] = { { 0.018570, 0.022270, 0.026710, 0.032010
937 , 0.038350, 0.045920, 0.054930, 0.065650
938 , 0.078370, 0.093420, 0.111150, 0.131940
939 , 0.156160, 0.184160, 0.216220, 0.252470
940 , 0.292860, 0.337030, 0.384330, 0.433750
941 , 0.484010, 0.533630, 0.581150, 0.625200
942 , 0.664710, 0.698860, 0.727130, 0.749230
943 , 0.765050, 0.774540, 0.777700, 0.774540
944 , 0.765050, 0.749230, 0.727130, 0.698860
945 , 0.664710, 0.625200, 0.581150, 0.533630
946 , 0.484010, 0.433750, 0.384330, 0.337030
947 , 0.292860, 0.252470, 0.216220, 0.184160
948 , 0.156160, 0.131940, 0.111150, 0.093420
949 , 0.078370, 0.065650, 0.054930, 0.045920
950 , 0.038350, 0.032010, 0.026710, 0.022270
951 , 0.018570 }
952 , { 0.015730, 0.019040, 0.023030, 0.027840
953 , 0.033650, 0.040650, 0.049060, 0.059160
954 , 0.071260, 0.085710, 0.102910, 0.123270
955 , 0.147240, 0.175220, 0.207590, 0.244540
956 , 0.286090, 0.331920, 0.381350, 0.433290
957 , 0.486290, 0.538710, 0.588870, 0.635280
958 , 0.676760, 0.712460, 0.741890, 0.764810
959 , 0.781150, 0.790930, 0.794180, 0.790930
960 , 0.781150, 0.764810, 0.741890, 0.712460
961 , 0.676760, 0.635280, 0.588870, 0.538710
962 , 0.486290, 0.433290, 0.381350, 0.331920
963 , 0.286090, 0.244540, 0.207590, 0.175220
964 , 0.147240, 0.123270, 0.102910, 0.085710
965 , 0.071260, 0.059160, 0.049060, 0.040650
966 , 0.033650, 0.027840, 0.023030, 0.019040
967 , 0.015730 }
968 , { 0.013330, 0.016270, 0.019850, 0.024210
969 , 0.029510, 0.035960, 0.043790, 0.053280
970 , 0.064740, 0.078580, 0.095190, 0.115070
971 , 0.138700, 0.166570, 0.199120, 0.236660
972 , 0.279260, 0.326660, 0.378140, 0.432540
973 , 0.488260, 0.543440, 0.596200, 0.644900
974 , 0.688240, 0.725380, 0.755840, 0.779470
975 , 0.796260, 0.806280, 0.809610, 0.806280
976 , 0.796260, 0.779470, 0.755840, 0.725380
977 , 0.688240, 0.644900, 0.596200, 0.543440
978 , 0.488260, 0.432540, 0.378140, 0.326660
979 , 0.279260, 0.236660, 0.199120, 0.166570
980 , 0.138700, 0.115070, 0.095190, 0.078580
981 , 0.064740, 0.053280, 0.043790, 0.035960
982 , 0.029510, 0.024210, 0.019850, 0.016270
983 , 0.013330 }
984 , { 0.011280, 0.013890, 0.017090, 0.021030
985 , 0.025870, 0.031800, 0.039060, 0.047940
986 , 0.058790, 0.071980, 0.087990, 0.107330
987 , 0.130550, 0.158220, 0.190850, 0.228870
988 , 0.272410, 0.321270, 0.374740, 0.431560
989 , 0.489960, 0.547870, 0.603180, 0.654080
990 , 0.699190, 0.737640, 0.769030, 0.793260
991 , 0.810410, 0.820620, 0.824010, 0.820620
992 , 0.810410, 0.793260, 0.769030, 0.737640
993 , 0.699190, 0.654080, 0.603180, 0.547870
994 , 0.489960, 0.431560, 0.374740, 0.321270
995 , 0.272410, 0.228870, 0.190850, 0.158220
996 , 0.130550, 0.107330, 0.087990, 0.071980
997 , 0.058790, 0.047940, 0.039060, 0.031800
998 , 0.025870, 0.021030, 0.017090, 0.013890
999 , 0.011280 }
1000 , { 0.009550, 0.011860, 0.014720, 0.018270
1001 , 0.022660, 0.028100, 0.034820, 0.043120
1002 , 0.053340, 0.065900, 0.081280, 0.100040
1003 , 0.122800, 0.150180, 0.182800, 0.221170
1004 , 0.265550, 0.315790, 0.371180, 0.430370
1005 , 0.491430, 0.552030, 0.609840, 0.662860
1006 , 0.709630, 0.749290, 0.781490, 0.806220
1007 , 0.823650, 0.834000, 0.837430, 0.834000
1008 , 0.823650, 0.806220, 0.781490, 0.749290
1009 , 0.709630, 0.662860, 0.609840, 0.552030
1010 , 0.491430, 0.430370, 0.371180, 0.315790
1011 , 0.265550, 0.221170, 0.182800, 0.150180
1012 , 0.122800, 0.100040, 0.081280, 0.065900
1013 , 0.053340, 0.043120, 0.034820, 0.028100
1014 , 0.022660, 0.018270, 0.014720, 0.011860
1015 , 0.009550 }
1016 , { 0.008080, 0.010120, 0.012670, 0.015860
1017 , 0.019840, 0.024820, 0.031030, 0.038760
1018 , 0.048370, 0.060300, 0.075040, 0.093200
1019 , 0.115430, 0.142450, 0.174980, 0.213610
1020 , 0.258720, 0.310250, 0.367480, 0.429010
1021 , 0.492690, 0.555950, 0.616210, 0.671280
1022 , 0.719600, 0.760350, 0.793250, 0.818380
1023 , 0.836020, 0.846460, 0.849920, 0.846460
1024 , 0.836020, 0.818380, 0.793250, 0.760350
1025 , 0.719600, 0.671280, 0.616210, 0.555950
1026 , 0.492690, 0.429010, 0.367480, 0.310250
1027 , 0.258720, 0.213610, 0.174980, 0.142450
1028 , 0.115430, 0.093200, 0.075040, 0.060300
1029 , 0.048370, 0.038760, 0.031030, 0.024820
1030 , 0.019840, 0.015860, 0.012670, 0.010120
1031 , 0.008080 } };
1032
1033 // More sampling precision with linear interpolation
1034 fPRFlo = -1.5;
1035 fPRFhi = 1.5;
1036 Float_t pad[kPRFbin];
1037 Int_t sPRFbin = kPRFbin;
1038 Float_t sPRFwid = (fPRFhi - fPRFlo) / ((Float_t) sPRFbin);
1039 for (Int_t iPad = 0; iPad < sPRFbin; iPad++) {
1040 pad[iPad] = ((Float_t) iPad + 0.5) * sPRFwid + fPRFlo;
1041 }
1042 fPRFbin = 500;
1043 fPRFwid = (fPRFhi - fPRFlo) / ((Float_t) fPRFbin);
1044 fPRFpad = ((Int_t) (1.0 / fPRFwid));
1045
1046 if (fPRFsmp) delete [] fPRFsmp;
1047 fPRFsmp = new Float_t[kNplan*fPRFbin];
1048
1049 Int_t ipos1;
1050 Int_t ipos2;
1051 Float_t diff;
1052
1053 for (Int_t iPla = 0; iPla < kNplan; iPla++) {
1054
1055 for (Int_t iBin = 0; iBin < fPRFbin; iBin++) {
1056
1057 Float_t bin = (((Float_t) iBin) + 0.5) * fPRFwid + fPRFlo;
1058 ipos1 = ipos2 = 0;
1059 diff = 0;
1060 do {
1061 diff = bin - pad[ipos2++];
8e64dd77 1062 } while ((diff > 0) && (ipos2 < kPRFbin));
1063 if (ipos2 == kPRFbin) {
1064 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2-1];
1065 }
1066 else if (ipos2 == 1) {
1067 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2-1];
1068 }
1069 else {
1070 ipos2--;
1071 if (ipos2 >= kPRFbin) ipos2 = kPRFbin - 1;
1072 ipos1 = ipos2 - 1;
1073 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2]
1074 + diff * (prf[iPla][ipos2] - prf[iPla][ipos1])
1075 / sPRFwid;
1076 }
17b26de4 1077
1078 }
1079 }
1080
1081}
1082
1083//_____________________________________________________________________________
1084void AliTRDparameter::FillLUT()
1085{
1086 //
1087 // Create the LUT
1088 //
1089
1090 const Int_t kNplan = AliTRDgeometry::kNplan;
1091 const Int_t kNlut = 128;
1092
1093 fLUTbin = kNplan * kNlut;
1094
1095 // The lookup table from Bogdan
1096 Float_t lut[kNplan][kNlut] = {
1097 {
1098 0.0070, 0.0150, 0.0224, 0.0298, 0.0374, 0.0454, 0.0533, 0.0611,
1099 0.0684, 0.0755, 0.0827, 0.0900, 0.0975, 0.1049, 0.1120, 0.1187,
1100 0.1253, 0.1318, 0.1385, 0.1453, 0.1519, 0.1584, 0.1646, 0.1704,
1101 0.1762, 0.1821, 0.1879, 0.1938, 0.1996, 0.2053, 0.2108, 0.2160,
1102 0.2210, 0.2260, 0.2310, 0.2361, 0.2411, 0.2461, 0.2509, 0.2557,
1103 0.2602, 0.2646, 0.2689, 0.2732, 0.2774, 0.2816, 0.2859, 0.2901,
1104 0.2942, 0.2983, 0.3022, 0.3061, 0.3099, 0.3136, 0.3172, 0.3207,
1105 0.3242, 0.3278, 0.3312, 0.3347, 0.3382, 0.3416, 0.3450, 0.3483,
1106 0.3515, 0.3547, 0.3579, 0.3609, 0.3639, 0.3669, 0.3698, 0.3727,
1107 0.3756, 0.3785, 0.3813, 0.3842, 0.3870, 0.3898, 0.3926, 0.3952,
1108 0.3979, 0.4005, 0.4032, 0.4057, 0.4082, 0.4108, 0.4132, 0.4157,
1109 0.4181, 0.4205, 0.4228, 0.4252, 0.4275, 0.4299, 0.4322, 0.4345,
1110 0.4367, 0.4390, 0.4412, 0.4434, 0.4456, 0.4478, 0.4499, 0.4520,
1111 0.4541, 0.4562, 0.4583, 0.4603, 0.4623, 0.4643, 0.4663, 0.4683,
1112 0.4702, 0.4722, 0.4741, 0.4758, 0.4774, 0.4790, 0.4805, 0.4824,
1113 0.4844, 0.4863, 0.4883, 0.4902, 0.4921, 0.4940, 0.4959, 0.4978
1114 },
1115 {
1116 0.0072, 0.0156, 0.0235, 0.0313, 0.0394, 0.0478, 0.0561, 0.0642,
1117 0.0718, 0.0792, 0.0868, 0.0947, 0.1025, 0.1101, 0.1172, 0.1241,
1118 0.1309, 0.1378, 0.1449, 0.1518, 0.1586, 0.1650, 0.1710, 0.1770,
1119 0.1830, 0.1891, 0.1952, 0.2011, 0.2070, 0.2125, 0.2177, 0.2229,
1120 0.2280, 0.2332, 0.2383, 0.2435, 0.2484, 0.2533, 0.2581, 0.2627,
1121 0.2670, 0.2714, 0.2757, 0.2799, 0.2842, 0.2884, 0.2927, 0.2968,
1122 0.3008, 0.3048, 0.3086, 0.3123, 0.3159, 0.3195, 0.3231, 0.3266,
1123 0.3301, 0.3335, 0.3370, 0.3404, 0.3438, 0.3471, 0.3504, 0.3536,
1124 0.3567, 0.3598, 0.3628, 0.3657, 0.3686, 0.3715, 0.3744, 0.3772,
1125 0.3800, 0.3828, 0.3856, 0.3884, 0.3911, 0.3938, 0.3965, 0.3991,
1126 0.4016, 0.4042, 0.4067, 0.4092, 0.4116, 0.4140, 0.4164, 0.4187,
1127 0.4211, 0.4234, 0.4257, 0.4280, 0.4302, 0.4325, 0.4347, 0.4369,
1128 0.4391, 0.4413, 0.4434, 0.4456, 0.4477, 0.4497, 0.4518, 0.4538,
1129 0.4558, 0.4578, 0.4598, 0.4618, 0.4637, 0.4656, 0.4675, 0.4694,
1130 0.4713, 0.4732, 0.4750, 0.4766, 0.4781, 0.4797, 0.4813, 0.4832,
1131 0.4851, 0.4870, 0.4888, 0.4906, 0.4925, 0.4942, 0.4960, 0.4978
1132 },
1133 {
1134 0.0075, 0.0163, 0.0246, 0.0328, 0.0415, 0.0504, 0.0592, 0.0674,
1135 0.0753, 0.0832, 0.0914, 0.0996, 0.1077, 0.1154, 0.1225, 0.1296,
1136 0.1369, 0.1442, 0.1515, 0.1585, 0.1652, 0.1714, 0.1776, 0.1839,
1137 0.1902, 0.1965, 0.2025, 0.2085, 0.2141, 0.2194, 0.2247, 0.2299,
1138 0.2352, 0.2405, 0.2457, 0.2507, 0.2557, 0.2604, 0.2649, 0.2693,
1139 0.2737, 0.2780, 0.2823, 0.2867, 0.2909, 0.2951, 0.2992, 0.3033,
1140 0.3072, 0.3110, 0.3146, 0.3182, 0.3218, 0.3253, 0.3288, 0.3323,
1141 0.3357, 0.3392, 0.3426, 0.3459, 0.3492, 0.3524, 0.3555, 0.3586,
1142 0.3616, 0.3645, 0.3674, 0.3703, 0.3731, 0.3759, 0.3787, 0.3815,
1143 0.3843, 0.3870, 0.3897, 0.3925, 0.3950, 0.3976, 0.4002, 0.4027,
1144 0.4052, 0.4076, 0.4101, 0.4124, 0.4148, 0.4171, 0.4194, 0.4217,
1145 0.4239, 0.4262, 0.4284, 0.4306, 0.4328, 0.4350, 0.4371, 0.4393,
1146 0.4414, 0.4435, 0.4455, 0.4476, 0.4496, 0.4516, 0.4536, 0.4555,
1147 0.4575, 0.4594, 0.4613, 0.4632, 0.4650, 0.4669, 0.4687, 0.4705,
1148 0.4723, 0.4741, 0.4758, 0.4773, 0.4789, 0.4804, 0.4821, 0.4839,
1149 0.4857, 0.4875, 0.4893, 0.4910, 0.4928, 0.4945, 0.4961, 0.4978
1150 },
1151 {
1152 0.0078, 0.0171, 0.0258, 0.0345, 0.0438, 0.0532, 0.0624, 0.0708,
1153 0.0791, 0.0875, 0.0962, 0.1048, 0.1130, 0.1206, 0.1281, 0.1356,
1154 0.1432, 0.1508, 0.1582, 0.1651, 0.1716, 0.1780, 0.1845, 0.1910,
1155 0.1975, 0.2038, 0.2099, 0.2155, 0.2210, 0.2263, 0.2317, 0.2371,
1156 0.2425, 0.2477, 0.2528, 0.2578, 0.2626, 0.2671, 0.2715, 0.2759,
1157 0.2803, 0.2846, 0.2890, 0.2933, 0.2975, 0.3016, 0.3056, 0.3095,
1158 0.3132, 0.3168, 0.3204, 0.3239, 0.3274, 0.3309, 0.3344, 0.3378,
1159 0.3412, 0.3446, 0.3479, 0.3511, 0.3543, 0.3574, 0.3603, 0.3633,
1160 0.3662, 0.3690, 0.3718, 0.3747, 0.3774, 0.3802, 0.3829, 0.3857,
1161 0.3883, 0.3910, 0.3936, 0.3962, 0.3987, 0.4012, 0.4037, 0.4061,
1162 0.4085, 0.4109, 0.4132, 0.4155, 0.4177, 0.4200, 0.4222, 0.4244,
1163 0.4266, 0.4288, 0.4309, 0.4331, 0.4352, 0.4373, 0.4394, 0.4414,
1164 0.4435, 0.4455, 0.4475, 0.4494, 0.4514, 0.4533, 0.4552, 0.4571,
1165 0.4590, 0.4608, 0.4626, 0.4645, 0.4662, 0.4680, 0.4698, 0.4715,
1166 0.4733, 0.4750, 0.4766, 0.4781, 0.4796, 0.4812, 0.4829, 0.4846,
1167 0.4863, 0.4880, 0.4897, 0.4914, 0.4930, 0.4946, 0.4963, 0.4979
1168 },
1169 {
1170 0.0081, 0.0178, 0.0270, 0.0364, 0.0463, 0.0562, 0.0656, 0.0744,
1171 0.0831, 0.0921, 0.1013, 0.1102, 0.1183, 0.1261, 0.1339, 0.1419,
1172 0.1499, 0.1576, 0.1648, 0.1715, 0.1782, 0.1849, 0.1917, 0.1984,
1173 0.2048, 0.2110, 0.2167, 0.2223, 0.2278, 0.2333, 0.2389, 0.2444,
1174 0.2497, 0.2548, 0.2598, 0.2645, 0.2691, 0.2735, 0.2780, 0.2824,
1175 0.2868, 0.2912, 0.2955, 0.2997, 0.3038, 0.3078, 0.3116, 0.3152,
1176 0.3188, 0.3224, 0.3259, 0.3294, 0.3329, 0.3364, 0.3398, 0.3432,
1177 0.3465, 0.3497, 0.3529, 0.3561, 0.3591, 0.3620, 0.3649, 0.3677,
1178 0.3705, 0.3733, 0.3761, 0.3788, 0.3816, 0.3843, 0.3869, 0.3896,
1179 0.3922, 0.3948, 0.3973, 0.3998, 0.4022, 0.4047, 0.4070, 0.4094,
1180 0.4117, 0.4139, 0.4162, 0.4184, 0.4206, 0.4227, 0.4249, 0.4270,
1181 0.4291, 0.4313, 0.4334, 0.4354, 0.4375, 0.4395, 0.4415, 0.4435,
1182 0.4455, 0.4474, 0.4493, 0.4512, 0.4531, 0.4550, 0.4568, 0.4586,
1183 0.4604, 0.4622, 0.4639, 0.4657, 0.4674, 0.4691, 0.4708, 0.4725,
1184 0.4742, 0.4758, 0.4773, 0.4788, 0.4803, 0.4819, 0.4836, 0.4852,
1185 0.4869, 0.4885, 0.4901, 0.4917, 0.4933, 0.4948, 0.4964, 0.4979
1186 },
1187 {
1188 0.0085, 0.0189, 0.0288, 0.0389, 0.0497, 0.0603, 0.0699, 0.0792,
1189 0.0887, 0.0985, 0.1082, 0.1170, 0.1253, 0.1336, 0.1421, 0.1505,
1190 0.1587, 0.1662, 0.1733, 0.1803, 0.1874, 0.1945, 0.2014, 0.2081,
1191 0.2143, 0.2201, 0.2259, 0.2316, 0.2374, 0.2431, 0.2487, 0.2541,
1192 0.2593, 0.2642, 0.2689, 0.2735, 0.2781, 0.2826, 0.2872, 0.2917,
1193 0.2961, 0.3003, 0.3045, 0.3086, 0.3125, 0.3162, 0.3198, 0.3235,
1194 0.3270, 0.3306, 0.3342, 0.3377, 0.3411, 0.3446, 0.3479, 0.3511,
1195 0.3543, 0.3575, 0.3605, 0.3634, 0.3663, 0.3691, 0.3720, 0.3748,
1196 0.3775, 0.3803, 0.3830, 0.3857, 0.3884, 0.3911, 0.3937, 0.3962,
1197 0.3987, 0.4012, 0.4036, 0.4060, 0.4084, 0.4107, 0.4129, 0.4152,
1198 0.4174, 0.4196, 0.4218, 0.4239, 0.4261, 0.4282, 0.4303, 0.4324,
1199 0.4344, 0.4365, 0.4385, 0.4405, 0.4425, 0.4445, 0.4464, 0.4483,
1200 0.4502, 0.4521, 0.4539, 0.4558, 0.4576, 0.4593, 0.4611, 0.4629,
1201 0.4646, 0.4663, 0.4680, 0.4697, 0.4714, 0.4730, 0.4747, 0.4759,
1202 0.4769, 0.4780, 0.4790, 0.4800, 0.4811, 0.4827, 0.4843, 0.4859,
1203 0.4874, 0.4889, 0.4905, 0.4920, 0.4935, 0.4950, 0.4965, 0.4979
1204 }
1205 };
1206
1207 if (fLUT) delete [] fLUT;
1208 fLUT = new Float_t[fLUTbin];
1209
1210 for (Int_t iplan = 0; iplan < kNplan; iplan++) {
1211 for (Int_t ilut = 0; ilut < kNlut; ilut++) {
1212 fLUT[iplan*kNplan+ilut] = lut[iplan][ilut];
1213 }
1214 }
1215
1216}
1217
1218//_____________________________________________________________________________
1219void AliTRDparameter::SetTiltingAngle(Float_t v)
1220{
1221 //
1222 // Set the tilting angle for the readout pads
1223 //
1224
1225 fTiltingAngle = TMath::Tan(TMath::Pi()/180.0 * v);
1226
1227}
1228
1229//_____________________________________________________________________________
1230Float_t AliTRDparameter::GetTiltingAngle() const
1231{
1232 //
1233 // Get the tilting angle for the readout pads
1234 //
1235
1236 return TMath::ATan(180.0/TMath::Pi() * fTiltingAngle);
1237
1238}
1239
1240//_____________________________________________________________________________
1241Float_t AliTRDparameter::GetDiffusionL(Float_t vd, Float_t b)
1242{
1243 //
1244 // Returns the longitudinal diffusion coefficient for a given drift
1245 // velocity <vd> and a B-field <b> for Xe/CO2 (15%).
1246 // The values are according to a GARFIELD simulation.
1247 //
1248
1249 const Int_t kNb = 5;
1250 Float_t p0[kNb] = { 0.007440, 0.007493, 0.007513, 0.007672, 0.007831 };
1251 Float_t p1[kNb] = { 0.019252, 0.018912, 0.018636, 0.018012, 0.017343 };
1252 Float_t p2[kNb] = { -0.005042, -0.004926, -0.004867, -0.004650, -0.004424 };
1253 Float_t p3[kNb] = { 0.000195, 0.000189, 0.000195, 0.000182, 0.000169 };
1254
1255 Int_t ib = ((Int_t) (10 * (b - 0.15)));
1256 ib = TMath::Max( 0,ib);
1257 ib = TMath::Min(kNb,ib);
1258
1259 Float_t diff = p0[ib]
1260 + p1[ib] * vd
1261 + p2[ib] * vd*vd
1262 + p3[ib] * vd*vd*vd;
1263
1264 return diff;
1265
1266}
1267
1268//_____________________________________________________________________________
1269Float_t AliTRDparameter::GetDiffusionT(Float_t vd, Float_t b)
1270{
1271 //
1272 // Returns the transverse diffusion coefficient for a given drift
1273 // velocity <vd> and a B-field <b> for Xe/CO2 (15%).
1274 // The values are according to a GARFIELD simulation.
1275 //
1276
1277 const Int_t kNb = 5;
1278 Float_t p0[kNb] = { 0.009550, 0.009599, 0.009674, 0.009757, 0.009850 };
1279 Float_t p1[kNb] = { 0.006667, 0.006539, 0.006359, 0.006153, 0.005925 };
1280 Float_t p2[kNb] = { -0.000853, -0.000798, -0.000721, -0.000635, -0.000541 };
1281 Float_t p3[kNb] = { 0.000131, 0.000122, 0.000111, 0.000098, 0.000085 };
1282
1283 Int_t ib = ((Int_t) (10 * (b - 0.15)));
1284 ib = TMath::Max( 0,ib);
1285 ib = TMath::Min(kNb,ib);
1286
1287 Float_t diff = p0[ib]
1288 + p1[ib] * vd
1289 + p2[ib] * vd*vd
1290 + p3[ib] * vd*vd*vd;
1291
1292 return diff;
1293
1294}
1295
1296//_____________________________________________________________________________
1297Float_t AliTRDparameter::GetOmegaTau(Float_t vd, Float_t b)
1298{
1299 //
1300 // Returns omega*tau (tan(Lorentz-angle)) for a given drift velocity <vd>
1301 // and a B-field <b> for Xe/CO2 (15%).
1302 // The values are according to a GARFIELD simulation.
1303 //
1304
1305 const Int_t kNb = 5;
1306 Float_t p0[kNb] = { 0.004810, 0.007412, 0.010252, 0.013409, 0.016888 };
1307 Float_t p1[kNb] = { 0.054875, 0.081534, 0.107333, 0.131983, 0.155455 };
1308 Float_t p2[kNb] = { -0.008682, -0.012896, -0.016987, -0.020880, -0.024623 };
1309 Float_t p3[kNb] = { 0.000155, 0.000238, 0.000330, 0.000428, 0.000541 };
1310
1311 Int_t ib = ((Int_t) (10 * (b - 0.15)));
1312 ib = TMath::Max( 0,ib);
1313 ib = TMath::Min(kNb,ib);
1314
1315 Float_t alphaL = p0[ib]
1316 + p1[ib] * vd
1317 + p2[ib] * vd*vd
1318 + p3[ib] * vd*vd*vd;
1319
1320 return TMath::Tan(alphaL);
1321
1322}
1323
1324//_____________________________________________________________________________
0a29d0f1 1325Double_t AliTRDparameter::LUTposition(Int_t iplane, Double_t ampL
1326 , Double_t ampC
1327 , Double_t ampR) const
17b26de4 1328{
1329 //
1330 // Calculates the cluster position using the lookup table.
1331 // Method provided by Bogdan Vulpescu.
1332 //
1333
1334 const Int_t kNplan = AliTRDgeometry::kNplan;
1335
1336 Double_t pos;
1337 Double_t x = 0.0;
1338 Double_t xmin;
1339 Double_t xmax;
1340 Double_t xwid;
1341
1342 Int_t side = 0;
1343 Int_t ix;
1344
1345 Double_t xMin[kNplan] = { 0.006492, 0.006377, 0.006258, 0.006144, 0.006030, 0.005980 };
1346 Double_t xMax[kNplan] = { 0.960351, 0.965870, 0.970445, 0.974352, 0.977667, 0.996101 };
1347
1348 if (ampL > ampR) {
1349 x = (ampL - ampR) / ampC;
1350 side = -1;
1351 }
1352 else if (ampL < ampR) {
1353 x = (ampR - ampL) / ampC;
1354 side = +1;
1355 }
1356
1357 if (ampL != ampR) {
1358
1359 xmin = xMin[iplane] + 0.000005;
1360 xmax = xMax[iplane] - 0.000005;
1361 xwid = (xmax - xmin) / 127.0;
1362
1363 if (x < xmin) {
1364 pos = 0.0000;
1365 }
1366 else if (x > xmax) {
1367 pos = side * 0.5000;
1368 }
1369 else {
1370 ix = (Int_t) ((x - xmin) / xwid);
1371 pos = side * fLUT[iplane*kNplan+ix];
1372 }
1373
1374 }
1375 else {
1376
1377 pos = 0.0;
1378
1379 }
1380
1381 return pos;
1382
1383}