]> 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
Remove compiler warnings
[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
88cb7938 16/* $Id$ */
17b26de4 17
18///////////////////////////////////////////////////////////////////////////////
19// //
20// TRD parameter class //
21// //
22///////////////////////////////////////////////////////////////////////////////
23
bdbb05bb 24#include <TRandom.h>
25
17b26de4 26#include "AliRun.h"
27#include "AliMagF.h"
28
17b26de4 29#include "AliTRDparameter.h"
30#include "AliTRDgeometryFull.h"
31
32ClassImp(AliTRDparameter)
33
34//_____________________________________________________________________________
35AliTRDparameter::AliTRDparameter():TNamed()
36{
37 //
38 // AliTRDparameter default constructor
39 //
40
41 fGeo = 0;
42 fPRFsmp = 0;
43 fTRFsmp = 0;
44 fCTsmp = 0;
45 fGasGain = 0.0;
46 fNoise = 0.0;
47 fChipGain = 0.0;
48 fADCoutRange = 0.0;
49 fADCinRange = 0.0;
50 fADCthreshold = 0;
7e582e9f 51 fADCbaseline = 0;
17b26de4 52 fDiffusionOn = 0;
53 fDiffusionT = 0.0;
54 fDiffusionL = 0.0;
55 fElAttachOn = 0;
56 fElAttachProp = 0.0;
57 fExBOn = 0;
58 fOmegaTau = 0.0;
59 fPRFOn = 0;
60 fTRFOn = 0;
61 fCTOn = 0;
62 fTCOn = 0;
63 fDriftVelocity = 0.0;
64 fPadCoupling = 0.0;
65 fTimeCoupling = 0.0;
66 fTimeBinWidth = 0.0;
67 fField = 0.0;
68 fTiltingAngle = 0.0;
69 fPRFbin = 0;
70 fPRFlo = 0.0;
71 fPRFhi = 0.0;
72 fPRFwid = 0.0;
73 fPRFpad = 0;
74 fTRFbin = 0;
75 fTRFlo = 0.0;
76 fTRFhi = 0.0;
77 fTRFwid = 0.0;
78 fTCnexp = 0;
79
80 fLUTOn = 0;
81 fLUT = 0;
82 fClusMaxThresh = 0;
83 fClusSigThresh = 0;
84
a328fff9 85 fTimeStructOn = 0;
86 fTimeStruct = 0;
a328fff9 87 fAnodeWireOffset = 0.0;
88
17b26de4 89}
90
91//_____________________________________________________________________________
92AliTRDparameter::AliTRDparameter(const Text_t *name, const Text_t *title)
93 :TNamed(name,title)
94{
95 //
96 // AliTRDparameter constructor
97 //
98
99 fGeo = new AliTRDgeometryFull();
67689223 100 fPRFsmp = 0;
101 fTRFsmp = 0;
102 fCTsmp = 0;
103 fGasGain = 0.0;
104 fNoise = 0.0;
105 fChipGain = 0.0;
106 fADCoutRange = 0.0;
107 fADCinRange = 0.0;
108 fADCthreshold = 0;
7e582e9f 109 fADCbaseline = 0;
67689223 110 fDiffusionOn = 0;
111 fDiffusionT = 0.0;
112 fDiffusionL = 0.0;
113 fElAttachOn = 0;
114 fElAttachProp = 0.0;
115 fExBOn = 0;
116 fOmegaTau = 0.0;
117 fPRFOn = 0;
118 fTRFOn = 0;
119 fCTOn = 0;
120 fTCOn = 0;
121 fDriftVelocity = 0.0;
122 fPadCoupling = 0.0;
123 fTimeCoupling = 0.0;
124 fTimeBinWidth = 0.0;
125 fField = 0.0;
126 fTiltingAngle = 0.0;
127 fPRFbin = 0;
128 fPRFlo = 0.0;
129 fPRFhi = 0.0;
130 fPRFwid = 0.0;
131 fPRFpad = 0;
132 fTRFbin = 0;
133 fTRFlo = 0.0;
134 fTRFhi = 0.0;
135 fTRFwid = 0.0;
136 fTCnexp = 0;
137
138 fLUTOn = 0;
139 fLUT = 0;
140 fClusMaxThresh = 0;
141 fClusSigThresh = 0;
17b26de4 142
a328fff9 143 fTimeStructOn = 0;
144 fTimeStruct = 0;
a328fff9 145 fAnodeWireOffset = 0.0;
146
17b26de4 147 Init();
148
149}
150
151
152//_____________________________________________________________________________
73ae7b59 153AliTRDparameter::AliTRDparameter(const AliTRDparameter &p):TNamed(p)
17b26de4 154{
155 //
156 // AliTRDparameter copy constructor
157 //
158
159 ((AliTRDparameter &) p).Copy(*this);
160
161}
162
163///_____________________________________________________________________________
164AliTRDparameter::~AliTRDparameter()
165{
166 //
167 // AliTRDparameter destructor
168 //
169
170 if (fTRFsmp) {
171 delete [] fTRFsmp;
172 fTRFsmp = 0;
173 }
174
175 if (fPRFsmp) {
176 delete [] fPRFsmp;
177 fPRFsmp = 0;
178 }
179
180 if (fCTsmp) {
181 delete [] fCTsmp;
182 fCTsmp = 0;
183 }
184
185 if (fLUT) {
186 delete [] fLUT;
187 fLUT = 0;
188 }
189
190 if (fGeo) {
191 delete fGeo;
192 fGeo = 0;
193 }
194
a328fff9 195 if (fTimeStruct) {
196 delete [] fTimeStruct;
197 fTimeStruct = 0;
198 }
199
17b26de4 200}
201
202//_____________________________________________________________________________
203AliTRDparameter &AliTRDparameter::operator=(const AliTRDparameter &p)
204{
205 //
206 // Assignment operator
207 //
208
209 if (this != &p) ((AliTRDparameter &) p).Copy(*this);
210 return *this;
211
212}
213
214//_____________________________________________________________________________
e0d47c25 215void AliTRDparameter::Copy(TObject &p) const
17b26de4 216{
217 //
218 // Copy function
219 //
220
221 Int_t iBin;
222
223 ((AliTRDparameter &) p).fGasGain = fGasGain;
224 ((AliTRDparameter &) p).fNoise = fNoise;
225 ((AliTRDparameter &) p).fChipGain = fChipGain;
226 ((AliTRDparameter &) p).fADCoutRange = fADCoutRange;
227 ((AliTRDparameter &) p).fADCinRange = fADCinRange;
228 ((AliTRDparameter &) p).fADCthreshold = fADCthreshold;
7e582e9f 229 ((AliTRDparameter &) p).fADCbaseline = fADCbaseline;
17b26de4 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
a328fff9 284 ((AliTRDparameter &) p).fAnodeWireOffset = fAnodeWireOffset;
a328fff9 285 ((AliTRDparameter &) p).fTimeStructOn = fTimeStructOn;
286 if (((AliTRDparameter &) p).fTimeStruct)
287 delete [] ((AliTRDparameter &) p).fTimeStruct;
682619b9 288 ((AliTRDparameter &) p).fTimeStruct = new Float_t[38*11];
289 for (Int_t i = 0; i < 38*11; i++) {
a328fff9 290 ((AliTRDparameter &) p).fTimeStruct[i] = fTimeStruct[i];
291 }
292
17b26de4 293}
294
295//_____________________________________________________________________________
296void AliTRDparameter::Init()
297{
298 //
299 // Initializes the parameter
300 //
301 // The maximum number of pads
302 // and the position of pad 0,0,0
303 //
304 // chambers seen from the top:
305 // +----------------------------+
306 // | |
307 // | | ^
308 // | | rphi|
309 // | | |
310 // |0 | |
311 // +----------------------------+ +------>
312 // z
313 // chambers seen from the side: ^
314 // +----------------------------+ drift|
315 // |0 | |
316 // | | |
317 // +----------------------------+ +------>
318 // z
319 //
320 // IMPORTANT: time bin 0 is now the first one in the drift region
321 // closest to the readout !!!
322 //
323
324 //
325 // ----------------------------------------------------------------------------
326 // The pad definition
327 // ----------------------------------------------------------------------------
328 //
329
330 // The pad size in column direction (rphi-direction)
a328fff9 331 SetColPadSize(0,0.664);
332 SetColPadSize(1,0.695);
333 SetColPadSize(2,0.726);
334 SetColPadSize(3,0.756);
335 SetColPadSize(4,0.788);
336 SetColPadSize(5,0.818);
73ae7b59 337
17b26de4 338 // The pad row (z-direction)
339 SetNRowPad();
340
341 // The number of time bins. Default is 100 ns timbin size
81a6930e 342 SetNTimeBin(18);
17b26de4 343
344 // Additional time bins before and after the drift region.
345 // Default is to only sample the drift region
346 SetExpandTimeBin(0,0);
347
348 //
349 // ----------------------------------------------------------------------------
350 // The digitization parameter
351 // ----------------------------------------------------------------------------
352 //
353
354 // The default parameter for the digitization
a328fff9 355 fGasGain = 4000.;
17b26de4 356 fChipGain = 12.4;
357 fNoise = 1000.;
358 fADCoutRange = 1023.; // 10-bit ADC
7e582e9f 359 fADCinRange = 2000.; // 2V input range
17b26de4 360 fADCthreshold = 1;
7e582e9f 361 fADCbaseline = 0;
17b26de4 362
363 // The drift velocity (cm / mus)
81a6930e 364 fDriftVelocity = 1.62;
17b26de4 365
366 // Diffusion on
367 fDiffusionOn = 1;
368
369 // E x B effects
7e582e9f 370 fExBOn = 1;
17b26de4 371
372 // Propability for electron attachment
373 fElAttachOn = 0;
374 fElAttachProp = 0.0;
375
376 // The pad response function
377 fPRFOn = 1;
378
379 // The time response function
380 fTRFOn = 1;
381
382 // The cross talk
7e582e9f 383 fCTOn = 1;
17b26de4 384
385 // The tail cancelation
386 fTCOn = 1;
387
388 // The number of exponentials
81a6930e 389 fTCnexp = 1;
17b26de4 390
b258a7c7 391 // The pad coupling factor
392 //fPadCoupling = 0.3;
682619b9 393 // Use 0.46 instead which reproduces better the test beam
b258a7c7 394 // data, even tough it is not understood why.
682619b9 395 fPadCoupling = 0.46;
17b26de4 396
397 // The time coupling factor (same number as for the TPC)
398 fTimeCoupling = 0.4;
399
a328fff9 400 // Drift time non-isochronity on
401 fTimeStructOn = 1;
402
a328fff9 403 // Distance of first Anode wire from first pad edge
404 fAnodeWireOffset = 0.25;
405
17b26de4 406 // The tilting angle for the readout pads
73ae7b59 407 SetTiltingAngle(2.0);
17b26de4 408
409 // The magnetic field strength in Tesla
81a6930e 410 Double_t x[3] = { 0.0, 0.0, 0.0 };
411 Double_t b[3];
412 gAlice->Field(x,b); // b[] is in kilo Gauss
413 fField = b[2] * 0.1; // Tesla
17b26de4 414
415 //
416 // ----------------------------------------------------------------------------
417 // The clusterization parameter
418 // ----------------------------------------------------------------------------
419 //
420
421 // The default parameter for the clustering
422 fClusMaxThresh = 3;
423 fClusSigThresh = 1;
424
425 // Use the LUT
426 fLUTOn = 1;
427
428 ReInit();
429
430}
431
432//_____________________________________________________________________________
433void AliTRDparameter::ReInit()
434{
435 //
436 // Reinitializes the parameter class after a change
437 //
438
a328fff9 439 // Create the sampled timing structure
440 SampleTimeStruct();
441
17b26de4 442 // Calculate the time bin width in ns
443 fTimeBinWidth = fTimeBinSize / fDriftVelocity * 1000.0;
444
445 // The range and the binwidth for the sampled TRF
446 fTRFbin = 100;
447 // Start 0.2 mus before the signal
448 fTRFlo = -0.2 * fDriftVelocity;
449 // End the maximum driftlength after the signal
450 fTRFhi = AliTRDgeometry::DrThick()
451 + fTimeAfter * fTimeBinSize;
452 fTRFwid = (fTRFhi - fTRFlo) / ((Float_t) fTRFbin);
453
454 // Transverse and longitudinal diffusion coefficients (Xe/CO2)
455 fDiffusionT = GetDiffusionT(fDriftVelocity,fField);
456 fDiffusionL = GetDiffusionL(fDriftVelocity,fField);
457
458 // omega * tau.= tan(Lorentz-angle)
459 fOmegaTau = GetOmegaTau(fDriftVelocity,fField);
460
461 // The Lorentz factor
462 if (fExBOn) {
463 fLorentzFactor = 1.0 / (1.0 + fOmegaTau*fOmegaTau);
464 }
465 else {
466 fLorentzFactor = 1.0;
467 }
468
469 // Create the sampled PRF
470 SamplePRF();
471
472 // Create the sampled TRF
473 SampleTRF();
474
475 // Create the LUT
476 FillLUT();
477
478}
479
480//_____________________________________________________________________________
afc51ac2 481void AliTRDparameter::SetNRowPad(Int_t p, Int_t c, Int_t npad)
17b26de4 482{
483 //
484 // Redefines the number of pads in raw direction for
485 // a given plane and chamber number
486 //
487
488 for (Int_t isect = 0; isect < AliTRDgeometry::Nsect(); isect++) {
489
490 fRowMax[p][c][isect] = npad;
491
492 fRowPadSize[p][c][isect] = (fGeo->GetChamberLength(p,c)
493 - 2.* AliTRDgeometry::RpadW())
494 / ((Float_t) npad);
495
496 }
497
498}
499
500//_____________________________________________________________________________
501void AliTRDparameter::SetNRowPad()
502{
503 //
504 // Defines the number of pads in row direction
505 //
506
507 Int_t isect;
508 Int_t icham;
509 Int_t iplan;
510
511 Int_t rowMax[kNplan][kNcham] = { { 16, 16, 12, 16, 16 }
512 , { 16, 16, 12, 16, 16 }
513 , { 16, 16, 12, 16, 16 }
514 , { 16, 16, 12, 16, 16 }
ad973dd4 515 , { 16, 16, 12, 16, 16 }
516 , { 16, 16, 12, 16, 16 } };
17b26de4 517
518 Float_t rpadW = AliTRDgeometry::RpadW();
519
520 for (isect = 0; isect < kNsect; isect++) {
521 for (icham = 0; icham < kNcham; icham++) {
522 for (iplan = 0; iplan < kNplan; iplan++) {
523
524 fRowMax[iplan][icham][isect] = rowMax[iplan][icham];
525
526 fRowPadSize[iplan][icham][isect] = (fGeo->GetChamberLength(iplan,icham)
527 - 2.*rpadW)
528 / ((Float_t) rowMax[iplan][icham]);
529
530 Float_t row0 = rpadW - fGeo->GetChamberLength(iplan,0)
531 - fGeo->GetChamberLength(iplan,1)
532 - fGeo->GetChamberLength(iplan,2) / 2.;
533 for (Int_t ic = 0; ic < icham; ic++) {
4bd34b78 534 row0 += fGeo->GetChamberLength(iplan,ic);
17b26de4 535 }
81a6930e 536
537 fRow0[iplan][icham][isect] = row0;
538 // For new chamber ordering
539 //fRow0[iplan][kNcham-icham-1][isect] = row0;
17b26de4 540
541 }
542 }
543 }
544
545}
546
547//_____________________________________________________________________________
afc51ac2 548void AliTRDparameter::SetColPadSize(Int_t p, Float_t s)
17b26de4 549{
550 //
551 // Redefines the pad size in column direction
552 //
553
554 Float_t cpadW = AliTRDgeometry::CpadW();
555
556 fColPadSize[p] = s;
557 fCol0[p] = - fGeo->GetChamberWidth(p)/2. + cpadW;
558 fColMax[p] = ((Int_t) ((fGeo->GetChamberWidth(p) - 2.*cpadW) / s));
559
560}
561
562//_____________________________________________________________________________
afc51ac2 563void AliTRDparameter::SetNTimeBin(Int_t nbin)
17b26de4 564{
565 //
566 // Redefines the number of time bins in the drift region.
567 // The time bin width is defined by the length of the
568 // drift region divided by <nbin>.
569 //
570
571 fTimeMax = nbin;
572 fTimeBinSize = AliTRDgeometry::DrThick() / ((Float_t) fTimeMax);
573 for (Int_t iplan = 0; iplan < AliTRDgeometry::Nplan(); iplan++) {
574 fTime0[iplan] = AliTRDgeometry::Rmin()
575 + AliTRDgeometry::CraHght()
576 + AliTRDgeometry::CdrHght()
577 + iplan * (AliTRDgeometry::Cheight()
578 + AliTRDgeometry::Cspace());
579 }
580
581}
582
583//_____________________________________________________________________________
0a29d0f1 584Float_t AliTRDparameter::CrossTalk(Float_t time) const
17b26de4 585{
586 //
587 // Applies the pad-pad capacitive cross talk
588 //
589
590 Int_t iBin = ((Int_t) ((time - fTRFlo) / fTRFwid));
591 if ((iBin >= 0) && (iBin < fTRFbin)) {
592 return fCTsmp[iBin];
593 }
594 else {
595 return 0.0;
596 }
597
598}
599
a328fff9 600//_____________________________________________________________________________
682619b9 601Int_t AliTRDparameter::TimeStruct(Float_t dist, Float_t z, Float_t *xyz) const
a328fff9 602{
603 //
604 // Applies the time structure of the drift cells (by C.Lippmann).
605 // The drift time of electrons to the anode wires depends on the
606 // distance to the wire (z) and on the position in the drift region.
682619b9 607 // Here we omnly reposition each electron, since the drift velocity
608 // in the digitisation is still constant. If the electrons at a given
609 // position has a long drift time, it is just moved further away from
610 // the anode wire plane.
611 //
612 // input :
613 // dist = radial distance from anode wire plane [cm]
614 // z = distance from anode wire (parallel to cathode planes) [cm]
a328fff9 615 //
616
a328fff9 617
682619b9 618 // indices:
619 Int_t r1 = (Int_t)(10*dist);
a328fff9 620 Int_t r2 = r1+1;
682619b9 621 Int_t z1 = (Int_t)(100*z/2.5);
a328fff9 622 Int_t z2 = z1+1;
623
682619b9 624 if (r1<0 || r1>37 || z1<0 || z1>10) {
a328fff9 625 printf("<AliTRDparameter::TimeStruct> Warning. Indices out of range: ");
682619b9 626 printf("dist=%.2f, z=%.2f, r1=%d, z1=%d\n",dist,z,r1,z1);
a328fff9 627 return kFALSE;
a328fff9 628 }
629
630 // 2D Interpolation:
682619b9 631 Float_t y11 = fTimeStruct[r1+38*z1];
632 Float_t y22 = (r2 <= 37 && z2 <= 10) ? fTimeStruct[r2+38*z2] : fTimeStruct[37+38*10];
633 Float_t y12 = (z2 <= 10) ? fTimeStruct[r1+38*z2] : fTimeStruct[r1+38*10];
634 Float_t y21 = (r2 <= 37) ? fTimeStruct[r2+38*z1] : fTimeStruct[37+38*z1];
635
636 Float_t y1 = (y21-y11)*10*dist + y11 - (y21-y11)*r1;
637 Float_t y2 = (y22-y12)*10*dist + y12 - (y22-y12)*r1;
638
639 Float_t AmTh = AliTRDgeometry::AmThick()/2.0;
640
641 // dist now is the drift distance to anode wires (negative if electrons are
642 // between anode wire plane and cathode pad plane)
643 dist -= AmTh;
644
645 // Get the drift time from the interpolation:
646 Float_t tdrift =
647 ( TMath::Abs(dist)>0.005 || z>0.005 ) ? tdrift = (y2-y1)*100*z/2.5+y1-(y2-y1)*z1 : 0.0;
648
649 // We move electrons further away from the anode wire plane, if the drift time is
650 // larger than the one expected for that distance for a constant drift velocity. We
651 // increase (decrease) the drift distance by a factor
652 //
653 // c = t(x0)*vd/x0 ,
654 //
655 // where x0 is the distance from the anode wire plane (dist), vd is the constant
656 // drift velocity (fDriftVelocity) and t(x0) is the real drift time (tdrift, as
657 // calculated with GARFIELD).
658 //
659 // The factor is negative for electrons between anode wire plane and cathode pads.
660 //
661 // The new position of the electron is then given by: t(x0)*vd:
662
663 Float_t zdrift = tdrift * fDriftVelocity;
664 if (dist < 0.) zdrift *= -1.;
665
666 xyz[0] = xyz[0] + dist - zdrift;
a328fff9 667 xyz[1] = xyz[1];
668 xyz[2] = xyz[2];
669
670 return 1;
671
672}
673
17b26de4 674//_____________________________________________________________________________
675Int_t AliTRDparameter::Diffusion(Float_t driftlength, Float_t *xyz)
676{
677 //
678 // Applies the diffusion smearing to the position of a single electron
679 //
680
681 Float_t driftSqrt = TMath::Sqrt(driftlength);
682 Float_t sigmaT = driftSqrt * fDiffusionT;
683 Float_t sigmaL = driftSqrt * fDiffusionL;
684 xyz[0] = gRandom->Gaus(xyz[0], sigmaL * fLorentzFactor);
685 xyz[1] = gRandom->Gaus(xyz[1], sigmaT * fLorentzFactor);
686 xyz[2] = gRandom->Gaus(xyz[2], sigmaT);
687
688 return 1;
689
690}
691
692//_____________________________________________________________________________
0a29d0f1 693Int_t AliTRDparameter::ExB(Float_t driftlength, Float_t *xyz) const
17b26de4 694{
695 //
696 // Applies E x B effects to the position of a single electron
697 //
698
699 xyz[0] = xyz[0];
700 xyz[1] = xyz[1] + fOmegaTau * driftlength;
701 xyz[2] = xyz[2];
702
703 return 1;
704
705}
706
707//_____________________________________________________________________________
708Int_t AliTRDparameter::PadResponse(Float_t signal, Float_t dist
0a29d0f1 709 , Int_t plane, Float_t *pad) const
17b26de4 710{
711 //
712 // Applies the pad response
713 //
714
715 const Int_t kNplan = AliTRDgeometry::kNplan;
716
717 Int_t iBin = ((Int_t) (( - dist - fPRFlo) / fPRFwid));
718 Int_t iOff = plane * fPRFbin;
719
720 Int_t iBin0 = iBin - fPRFpad + iOff;
721 Int_t iBin1 = iBin + iOff;
722 Int_t iBin2 = iBin + fPRFpad + iOff;
723
724 pad[0] = 0.0;
725 pad[1] = 0.0;
726 pad[2] = 0.0;
727 if ((iBin1 >= 0) && (iBin1 < (fPRFbin*kNplan))) {
728
729 if (iBin0 >= 0) {
730 pad[0] = signal * fPRFsmp[iBin0];
731 }
732 pad[1] = signal * fPRFsmp[iBin1];
733 if (iBin2 < (fPRFbin*kNplan)) {
734 pad[2] = signal * fPRFsmp[iBin2];
735 }
736
737 return 1;
738
739 }
740 else {
741
742 return 0;
743
744 }
745
746}
747
748//_____________________________________________________________________________
0a29d0f1 749Float_t AliTRDparameter::TimeResponse(Float_t time) const
17b26de4 750{
751 //
752 // Applies the preamp shaper time response
753 //
754
81a6930e 755 Int_t iBin = ((Int_t) ((time) / fTRFwid));
756 //Int_t iBin = ((Int_t) ((time - fTRFlo) / fTRFwid));
17b26de4 757 if ((iBin >= 0) && (iBin < fTRFbin)) {
758 return fTRFsmp[iBin];
759 }
760 else {
761 return 0.0;
762 }
763
764}
765
766//_____________________________________________________________________________
767Float_t AliTRDparameter::Col0Tilted(Float_t col0, Float_t rowOffset
768 , Int_t plane)
769{
770 //
771 // Calculates col0 for tilted pads
772 //
773
774 Float_t diff = fTiltingAngle * rowOffset;
a328fff9 775 return (col0 + TMath::Power(-1.0,(plane+1)) * diff);
17b26de4 776
777}
778
779//_____________________________________________________________________________
780void AliTRDparameter::SampleTRF()
781{
782 //
783 // Samples the time response function
784 //
785 // New TRF from Venelin Angelov, simulated with CADENCE
786 // Pad-ground capacitance = 25 pF
787 // Pad-pad cross talk capacitance = 6 pF
788 //
789
790 Int_t ipos1;
791 Int_t ipos2;
792 Float_t diff;
793
794 const Int_t kNpasa = 252;
795
796 Float_t time[kNpasa] = { -0.220000, -0.210000, -0.200000, -0.190000
797 , -0.180000, -0.170000, -0.160000, -0.150000
798 , -0.140000, -0.130000, -0.120000, -0.110000
799 , -0.100000, -0.090000, -0.080000, -0.070000
800 , -0.060000, -0.050000, -0.040000, -0.030000
801 , -0.020000, -0.010000, -0.000000, 0.010000
802 , 0.020000, 0.030000, 0.040000, 0.050000
803 , 0.060000, 0.070000, 0.080000, 0.090000
804 , 0.100000, 0.110000, 0.120000, 0.130000
805 , 0.140000, 0.150000, 0.160000, 0.170000
806 , 0.180000, 0.190000, 0.200000, 0.210000
807 , 0.220000, 0.230000, 0.240000, 0.250000
808 , 0.260000, 0.270000, 0.280000, 0.290000
809 , 0.300000, 0.310000, 0.320000, 0.330000
810 , 0.340000, 0.350000, 0.360000, 0.370000
811 , 0.380000, 0.390000, 0.400000, 0.410000
812 , 0.420000, 0.430000, 0.440000, 0.450000
813 , 0.460000, 0.470000, 0.480000, 0.490000
814 , 0.500000, 0.510000, 0.520000, 0.530000
815 , 0.540000, 0.550000, 0.560000, 0.570000
816 , 0.580000, 0.590000, 0.600000, 0.610000
817 , 0.620000, 0.630000, 0.640000, 0.650000
818 , 0.660000, 0.670000, 0.680000, 0.690000
819 , 0.700000, 0.710000, 0.720000, 0.730000
820 , 0.740000, 0.750000, 0.760000, 0.770000
821 , 0.780000, 0.790000, 0.800000, 0.810000
822 , 0.820000, 0.830000, 0.840000, 0.850000
823 , 0.860000, 0.870000, 0.880000, 0.890000
824 , 0.900000, 0.910000, 0.920000, 0.930000
825 , 0.940000, 0.950000, 0.960000, 0.970000
826 , 0.980000, 0.990000, 1.000000, 1.010000
827 , 1.020000, 1.030000, 1.040000, 1.050000
828 , 1.060000, 1.070000, 1.080000, 1.090000
829 , 1.100000, 1.110000, 1.120000, 1.130000
830 , 1.140000, 1.150000, 1.160000, 1.170000
831 , 1.180000, 1.190000, 1.200000, 1.210000
832 , 1.220000, 1.230000, 1.240000, 1.250000
833 , 1.260000, 1.270000, 1.280000, 1.290000
834 , 1.300000, 1.310000, 1.320000, 1.330000
835 , 1.340000, 1.350000, 1.360000, 1.370000
836 , 1.380000, 1.390000, 1.400000, 1.410000
837 , 1.420000, 1.430000, 1.440000, 1.450000
838 , 1.460000, 1.470000, 1.480000, 1.490000
839 , 1.500000, 1.510000, 1.520000, 1.530000
840 , 1.540000, 1.550000, 1.560000, 1.570000
841 , 1.580000, 1.590000, 1.600000, 1.610000
842 , 1.620000, 1.630000, 1.640000, 1.650000
843 , 1.660000, 1.670000, 1.680000, 1.690000
844 , 1.700000, 1.710000, 1.720000, 1.730000
845 , 1.740000, 1.750000, 1.760000, 1.770000
846 , 1.780000, 1.790000, 1.800000, 1.810000
847 , 1.820000, 1.830000, 1.840000, 1.850000
848 , 1.860000, 1.870000, 1.880000, 1.890000
849 , 1.900000, 1.910000, 1.920000, 1.930000
850 , 1.940000, 1.950000, 1.960000, 1.970000
851 , 1.980000, 1.990000, 2.000000, 2.010000
852 , 2.020000, 2.030000, 2.040000, 2.050000
853 , 2.060000, 2.070000, 2.080000, 2.090000
854 , 2.100000, 2.110000, 2.120000, 2.130000
855 , 2.140000, 2.150000, 2.160000, 2.170000
856 , 2.180000, 2.190000, 2.200000, 2.210000
857 , 2.220000, 2.230000, 2.240000, 2.250000
858 , 2.260000, 2.270000, 2.280000, 2.290000 };
859
860 Float_t signal[kNpasa] = { 0.000000, 0.000000, 0.000000, 0.000000
861 , 0.000000, 0.000000, 0.000000, 0.000396
862 , 0.005096, 0.022877, 0.061891, 0.126614
863 , 0.215798, 0.324406, 0.444507, 0.566817
864 , 0.683465, 0.787089, 0.873159, 0.937146
865 , 0.979049, 0.999434, 1.000000, 0.983579
866 , 0.954134, 0.913364, 0.866365, 0.813703
867 , 0.759910, 0.706116, 0.653454, 0.603624
868 , 0.556625, 0.514156, 0.475085, 0.439977
869 , 0.408834, 0.380578, 0.355549, 0.333352
870 , 0.313647, 0.296093, 0.280351, 0.266195
871 , 0.253397, 0.241789, 0.231257, 0.221574
872 , 0.212627, 0.204417, 0.196772, 0.189581
873 , 0.182956, 0.176784, 0.171008, 0.165515
874 , 0.160419, 0.155606, 0.151076, 0.146716
875 , 0.142639, 0.138845, 0.135221, 0.131767
876 , 0.128482, 0.125368, 0.122424, 0.119592
877 , 0.116931, 0.114326, 0.111891, 0.109513
878 , 0.107248, 0.105096, 0.103058, 0.101019
879 , 0.099151, 0.097282, 0.095527, 0.093715
880 , 0.092129, 0.090544, 0.088958, 0.087429
881 , 0.086014, 0.084598, 0.083239, 0.081880
882 , 0.080634, 0.079388, 0.078143, 0.077010
883 , 0.075878, 0.074745, 0.073669, 0.072593
884 , 0.071574, 0.070612, 0.069649, 0.068686
885 , 0.067780, 0.066874, 0.066025, 0.065176
886 , 0.064326, 0.063533, 0.062684, 0.061948
887 , 0.061212, 0.060419, 0.059740, 0.059003
888 , 0.058324, 0.057644, 0.057022, 0.056342
889 , 0.055663, 0.055096, 0.054473, 0.053851
890 , 0.053284, 0.052718, 0.052152, 0.051585
891 , 0.051019, 0.050566, 0.050000, 0.049490
892 , 0.048981, 0.048528, 0.048018, 0.047508
893 , 0.047055, 0.046602, 0.046149, 0.045696
894 , 0.045300, 0.044904, 0.044451, 0.044054
895 , 0.043658, 0.043205, 0.042865, 0.042469
896 , 0.042072, 0.041733, 0.041336, 0.040997
897 , 0.040657, 0.040260, 0.039921, 0.039581
898 , 0.039241, 0.038958, 0.038618, 0.038335
899 , 0.037995, 0.037656, 0.037373, 0.037089
900 , 0.036806, 0.036467, 0.036183, 0.035900
901 , 0.035617, 0.035334, 0.035108, 0.034824
902 , 0.034541, 0.034315, 0.034032, 0.033805
903 , 0.033522, 0.033296, 0.033069, 0.032786
904 , 0.032559, 0.032333, 0.032106, 0.031880
905 , 0.031653, 0.031427, 0.031200, 0.030974
906 , 0.030804, 0.030578, 0.030351, 0.030125
907 , 0.029955, 0.029785, 0.029558, 0.029332
908 , 0.029162, 0.028992, 0.028766, 0.028596
909 , 0.028426, 0.028199, 0.028086, 0.027860
910 , 0.027746, 0.027633, 0.027463, 0.027293
911 , 0.027180, 0.027067, 0.026954, 0.026954
912 , 0.026840, 0.026727, 0.026727, 0.026614
913 , 0.026614, 0.026614, 0.026557, 0.026501
914 , 0.026501, 0.026501, 0.026501, 0.026501
915 , 0.026501, 0.026501, 0.026501, 0.026387
916 , 0.026387, 0.026387, 0.026387, 0.026387
917 , 0.026387, 0.026387, 0.026387, 0.026387
918 , 0.026387, 0.026387, 0.026387, 0.026387
919 , 0.026387, 0.026274, 0.026274, 0.026274
920 , 0.026274, 0.026274, 0.026274, 0.026274
921 , 0.026274, 0.026274, 0.026274, 0.026274
922 , 0.026274, 0.026274, 0.026274, 0.026161 };
923
924 Float_t xtalk[kNpasa] = { 0.000000, 0.000000, 0.000000, 0.000000
925 , 0.000000, 0.000000, 0.000000, 0.000113
926 , 0.000793, 0.003058, 0.007305, 0.013194
927 , 0.019706, 0.025821, 0.030634, 0.033465
928 , 0.034145, 0.032729, 0.029615, 0.025198
929 , 0.019989, 0.014496, 0.009003, 0.003964
930 , -0.000510, -0.004190, -0.007191, -0.009400
931 , -0.010872, -0.011835, -0.012288, -0.012288
932 , -0.012005, -0.011495, -0.010872, -0.010136
933 , -0.009343, -0.008607, -0.007871, -0.007191
934 , -0.006512, -0.005946, -0.005379, -0.004926
935 , -0.004473, -0.004077, -0.003737, -0.003398
936 , -0.003114, -0.002831, -0.002605, -0.002378
937 , -0.002208, -0.002039, -0.001869, -0.001699
938 , -0.001585, -0.001472, -0.001359, -0.001246
939 , -0.001132, -0.001019, -0.001019, -0.000906
940 , -0.000906, -0.000793, -0.000793, -0.000680
941 , -0.000680, -0.000680, -0.000566, -0.000566
942 , -0.000566, -0.000566, -0.000453, -0.000453
943 , -0.000453, -0.000453, -0.000453, -0.000453
944 , -0.000340, -0.000340, -0.000340, -0.000340
945 , -0.000340, -0.000340, -0.000340, -0.000340
946 , -0.000340, -0.000340, -0.000340, -0.000340
947 , -0.000340, -0.000227, -0.000227, -0.000227
948 , -0.000227, -0.000227, -0.000227, -0.000227
949 , -0.000227, -0.000227, -0.000227, -0.000227
950 , -0.000227, -0.000227, -0.000227, -0.000227
951 , -0.000227, -0.000227, -0.000227, -0.000227
952 , -0.000227, -0.000227, -0.000227, -0.000227
953 , -0.000227, -0.000227, -0.000227, -0.000227
954 , -0.000227, -0.000227, -0.000227, -0.000227
955 , -0.000227, -0.000227, -0.000227, -0.000227
956 , -0.000227, -0.000227, -0.000227, -0.000113
957 , -0.000113, -0.000113, -0.000113, -0.000113
958 , -0.000113, -0.000113, -0.000113, -0.000113
959 , -0.000113, -0.000113, -0.000113, -0.000113
960 , -0.000113, -0.000113, -0.000113, -0.000113
961 , -0.000113, -0.000113, -0.000113, -0.000113
962 , -0.000113, -0.000113, -0.000113, -0.000113
963 , -0.000113, -0.000113, -0.000113, -0.000113
964 , -0.000113, -0.000113, -0.000113, -0.000113
965 , -0.000113, -0.000113, -0.000113, -0.000113
966 , -0.000113, -0.000113, -0.000113, -0.000113
967 , -0.000113, -0.000113, -0.000113, -0.000113
968 , -0.000113, -0.000113, -0.000113, -0.000113
969 , -0.000113, -0.000113, -0.000113, -0.000113
970 , -0.000113, -0.000113, -0.000113, -0.000113
971 , -0.000113, -0.000113, -0.000113, -0.000113
972 , -0.000113, -0.000113, -0.000113, -0.000113
973 , -0.000113, -0.000113, -0.000113, -0.000113
974 , -0.000113, -0.000113, -0.000113, -0.000113
975 , -0.000113, -0.000113, -0.000113, -0.000113
976 , -0.000113, -0.000113, -0.000113, 0.000000
977 , 0.000000, 0.000000, 0.000000, 0.000000
978 , 0.000000, 0.000000, 0.000000, 0.000000
979 , 0.000000, 0.000000, 0.000000, 0.000000
980 , 0.000000, 0.000000, 0.000000, 0.000000
981 , 0.000000, 0.000000, 0.000000, 0.000000
982 , 0.000000, 0.000000, 0.000000, 0.000000
983 , 0.000000, 0.000000, 0.000000, 0.000000
984 , 0.000000, 0.000000, 0.000000, 0.000000
985 , 0.000000, 0.000000, 0.000000, 0.000000
986 , 0.000000, 0.000000, 0.000000, 0.000000 };
987
988 // increase CrossTalk to measurements
989 for (Int_t ipasa = 0; ipasa < kNpasa; ipasa++) {
990 xtalk[ipasa] *= 1.75;
991 }
992
993 if (fTRFsmp) delete [] fTRFsmp;
994 fTRFsmp = new Float_t[fTRFbin];
995 if (fCTsmp) delete [] fCTsmp;
996 fCTsmp = new Float_t[fTRFbin];
997
998 Float_t loTRF = TMath::Max(fTRFlo / fDriftVelocity,time[0]);
999 Float_t hiTRF = TMath::Min(fTRFhi / fDriftVelocity,time[kNpasa-1]);
1000 Float_t binWidth = (hiTRF - loTRF) / ((Float_t) fTRFbin);
1001
1002 // Take the linear interpolation
1003 for (Int_t iBin = 0; iBin < fTRFbin; iBin++) {
1004
1005 Float_t bin = (((Float_t) iBin) + 0.5) * binWidth + loTRF;
1006 ipos1 = ipos2 = 0;
1007 diff = 0;
1008 do {
1009 diff = bin - time[ipos2++];
1010 } while (diff > 0);
1011 ipos2--;
1012 if (ipos2 >= kNpasa) ipos2 = kNpasa - 1;
1013 ipos1 = ipos2 - 1;
1014
1015 fTRFsmp[iBin] = signal[ipos2]
1016 + diff * (signal[ipos2] - signal[ipos1])
1017 / ( time[ipos2] - time[ipos1]);
1018
1019 fCTsmp[iBin] = xtalk[ipos2]
1020 + diff * (xtalk[ipos2] - xtalk[ipos1])
1021 / ( time[ipos2] - time[ipos1]);
1022
1023 }
1024
1025}
1026
a328fff9 1027//_____________________________________________________________________________
1028void AliTRDparameter::SampleTimeStruct()
1029{
1030 //
1031 // Samples the timing structure of a drift cell
1032 // Drift Time data calculated with Garfield (by C.Lippmann)
1033 //
1034
1035 const Int_t ktimebin = 38;
682619b9 1036 const Int_t kZbin = 11;
1037
1038 Float_t time[ktimebin][kZbin] =
1039 {{0.10176, 0.10197, 0.10356, 0.10479, 0.10497, 0.10500, 0.10500,
1040 0.10698, 0.12213, 0.16118, 0.22154},
1041 {0.07500, 0.07500, 0.07500, 0.07500, 0.07500, 0.07500, 0.07509,
1042 0.08010, 0.10176, 0.13244, 0.22848},
1043 {0.04500, 0.04500, 0.04500, 0.04500, 0.04500, 0.04503, 0.05130,
1044 0.07398, 0.07557, 0.11064, 0.26833},
1045 {0.01500, 0.01500, 0.01500, 0.01500, 0.03219, 0.04500, 0.04500,
1046 0.04500, 0.07482, 0.11628, 0.40846},
1047 {0.01500, 0.01500, 0.01500, 0.01500, 0.02886, 0.04500, 0.04500,
1048 0.04512, 0.07503, 0.13230, 0.43602},
1049 {0.04500, 0.04500, 0.04500, 0.04500, 0.04500, 0.04500, 0.04872,
1050 0.07437, 0.07635, 0.11925, 0.30372},
1051 {0.07500, 0.07497, 0.07500, 0.07500, 0.07500, 0.07500, 0.07500,
1052 0.08232, 0.10458, 0.14100, 0.24799},
1053 {0.10098, 0.10125, 0.10218, 0.10299, 0.10467, 0.15000, 0.10959,
1054 0.12207, 0.14241, 0.19485, 0.25688},
1055 {0.13500, 0.13503, 0.13509, 0.13782, 0.16182, 0.25076, 0.19593,
1056 0.19254, 0.21060, 0.25561, 0.29469},
1057 {0.19305, 0.19374, 0.19512, 0.19974, 0.22020, 0.29420, 0.25992,
1058 0.25461, 0.27515, 0.31938, 0.35238},
1059 {0.25143, 0.25320, 0.25476, 0.26166, 0.27993, 0.34470, 0.32160,
1060 0.31623, 0.33687, 0.37420, 0.40797},
1061 {0.31044, 0.31197, 0.31491, 0.32280, 0.34125, 0.39702, 0.38401,
1062 0.37794, 0.39545, 0.43494, 0.46956},
1063 {0.37128, 0.37146, 0.37668, 0.38424, 0.40326, 0.45271, 0.44643,
1064 0.43783, 0.45882, 0.49726, 0.52348},
1065 {0.43065, 0.43242, 0.43722, 0.44379, 0.46551, 0.51167, 0.50685,
1066 0.49924, 0.52016, 0.55609, 0.57829},
1067 {0.49161, 0.49329, 0.49725, 0.50629, 0.52803, 0.56777, 0.56755,
1068 0.56057, 0.58062, 0.61293, 0.63839},
1069 {0.55284, 0.55311, 0.55743, 0.56703, 0.58824, 0.62747, 0.62833,
1070 0.62083, 0.63899, 0.67412, 0.70006},
1071 {0.61323, 0.61494, 0.61782, 0.62802, 0.65025, 0.68675, 0.68847,
1072 0.68248, 0.70177, 0.73588, 0.76314},
1073 {0.67293, 0.67482, 0.68034, 0.68814, 0.71169, 0.74695, 0.75024,
1074 0.74462, 0.76083, 0.79194, 0.81463},
1075 {0.73470, 0.73665, 0.74013, 0.75054, 0.77107, 0.80588, 0.81071,
1076 0.80508, 0.82337, 0.85380, 0.87045},
1077 {0.79509, 0.79677, 0.80115, 0.81150, 0.83400, 0.86360, 0.87159,
1078 0.86608, 0.88323, 0.91756, 0.93721},
1079 {0.85695, 0.85800, 0.86307, 0.87276, 0.89599, 0.92639, 0.93266,
1080 0.92962, 0.94291, 0.97445, 0.99502},
1081 {0.91680, 0.91875, 0.92400, 0.93219, 0.95555, 0.98729, 0.99281,
1082 0.99145, 1.00683, 1.03343, 1.05299},
1083 {0.97833, 0.98058, 0.98421, 0.99357, 1.01602, 1.04438, 1.05299,
1084 1.05168, 1.06528, 1.09644, 1.11557},
1085 {1.03938, 1.04061, 1.04568, 1.05617, 1.07923, 1.10648, 1.11466,
1086 1.11385, 1.12606, 1.15817, 1.17352},
1087 {1.10166, 1.10268, 1.10766, 1.11785, 1.14014, 1.16717, 1.17496,
1088 1.17644, 1.19115, 1.21622, 1.23605},
1089 {1.16334, 1.16364, 1.16985, 1.17996, 1.19935, 1.22911, 1.23561,
1090 1.23526, 1.25144, 1.27936, 1.29018},
1091 {1.22403, 1.22598, 1.23000, 1.24092, 1.26173, 1.28817, 1.29991,
1092 1.29698, 1.31238, 1.33793, 1.35595},
1093 {1.28556, 1.28898, 1.29300, 1.30356, 1.32617, 1.34912, 1.35925,
1094 1.35835, 1.37312, 1.40111, 1.41497},
1095 {1.34787, 1.34877, 1.35468, 1.36449, 1.38654, 1.41077, 1.41975,
1096 1.42136, 1.43843, 1.45897, 1.48025},
1097 {1.41045, 1.41171, 1.41555, 1.42656, 1.44698, 1.47026, 1.48324,
1098 1.48479, 1.50180, 1.52354, 1.53392},
1099 {1.47090, 1.47354, 1.47647, 1.48905, 1.51211, 1.53038, 1.54349,
1100 1.54803, 1.56136, 1.58662, 1.60480},
1101 {1.53357, 1.53429, 1.54017, 1.55116, 1.57247, 1.59498, 1.60576,
1102 1.60746, 1.62389, 1.64239, 1.66245},
1103 {1.59522, 1.59636, 1.60247, 1.61385, 1.63641, 1.65659, 1.66619,
1104 1.67044, 1.68213, 1.70379, 1.72372},
1105 {1.65867, 1.65711, 1.66320, 1.67434, 1.69603, 1.71788, 1.72927,
1106 1.73214, 1.74713, 1.76485, 1.77775},
1107 {1.71825, 1.72074, 1.72497, 1.73712, 1.75690, 1.77755, 1.78711,
1108 1.79849, 1.80587, 1.82576, 1.84050},
1109 {1.78194, 1.78320, 1.78875, 1.80047, 1.81716, 1.84044, 1.85328,
1110 1.85511, 1.86911, 1.88919, 1.89926},
1111 {1.84335, 1.84431, 1.84938, 1.86063, 1.88022, 1.90234, 1.91177,
1112 1.92173, 1.93202, 1.95097, 1.97105},
1113 {1.90491, 1.90665, 1.91303, 1.92575, 1.94432, 1.96400, 1.97574,
1114 1.98378, 1.99368, 2.01964, 2.02605}};
a328fff9 1115
1116 if (fTimeStruct) delete [] fTimeStruct;
1117 fTimeStruct = new Float_t[ktimebin*kZbin];
1118
1119 for (Int_t ctrt = 0; ctrt<ktimebin; ctrt++)
1120 for (Int_t ctrz = 0; ctrz<kZbin; ctrz++) {
682619b9 1121 fTimeStruct[ctrt+ctrz*ktimebin] = time[ctrt][ctrz];
a328fff9 1122 }
1123
1124}
1125
17b26de4 1126//_____________________________________________________________________________
1127void AliTRDparameter::SamplePRF()
1128{
1129 //
1130 // Samples the pad response function
1131 //
1132
1133 const Int_t kNplan = AliTRDgeometry::kNplan;
1134 const Int_t kPRFbin = 61;
1135
a328fff9 1136 Float_t prf[kNplan][kPRFbin] = { {0.0267681, 0.031056, 0.036111, 0.0419945, 0.0489578,
1137 0.0571508, 0.0666984, 0.078064, 0.0912508, 0.106579,
1138 0.124429, 0.144803, 0.168187, 0.194645, 0.224143,
1139 0.25717, 0.293146, 0.331943, 0.373128, 0.415798,
1140 0.459149, 0.502028, 0.54349, 0.582109, 0.617335,
1141 0.648306, 0.674189, 0.695082, 0.710222, 0.719711,
1142 0.723751, 0.719711, 0.710222, 0.695082, 0.674189,
1143 0.648306, 0.617335, 0.582109, 0.54349, 0.502028,
1144 0.459149, 0.415798, 0.373128, 0.331943, 0.293146,
1145 0.25717, 0.224143, 0.194645, 0.168187, 0.144803,
1146 0.124429, 0.106579, 0.0912508, 0.078064, 0.0666984,
1147 0.0571508, 0.0489578, 0.0419945, 0.036111, 0.031056,
1148 0.0267681},
1149 {0.0240831, 0.0281427, 0.0328987, 0.0384801, 0.0451249,
1150 0.0530119, 0.0623139, 0.0733268, 0.0864095, 0.101704,
1151 0.119473, 0.140308, 0.164191, 0.191333, 0.222141,
1152 0.256551, 0.29433, 0.335294, 0.37892, 0.424206,
1153 0.470202, 0.515621, 0.559255, 0.600068, 0.636821,
1154 0.668836, 0.695767, 0.717076, 0.73233, 0.741818,
1155 0.745491, 0.741818, 0.73233, 0.717076, 0.695767,
1156 0.668836, 0.636821, 0.600068, 0.559255, 0.515621,
1157 0.470202, 0.424206, 0.37892, 0.335294, 0.29433,
1158 0.256551, 0.222141, 0.191333, 0.164191, 0.140308,
1159 0.119473, 0.101704, 0.0864095, 0.0733268, 0.0623139,
1160 0.0530119, 0.0451249, 0.0384801, 0.0328987, 0.0281427,
1161 0.0240831},
1162 {0.0206855, 0.0243011, 0.0285603, 0.0335902, 0.0396015,
1163 0.0467843, 0.0553424, 0.0655316, 0.0776688, 0.0921782,
1164 0.109269, 0.129259, 0.152459, 0.179369, 0.210091,
1165 0.244681, 0.283099, 0.325176, 0.370396, 0.417732,
1166 0.466086, 0.514192, 0.560518, 0.603805, 0.643028,
1167 0.67738, 0.706095, 0.728728, 0.745307, 0.755731,
1168 0.759954, 0.755731, 0.745307, 0.728728, 0.706095,
1169 0.67738, 0.643028, 0.603805, 0.560518, 0.514192,
1170 0.466086, 0.417732, 0.370396, 0.325176, 0.283099,
1171 0.244681, 0.210091, 0.179369, 0.152459, 0.129259,
1172 0.109269, 0.0921782, 0.0776688, 0.0655316, 0.0553424,
1173 0.0467843, 0.0396015, 0.0335902, 0.0285603, 0.0243011,
1174 0.0206855},
1175 {0.0186168, 0.0219999, 0.0260102, 0.0307769, 0.0364947,
1176 0.0433655, 0.0516213, 0.0615466, 0.0734611, 0.0877121,
1177 0.104666, 0.124855, 0.14853, 0.176033, 0.207639,
1178 0.243511, 0.283633, 0.327786, 0.37537, 0.425281,
1179 0.476227, 0.526727, 0.575268, 0.620462, 0.66101,
1180 0.69611, 0.725313, 0.748317, 0.764969, 0.775206,
1181 0.779006, 0.775206, 0.764969, 0.748317, 0.725313,
1182 0.69611, 0.66101, 0.620462, 0.575268, 0.526727,
1183 0.476227, 0.425281, 0.37537, 0.327786, 0.283633,
1184 0.243511, 0.207639, 0.176033, 0.14853, 0.124855,
1185 0.104666, 0.0877121, 0.0734611, 0.0615466, 0.0516213,
1186 0.0433655, 0.0364947, 0.0307769, 0.0260102, 0.0219999,
1187 0.0186168, },
1188 {0.0159737, 0.0189921, 0.0225916, 0.0268927, 0.0320634,
1189 0.0382995, 0.0458393, 0.0549765, 0.0660512, 0.0794439,
1190 0.095565, 0.114844, 0.137713, 0.164586, 0.195824,
1191 0.231681, 0.272223, 0.31727, 0.366248, 0.418078,
1192 0.471358, 0.524425, 0.575561, 0.623193, 0.666055,
1193 0.703243, 0.734192, 0.7586, 0.776331, 0.787347,
1194 0.791646, 0.787347, 0.776331, 0.7586, 0.734192,
1195 0.703243, 0.666055, 0.623193, 0.575561, 0.524425,
1196 0.471358, 0.418078, 0.366248, 0.31727, 0.272223,
1197 0.231681, 0.195824, 0.164586, 0.137713, 0.114844,
1198 0.095565, 0.0794439, 0.0660512, 0.0549765, 0.0458393,
1199 0.0382995, 0.0320634, 0.0268927, 0.0225916, 0.0189921,
1200 0.0159737},
1201 {0.0143532, 0.0171745, 0.0205622, 0.024635, 0.0295461,
1202 0.0354915, 0.0427198, 0.0515368, 0.0623058, 0.0754416,
1203 0.0913994, 0.110662, 0.133721, 0.161057, 0.193093,
1204 0.230138, 0.27229, 0.319327, 0.370597, 0.424961,
1205 0.480832, 0.536333, 0.589555, 0.638805, 0.682781,
1206 0.720621, 0.751845, 0.776243, 0.793764, 0.804426,
1207 0.808259, 0.804426, 0.793764, 0.776243, 0.751845,
1208 0.720621, 0.682781, 0.638805, 0.589555, 0.536333,
1209 0.480832, 0.424961, 0.370597, 0.319327, 0.27229,
1210 0.230138, 0.193093, 0.161057, 0.133721, 0.110662,
1211 0.0913994, 0.0754416, 0.0623058, 0.0515368, 0.0427198,
1212 0.0354915, 0.0295461, 0.024635, 0.0205622, 0.0171745,
1213 0.0143532} };
17b26de4 1214
1215 // More sampling precision with linear interpolation
1216 fPRFlo = -1.5;
1217 fPRFhi = 1.5;
1218 Float_t pad[kPRFbin];
1219 Int_t sPRFbin = kPRFbin;
1220 Float_t sPRFwid = (fPRFhi - fPRFlo) / ((Float_t) sPRFbin);
1221 for (Int_t iPad = 0; iPad < sPRFbin; iPad++) {
1222 pad[iPad] = ((Float_t) iPad + 0.5) * sPRFwid + fPRFlo;
1223 }
1224 fPRFbin = 500;
1225 fPRFwid = (fPRFhi - fPRFlo) / ((Float_t) fPRFbin);
1226 fPRFpad = ((Int_t) (1.0 / fPRFwid));
1227
1228 if (fPRFsmp) delete [] fPRFsmp;
1229 fPRFsmp = new Float_t[kNplan*fPRFbin];
1230
1231 Int_t ipos1;
1232 Int_t ipos2;
1233 Float_t diff;
1234
1235 for (Int_t iPla = 0; iPla < kNplan; iPla++) {
1236
1237 for (Int_t iBin = 0; iBin < fPRFbin; iBin++) {
1238
1239 Float_t bin = (((Float_t) iBin) + 0.5) * fPRFwid + fPRFlo;
1240 ipos1 = ipos2 = 0;
1241 diff = 0;
1242 do {
1243 diff = bin - pad[ipos2++];
8e64dd77 1244 } while ((diff > 0) && (ipos2 < kPRFbin));
1245 if (ipos2 == kPRFbin) {
1246 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2-1];
1247 }
1248 else if (ipos2 == 1) {
1249 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2-1];
1250 }
1251 else {
1252 ipos2--;
1253 if (ipos2 >= kPRFbin) ipos2 = kPRFbin - 1;
1254 ipos1 = ipos2 - 1;
1255 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2]
1256 + diff * (prf[iPla][ipos2] - prf[iPla][ipos1])
1257 / sPRFwid;
1258 }
17b26de4 1259
1260 }
1261 }
1262
1263}
1264
1265//_____________________________________________________________________________
1266void AliTRDparameter::FillLUT()
1267{
1268 //
1269 // Create the LUT
1270 //
1271
1272 const Int_t kNplan = AliTRDgeometry::kNplan;
1273 const Int_t kNlut = 128;
1274
1275 fLUTbin = kNplan * kNlut;
1276
1277 // The lookup table from Bogdan
1278 Float_t lut[kNplan][kNlut] = {
1279 {
1280 0.0070, 0.0150, 0.0224, 0.0298, 0.0374, 0.0454, 0.0533, 0.0611,
1281 0.0684, 0.0755, 0.0827, 0.0900, 0.0975, 0.1049, 0.1120, 0.1187,
1282 0.1253, 0.1318, 0.1385, 0.1453, 0.1519, 0.1584, 0.1646, 0.1704,
1283 0.1762, 0.1821, 0.1879, 0.1938, 0.1996, 0.2053, 0.2108, 0.2160,
1284 0.2210, 0.2260, 0.2310, 0.2361, 0.2411, 0.2461, 0.2509, 0.2557,
1285 0.2602, 0.2646, 0.2689, 0.2732, 0.2774, 0.2816, 0.2859, 0.2901,
1286 0.2942, 0.2983, 0.3022, 0.3061, 0.3099, 0.3136, 0.3172, 0.3207,
1287 0.3242, 0.3278, 0.3312, 0.3347, 0.3382, 0.3416, 0.3450, 0.3483,
1288 0.3515, 0.3547, 0.3579, 0.3609, 0.3639, 0.3669, 0.3698, 0.3727,
1289 0.3756, 0.3785, 0.3813, 0.3842, 0.3870, 0.3898, 0.3926, 0.3952,
1290 0.3979, 0.4005, 0.4032, 0.4057, 0.4082, 0.4108, 0.4132, 0.4157,
1291 0.4181, 0.4205, 0.4228, 0.4252, 0.4275, 0.4299, 0.4322, 0.4345,
1292 0.4367, 0.4390, 0.4412, 0.4434, 0.4456, 0.4478, 0.4499, 0.4520,
1293 0.4541, 0.4562, 0.4583, 0.4603, 0.4623, 0.4643, 0.4663, 0.4683,
1294 0.4702, 0.4722, 0.4741, 0.4758, 0.4774, 0.4790, 0.4805, 0.4824,
1295 0.4844, 0.4863, 0.4883, 0.4902, 0.4921, 0.4940, 0.4959, 0.4978
1296 },
1297 {
1298 0.0072, 0.0156, 0.0235, 0.0313, 0.0394, 0.0478, 0.0561, 0.0642,
1299 0.0718, 0.0792, 0.0868, 0.0947, 0.1025, 0.1101, 0.1172, 0.1241,
1300 0.1309, 0.1378, 0.1449, 0.1518, 0.1586, 0.1650, 0.1710, 0.1770,
1301 0.1830, 0.1891, 0.1952, 0.2011, 0.2070, 0.2125, 0.2177, 0.2229,
1302 0.2280, 0.2332, 0.2383, 0.2435, 0.2484, 0.2533, 0.2581, 0.2627,
1303 0.2670, 0.2714, 0.2757, 0.2799, 0.2842, 0.2884, 0.2927, 0.2968,
1304 0.3008, 0.3048, 0.3086, 0.3123, 0.3159, 0.3195, 0.3231, 0.3266,
1305 0.3301, 0.3335, 0.3370, 0.3404, 0.3438, 0.3471, 0.3504, 0.3536,
1306 0.3567, 0.3598, 0.3628, 0.3657, 0.3686, 0.3715, 0.3744, 0.3772,
1307 0.3800, 0.3828, 0.3856, 0.3884, 0.3911, 0.3938, 0.3965, 0.3991,
1308 0.4016, 0.4042, 0.4067, 0.4092, 0.4116, 0.4140, 0.4164, 0.4187,
1309 0.4211, 0.4234, 0.4257, 0.4280, 0.4302, 0.4325, 0.4347, 0.4369,
1310 0.4391, 0.4413, 0.4434, 0.4456, 0.4477, 0.4497, 0.4518, 0.4538,
1311 0.4558, 0.4578, 0.4598, 0.4618, 0.4637, 0.4656, 0.4675, 0.4694,
1312 0.4713, 0.4732, 0.4750, 0.4766, 0.4781, 0.4797, 0.4813, 0.4832,
1313 0.4851, 0.4870, 0.4888, 0.4906, 0.4925, 0.4942, 0.4960, 0.4978
1314 },
1315 {
1316 0.0075, 0.0163, 0.0246, 0.0328, 0.0415, 0.0504, 0.0592, 0.0674,
1317 0.0753, 0.0832, 0.0914, 0.0996, 0.1077, 0.1154, 0.1225, 0.1296,
1318 0.1369, 0.1442, 0.1515, 0.1585, 0.1652, 0.1714, 0.1776, 0.1839,
1319 0.1902, 0.1965, 0.2025, 0.2085, 0.2141, 0.2194, 0.2247, 0.2299,
1320 0.2352, 0.2405, 0.2457, 0.2507, 0.2557, 0.2604, 0.2649, 0.2693,
1321 0.2737, 0.2780, 0.2823, 0.2867, 0.2909, 0.2951, 0.2992, 0.3033,
1322 0.3072, 0.3110, 0.3146, 0.3182, 0.3218, 0.3253, 0.3288, 0.3323,
1323 0.3357, 0.3392, 0.3426, 0.3459, 0.3492, 0.3524, 0.3555, 0.3586,
1324 0.3616, 0.3645, 0.3674, 0.3703, 0.3731, 0.3759, 0.3787, 0.3815,
1325 0.3843, 0.3870, 0.3897, 0.3925, 0.3950, 0.3976, 0.4002, 0.4027,
1326 0.4052, 0.4076, 0.4101, 0.4124, 0.4148, 0.4171, 0.4194, 0.4217,
1327 0.4239, 0.4262, 0.4284, 0.4306, 0.4328, 0.4350, 0.4371, 0.4393,
1328 0.4414, 0.4435, 0.4455, 0.4476, 0.4496, 0.4516, 0.4536, 0.4555,
1329 0.4575, 0.4594, 0.4613, 0.4632, 0.4650, 0.4669, 0.4687, 0.4705,
1330 0.4723, 0.4741, 0.4758, 0.4773, 0.4789, 0.4804, 0.4821, 0.4839,
1331 0.4857, 0.4875, 0.4893, 0.4910, 0.4928, 0.4945, 0.4961, 0.4978
1332 },
1333 {
1334 0.0078, 0.0171, 0.0258, 0.0345, 0.0438, 0.0532, 0.0624, 0.0708,
1335 0.0791, 0.0875, 0.0962, 0.1048, 0.1130, 0.1206, 0.1281, 0.1356,
1336 0.1432, 0.1508, 0.1582, 0.1651, 0.1716, 0.1780, 0.1845, 0.1910,
1337 0.1975, 0.2038, 0.2099, 0.2155, 0.2210, 0.2263, 0.2317, 0.2371,
1338 0.2425, 0.2477, 0.2528, 0.2578, 0.2626, 0.2671, 0.2715, 0.2759,
1339 0.2803, 0.2846, 0.2890, 0.2933, 0.2975, 0.3016, 0.3056, 0.3095,
1340 0.3132, 0.3168, 0.3204, 0.3239, 0.3274, 0.3309, 0.3344, 0.3378,
1341 0.3412, 0.3446, 0.3479, 0.3511, 0.3543, 0.3574, 0.3603, 0.3633,
1342 0.3662, 0.3690, 0.3718, 0.3747, 0.3774, 0.3802, 0.3829, 0.3857,
1343 0.3883, 0.3910, 0.3936, 0.3962, 0.3987, 0.4012, 0.4037, 0.4061,
1344 0.4085, 0.4109, 0.4132, 0.4155, 0.4177, 0.4200, 0.4222, 0.4244,
1345 0.4266, 0.4288, 0.4309, 0.4331, 0.4352, 0.4373, 0.4394, 0.4414,
1346 0.4435, 0.4455, 0.4475, 0.4494, 0.4514, 0.4533, 0.4552, 0.4571,
1347 0.4590, 0.4608, 0.4626, 0.4645, 0.4662, 0.4680, 0.4698, 0.4715,
1348 0.4733, 0.4750, 0.4766, 0.4781, 0.4796, 0.4812, 0.4829, 0.4846,
1349 0.4863, 0.4880, 0.4897, 0.4914, 0.4930, 0.4946, 0.4963, 0.4979
1350 },
1351 {
1352 0.0081, 0.0178, 0.0270, 0.0364, 0.0463, 0.0562, 0.0656, 0.0744,
1353 0.0831, 0.0921, 0.1013, 0.1102, 0.1183, 0.1261, 0.1339, 0.1419,
1354 0.1499, 0.1576, 0.1648, 0.1715, 0.1782, 0.1849, 0.1917, 0.1984,
1355 0.2048, 0.2110, 0.2167, 0.2223, 0.2278, 0.2333, 0.2389, 0.2444,
1356 0.2497, 0.2548, 0.2598, 0.2645, 0.2691, 0.2735, 0.2780, 0.2824,
1357 0.2868, 0.2912, 0.2955, 0.2997, 0.3038, 0.3078, 0.3116, 0.3152,
1358 0.3188, 0.3224, 0.3259, 0.3294, 0.3329, 0.3364, 0.3398, 0.3432,
1359 0.3465, 0.3497, 0.3529, 0.3561, 0.3591, 0.3620, 0.3649, 0.3677,
1360 0.3705, 0.3733, 0.3761, 0.3788, 0.3816, 0.3843, 0.3869, 0.3896,
1361 0.3922, 0.3948, 0.3973, 0.3998, 0.4022, 0.4047, 0.4070, 0.4094,
1362 0.4117, 0.4139, 0.4162, 0.4184, 0.4206, 0.4227, 0.4249, 0.4270,
1363 0.4291, 0.4313, 0.4334, 0.4354, 0.4375, 0.4395, 0.4415, 0.4435,
1364 0.4455, 0.4474, 0.4493, 0.4512, 0.4531, 0.4550, 0.4568, 0.4586,
1365 0.4604, 0.4622, 0.4639, 0.4657, 0.4674, 0.4691, 0.4708, 0.4725,
1366 0.4742, 0.4758, 0.4773, 0.4788, 0.4803, 0.4819, 0.4836, 0.4852,
1367 0.4869, 0.4885, 0.4901, 0.4917, 0.4933, 0.4948, 0.4964, 0.4979
1368 },
1369 {
1370 0.0085, 0.0189, 0.0288, 0.0389, 0.0497, 0.0603, 0.0699, 0.0792,
1371 0.0887, 0.0985, 0.1082, 0.1170, 0.1253, 0.1336, 0.1421, 0.1505,
1372 0.1587, 0.1662, 0.1733, 0.1803, 0.1874, 0.1945, 0.2014, 0.2081,
1373 0.2143, 0.2201, 0.2259, 0.2316, 0.2374, 0.2431, 0.2487, 0.2541,
1374 0.2593, 0.2642, 0.2689, 0.2735, 0.2781, 0.2826, 0.2872, 0.2917,
1375 0.2961, 0.3003, 0.3045, 0.3086, 0.3125, 0.3162, 0.3198, 0.3235,
1376 0.3270, 0.3306, 0.3342, 0.3377, 0.3411, 0.3446, 0.3479, 0.3511,
1377 0.3543, 0.3575, 0.3605, 0.3634, 0.3663, 0.3691, 0.3720, 0.3748,
1378 0.3775, 0.3803, 0.3830, 0.3857, 0.3884, 0.3911, 0.3937, 0.3962,
1379 0.3987, 0.4012, 0.4036, 0.4060, 0.4084, 0.4107, 0.4129, 0.4152,
1380 0.4174, 0.4196, 0.4218, 0.4239, 0.4261, 0.4282, 0.4303, 0.4324,
1381 0.4344, 0.4365, 0.4385, 0.4405, 0.4425, 0.4445, 0.4464, 0.4483,
1382 0.4502, 0.4521, 0.4539, 0.4558, 0.4576, 0.4593, 0.4611, 0.4629,
1383 0.4646, 0.4663, 0.4680, 0.4697, 0.4714, 0.4730, 0.4747, 0.4759,
1384 0.4769, 0.4780, 0.4790, 0.4800, 0.4811, 0.4827, 0.4843, 0.4859,
1385 0.4874, 0.4889, 0.4905, 0.4920, 0.4935, 0.4950, 0.4965, 0.4979
1386 }
1387 };
1388
1389 if (fLUT) delete [] fLUT;
1390 fLUT = new Float_t[fLUTbin];
1391
1392 for (Int_t iplan = 0; iplan < kNplan; iplan++) {
1393 for (Int_t ilut = 0; ilut < kNlut; ilut++) {
cfed70d4 1394 fLUT[iplan*kNlut+ilut] = lut[iplan][ilut];
17b26de4 1395 }
1396 }
1397
1398}
1399
1400//_____________________________________________________________________________
1401void AliTRDparameter::SetTiltingAngle(Float_t v)
1402{
1403 //
1404 // Set the tilting angle for the readout pads
1405 //
1406
1407 fTiltingAngle = TMath::Tan(TMath::Pi()/180.0 * v);
1408
1409}
1410
1411//_____________________________________________________________________________
1412Float_t AliTRDparameter::GetTiltingAngle() const
1413{
1414 //
1415 // Get the tilting angle for the readout pads
1416 //
1417
ab35e265 1418 return 180.0 / TMath::Pi() * TMath::ATan(fTiltingAngle);
17b26de4 1419
1420}
1421
1422//_____________________________________________________________________________
1423Float_t AliTRDparameter::GetDiffusionL(Float_t vd, Float_t b)
1424{
1425 //
1426 // Returns the longitudinal diffusion coefficient for a given drift
1427 // velocity <vd> and a B-field <b> for Xe/CO2 (15%).
1428 // The values are according to a GARFIELD simulation.
1429 //
1430
1431 const Int_t kNb = 5;
1432 Float_t p0[kNb] = { 0.007440, 0.007493, 0.007513, 0.007672, 0.007831 };
1433 Float_t p1[kNb] = { 0.019252, 0.018912, 0.018636, 0.018012, 0.017343 };
1434 Float_t p2[kNb] = { -0.005042, -0.004926, -0.004867, -0.004650, -0.004424 };
1435 Float_t p3[kNb] = { 0.000195, 0.000189, 0.000195, 0.000182, 0.000169 };
1436
1437 Int_t ib = ((Int_t) (10 * (b - 0.15)));
1438 ib = TMath::Max( 0,ib);
1439 ib = TMath::Min(kNb,ib);
1440
1441 Float_t diff = p0[ib]
1442 + p1[ib] * vd
1443 + p2[ib] * vd*vd
1444 + p3[ib] * vd*vd*vd;
1445
1446 return diff;
1447
1448}
1449
1450//_____________________________________________________________________________
1451Float_t AliTRDparameter::GetDiffusionT(Float_t vd, Float_t b)
1452{
1453 //
1454 // Returns the transverse diffusion coefficient for a given drift
1455 // velocity <vd> and a B-field <b> for Xe/CO2 (15%).
1456 // The values are according to a GARFIELD simulation.
1457 //
1458
1459 const Int_t kNb = 5;
1460 Float_t p0[kNb] = { 0.009550, 0.009599, 0.009674, 0.009757, 0.009850 };
1461 Float_t p1[kNb] = { 0.006667, 0.006539, 0.006359, 0.006153, 0.005925 };
1462 Float_t p2[kNb] = { -0.000853, -0.000798, -0.000721, -0.000635, -0.000541 };
1463 Float_t p3[kNb] = { 0.000131, 0.000122, 0.000111, 0.000098, 0.000085 };
1464
1465 Int_t ib = ((Int_t) (10 * (b - 0.15)));
1466 ib = TMath::Max( 0,ib);
1467 ib = TMath::Min(kNb,ib);
1468
1469 Float_t diff = p0[ib]
1470 + p1[ib] * vd
1471 + p2[ib] * vd*vd
1472 + p3[ib] * vd*vd*vd;
1473
1474 return diff;
1475
1476}
1477
1478//_____________________________________________________________________________
1479Float_t AliTRDparameter::GetOmegaTau(Float_t vd, Float_t b)
1480{
1481 //
1482 // Returns omega*tau (tan(Lorentz-angle)) for a given drift velocity <vd>
1483 // and a B-field <b> for Xe/CO2 (15%).
1484 // The values are according to a GARFIELD simulation.
1485 //
1486
1487 const Int_t kNb = 5;
1488 Float_t p0[kNb] = { 0.004810, 0.007412, 0.010252, 0.013409, 0.016888 };
1489 Float_t p1[kNb] = { 0.054875, 0.081534, 0.107333, 0.131983, 0.155455 };
1490 Float_t p2[kNb] = { -0.008682, -0.012896, -0.016987, -0.020880, -0.024623 };
1491 Float_t p3[kNb] = { 0.000155, 0.000238, 0.000330, 0.000428, 0.000541 };
1492
1493 Int_t ib = ((Int_t) (10 * (b - 0.15)));
1494 ib = TMath::Max( 0,ib);
1495 ib = TMath::Min(kNb,ib);
1496
1497 Float_t alphaL = p0[ib]
1498 + p1[ib] * vd
1499 + p2[ib] * vd*vd
1500 + p3[ib] * vd*vd*vd;
1501
1502 return TMath::Tan(alphaL);
1503
1504}
1505
1506//_____________________________________________________________________________
0a29d0f1 1507Double_t AliTRDparameter::LUTposition(Int_t iplane, Double_t ampL
1508 , Double_t ampC
1509 , Double_t ampR) const
17b26de4 1510{
1511 //
1512 // Calculates the cluster position using the lookup table.
1513 // Method provided by Bogdan Vulpescu.
1514 //
1515
1516 const Int_t kNplan = AliTRDgeometry::kNplan;
cfed70d4 1517 const Int_t kNlut = 128;
17b26de4 1518
1519 Double_t pos;
1520 Double_t x = 0.0;
1521 Double_t xmin;
1522 Double_t xmax;
1523 Double_t xwid;
1524
1525 Int_t side = 0;
1526 Int_t ix;
1527
1528 Double_t xMin[kNplan] = { 0.006492, 0.006377, 0.006258, 0.006144, 0.006030, 0.005980 };
1529 Double_t xMax[kNplan] = { 0.960351, 0.965870, 0.970445, 0.974352, 0.977667, 0.996101 };
1530
1531 if (ampL > ampR) {
1532 x = (ampL - ampR) / ampC;
1533 side = -1;
1534 }
1535 else if (ampL < ampR) {
1536 x = (ampR - ampL) / ampC;
1537 side = +1;
1538 }
1539
1540 if (ampL != ampR) {
1541
1542 xmin = xMin[iplane] + 0.000005;
1543 xmax = xMax[iplane] - 0.000005;
1544 xwid = (xmax - xmin) / 127.0;
1545
1546 if (x < xmin) {
1547 pos = 0.0000;
1548 }
1549 else if (x > xmax) {
1550 pos = side * 0.5000;
1551 }
1552 else {
1553 ix = (Int_t) ((x - xmin) / xwid);
cfed70d4 1554 pos = side * fLUT[iplane*kNlut+ix];
17b26de4 1555 }
1556
1557 }
1558 else {
1559
1560 pos = 0.0;
1561
1562 }
1563
1564 return pos;
1565
1566}