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