]> 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
set magnetic field for tracks
[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;
87 fStaggeringOn = 0;
88 fAnodeWireOffset = 0.0;
89
17b26de4 90}
91
92//_____________________________________________________________________________
93AliTRDparameter::AliTRDparameter(const Text_t *name, const Text_t *title)
94 :TNamed(name,title)
95{
96 //
97 // AliTRDparameter constructor
98 //
99
100 fGeo = new AliTRDgeometryFull();
67689223 101 fPRFsmp = 0;
102 fTRFsmp = 0;
103 fCTsmp = 0;
104 fGasGain = 0.0;
105 fNoise = 0.0;
106 fChipGain = 0.0;
107 fADCoutRange = 0.0;
108 fADCinRange = 0.0;
109 fADCthreshold = 0;
7e582e9f 110 fADCbaseline = 0;
67689223 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
a328fff9 144 fTimeStructOn = 0;
145 fTimeStruct = 0;
146 fStaggeringOn = 0;
147 fAnodeWireOffset = 0.0;
148
17b26de4 149 Init();
150
151}
152
153
154//_____________________________________________________________________________
73ae7b59 155AliTRDparameter::AliTRDparameter(const AliTRDparameter &p):TNamed(p)
17b26de4 156{
157 //
158 // AliTRDparameter copy constructor
159 //
160
161 ((AliTRDparameter &) p).Copy(*this);
162
163}
164
165///_____________________________________________________________________________
166AliTRDparameter::~AliTRDparameter()
167{
168 //
169 // AliTRDparameter destructor
170 //
171
172 if (fTRFsmp) {
173 delete [] fTRFsmp;
174 fTRFsmp = 0;
175 }
176
177 if (fPRFsmp) {
178 delete [] fPRFsmp;
179 fPRFsmp = 0;
180 }
181
182 if (fCTsmp) {
183 delete [] fCTsmp;
184 fCTsmp = 0;
185 }
186
187 if (fLUT) {
188 delete [] fLUT;
189 fLUT = 0;
190 }
191
192 if (fGeo) {
193 delete fGeo;
194 fGeo = 0;
195 }
196
a328fff9 197 if (fTimeStruct) {
198 delete [] fTimeStruct;
199 fTimeStruct = 0;
200 }
201
17b26de4 202}
203
204//_____________________________________________________________________________
205AliTRDparameter &AliTRDparameter::operator=(const AliTRDparameter &p)
206{
207 //
208 // Assignment operator
209 //
210
211 if (this != &p) ((AliTRDparameter &) p).Copy(*this);
212 return *this;
213
214}
215
216//_____________________________________________________________________________
217void AliTRDparameter::Copy(TObject &p)
218{
219 //
220 // Copy function
221 //
222
223 Int_t iBin;
224
225 ((AliTRDparameter &) p).fGasGain = fGasGain;
226 ((AliTRDparameter &) p).fNoise = fNoise;
227 ((AliTRDparameter &) p).fChipGain = fChipGain;
228 ((AliTRDparameter &) p).fADCoutRange = fADCoutRange;
229 ((AliTRDparameter &) p).fADCinRange = fADCinRange;
230 ((AliTRDparameter &) p).fADCthreshold = fADCthreshold;
7e582e9f 231 ((AliTRDparameter &) p).fADCbaseline = fADCbaseline;
17b26de4 232 ((AliTRDparameter &) p).fDiffusionOn = fDiffusionOn;
233 ((AliTRDparameter &) p).fDiffusionT = fDiffusionT;
234 ((AliTRDparameter &) p).fDiffusionL = fDiffusionL;
235 ((AliTRDparameter &) p).fElAttachOn = fElAttachOn;
236 ((AliTRDparameter &) p).fElAttachProp = fElAttachProp;
237 ((AliTRDparameter &) p).fExBOn = fExBOn;
238 ((AliTRDparameter &) p).fOmegaTau = fOmegaTau;
239 ((AliTRDparameter &) p).fLorentzFactor = fLorentzFactor;
240 ((AliTRDparameter &) p).fDriftVelocity = fDriftVelocity;
241 ((AliTRDparameter &) p).fPadCoupling = fPadCoupling;
242 ((AliTRDparameter &) p).fTimeCoupling = fTimeCoupling;
243 ((AliTRDparameter &) p).fTimeBinWidth = fTimeBinWidth;
244 ((AliTRDparameter &) p).fField = fField;
245 ((AliTRDparameter &) p).fPRFOn = fPRFOn;
246 ((AliTRDparameter &) p).fTRFOn = fTRFOn;
247 ((AliTRDparameter &) p).fCTOn = fCTOn;
248 ((AliTRDparameter &) p).fTCOn = fTCOn;
249 ((AliTRDparameter &) p).fTiltingAngle = fTiltingAngle;
250 ((AliTRDparameter &) p).fPRFbin = fPRFbin;
251 ((AliTRDparameter &) p).fPRFlo = fPRFlo;
252 ((AliTRDparameter &) p).fPRFhi = fPRFhi;
253 ((AliTRDparameter &) p).fPRFwid = fPRFwid;
254 ((AliTRDparameter &) p).fPRFpad = fPRFpad;
255 if (((AliTRDparameter &) p).fPRFsmp) delete [] ((AliTRDparameter &) p).fPRFsmp;
256 ((AliTRDparameter &) p).fPRFsmp = new Float_t[fPRFbin];
257 for (iBin = 0; iBin < fPRFbin; iBin++) {
258 ((AliTRDparameter &) p).fPRFsmp[iBin] = fPRFsmp[iBin];
259 }
260 ((AliTRDparameter &) p).fTRFbin = fTRFbin;
261 ((AliTRDparameter &) p).fTRFlo = fTRFlo;
262 ((AliTRDparameter &) p).fTRFhi = fTRFhi;
263 ((AliTRDparameter &) p).fTRFwid = fTRFwid;
264 if (((AliTRDparameter &) p).fTRFsmp) delete [] ((AliTRDparameter &) p).fTRFsmp;
265 ((AliTRDparameter &) p).fTRFsmp = new Float_t[fTRFbin];
266 for (iBin = 0; iBin < fTRFbin; iBin++) {
267 ((AliTRDparameter &) p).fTRFsmp[iBin] = fTRFsmp[iBin];
268 }
269 if (((AliTRDparameter &) p).fCTsmp) delete [] ((AliTRDparameter &) p).fCTsmp;
270 ((AliTRDparameter &) p).fCTsmp = new Float_t[fTRFbin];
271 for (iBin = 0; iBin < fTRFbin; iBin++) {
272 ((AliTRDparameter &) p).fCTsmp[iBin] = fCTsmp[iBin];
273 }
274 ((AliTRDparameter &) p).fTCnexp = fTCnexp;
275
276 ((AliTRDparameter &) p).fLUTOn = fLUTOn;
277 ((AliTRDparameter &) p).fLUTbin = fLUTbin;
278 if (((AliTRDparameter &) p).fLUT) delete [] ((AliTRDparameter &) p).fLUT;
279 ((AliTRDparameter &) p).fLUT = new Float_t[fLUTbin];
280 for (iBin = 0; iBin < fLUTbin; iBin++) {
281 ((AliTRDparameter &) p).fLUT[iBin] = fLUT[iBin];
282 }
283 ((AliTRDparameter &) p).fClusMaxThresh = fClusMaxThresh;
284 ((AliTRDparameter &) p).fClusSigThresh = fClusSigThresh;
285
a328fff9 286 ((AliTRDparameter &) p).fAnodeWireOffset = fAnodeWireOffset;
287 ((AliTRDparameter &) p).fStaggeringOn = fStaggeringOn;
288 ((AliTRDparameter &) p).fTimeStructOn = fTimeStructOn;
289 if (((AliTRDparameter &) p).fTimeStruct)
290 delete [] ((AliTRDparameter &) p).fTimeStruct;
291 ((AliTRDparameter &) p).fTimeStruct = new Float_t[38*26];
292 for (Int_t i = 0; i < 38*26; i++) {
293 ((AliTRDparameter &) p).fTimeStruct[i] = fTimeStruct[i];
294 }
295
17b26de4 296}
297
298//_____________________________________________________________________________
299void AliTRDparameter::Init()
300{
301 //
302 // Initializes the parameter
303 //
304 // The maximum number of pads
305 // and the position of pad 0,0,0
306 //
307 // chambers seen from the top:
308 // +----------------------------+
309 // | |
310 // | | ^
311 // | | rphi|
312 // | | |
313 // |0 | |
314 // +----------------------------+ +------>
315 // z
316 // chambers seen from the side: ^
317 // +----------------------------+ drift|
318 // |0 | |
319 // | | |
320 // +----------------------------+ +------>
321 // z
322 //
323 // IMPORTANT: time bin 0 is now the first one in the drift region
324 // closest to the readout !!!
325 //
326
327 //
328 // ----------------------------------------------------------------------------
329 // The pad definition
330 // ----------------------------------------------------------------------------
331 //
332
333 // The pad size in column direction (rphi-direction)
a328fff9 334 //SetColPadSize(0,0.65);
335 //SetColPadSize(1,0.68);
336 //SetColPadSize(2,0.71);
337 //SetColPadSize(3,0.74);
338 //SetColPadSize(4,0.77);
339 //SetColPadSize(5,0.80);
340
341 SetColPadSize(0,0.664);
342 SetColPadSize(1,0.695);
343 SetColPadSize(2,0.726);
344 SetColPadSize(3,0.756);
345 SetColPadSize(4,0.788);
346 SetColPadSize(5,0.818);
73ae7b59 347
17b26de4 348 // The pad row (z-direction)
349 SetNRowPad();
350
351 // The number of time bins. Default is 100 ns timbin size
a328fff9 352 SetNTimeBin(20);
17b26de4 353
354 // Additional time bins before and after the drift region.
355 // Default is to only sample the drift region
356 SetExpandTimeBin(0,0);
357
358 //
359 // ----------------------------------------------------------------------------
360 // The digitization parameter
361 // ----------------------------------------------------------------------------
362 //
363
364 // The default parameter for the digitization
a328fff9 365 fGasGain = 4000.;
17b26de4 366 fChipGain = 12.4;
367 fNoise = 1000.;
368 fADCoutRange = 1023.; // 10-bit ADC
7e582e9f 369 fADCinRange = 2000.; // 2V input range
17b26de4 370 fADCthreshold = 1;
7e582e9f 371 fADCbaseline = 0;
17b26de4 372
373 // The drift velocity (cm / mus)
374 fDriftVelocity = 1.5;
375
376 // Diffusion on
377 fDiffusionOn = 1;
378
379 // E x B effects
7e582e9f 380 fExBOn = 1;
17b26de4 381
382 // Propability for electron attachment
383 fElAttachOn = 0;
384 fElAttachProp = 0.0;
385
386 // The pad response function
387 fPRFOn = 1;
388
389 // The time response function
390 fTRFOn = 1;
391
392 // The cross talk
7e582e9f 393 fCTOn = 1;
17b26de4 394
395 // The tail cancelation
396 fTCOn = 1;
397
398 // The number of exponentials
399 fTCnexp = 2;
400
b258a7c7 401 // The pad coupling factor
402 //fPadCoupling = 0.3;
403 // Use 0.49 instead which reproduces better the test beam
404 // data, even tough it is not understood why.
405 fPadCoupling = 0.49;
17b26de4 406
407 // The time coupling factor (same number as for the TPC)
408 fTimeCoupling = 0.4;
409
a328fff9 410 // Drift time non-isochronity on
411 fTimeStructOn = 1;
412
413 // Wire planes staggered
414 fStaggeringOn = 1;
415
416 // Distance of first Anode wire from first pad edge
417 fAnodeWireOffset = 0.25;
418
17b26de4 419 // The tilting angle for the readout pads
73ae7b59 420 SetTiltingAngle(2.0);
17b26de4 421
422 // The magnetic field strength in Tesla
c57e2264 423 fField = 0.4;
17b26de4 424
425 //
426 // ----------------------------------------------------------------------------
427 // The clusterization parameter
428 // ----------------------------------------------------------------------------
429 //
430
431 // The default parameter for the clustering
432 fClusMaxThresh = 3;
433 fClusSigThresh = 1;
434
435 // Use the LUT
436 fLUTOn = 1;
437
438 ReInit();
439
440}
441
442//_____________________________________________________________________________
443void AliTRDparameter::ReInit()
444{
445 //
446 // Reinitializes the parameter class after a change
447 //
448
a328fff9 449 // Create the sampled timing structure
450 SampleTimeStruct();
451
17b26de4 452 // Calculate the time bin width in ns
453 fTimeBinWidth = fTimeBinSize / fDriftVelocity * 1000.0;
454
455 // The range and the binwidth for the sampled TRF
456 fTRFbin = 100;
457 // Start 0.2 mus before the signal
458 fTRFlo = -0.2 * fDriftVelocity;
459 // End the maximum driftlength after the signal
460 fTRFhi = AliTRDgeometry::DrThick()
461 + fTimeAfter * fTimeBinSize;
462 fTRFwid = (fTRFhi - fTRFlo) / ((Float_t) fTRFbin);
463
464 // Transverse and longitudinal diffusion coefficients (Xe/CO2)
465 fDiffusionT = GetDiffusionT(fDriftVelocity,fField);
466 fDiffusionL = GetDiffusionL(fDriftVelocity,fField);
467
468 // omega * tau.= tan(Lorentz-angle)
469 fOmegaTau = GetOmegaTau(fDriftVelocity,fField);
470
471 // The Lorentz factor
472 if (fExBOn) {
473 fLorentzFactor = 1.0 / (1.0 + fOmegaTau*fOmegaTau);
474 }
475 else {
476 fLorentzFactor = 1.0;
477 }
478
479 // Create the sampled PRF
480 SamplePRF();
481
482 // Create the sampled TRF
483 SampleTRF();
484
485 // Create the LUT
486 FillLUT();
487
488}
489
490//_____________________________________________________________________________
afc51ac2 491void AliTRDparameter::SetNRowPad(Int_t p, Int_t c, Int_t npad)
17b26de4 492{
493 //
494 // Redefines the number of pads in raw direction for
495 // a given plane and chamber number
496 //
497
498 for (Int_t isect = 0; isect < AliTRDgeometry::Nsect(); isect++) {
499
500 fRowMax[p][c][isect] = npad;
501
502 fRowPadSize[p][c][isect] = (fGeo->GetChamberLength(p,c)
503 - 2.* AliTRDgeometry::RpadW())
504 / ((Float_t) npad);
505
506 }
507
508}
509
510//_____________________________________________________________________________
511void AliTRDparameter::SetNRowPad()
512{
513 //
514 // Defines the number of pads in row direction
515 //
516
517 Int_t isect;
518 Int_t icham;
519 Int_t iplan;
520
521 Int_t rowMax[kNplan][kNcham] = { { 16, 16, 12, 16, 16 }
522 , { 16, 16, 12, 16, 16 }
523 , { 16, 16, 12, 16, 16 }
524 , { 16, 16, 12, 16, 16 }
ad973dd4 525 , { 16, 16, 12, 16, 16 }
526 , { 16, 16, 12, 16, 16 } };
17b26de4 527
528 Float_t rpadW = AliTRDgeometry::RpadW();
529
530 for (isect = 0; isect < kNsect; isect++) {
531 for (icham = 0; icham < kNcham; icham++) {
532 for (iplan = 0; iplan < kNplan; iplan++) {
533
534 fRowMax[iplan][icham][isect] = rowMax[iplan][icham];
535
536 fRowPadSize[iplan][icham][isect] = (fGeo->GetChamberLength(iplan,icham)
537 - 2.*rpadW)
538 / ((Float_t) rowMax[iplan][icham]);
539
540 Float_t row0 = rpadW - fGeo->GetChamberLength(iplan,0)
541 - fGeo->GetChamberLength(iplan,1)
542 - fGeo->GetChamberLength(iplan,2) / 2.;
543 for (Int_t ic = 0; ic < icham; ic++) {
4bd34b78 544 row0 += fGeo->GetChamberLength(iplan,ic);
17b26de4 545 }
546 fRow0[iplan][icham][isect] = row0;
547
548 }
549 }
550 }
551
552}
553
554//_____________________________________________________________________________
afc51ac2 555void AliTRDparameter::SetColPadSize(Int_t p, Float_t s)
17b26de4 556{
557 //
558 // Redefines the pad size in column direction
559 //
560
561 Float_t cpadW = AliTRDgeometry::CpadW();
562
563 fColPadSize[p] = s;
564 fCol0[p] = - fGeo->GetChamberWidth(p)/2. + cpadW;
565 fColMax[p] = ((Int_t) ((fGeo->GetChamberWidth(p) - 2.*cpadW) / s));
566
567}
568
569//_____________________________________________________________________________
afc51ac2 570void AliTRDparameter::SetNTimeBin(Int_t nbin)
17b26de4 571{
572 //
573 // Redefines the number of time bins in the drift region.
574 // The time bin width is defined by the length of the
575 // drift region divided by <nbin>.
576 //
577
578 fTimeMax = nbin;
579 fTimeBinSize = AliTRDgeometry::DrThick() / ((Float_t) fTimeMax);
580 for (Int_t iplan = 0; iplan < AliTRDgeometry::Nplan(); iplan++) {
581 fTime0[iplan] = AliTRDgeometry::Rmin()
582 + AliTRDgeometry::CraHght()
583 + AliTRDgeometry::CdrHght()
584 + iplan * (AliTRDgeometry::Cheight()
585 + AliTRDgeometry::Cspace());
586 }
587
588}
589
590//_____________________________________________________________________________
0a29d0f1 591Float_t AliTRDparameter::CrossTalk(Float_t time) const
17b26de4 592{
593 //
594 // Applies the pad-pad capacitive cross talk
595 //
596
597 Int_t iBin = ((Int_t) ((time - fTRFlo) / fTRFwid));
598 if ((iBin >= 0) && (iBin < fTRFbin)) {
599 return fCTsmp[iBin];
600 }
601 else {
602 return 0.0;
603 }
604
605}
606
a328fff9 607//_____________________________________________________________________________
608Int_t AliTRDparameter::TimeStruct(Float_t time, Float_t z, Float_t *xyz) const
609{
610 //
611 // Applies the time structure of the drift cells (by C.Lippmann).
612 // The drift time of electrons to the anode wires depends on the
613 // distance to the wire (z) and on the position in the drift region.
614 //
615 // input : time (radial distance from pad cathode) in [cm] and z (distance
616 // from anode wire parallel to cathode planes) in [cm]
617 //
618
619 if ( z > 0.24 ) z = 0.24;
620 if ( !fStaggeringOn && z < 0.01 ) z = 0.01;
621
622 Int_t r1 = (Int_t)(10*time);
623 Int_t r2 = r1+1;
624 Int_t z1 = (Int_t)(100*z);
625 Int_t z2 = z1+1;
626
627 if (r1<0 || r1>37 || z1<0 || z1>25) {
628 printf("<AliTRDparameter::TimeStruct> Warning. Indices out of range: ");
629 printf("time=%.2f, z=%.2f, r1=%d, z1=%d\n",time,z,r1,z1);
630 return kFALSE;
631 if (r1<0) r1 = 0;
632 if (r1>37) r1 = 37;
633 if (z1<0) z1 = 0;
634 if (z1>25) z1 = 25;
635 }
636
637 // 2D Interpolation:
638 Float_t y11, y12, y21, y22, y1, y2;
639 y11 = fTimeStruct[r1+38*z1];
640 y22 = (r2 <= 37 && z2 <= 25) ? fTimeStruct[r2+38*z2] : fTimeStruct[37+38*25];
641 y12 = (z2 <= 25) ? fTimeStruct[r1+38*z2] : fTimeStruct[r1+38*25];
642 y21 = (r2 <= 37) ? fTimeStruct[r2+38*z1] : fTimeStruct[37+38*z1];
643
644 y1 = (y21-y11)*10*time + y11 - (y21-y11)*r1;
645 y2 = (y22-y12)*10*time + y12 - (y22-y12)*r1;
646
647 Float_t AmTh = AliTRDgeometry::AmThick() / 2.;
648 Float_t tdrift = 0.0; // drift time (GARFIELD)
649 if (TMath::Abs(time-AmTh)>0.005 || z>0.005)
650 tdrift = (y2-y1)*100*z + y1 - (y2-y1)*z1;
651 if (time < AmTh) tdrift *= -1.;
652 Float_t xdrift = tdrift*(AmTh + AliTRDgeometry::DrThick())/fTimeStruct[37] + AmTh;
653 Float_t offset = time - xdrift;
654
655 xyz[0] = xyz[0] + offset;
656 xyz[1] = xyz[1];
657 xyz[2] = xyz[2];
658
659 return 1;
660
661}
662
17b26de4 663//_____________________________________________________________________________
664Int_t AliTRDparameter::Diffusion(Float_t driftlength, Float_t *xyz)
665{
666 //
667 // Applies the diffusion smearing to the position of a single electron
668 //
669
670 Float_t driftSqrt = TMath::Sqrt(driftlength);
671 Float_t sigmaT = driftSqrt * fDiffusionT;
672 Float_t sigmaL = driftSqrt * fDiffusionL;
673 xyz[0] = gRandom->Gaus(xyz[0], sigmaL * fLorentzFactor);
674 xyz[1] = gRandom->Gaus(xyz[1], sigmaT * fLorentzFactor);
675 xyz[2] = gRandom->Gaus(xyz[2], sigmaT);
676
677 return 1;
678
679}
680
681//_____________________________________________________________________________
0a29d0f1 682Int_t AliTRDparameter::ExB(Float_t driftlength, Float_t *xyz) const
17b26de4 683{
684 //
685 // Applies E x B effects to the position of a single electron
686 //
687
688 xyz[0] = xyz[0];
689 xyz[1] = xyz[1] + fOmegaTau * driftlength;
690 xyz[2] = xyz[2];
691
692 return 1;
693
694}
695
696//_____________________________________________________________________________
697Int_t AliTRDparameter::PadResponse(Float_t signal, Float_t dist
0a29d0f1 698 , Int_t plane, Float_t *pad) const
17b26de4 699{
700 //
701 // Applies the pad response
702 //
703
704 const Int_t kNplan = AliTRDgeometry::kNplan;
705
706 Int_t iBin = ((Int_t) (( - dist - fPRFlo) / fPRFwid));
707 Int_t iOff = plane * fPRFbin;
708
709 Int_t iBin0 = iBin - fPRFpad + iOff;
710 Int_t iBin1 = iBin + iOff;
711 Int_t iBin2 = iBin + fPRFpad + iOff;
712
713 pad[0] = 0.0;
714 pad[1] = 0.0;
715 pad[2] = 0.0;
716 if ((iBin1 >= 0) && (iBin1 < (fPRFbin*kNplan))) {
717
718 if (iBin0 >= 0) {
719 pad[0] = signal * fPRFsmp[iBin0];
720 }
721 pad[1] = signal * fPRFsmp[iBin1];
722 if (iBin2 < (fPRFbin*kNplan)) {
723 pad[2] = signal * fPRFsmp[iBin2];
724 }
725
726 return 1;
727
728 }
729 else {
730
731 return 0;
732
733 }
734
735}
736
737//_____________________________________________________________________________
0a29d0f1 738Float_t AliTRDparameter::TimeResponse(Float_t time) const
17b26de4 739{
740 //
741 // Applies the preamp shaper time response
742 //
743
744 Int_t iBin = ((Int_t) ((time - fTRFlo) / fTRFwid));
745 if ((iBin >= 0) && (iBin < fTRFbin)) {
746 return fTRFsmp[iBin];
747 }
748 else {
749 return 0.0;
750 }
751
752}
753
754//_____________________________________________________________________________
755Float_t AliTRDparameter::Col0Tilted(Float_t col0, Float_t rowOffset
756 , Int_t plane)
757{
758 //
759 // Calculates col0 for tilted pads
760 //
761
762 Float_t diff = fTiltingAngle * rowOffset;
a328fff9 763 return (col0 + TMath::Power(-1.0,(plane+1)) * diff);
17b26de4 764
765}
766
767//_____________________________________________________________________________
768void AliTRDparameter::SampleTRF()
769{
770 //
771 // Samples the time response function
772 //
773 // New TRF from Venelin Angelov, simulated with CADENCE
774 // Pad-ground capacitance = 25 pF
775 // Pad-pad cross talk capacitance = 6 pF
776 //
777
778 Int_t ipos1;
779 Int_t ipos2;
780 Float_t diff;
781
782 const Int_t kNpasa = 252;
783
784 Float_t time[kNpasa] = { -0.220000, -0.210000, -0.200000, -0.190000
785 , -0.180000, -0.170000, -0.160000, -0.150000
786 , -0.140000, -0.130000, -0.120000, -0.110000
787 , -0.100000, -0.090000, -0.080000, -0.070000
788 , -0.060000, -0.050000, -0.040000, -0.030000
789 , -0.020000, -0.010000, -0.000000, 0.010000
790 , 0.020000, 0.030000, 0.040000, 0.050000
791 , 0.060000, 0.070000, 0.080000, 0.090000
792 , 0.100000, 0.110000, 0.120000, 0.130000
793 , 0.140000, 0.150000, 0.160000, 0.170000
794 , 0.180000, 0.190000, 0.200000, 0.210000
795 , 0.220000, 0.230000, 0.240000, 0.250000
796 , 0.260000, 0.270000, 0.280000, 0.290000
797 , 0.300000, 0.310000, 0.320000, 0.330000
798 , 0.340000, 0.350000, 0.360000, 0.370000
799 , 0.380000, 0.390000, 0.400000, 0.410000
800 , 0.420000, 0.430000, 0.440000, 0.450000
801 , 0.460000, 0.470000, 0.480000, 0.490000
802 , 0.500000, 0.510000, 0.520000, 0.530000
803 , 0.540000, 0.550000, 0.560000, 0.570000
804 , 0.580000, 0.590000, 0.600000, 0.610000
805 , 0.620000, 0.630000, 0.640000, 0.650000
806 , 0.660000, 0.670000, 0.680000, 0.690000
807 , 0.700000, 0.710000, 0.720000, 0.730000
808 , 0.740000, 0.750000, 0.760000, 0.770000
809 , 0.780000, 0.790000, 0.800000, 0.810000
810 , 0.820000, 0.830000, 0.840000, 0.850000
811 , 0.860000, 0.870000, 0.880000, 0.890000
812 , 0.900000, 0.910000, 0.920000, 0.930000
813 , 0.940000, 0.950000, 0.960000, 0.970000
814 , 0.980000, 0.990000, 1.000000, 1.010000
815 , 1.020000, 1.030000, 1.040000, 1.050000
816 , 1.060000, 1.070000, 1.080000, 1.090000
817 , 1.100000, 1.110000, 1.120000, 1.130000
818 , 1.140000, 1.150000, 1.160000, 1.170000
819 , 1.180000, 1.190000, 1.200000, 1.210000
820 , 1.220000, 1.230000, 1.240000, 1.250000
821 , 1.260000, 1.270000, 1.280000, 1.290000
822 , 1.300000, 1.310000, 1.320000, 1.330000
823 , 1.340000, 1.350000, 1.360000, 1.370000
824 , 1.380000, 1.390000, 1.400000, 1.410000
825 , 1.420000, 1.430000, 1.440000, 1.450000
826 , 1.460000, 1.470000, 1.480000, 1.490000
827 , 1.500000, 1.510000, 1.520000, 1.530000
828 , 1.540000, 1.550000, 1.560000, 1.570000
829 , 1.580000, 1.590000, 1.600000, 1.610000
830 , 1.620000, 1.630000, 1.640000, 1.650000
831 , 1.660000, 1.670000, 1.680000, 1.690000
832 , 1.700000, 1.710000, 1.720000, 1.730000
833 , 1.740000, 1.750000, 1.760000, 1.770000
834 , 1.780000, 1.790000, 1.800000, 1.810000
835 , 1.820000, 1.830000, 1.840000, 1.850000
836 , 1.860000, 1.870000, 1.880000, 1.890000
837 , 1.900000, 1.910000, 1.920000, 1.930000
838 , 1.940000, 1.950000, 1.960000, 1.970000
839 , 1.980000, 1.990000, 2.000000, 2.010000
840 , 2.020000, 2.030000, 2.040000, 2.050000
841 , 2.060000, 2.070000, 2.080000, 2.090000
842 , 2.100000, 2.110000, 2.120000, 2.130000
843 , 2.140000, 2.150000, 2.160000, 2.170000
844 , 2.180000, 2.190000, 2.200000, 2.210000
845 , 2.220000, 2.230000, 2.240000, 2.250000
846 , 2.260000, 2.270000, 2.280000, 2.290000 };
847
848 Float_t signal[kNpasa] = { 0.000000, 0.000000, 0.000000, 0.000000
849 , 0.000000, 0.000000, 0.000000, 0.000396
850 , 0.005096, 0.022877, 0.061891, 0.126614
851 , 0.215798, 0.324406, 0.444507, 0.566817
852 , 0.683465, 0.787089, 0.873159, 0.937146
853 , 0.979049, 0.999434, 1.000000, 0.983579
854 , 0.954134, 0.913364, 0.866365, 0.813703
855 , 0.759910, 0.706116, 0.653454, 0.603624
856 , 0.556625, 0.514156, 0.475085, 0.439977
857 , 0.408834, 0.380578, 0.355549, 0.333352
858 , 0.313647, 0.296093, 0.280351, 0.266195
859 , 0.253397, 0.241789, 0.231257, 0.221574
860 , 0.212627, 0.204417, 0.196772, 0.189581
861 , 0.182956, 0.176784, 0.171008, 0.165515
862 , 0.160419, 0.155606, 0.151076, 0.146716
863 , 0.142639, 0.138845, 0.135221, 0.131767
864 , 0.128482, 0.125368, 0.122424, 0.119592
865 , 0.116931, 0.114326, 0.111891, 0.109513
866 , 0.107248, 0.105096, 0.103058, 0.101019
867 , 0.099151, 0.097282, 0.095527, 0.093715
868 , 0.092129, 0.090544, 0.088958, 0.087429
869 , 0.086014, 0.084598, 0.083239, 0.081880
870 , 0.080634, 0.079388, 0.078143, 0.077010
871 , 0.075878, 0.074745, 0.073669, 0.072593
872 , 0.071574, 0.070612, 0.069649, 0.068686
873 , 0.067780, 0.066874, 0.066025, 0.065176
874 , 0.064326, 0.063533, 0.062684, 0.061948
875 , 0.061212, 0.060419, 0.059740, 0.059003
876 , 0.058324, 0.057644, 0.057022, 0.056342
877 , 0.055663, 0.055096, 0.054473, 0.053851
878 , 0.053284, 0.052718, 0.052152, 0.051585
879 , 0.051019, 0.050566, 0.050000, 0.049490
880 , 0.048981, 0.048528, 0.048018, 0.047508
881 , 0.047055, 0.046602, 0.046149, 0.045696
882 , 0.045300, 0.044904, 0.044451, 0.044054
883 , 0.043658, 0.043205, 0.042865, 0.042469
884 , 0.042072, 0.041733, 0.041336, 0.040997
885 , 0.040657, 0.040260, 0.039921, 0.039581
886 , 0.039241, 0.038958, 0.038618, 0.038335
887 , 0.037995, 0.037656, 0.037373, 0.037089
888 , 0.036806, 0.036467, 0.036183, 0.035900
889 , 0.035617, 0.035334, 0.035108, 0.034824
890 , 0.034541, 0.034315, 0.034032, 0.033805
891 , 0.033522, 0.033296, 0.033069, 0.032786
892 , 0.032559, 0.032333, 0.032106, 0.031880
893 , 0.031653, 0.031427, 0.031200, 0.030974
894 , 0.030804, 0.030578, 0.030351, 0.030125
895 , 0.029955, 0.029785, 0.029558, 0.029332
896 , 0.029162, 0.028992, 0.028766, 0.028596
897 , 0.028426, 0.028199, 0.028086, 0.027860
898 , 0.027746, 0.027633, 0.027463, 0.027293
899 , 0.027180, 0.027067, 0.026954, 0.026954
900 , 0.026840, 0.026727, 0.026727, 0.026614
901 , 0.026614, 0.026614, 0.026557, 0.026501
902 , 0.026501, 0.026501, 0.026501, 0.026501
903 , 0.026501, 0.026501, 0.026501, 0.026387
904 , 0.026387, 0.026387, 0.026387, 0.026387
905 , 0.026387, 0.026387, 0.026387, 0.026387
906 , 0.026387, 0.026387, 0.026387, 0.026387
907 , 0.026387, 0.026274, 0.026274, 0.026274
908 , 0.026274, 0.026274, 0.026274, 0.026274
909 , 0.026274, 0.026274, 0.026274, 0.026274
910 , 0.026274, 0.026274, 0.026274, 0.026161 };
911
912 Float_t xtalk[kNpasa] = { 0.000000, 0.000000, 0.000000, 0.000000
913 , 0.000000, 0.000000, 0.000000, 0.000113
914 , 0.000793, 0.003058, 0.007305, 0.013194
915 , 0.019706, 0.025821, 0.030634, 0.033465
916 , 0.034145, 0.032729, 0.029615, 0.025198
917 , 0.019989, 0.014496, 0.009003, 0.003964
918 , -0.000510, -0.004190, -0.007191, -0.009400
919 , -0.010872, -0.011835, -0.012288, -0.012288
920 , -0.012005, -0.011495, -0.010872, -0.010136
921 , -0.009343, -0.008607, -0.007871, -0.007191
922 , -0.006512, -0.005946, -0.005379, -0.004926
923 , -0.004473, -0.004077, -0.003737, -0.003398
924 , -0.003114, -0.002831, -0.002605, -0.002378
925 , -0.002208, -0.002039, -0.001869, -0.001699
926 , -0.001585, -0.001472, -0.001359, -0.001246
927 , -0.001132, -0.001019, -0.001019, -0.000906
928 , -0.000906, -0.000793, -0.000793, -0.000680
929 , -0.000680, -0.000680, -0.000566, -0.000566
930 , -0.000566, -0.000566, -0.000453, -0.000453
931 , -0.000453, -0.000453, -0.000453, -0.000453
932 , -0.000340, -0.000340, -0.000340, -0.000340
933 , -0.000340, -0.000340, -0.000340, -0.000340
934 , -0.000340, -0.000340, -0.000340, -0.000340
935 , -0.000340, -0.000227, -0.000227, -0.000227
936 , -0.000227, -0.000227, -0.000227, -0.000227
937 , -0.000227, -0.000227, -0.000227, -0.000227
938 , -0.000227, -0.000227, -0.000227, -0.000227
939 , -0.000227, -0.000227, -0.000227, -0.000227
940 , -0.000227, -0.000227, -0.000227, -0.000227
941 , -0.000227, -0.000227, -0.000227, -0.000227
942 , -0.000227, -0.000227, -0.000227, -0.000227
943 , -0.000227, -0.000227, -0.000227, -0.000227
944 , -0.000227, -0.000227, -0.000227, -0.000113
945 , -0.000113, -0.000113, -0.000113, -0.000113
946 , -0.000113, -0.000113, -0.000113, -0.000113
947 , -0.000113, -0.000113, -0.000113, -0.000113
948 , -0.000113, -0.000113, -0.000113, -0.000113
949 , -0.000113, -0.000113, -0.000113, -0.000113
950 , -0.000113, -0.000113, -0.000113, -0.000113
951 , -0.000113, -0.000113, -0.000113, -0.000113
952 , -0.000113, -0.000113, -0.000113, -0.000113
953 , -0.000113, -0.000113, -0.000113, -0.000113
954 , -0.000113, -0.000113, -0.000113, -0.000113
955 , -0.000113, -0.000113, -0.000113, -0.000113
956 , -0.000113, -0.000113, -0.000113, -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.000000
965 , 0.000000, 0.000000, 0.000000, 0.000000
966 , 0.000000, 0.000000, 0.000000, 0.000000
967 , 0.000000, 0.000000, 0.000000, 0.000000
968 , 0.000000, 0.000000, 0.000000, 0.000000
969 , 0.000000, 0.000000, 0.000000, 0.000000
970 , 0.000000, 0.000000, 0.000000, 0.000000
971 , 0.000000, 0.000000, 0.000000, 0.000000
972 , 0.000000, 0.000000, 0.000000, 0.000000
973 , 0.000000, 0.000000, 0.000000, 0.000000
974 , 0.000000, 0.000000, 0.000000, 0.000000 };
975
976 // increase CrossTalk to measurements
977 for (Int_t ipasa = 0; ipasa < kNpasa; ipasa++) {
978 xtalk[ipasa] *= 1.75;
979 }
980
981 if (fTRFsmp) delete [] fTRFsmp;
982 fTRFsmp = new Float_t[fTRFbin];
983 if (fCTsmp) delete [] fCTsmp;
984 fCTsmp = new Float_t[fTRFbin];
985
986 Float_t loTRF = TMath::Max(fTRFlo / fDriftVelocity,time[0]);
987 Float_t hiTRF = TMath::Min(fTRFhi / fDriftVelocity,time[kNpasa-1]);
988 Float_t binWidth = (hiTRF - loTRF) / ((Float_t) fTRFbin);
989
990 // Take the linear interpolation
991 for (Int_t iBin = 0; iBin < fTRFbin; iBin++) {
992
993 Float_t bin = (((Float_t) iBin) + 0.5) * binWidth + loTRF;
994 ipos1 = ipos2 = 0;
995 diff = 0;
996 do {
997 diff = bin - time[ipos2++];
998 } while (diff > 0);
999 ipos2--;
1000 if (ipos2 >= kNpasa) ipos2 = kNpasa - 1;
1001 ipos1 = ipos2 - 1;
1002
1003 fTRFsmp[iBin] = signal[ipos2]
1004 + diff * (signal[ipos2] - signal[ipos1])
1005 / ( time[ipos2] - time[ipos1]);
1006
1007 fCTsmp[iBin] = xtalk[ipos2]
1008 + diff * (xtalk[ipos2] - xtalk[ipos1])
1009 / ( time[ipos2] - time[ipos1]);
1010
1011 }
1012
1013}
1014
a328fff9 1015//_____________________________________________________________________________
1016void AliTRDparameter::SampleTimeStruct()
1017{
1018 //
1019 // Samples the timing structure of a drift cell
1020 // Drift Time data calculated with Garfield (by C.Lippmann)
1021 //
1022
1023 const Int_t ktimebin = 38;
1024 const Int_t kZbin = 26;
1025
1026 Float_t stagg[ktimebin][kZbin] =
1027 {{9.1221e-02, 9.1269e-02, 9.1426e-02, 9.1678e-02, 9.2041e-02, 9.2500e-02,
1028 9.3081e-02, 9.3760e-02, 9.4540e-02, 9.5434e-02, 9.6436e-02, 9.7554e-02,
1029 9.8796e-02, 1.0016e-01, 1.0167e-01, 1.0336e-01, 1.0528e-01, 1.0748e-01,
1030 1.1012e-01, 1.1348e-01, 1.1790e-01, 1.2395e-01, 1.3317e-01, 1.5036e-01,
1031 1.9348e-01, 5.7191e-01},
1032 {6.5708e-02, 6.5787e-02, 6.5954e-02, 6.6237e-02, 6.6646e-02, 6.7170e-02,
1033 6.7804e-02, 6.8554e-02, 6.9424e-02, 7.0392e-02, 7.1492e-02, 7.2693e-02,
1034 7.4006e-02, 7.5459e-02, 7.7042e-02, 7.8795e-02, 8.0786e-02, 8.3081e-02,
1035 8.5798e-02, 8.9290e-02, 9.3841e-02, 1.0007e-01, 1.0959e-01, 1.2741e-01,
1036 1.7152e-01, 5.4720e-01},
1037 {3.9602e-02, 3.9699e-02, 3.9969e-02, 4.0417e-02, 4.1015e-02, 4.1802e-02,
1038 4.2726e-02, 4.3804e-02, 4.4992e-02, 4.6309e-02, 4.7720e-02, 4.9244e-02,
1039 5.0874e-02, 5.2592e-02, 5.4437e-02, 5.6461e-02, 5.8704e-02, 6.1340e-02,
1040 6.4533e-02, 6.8509e-02, 7.3620e-02, 8.0737e-02, 9.1920e-02, 1.1276e-01,
1041 1.6078e-01, 5.2455e-01},
1042 {1.1236e-02, 1.1514e-02, 1.2355e-02, 1.3666e-02, 1.5354e-02, 1.7292e-02,
1043 1.9469e-02, 2.1756e-02, 2.4118e-02, 2.6505e-02, 2.8902e-02, 3.1273e-02,
1044 3.3618e-02, 3.5960e-02, 3.8367e-02, 4.0957e-02, 4.3859e-02, 4.7286e-02,
1045 5.1516e-02, 5.7160e-02, 6.5086e-02, 7.6740e-02, 9.4905e-02, 1.2426e-01,
1046 1.8267e-01, 4.8444e-01},
1047 {1.1285e-02, 1.1581e-02, 1.2430e-02, 1.3739e-02, 1.5383e-02, 1.7329e-02,
1048 1.9477e-02, 2.1732e-02, 2.4063e-02, 2.6412e-02, 2.8746e-02, 3.1073e-02,
1049 3.3329e-02, 3.5733e-02, 3.8146e-02, 4.0847e-02, 4.3908e-02, 4.7650e-02,
1050 5.2475e-02, 5.9069e-02, 6.8420e-02, 8.2181e-02, 1.0199e-01, 1.3525e-01,
1051 1.9150e-01, 4.6381e-01},
1052 {3.9345e-02, 3.9430e-02, 3.9698e-02, 4.0124e-02, 4.0714e-02, 4.1462e-02,
1053 4.2325e-02, 4.3327e-02, 4.4476e-02, 4.5727e-02, 4.7075e-02, 4.8532e-02,
1054 5.0118e-02, 5.1838e-02, 5.3758e-02, 5.5958e-02, 5.8528e-02, 6.1589e-02,
1055 6.5285e-02, 6.9903e-02, 7.6022e-02, 8.4799e-02, 9.8707e-02, 1.2316e-01,
1056 1.7495e-01, 5.2371e-01},
1057 {6.4837e-02, 6.4889e-02, 6.5062e-02, 6.5329e-02, 6.5688e-02, 6.6204e-02,
1058 6.6770e-02, 6.7494e-02, 6.8306e-02, 6.9234e-02, 7.0297e-02, 7.1484e-02,
1059 7.2842e-02, 7.4318e-02, 7.6041e-02, 7.8040e-02, 8.0428e-02, 8.3315e-02,
1060 8.6827e-02, 9.1171e-02, 9.6759e-02, 1.0445e-01, 1.1651e-01, 1.3841e-01,
1061 1.8746e-01, 5.4644e-01},
1062 {9.0789e-02, 9.0824e-02, 9.0908e-02, 9.1065e-02, 9.1268e-02, 9.1522e-02,
1063 9.1838e-02, 9.2207e-02, 9.2643e-02, 9.3230e-02, 9.4073e-02, 9.5341e-02,
1064 9.7229e-02, 1.0633e-01, 1.0703e-01, 1.0868e-01, 1.1124e-01, 1.1451e-01,
1065 1.1853e-01, 1.2342e-01, 1.2981e-01, 1.3868e-01, 1.5245e-01, 1.7661e-01,
1066 2.2796e-01, 5.6931e-01},
1067 {1.2735e-01, 1.2744e-01, 1.2775e-01, 1.2828e-01, 1.2909e-01, 1.3021e-01,
1068 1.3172e-01, 1.3377e-01, 1.3658e-01, 1.4055e-01, 1.4624e-01, 1.5549e-01,
1069 1.7490e-01, 2.1347e-01, 1.8445e-01, 1.7581e-01, 1.7263e-01, 1.7216e-01,
1070 1.7377e-01, 1.7731e-01, 1.8337e-01, 1.9306e-01, 2.0875e-01, 2.3610e-01,
1071 2.9100e-01, 6.0131e-01},
1072 {1.7654e-01, 1.7663e-01, 1.7711e-01, 1.7789e-01, 1.7901e-01, 1.8055e-01,
1073 1.8252e-01, 1.8511e-01, 1.8844e-01, 1.9288e-01, 1.9889e-01, 2.0816e-01,
1074 2.2673e-01, 2.7253e-01, 2.3949e-01, 2.3043e-01, 2.2695e-01, 2.2620e-01,
1075 2.2751e-01, 2.3091e-01, 2.3691e-01, 2.4671e-01, 2.6258e-01, 2.9036e-01,
1076 3.4575e-01, 6.4208e-01},
1077 {2.2845e-01, 2.2860e-01, 2.2907e-01, 2.2987e-01, 2.3103e-01, 2.3262e-01,
1078 2.3466e-01, 2.3727e-01, 2.4063e-01, 2.4510e-01, 2.5116e-01, 2.6032e-01,
1079 2.7853e-01, 3.2713e-01, 2.9240e-01, 2.8314e-01, 2.7962e-01, 2.7878e-01,
1080 2.8008e-01, 2.8341e-01, 2.8944e-01, 2.9918e-01, 3.1503e-01, 3.4275e-01,
1081 3.9786e-01, 6.8482e-01},
1082 {2.8077e-01, 2.8094e-01, 2.8142e-01, 2.8222e-01, 2.8341e-01, 2.8496e-01,
1083 2.8702e-01, 2.8962e-01, 2.9301e-01, 2.9745e-01, 3.0350e-01, 3.1263e-01,
1084 3.3069e-01, 3.8049e-01, 3.4500e-01, 3.3562e-01, 3.3206e-01, 3.3126e-01,
1085 3.3253e-01, 3.3588e-01, 3.4183e-01, 3.5155e-01, 3.6729e-01, 3.9494e-01,
1086 4.4955e-01, 7.2458e-01},
1087 {3.3321e-01, 3.3336e-01, 3.3385e-01, 3.3466e-01, 3.3584e-01, 3.3744e-01,
1088 3.3946e-01, 3.4213e-01, 3.4544e-01, 3.4987e-01, 3.5592e-01, 3.6502e-01,
1089 3.8298e-01, 4.3361e-01, 3.9758e-01, 3.8814e-01, 3.8456e-01, 3.8373e-01,
1090 3.8498e-01, 3.8831e-01, 3.9424e-01, 4.0392e-01, 4.1961e-01, 4.4703e-01,
1091 5.0121e-01, 7.6594e-01},
1092 {3.8572e-01, 3.8588e-01, 3.8634e-01, 3.8716e-01, 3.8834e-01, 3.8993e-01,
1093 3.9196e-01, 3.9456e-01, 3.9792e-01, 4.0236e-01, 4.0848e-01, 4.1745e-01,
1094 4.3532e-01, 4.8671e-01, 4.5019e-01, 4.4067e-01, 4.3708e-01, 4.3625e-01,
1095 4.3746e-01, 4.4078e-01, 4.4663e-01, 4.5626e-01, 4.7192e-01, 4.9913e-01,
1096 5.5270e-01, 8.0632e-01},
1097 {4.3847e-01, 4.3855e-01, 4.3892e-01, 4.3976e-01, 4.4089e-01, 4.4247e-01,
1098 4.4451e-01, 4.4711e-01, 4.5048e-01, 4.5489e-01, 4.6092e-01, 4.6997e-01,
1099 4.8773e-01, 5.3996e-01, 5.0285e-01, 4.9328e-01, 4.8966e-01, 4.8878e-01,
1100 4.9006e-01, 4.9327e-01, 4.9912e-01, 5.0870e-01, 5.2424e-01, 5.5126e-01,
1101 6.0420e-01, 8.4589e-01},
1102 {4.9089e-01, 4.9105e-01, 4.9153e-01, 4.9234e-01, 4.9351e-01, 4.9509e-01,
1103 4.9712e-01, 4.9970e-01, 5.0306e-01, 5.0749e-01, 5.1348e-01, 5.2250e-01,
1104 5.4016e-01, 5.9329e-01, 5.5559e-01, 5.4594e-01, 5.4229e-01, 5.4138e-01,
1105 5.4257e-01, 5.4584e-01, 5.5162e-01, 5.6114e-01, 5.7660e-01, 6.0338e-01,
1106 6.5560e-01, 8.8395e-01},
1107 {5.4358e-01, 5.4372e-01, 5.4422e-01, 5.4500e-01, 5.4618e-01, 5.4776e-01,
1108 5.4978e-01, 5.5238e-01, 5.5572e-01, 5.6014e-01, 5.6613e-01, 5.7508e-01,
1109 5.9263e-01, 6.4680e-01, 6.0840e-01, 5.9867e-01, 5.9497e-01, 5.9416e-01,
1110 5.9522e-01, 5.9846e-01, 6.0424e-01, 6.1369e-01, 6.2902e-01, 6.5556e-01,
1111 7.0707e-01, 9.1967e-01},
1112 {5.9631e-01, 5.9648e-01, 5.9694e-01, 5.9776e-01, 5.9892e-01, 6.0049e-01,
1113 6.0252e-01, 6.0511e-01, 6.0845e-01, 6.1287e-01, 6.1881e-01, 6.2775e-01,
1114 6.4519e-01, 7.0051e-01, 6.6128e-01, 6.5147e-01, 6.4773e-01, 6.4680e-01,
1115 6.4793e-01, 6.5114e-01, 6.5690e-01, 6.6624e-01, 6.8147e-01, 7.0775e-01,
1116 7.5854e-01, 9.5826e-01},
1117 {6.4916e-01, 6.4928e-01, 6.4976e-01, 6.5057e-01, 6.5173e-01, 6.5329e-01,
1118 6.5533e-01, 6.5792e-01, 6.6124e-01, 6.6563e-01, 6.7159e-01, 6.8051e-01,
1119 6.9779e-01, 7.5441e-01, 7.1424e-01, 7.0434e-01, 7.0057e-01, 6.9962e-01,
1120 7.0071e-01, 7.0392e-01, 7.0959e-01, 7.1891e-01, 7.3396e-01, 7.6003e-01,
1121 8.1008e-01, 9.9861e-01},
1122 {7.0207e-01, 7.0221e-01, 7.0264e-01, 7.0345e-01, 7.0461e-01, 7.0621e-01,
1123 7.0820e-01, 7.1079e-01, 7.1410e-01, 7.1849e-01, 7.2442e-01, 7.3329e-01,
1124 7.5047e-01, 8.0863e-01, 7.6728e-01, 7.5729e-01, 7.5348e-01, 7.5250e-01,
1125 7.5358e-01, 7.5674e-01, 7.6239e-01, 7.7162e-01, 7.8659e-01, 8.1236e-01,
1126 8.6164e-01, 1.0402},
1127 {7.5498e-01, 7.5573e-01, 7.5560e-01, 7.5642e-01, 7.5759e-01, 7.5913e-01,
1128 7.6116e-01, 7.6373e-01, 7.6704e-01, 7.7142e-01, 7.7734e-01, 7.8617e-01,
1129 8.0322e-01, 8.6277e-01, 8.2039e-01, 8.1030e-01, 8.0646e-01, 8.0546e-01,
1130 8.0650e-01, 8.0964e-01, 8.1524e-01, 8.2443e-01, 8.3926e-01, 8.6476e-01,
1131 9.1328e-01, 1.0826},
1132 {8.0804e-01, 8.0817e-01, 8.0865e-01, 8.0945e-01, 8.1060e-01, 8.1219e-01,
1133 8.1420e-01, 8.1674e-01, 8.2006e-01, 8.2444e-01, 8.3032e-01, 8.3912e-01,
1134 8.5604e-01, 9.1704e-01, 8.7360e-01, 8.6340e-01, 8.5952e-01, 8.5850e-01,
1135 8.5951e-01, 8.6262e-01, 8.6815e-01, 8.7727e-01, 8.9200e-01, 9.1725e-01,
1136 9.6506e-01, 1.1263},
1137 {8.6112e-01, 8.6129e-01, 8.6175e-01, 8.6255e-01, 8.6370e-01, 8.6526e-01,
1138 8.6729e-01, 8.6986e-01, 8.7316e-01, 8.7749e-01, 8.8338e-01, 8.9214e-01,
1139 9.0895e-01, 9.7166e-01, 9.2687e-01, 9.1657e-01, 9.1266e-01, 9.1157e-01,
1140 9.1260e-01, 9.1569e-01, 9.2125e-01, 9.3021e-01, 9.4481e-01, 9.6985e-01,
1141 1.0169, 1.1709},
1142 {9.1435e-01, 9.1448e-01, 9.1492e-01, 9.1575e-01, 9.1689e-01, 9.1846e-01,
1143 9.2046e-01, 9.2304e-01, 9.2633e-01, 9.3068e-01, 9.3652e-01, 9.4523e-01,
1144 9.6191e-01, 1.0265, 9.8022e-01, 9.6983e-01, 9.6588e-01, 9.6480e-01,
1145 9.6576e-01, 9.6877e-01, 9.7426e-01, 9.8323e-01, 9.9774e-01, 1.0225,
1146 1.0689, 1.2167},
1147 {9.6759e-01, 9.6773e-01, 9.6823e-01, 9.6900e-01, 9.7017e-01, 9.7170e-01,
1148 9.7372e-01, 9.7628e-01, 9.7957e-01, 9.8389e-01, 9.8975e-01, 9.9845e-01,
1149 1.0150, 1.0814, 1.0336, 1.0232, 1.0192, 1.0181, 1.0190, 1.0220, 1.0274,
1150 1.0363, 1.0507, 1.0752, 1.1210, 1.2635},
1151 {1.0219, 1.0211, 1.0215, 1.0223, 1.0235, 1.0250, 1.0270, 1.0296, 1.0329,
1152 1.0372, 1.0432, 1.0517, 1.0681, 1.1366, 1.0872, 1.0766, 1.0725, 1.0714,
1153 1.0723, 1.0753, 1.0807, 1.0895, 1.1038, 1.1281, 1.1732, 1.3110},
1154 {1.0743, 1.0748, 1.0749, 1.0757, 1.0770, 1.0784, 1.0804, 1.0830, 1.0863,
1155 1.0906, 1.0964, 1.1050, 1.1214, 1.1919, 1.1407, 1.1300, 1.1260, 1.1248,
1156 1.1257, 1.1287, 1.1340, 1.1428, 1.1570, 1.1810, 1.2256, 1.3591},
1157 {1.1278, 1.1279, 1.1284, 1.1292, 1.1304, 1.1319, 1.1339, 1.1364, 1.1397,
1158 1.1440, 1.1498, 1.1584, 1.1747, 1.2472, 1.1943, 1.1836, 1.1795, 1.1783,
1159 1.1792, 1.1821, 1.1874, 1.1961, 1.2102, 1.2340, 1.2781, 1.4079},
1160 {1.1912, 1.1815, 1.1820, 1.1827, 1.1839, 1.1854, 1.1875, 1.1900, 1.1933,
1161 1.1975, 1.2033, 1.2119, 1.2281, 1.3027, 1.2480, 1.2371, 1.2330, 1.2318,
1162 1.2327, 1.2356, 1.2409, 1.2495, 1.2635, 1.2872, 1.3308, 1.4571},
1163 {1.2349, 1.2351, 1.2356, 1.2363, 1.2375, 1.2390, 1.2410, 1.2436, 1.2468,
1164 1.2511, 1.2569, 1.2654, 1.2815, 1.3583, 1.3017, 1.2908, 1.2866, 1.2854,
1165 1.2863, 1.2892, 1.2944, 1.3030, 1.3170, 1.3404, 1.3836, 1.5071},
1166 {1.3170, 1.2887, 1.2892, 1.2900, 1.2911, 1.2927, 1.2947, 1.2972, 1.3005,
1167 1.3047, 1.3105, 1.3190, 1.3350, 1.4139, 1.3555, 1.3445, 1.3403, 1.3391,
1168 1.3399, 1.3428, 1.3480, 1.3566, 1.3704, 1.3937, 1.4366, 1.5575},
1169 {1.3423, 1.3424, 1.3429, 1.3437, 1.3449, 1.3464, 1.3484, 1.3509, 1.3541,
1170 1.3584, 1.3642, 1.3727, 1.3886, 1.4695, 1.4093, 1.3982, 1.3940, 1.3928,
1171 1.3936, 1.3966, 1.4017, 1.4103, 1.4240, 1.4472, 1.4896, 1.6085},
1172 {1.3960, 1.3962, 1.3966, 1.3974, 1.3985, 1.4001, 1.4021, 1.4046, 1.4079,
1173 1.4122, 1.4179, 1.4264, 1.4423, 1.5251, 1.4632, 1.4520, 1.4477, 1.4465,
1174 1.4474, 1.4502, 1.4554, 1.4639, 1.4776, 1.5007, 1.5429, 1.6600},
1175 {1.4498, 1.4499, 1.4504, 1.4512, 1.4524, 1.4539, 1.4559, 1.4584, 1.4616,
1176 1.4659, 1.4717, 1.4801, 1.4959, 1.5806, 1.5170, 1.5058, 1.5016, 1.5003,
1177 1.5011, 1.5040, 1.5091, 1.5177, 1.5313, 1.5542, 1.5962, 1.7119},
1178 {1.5036, 1.5037, 1.5042, 1.5050, 1.5061, 1.5077, 1.5097, 1.5122, 1.5154,
1179 1.5197, 1.5254, 1.5339, 1.5497, 1.6359, 1.5708, 1.5596, 1.5554, 1.5541,
1180 1.5549, 1.5577, 1.5629, 1.5714, 1.5850, 1.6079, 1.6497, 1.7643},
1181 {1.5574, 1.5575, 1.5580, 1.5588, 1.5600, 1.5615, 1.5635, 1.5660, 1.5693,
1182 1.5735, 1.5792, 1.5877, 1.6035, 1.6909, 1.6247, 1.6135, 1.6092, 1.6080,
1183 1.6088, 1.6117, 1.6167, 1.6252, 1.6388, 1.6616, 1.7033, 1.8171},
1184 {1.6173, 1.6114, 1.6119, 1.6127, 1.6138, 1.6153, 1.6173, 1.6198, 1.6231,
1185 1.6274, 1.6331, 1.6415, 1.6573, 1.7454, 1.6786, 1.6674, 1.6631, 1.6618,
1186 1.6626, 1.6654, 1.6706, 1.6791, 1.6926, 1.7154, 1.7570, 1.8704},
1187 {1.6651, 1.6652, 1.6657, 1.6665, 1.6676, 1.6692, 1.6712, 1.6737, 1.6769,
1188 1.6812, 1.6869, 1.6953, 1.7111, 1.7995, 1.7324, 1.7212, 1.7169, 1.7157,
1189 1.7164, 1.7193, 1.7244, 1.7328, 1.7465, 1.7693, 1.8108, 1.9241}};
1190
1191 Float_t nonstagg[ktimebin][kZbin] =
1192 {{0.0912, 0.0913, 0.0914, 0.0917, 0.0920, 0.0925, 0.0931, 0.0938, 0.0946,
1193 0.0954, 0.0964, 0.0975, 0.0988, 0.1001, 0.1017, 0.1034, 0.1053, 0.1075,
1194 0.1101, 0.1135, 0.1179, 0.1239, 0.1332, 0.1504, 0.1935, 0.5719},
1195 {0.0657, 0.0658, 0.0660, 0.0662, 0.0666, 0.0672, 0.0678, 0.0685, 0.0694,
1196 0.0704, 0.0715, 0.0727, 0.0740, 0.0755, 0.0770, 0.0788, 0.0808, 0.0831,
1197 0.0858, 0.0893, 0.0938, 0.1001, 0.1096, 0.1274, 0.1715, 0.5475},
1198 {0.0396, 0.0397, 0.0400, 0.0404, 0.0410, 0.0418, 0.0427, 0.0438, 0.0450,
1199 0.0463, 0.0477, 0.0492, 0.0509, 0.0526, 0.0544, 0.0565, 0.0587, 0.0613,
1200 0.0645, 0.0685, 0.0736, 0.0807, 0.0919, 0.1128, 0.1608, 0.5233},
1201 {0.0112, 0.0115, 0.0124, 0.0137, 0.0154, 0.0173, 0.0195, 0.0218, 0.0241,
1202 0.0265, 0.0289, 0.0313, 0.0336, 0.0360, 0.0384, 0.0410, 0.0439, 0.0473,
1203 0.0515, 0.0571, 0.0651, 0.0767, 0.0949, 0.1243, 0.1826, 0.4858},
1204 {0.0113, 0.0116, 0.0124, 0.0137, 0.0154, 0.0173, 0.0195, 0.0217, 0.0241,
1205 0.0264, 0.0287, 0.0311, 0.0334, 0.0357, 0.0381, 0.0408, 0.0439, 0.0476,
1206 0.0525, 0.0591, 0.0684, 0.0821, 0.1019, 0.1352, 0.1914, 0.4601},
1207 {0.0394, 0.0394, 0.0397, 0.0401, 0.0407, 0.0415, 0.0423, 0.0434, 0.0445,
1208 0.0457, 0.0471, 0.0485, 0.0501, 0.0518, 0.0538, 0.0560, 0.0585, 0.0616,
1209 0.0652, 0.0698, 0.0759, 0.0846, 0.0983, 0.1226, 0.1742, 0.5225},
1210 {0.0650, 0.0650, 0.0652, 0.0655, 0.0658, 0.0663, 0.0669, 0.0676, 0.0685,
1211 0.0694, 0.0704, 0.0716, 0.0729, 0.0744, 0.0761, 0.0780, 0.0802, 0.0828,
1212 0.0860, 0.0899, 0.0950, 0.1020, 0.1130, 0.1334, 0.1809, 0.5444},
1213 {0.0929, 0.0929, 0.0922, 0.0921, 0.0923, 0.0928, 0.0934, 0.0941, 0.0950,
1214 0.0958, 0.0968, 0.0979, 0.0990, 0.1002, 0.1015, 0.1029, 0.1044, 0.1062,
1215 0.1080, 0.1101, 0.1126, 0.1155, 0.1192, 0.1241, 0.1312, 0.0000},
1216 {0.1688, 0.1657, 0.1518, 0.1450, 0.1409, 0.1382, 0.1366, 0.1357, 0.1352,
1217 0.1352, 0.1354, 0.1359, 0.1367, 0.1377, 0.1389, 0.1403, 0.1420, 0.1441,
1218 0.1465, 0.1495, 0.1532, 0.1580, 0.1643, 0.1736, 0.1902, 0.2919},
1219 {0.3413, 0.2191, 0.2052, 0.1981, 0.1936, 0.1905, 0.1883, 0.1869, 0.1860,
1220 0.1856, 0.1856, 0.1859, 0.1865, 0.1875, 0.1888, 0.1904, 0.1924, 0.1948,
1221 0.1977, 0.2011, 0.2054, 0.2107, 0.2175, 0.2272, 0.2440, 0.3413},
1222 {0.4016, 0.2716, 0.2576, 0.2505, 0.2459, 0.2428, 0.2406, 0.2392, 0.2383,
1223 0.2378, 0.2378, 0.2381, 0.2387, 0.2396, 0.2409, 0.2426, 0.2446, 0.2470,
1224 0.2500, 0.2535, 0.2578, 0.2631, 0.2700, 0.2796, 0.2964, 0.3877},
1225 {0.4388, 0.3240, 0.3100, 0.3029, 0.2983, 0.2952, 0.2931, 0.2916, 0.2907,
1226 0.2902, 0.2901, 0.2904, 0.2911, 0.2920, 0.2933, 0.2950, 0.2970, 0.2994,
1227 0.3024, 0.3059, 0.3102, 0.3155, 0.3224, 0.3319, 0.3486, 0.4320},
1228 {0.3580, 0.3764, 0.3625, 0.3554, 0.3509, 0.3477, 0.3455, 0.3440, 0.3431,
1229 0.3427, 0.3426, 0.3429, 0.3435, 0.3445, 0.3458, 0.3474, 0.3494, 0.3519,
1230 0.3548, 0.3583, 0.3626, 0.3679, 0.3747, 0.3843, 0.4008, 0.4768},
1231 {0.4023, 0.4289, 0.4150, 0.4079, 0.4033, 0.4002, 0.3980, 0.3965, 0.3956,
1232 0.3952, 0.3951, 0.3954, 0.3960, 0.3970, 0.3983, 0.3999, 0.4019, 0.4043,
1233 0.4072, 0.4107, 0.4150, 0.4203, 0.4271, 0.4367, 0.4530, 0.5228},
1234 {0.4334, 0.4815, 0.4676, 0.4604, 0.4559, 0.4527, 0.4506, 0.4491, 0.4482,
1235 0.4477, 0.4476, 0.4479, 0.4486, 0.4495, 0.4508, 0.4524, 0.4544, 0.4568,
1236 0.4597, 0.4632, 0.4675, 0.4728, 0.4796, 0.4890, 0.5052, 0.5699},
1237 {0.4695, 0.5341, 0.5202, 0.5130, 0.5085, 0.5053, 0.5032, 0.5017, 0.5008,
1238 0.5003, 0.5002, 0.5005, 0.5012, 0.5021, 0.5034, 0.5050, 0.5070, 0.5095,
1239 0.5123, 0.5158, 0.5200, 0.5253, 0.5321, 0.5415, 0.5576, 0.6179},
1240 {0.5905, 0.5868, 0.5729, 0.5657, 0.5612, 0.5580, 0.5559, 0.5544, 0.5535,
1241 0.5530, 0.5529, 0.5532, 0.5538, 0.5548, 0.5560, 0.5577, 0.5597, 0.5621,
1242 0.5649, 0.5684, 0.5726, 0.5779, 0.5846, 0.5940, 0.6099, 0.6666},
1243 {0.5940, 0.6307, 0.6254, 0.6185, 0.6141, 0.6108, 0.6086, 0.6071, 0.6062,
1244 0.6057, 0.6057, 0.6059, 0.6066, 0.6075, 0.6088, 0.6104, 0.6124, 0.6148,
1245 0.6176, 0.6211, 0.6253, 0.6306, 0.6373, 0.6465, 0.6623, 0.7156},
1246 {0.6986, 0.6924, 0.6784, 0.6713, 0.6667, 0.6636, 0.6614, 0.6599, 0.6590,
1247 0.6585, 0.6585, 0.6587, 0.6594, 0.6603, 0.6616, 0.6632, 0.6651, 0.6675,
1248 0.6704, 0.6739, 0.6781, 0.6833, 0.6900, 0.6992, 0.7147, 0.7651},
1249 {0.7449, 0.7453, 0.7312, 0.7242, 0.7196, 0.7165, 0.7143, 0.7128, 0.7119,
1250 0.7117, 0.7113, 0.7116, 0.7122, 0.7132, 0.7144, 0.7160, 0.7180, 0.7204,
1251 0.7233, 0.7267, 0.7309, 0.7361, 0.7427, 0.7519, 0.7672, 0.8151},
1252 {0.8011, 0.7982, 0.7843, 0.7772, 0.7726, 0.7694, 0.7673, 0.7658, 0.7649,
1253 0.7644, 0.7643, 0.7646, 0.7652, 0.7661, 0.7674, 0.7690, 0.7709, 0.7733,
1254 0.7763, 0.7796, 0.7838, 0.7890, 0.7956, 0.8046, 0.8198, 0.8656},
1255 {0.8261, 0.8513, 0.8374, 0.8302, 0.8256, 0.8225, 0.8203, 0.8188, 0.8179,
1256 0.8174, 0.8173, 0.8176, 0.8182, 0.8191, 0.8204, 0.8220, 0.8240, 0.8263,
1257 0.8292, 0.8326, 0.8367, 0.8419, 0.8485, 0.8575, 0.8725, 0.9163},
1258 {0.9145, 0.9043, 0.8905, 0.8834, 0.8788, 0.8756, 0.8734, 0.8720, 0.8710,
1259 0.8706, 0.8705, 0.8707, 0.8713, 0.8722, 0.8735, 0.8751, 0.8770, 0.8794,
1260 0.8822, 0.8856, 0.8898, 0.8949, 0.9015, 0.9104, 0.9252, 0.9674},
1261 {0.8915, 0.9576, 0.9437, 0.9365, 0.9320, 0.9288, 0.9267, 0.9251, 0.9242,
1262 0.9237, 0.9236, 0.9239, 0.9245, 0.9254, 0.9267, 0.9282, 0.9302, 0.9326,
1263 0.9354, 0.9388, 0.9429, 0.9480, 0.9545, 0.9635, 0.9781, 1.0188},
1264 {1.0290, 1.0109, 0.9970, 0.9898, 0.9852, 0.9821, 0.9799, 0.9784, 0.9775,
1265 0.9770, 0.9769, 0.9772, 0.9778, 0.9787, 0.9800, 0.9815, 0.9835, 0.9858,
1266 0.9886, 0.9920, 0.9961, 1.0012, 1.0077, 1.0165, 1.0310, 1.0705},
1267 {1.0470, 1.1170, 1.0503, 1.0431, 1.0386, 1.0354, 1.0333, 1.0318, 1.0308,
1268 1.0303, 1.0303, 1.0305, 1.0311, 1.0320, 1.0333, 1.0348, 1.0368, 1.0391,
1269 1.0419, 1.0453, 1.0494, 1.0545, 1.0609, 1.0697, 1.0840, 1.1223},
1270 {1.2284, 1.1176, 1.1037, 1.0966, 1.0922, 1.0889, 1.0867, 1.0852, 1.0842,
1271 1.0838, 1.0837, 1.0839, 1.0845, 1.0854, 1.0867, 1.0882, 1.0902, 1.0925,
1272 1.0953, 1.0987, 1.1028, 1.1078, 1.1143, 1.1230, 1.1371, 1.1744},
1273 {1.1548, 1.1712, 1.1572, 1.1501, 1.1455, 1.1424, 1.1402, 1.1387, 1.1377,
1274 1.1373, 1.1371, 1.1374, 1.1380, 1.1389, 1.1401, 1.1417, 1.1436, 1.1460,
1275 1.1488, 1.1521, 1.1562, 1.1613, 1.1676, 1.1763, 1.1903, 1.2267},
1276 {1.2050, 1.2248, 1.2108, 1.2037, 1.1991, 1.1959, 1.1937, 1.1922, 1.1913,
1277 1.1908, 1.1907, 1.1909, 1.1915, 1.1925, 1.1937, 1.1952, 1.1972, 1.1995,
1278 1.2023, 1.2056, 1.2097, 1.2147, 1.2211, 1.2297, 1.2436, 1.2791},
1279 {1.2799, 1.2784, 1.2644, 1.2573, 1.2527, 1.2495, 1.2473, 1.2458, 1.2449,
1280 1.2444, 1.2443, 1.2446, 1.2451, 1.2460, 1.2473, 1.2488, 1.2508, 1.2531,
1281 1.2559, 1.2592, 1.2633, 1.2683, 1.2746, 1.2833, 1.2970, 1.3317},
1282 {1.2857, 1.3320, 1.3181, 1.3109, 1.3063, 1.3032, 1.3010, 1.2995, 1.2986,
1283 1.2981, 1.2980, 1.2982, 1.2988, 1.2997, 1.3009, 1.3025, 1.3044, 1.3067,
1284 1.3095, 1.3129, 1.3169, 1.3219, 1.3283, 1.3368, 1.3505, 1.3845},
1285 {1.3924, 1.3856, 1.3718, 1.3646, 1.3601, 1.3569, 1.3547, 1.3532, 1.3523,
1286 1.3518, 1.3517, 1.3519, 1.3525, 1.3534, 1.3546, 1.3562, 1.3581, 1.3604,
1287 1.3632, 1.3665, 1.3706, 1.3756, 1.3819, 1.3904, 1.4041, 1.4375},
1288 {1.5732, 1.4395, 1.4255, 1.4184, 1.4138, 1.4107, 1.4085, 1.4070, 1.4060,
1289 1.4055, 1.4054, 1.4057, 1.4063, 1.4071, 1.4086, 1.4099, 1.4119, 1.4142,
1290 1.4170, 1.4203, 1.4243, 1.4293, 1.4356, 1.4441, 1.4576, 1.4906},
1291 {1.4352, 1.4933, 1.4793, 1.4722, 1.4676, 1.4644, 1.4623, 1.4607, 1.4598,
1292 1.4593, 1.4592, 1.4595, 1.4600, 1.4609, 1.4622, 1.4637, 1.4656, 1.4680,
1293 1.4707, 1.4741, 1.4781, 1.4831, 1.4894, 1.4979, 1.5113, 1.5439},
1294 {1.4933, 1.5471, 1.5332, 1.5260, 1.5214, 1.5183, 1.5161, 1.5146, 1.5136,
1295 1.5131, 1.5131, 1.5133, 1.5139, 1.5148, 1.5160, 1.5175, 1.5195, 1.5218,
1296 1.5247, 1.5279, 1.5319, 1.5369, 1.5432, 1.5516, 1.5650, 1.5973},
1297 {1.5709, 1.6010, 1.5870, 1.5798, 1.5753, 1.5721, 1.5699, 1.5684, 1.5675,
1298 1.5670, 1.5669, 1.5671, 1.5677, 1.5686, 1.5698, 1.5714, 1.5733, 1.5756,
1299 1.5784, 1.5817, 1.5857, 1.5907, 1.5970, 1.6054, 1.6188, 1.6509},
1300 {1.6675, 1.6548, 1.6408, 1.6337, 1.6291, 1.6259, 1.6238, 1.6223, 1.6213,
1301 1.6208, 1.6207, 1.6210, 1.6216, 1.6225, 1.6237, 1.6252, 1.6272, 1.6295,
1302 1.6322, 1.6356, 1.6396, 1.6445, 1.6508, 1.6593, 1.6727, 1.7046},
1303 {1.7042, 1.7087, 1.6947, 1.6876, 1.6830, 1.6798, 1.6776, 1.6761, 1.6752,
1304 1.6747, 1.6746, 1.6749, 1.6754, 1.6763, 1.6775, 1.6791, 1.6810, 1.6833,
1305 1.6861, 1.6894, 1.6935, 1.6984, 1.7047, 1.7132, 1.7265, 1.7584}};
1306
1307 if (fTimeStruct) delete [] fTimeStruct;
1308 fTimeStruct = new Float_t[ktimebin*kZbin];
1309
1310 for (Int_t ctrt = 0; ctrt<ktimebin; ctrt++)
1311 for (Int_t ctrz = 0; ctrz<kZbin; ctrz++) {
1312 if (fStaggeringOn) fTimeStruct[ctrt+ctrz*ktimebin] = stagg[ctrt][ctrz];
1313 else fTimeStruct[ctrt+ctrz*ktimebin] = nonstagg[ctrt][ctrz];
1314 }
1315
1316}
1317
17b26de4 1318//_____________________________________________________________________________
1319void AliTRDparameter::SamplePRF()
1320{
1321 //
1322 // Samples the pad response function
1323 //
1324
1325 const Int_t kNplan = AliTRDgeometry::kNplan;
1326 const Int_t kPRFbin = 61;
1327
a328fff9 1328 Float_t prf[kNplan][kPRFbin] = { {0.0267681, 0.031056, 0.036111, 0.0419945, 0.0489578,
1329 0.0571508, 0.0666984, 0.078064, 0.0912508, 0.106579,
1330 0.124429, 0.144803, 0.168187, 0.194645, 0.224143,
1331 0.25717, 0.293146, 0.331943, 0.373128, 0.415798,
1332 0.459149, 0.502028, 0.54349, 0.582109, 0.617335,
1333 0.648306, 0.674189, 0.695082, 0.710222, 0.719711,
1334 0.723751, 0.719711, 0.710222, 0.695082, 0.674189,
1335 0.648306, 0.617335, 0.582109, 0.54349, 0.502028,
1336 0.459149, 0.415798, 0.373128, 0.331943, 0.293146,
1337 0.25717, 0.224143, 0.194645, 0.168187, 0.144803,
1338 0.124429, 0.106579, 0.0912508, 0.078064, 0.0666984,
1339 0.0571508, 0.0489578, 0.0419945, 0.036111, 0.031056,
1340 0.0267681},
1341 {0.0240831, 0.0281427, 0.0328987, 0.0384801, 0.0451249,
1342 0.0530119, 0.0623139, 0.0733268, 0.0864095, 0.101704,
1343 0.119473, 0.140308, 0.164191, 0.191333, 0.222141,
1344 0.256551, 0.29433, 0.335294, 0.37892, 0.424206,
1345 0.470202, 0.515621, 0.559255, 0.600068, 0.636821,
1346 0.668836, 0.695767, 0.717076, 0.73233, 0.741818,
1347 0.745491, 0.741818, 0.73233, 0.717076, 0.695767,
1348 0.668836, 0.636821, 0.600068, 0.559255, 0.515621,
1349 0.470202, 0.424206, 0.37892, 0.335294, 0.29433,
1350 0.256551, 0.222141, 0.191333, 0.164191, 0.140308,
1351 0.119473, 0.101704, 0.0864095, 0.0733268, 0.0623139,
1352 0.0530119, 0.0451249, 0.0384801, 0.0328987, 0.0281427,
1353 0.0240831},
1354 {0.0206855, 0.0243011, 0.0285603, 0.0335902, 0.0396015,
1355 0.0467843, 0.0553424, 0.0655316, 0.0776688, 0.0921782,
1356 0.109269, 0.129259, 0.152459, 0.179369, 0.210091,
1357 0.244681, 0.283099, 0.325176, 0.370396, 0.417732,
1358 0.466086, 0.514192, 0.560518, 0.603805, 0.643028,
1359 0.67738, 0.706095, 0.728728, 0.745307, 0.755731,
1360 0.759954, 0.755731, 0.745307, 0.728728, 0.706095,
1361 0.67738, 0.643028, 0.603805, 0.560518, 0.514192,
1362 0.466086, 0.417732, 0.370396, 0.325176, 0.283099,
1363 0.244681, 0.210091, 0.179369, 0.152459, 0.129259,
1364 0.109269, 0.0921782, 0.0776688, 0.0655316, 0.0553424,
1365 0.0467843, 0.0396015, 0.0335902, 0.0285603, 0.0243011,
1366 0.0206855},
1367 {0.0186168, 0.0219999, 0.0260102, 0.0307769, 0.0364947,
1368 0.0433655, 0.0516213, 0.0615466, 0.0734611, 0.0877121,
1369 0.104666, 0.124855, 0.14853, 0.176033, 0.207639,
1370 0.243511, 0.283633, 0.327786, 0.37537, 0.425281,
1371 0.476227, 0.526727, 0.575268, 0.620462, 0.66101,
1372 0.69611, 0.725313, 0.748317, 0.764969, 0.775206,
1373 0.779006, 0.775206, 0.764969, 0.748317, 0.725313,
1374 0.69611, 0.66101, 0.620462, 0.575268, 0.526727,
1375 0.476227, 0.425281, 0.37537, 0.327786, 0.283633,
1376 0.243511, 0.207639, 0.176033, 0.14853, 0.124855,
1377 0.104666, 0.0877121, 0.0734611, 0.0615466, 0.0516213,
1378 0.0433655, 0.0364947, 0.0307769, 0.0260102, 0.0219999,
1379 0.0186168, },
1380 {0.0159737, 0.0189921, 0.0225916, 0.0268927, 0.0320634,
1381 0.0382995, 0.0458393, 0.0549765, 0.0660512, 0.0794439,
1382 0.095565, 0.114844, 0.137713, 0.164586, 0.195824,
1383 0.231681, 0.272223, 0.31727, 0.366248, 0.418078,
1384 0.471358, 0.524425, 0.575561, 0.623193, 0.666055,
1385 0.703243, 0.734192, 0.7586, 0.776331, 0.787347,
1386 0.791646, 0.787347, 0.776331, 0.7586, 0.734192,
1387 0.703243, 0.666055, 0.623193, 0.575561, 0.524425,
1388 0.471358, 0.418078, 0.366248, 0.31727, 0.272223,
1389 0.231681, 0.195824, 0.164586, 0.137713, 0.114844,
1390 0.095565, 0.0794439, 0.0660512, 0.0549765, 0.0458393,
1391 0.0382995, 0.0320634, 0.0268927, 0.0225916, 0.0189921,
1392 0.0159737},
1393 {0.0143532, 0.0171745, 0.0205622, 0.024635, 0.0295461,
1394 0.0354915, 0.0427198, 0.0515368, 0.0623058, 0.0754416,
1395 0.0913994, 0.110662, 0.133721, 0.161057, 0.193093,
1396 0.230138, 0.27229, 0.319327, 0.370597, 0.424961,
1397 0.480832, 0.536333, 0.589555, 0.638805, 0.682781,
1398 0.720621, 0.751845, 0.776243, 0.793764, 0.804426,
1399 0.808259, 0.804426, 0.793764, 0.776243, 0.751845,
1400 0.720621, 0.682781, 0.638805, 0.589555, 0.536333,
1401 0.480832, 0.424961, 0.370597, 0.319327, 0.27229,
1402 0.230138, 0.193093, 0.161057, 0.133721, 0.110662,
1403 0.0913994, 0.0754416, 0.0623058, 0.0515368, 0.0427198,
1404 0.0354915, 0.0295461, 0.024635, 0.0205622, 0.0171745,
1405 0.0143532} };
17b26de4 1406
1407 // More sampling precision with linear interpolation
1408 fPRFlo = -1.5;
1409 fPRFhi = 1.5;
1410 Float_t pad[kPRFbin];
1411 Int_t sPRFbin = kPRFbin;
1412 Float_t sPRFwid = (fPRFhi - fPRFlo) / ((Float_t) sPRFbin);
1413 for (Int_t iPad = 0; iPad < sPRFbin; iPad++) {
1414 pad[iPad] = ((Float_t) iPad + 0.5) * sPRFwid + fPRFlo;
1415 }
1416 fPRFbin = 500;
1417 fPRFwid = (fPRFhi - fPRFlo) / ((Float_t) fPRFbin);
1418 fPRFpad = ((Int_t) (1.0 / fPRFwid));
1419
1420 if (fPRFsmp) delete [] fPRFsmp;
1421 fPRFsmp = new Float_t[kNplan*fPRFbin];
1422
1423 Int_t ipos1;
1424 Int_t ipos2;
1425 Float_t diff;
1426
1427 for (Int_t iPla = 0; iPla < kNplan; iPla++) {
1428
1429 for (Int_t iBin = 0; iBin < fPRFbin; iBin++) {
1430
1431 Float_t bin = (((Float_t) iBin) + 0.5) * fPRFwid + fPRFlo;
1432 ipos1 = ipos2 = 0;
1433 diff = 0;
1434 do {
1435 diff = bin - pad[ipos2++];
8e64dd77 1436 } while ((diff > 0) && (ipos2 < kPRFbin));
1437 if (ipos2 == kPRFbin) {
1438 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2-1];
1439 }
1440 else if (ipos2 == 1) {
1441 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2-1];
1442 }
1443 else {
1444 ipos2--;
1445 if (ipos2 >= kPRFbin) ipos2 = kPRFbin - 1;
1446 ipos1 = ipos2 - 1;
1447 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2]
1448 + diff * (prf[iPla][ipos2] - prf[iPla][ipos1])
1449 / sPRFwid;
1450 }
17b26de4 1451
1452 }
1453 }
1454
1455}
1456
1457//_____________________________________________________________________________
1458void AliTRDparameter::FillLUT()
1459{
1460 //
1461 // Create the LUT
1462 //
1463
1464 const Int_t kNplan = AliTRDgeometry::kNplan;
1465 const Int_t kNlut = 128;
1466
1467 fLUTbin = kNplan * kNlut;
1468
1469 // The lookup table from Bogdan
1470 Float_t lut[kNplan][kNlut] = {
1471 {
1472 0.0070, 0.0150, 0.0224, 0.0298, 0.0374, 0.0454, 0.0533, 0.0611,
1473 0.0684, 0.0755, 0.0827, 0.0900, 0.0975, 0.1049, 0.1120, 0.1187,
1474 0.1253, 0.1318, 0.1385, 0.1453, 0.1519, 0.1584, 0.1646, 0.1704,
1475 0.1762, 0.1821, 0.1879, 0.1938, 0.1996, 0.2053, 0.2108, 0.2160,
1476 0.2210, 0.2260, 0.2310, 0.2361, 0.2411, 0.2461, 0.2509, 0.2557,
1477 0.2602, 0.2646, 0.2689, 0.2732, 0.2774, 0.2816, 0.2859, 0.2901,
1478 0.2942, 0.2983, 0.3022, 0.3061, 0.3099, 0.3136, 0.3172, 0.3207,
1479 0.3242, 0.3278, 0.3312, 0.3347, 0.3382, 0.3416, 0.3450, 0.3483,
1480 0.3515, 0.3547, 0.3579, 0.3609, 0.3639, 0.3669, 0.3698, 0.3727,
1481 0.3756, 0.3785, 0.3813, 0.3842, 0.3870, 0.3898, 0.3926, 0.3952,
1482 0.3979, 0.4005, 0.4032, 0.4057, 0.4082, 0.4108, 0.4132, 0.4157,
1483 0.4181, 0.4205, 0.4228, 0.4252, 0.4275, 0.4299, 0.4322, 0.4345,
1484 0.4367, 0.4390, 0.4412, 0.4434, 0.4456, 0.4478, 0.4499, 0.4520,
1485 0.4541, 0.4562, 0.4583, 0.4603, 0.4623, 0.4643, 0.4663, 0.4683,
1486 0.4702, 0.4722, 0.4741, 0.4758, 0.4774, 0.4790, 0.4805, 0.4824,
1487 0.4844, 0.4863, 0.4883, 0.4902, 0.4921, 0.4940, 0.4959, 0.4978
1488 },
1489 {
1490 0.0072, 0.0156, 0.0235, 0.0313, 0.0394, 0.0478, 0.0561, 0.0642,
1491 0.0718, 0.0792, 0.0868, 0.0947, 0.1025, 0.1101, 0.1172, 0.1241,
1492 0.1309, 0.1378, 0.1449, 0.1518, 0.1586, 0.1650, 0.1710, 0.1770,
1493 0.1830, 0.1891, 0.1952, 0.2011, 0.2070, 0.2125, 0.2177, 0.2229,
1494 0.2280, 0.2332, 0.2383, 0.2435, 0.2484, 0.2533, 0.2581, 0.2627,
1495 0.2670, 0.2714, 0.2757, 0.2799, 0.2842, 0.2884, 0.2927, 0.2968,
1496 0.3008, 0.3048, 0.3086, 0.3123, 0.3159, 0.3195, 0.3231, 0.3266,
1497 0.3301, 0.3335, 0.3370, 0.3404, 0.3438, 0.3471, 0.3504, 0.3536,
1498 0.3567, 0.3598, 0.3628, 0.3657, 0.3686, 0.3715, 0.3744, 0.3772,
1499 0.3800, 0.3828, 0.3856, 0.3884, 0.3911, 0.3938, 0.3965, 0.3991,
1500 0.4016, 0.4042, 0.4067, 0.4092, 0.4116, 0.4140, 0.4164, 0.4187,
1501 0.4211, 0.4234, 0.4257, 0.4280, 0.4302, 0.4325, 0.4347, 0.4369,
1502 0.4391, 0.4413, 0.4434, 0.4456, 0.4477, 0.4497, 0.4518, 0.4538,
1503 0.4558, 0.4578, 0.4598, 0.4618, 0.4637, 0.4656, 0.4675, 0.4694,
1504 0.4713, 0.4732, 0.4750, 0.4766, 0.4781, 0.4797, 0.4813, 0.4832,
1505 0.4851, 0.4870, 0.4888, 0.4906, 0.4925, 0.4942, 0.4960, 0.4978
1506 },
1507 {
1508 0.0075, 0.0163, 0.0246, 0.0328, 0.0415, 0.0504, 0.0592, 0.0674,
1509 0.0753, 0.0832, 0.0914, 0.0996, 0.1077, 0.1154, 0.1225, 0.1296,
1510 0.1369, 0.1442, 0.1515, 0.1585, 0.1652, 0.1714, 0.1776, 0.1839,
1511 0.1902, 0.1965, 0.2025, 0.2085, 0.2141, 0.2194, 0.2247, 0.2299,
1512 0.2352, 0.2405, 0.2457, 0.2507, 0.2557, 0.2604, 0.2649, 0.2693,
1513 0.2737, 0.2780, 0.2823, 0.2867, 0.2909, 0.2951, 0.2992, 0.3033,
1514 0.3072, 0.3110, 0.3146, 0.3182, 0.3218, 0.3253, 0.3288, 0.3323,
1515 0.3357, 0.3392, 0.3426, 0.3459, 0.3492, 0.3524, 0.3555, 0.3586,
1516 0.3616, 0.3645, 0.3674, 0.3703, 0.3731, 0.3759, 0.3787, 0.3815,
1517 0.3843, 0.3870, 0.3897, 0.3925, 0.3950, 0.3976, 0.4002, 0.4027,
1518 0.4052, 0.4076, 0.4101, 0.4124, 0.4148, 0.4171, 0.4194, 0.4217,
1519 0.4239, 0.4262, 0.4284, 0.4306, 0.4328, 0.4350, 0.4371, 0.4393,
1520 0.4414, 0.4435, 0.4455, 0.4476, 0.4496, 0.4516, 0.4536, 0.4555,
1521 0.4575, 0.4594, 0.4613, 0.4632, 0.4650, 0.4669, 0.4687, 0.4705,
1522 0.4723, 0.4741, 0.4758, 0.4773, 0.4789, 0.4804, 0.4821, 0.4839,
1523 0.4857, 0.4875, 0.4893, 0.4910, 0.4928, 0.4945, 0.4961, 0.4978
1524 },
1525 {
1526 0.0078, 0.0171, 0.0258, 0.0345, 0.0438, 0.0532, 0.0624, 0.0708,
1527 0.0791, 0.0875, 0.0962, 0.1048, 0.1130, 0.1206, 0.1281, 0.1356,
1528 0.1432, 0.1508, 0.1582, 0.1651, 0.1716, 0.1780, 0.1845, 0.1910,
1529 0.1975, 0.2038, 0.2099, 0.2155, 0.2210, 0.2263, 0.2317, 0.2371,
1530 0.2425, 0.2477, 0.2528, 0.2578, 0.2626, 0.2671, 0.2715, 0.2759,
1531 0.2803, 0.2846, 0.2890, 0.2933, 0.2975, 0.3016, 0.3056, 0.3095,
1532 0.3132, 0.3168, 0.3204, 0.3239, 0.3274, 0.3309, 0.3344, 0.3378,
1533 0.3412, 0.3446, 0.3479, 0.3511, 0.3543, 0.3574, 0.3603, 0.3633,
1534 0.3662, 0.3690, 0.3718, 0.3747, 0.3774, 0.3802, 0.3829, 0.3857,
1535 0.3883, 0.3910, 0.3936, 0.3962, 0.3987, 0.4012, 0.4037, 0.4061,
1536 0.4085, 0.4109, 0.4132, 0.4155, 0.4177, 0.4200, 0.4222, 0.4244,
1537 0.4266, 0.4288, 0.4309, 0.4331, 0.4352, 0.4373, 0.4394, 0.4414,
1538 0.4435, 0.4455, 0.4475, 0.4494, 0.4514, 0.4533, 0.4552, 0.4571,
1539 0.4590, 0.4608, 0.4626, 0.4645, 0.4662, 0.4680, 0.4698, 0.4715,
1540 0.4733, 0.4750, 0.4766, 0.4781, 0.4796, 0.4812, 0.4829, 0.4846,
1541 0.4863, 0.4880, 0.4897, 0.4914, 0.4930, 0.4946, 0.4963, 0.4979
1542 },
1543 {
1544 0.0081, 0.0178, 0.0270, 0.0364, 0.0463, 0.0562, 0.0656, 0.0744,
1545 0.0831, 0.0921, 0.1013, 0.1102, 0.1183, 0.1261, 0.1339, 0.1419,
1546 0.1499, 0.1576, 0.1648, 0.1715, 0.1782, 0.1849, 0.1917, 0.1984,
1547 0.2048, 0.2110, 0.2167, 0.2223, 0.2278, 0.2333, 0.2389, 0.2444,
1548 0.2497, 0.2548, 0.2598, 0.2645, 0.2691, 0.2735, 0.2780, 0.2824,
1549 0.2868, 0.2912, 0.2955, 0.2997, 0.3038, 0.3078, 0.3116, 0.3152,
1550 0.3188, 0.3224, 0.3259, 0.3294, 0.3329, 0.3364, 0.3398, 0.3432,
1551 0.3465, 0.3497, 0.3529, 0.3561, 0.3591, 0.3620, 0.3649, 0.3677,
1552 0.3705, 0.3733, 0.3761, 0.3788, 0.3816, 0.3843, 0.3869, 0.3896,
1553 0.3922, 0.3948, 0.3973, 0.3998, 0.4022, 0.4047, 0.4070, 0.4094,
1554 0.4117, 0.4139, 0.4162, 0.4184, 0.4206, 0.4227, 0.4249, 0.4270,
1555 0.4291, 0.4313, 0.4334, 0.4354, 0.4375, 0.4395, 0.4415, 0.4435,
1556 0.4455, 0.4474, 0.4493, 0.4512, 0.4531, 0.4550, 0.4568, 0.4586,
1557 0.4604, 0.4622, 0.4639, 0.4657, 0.4674, 0.4691, 0.4708, 0.4725,
1558 0.4742, 0.4758, 0.4773, 0.4788, 0.4803, 0.4819, 0.4836, 0.4852,
1559 0.4869, 0.4885, 0.4901, 0.4917, 0.4933, 0.4948, 0.4964, 0.4979
1560 },
1561 {
1562 0.0085, 0.0189, 0.0288, 0.0389, 0.0497, 0.0603, 0.0699, 0.0792,
1563 0.0887, 0.0985, 0.1082, 0.1170, 0.1253, 0.1336, 0.1421, 0.1505,
1564 0.1587, 0.1662, 0.1733, 0.1803, 0.1874, 0.1945, 0.2014, 0.2081,
1565 0.2143, 0.2201, 0.2259, 0.2316, 0.2374, 0.2431, 0.2487, 0.2541,
1566 0.2593, 0.2642, 0.2689, 0.2735, 0.2781, 0.2826, 0.2872, 0.2917,
1567 0.2961, 0.3003, 0.3045, 0.3086, 0.3125, 0.3162, 0.3198, 0.3235,
1568 0.3270, 0.3306, 0.3342, 0.3377, 0.3411, 0.3446, 0.3479, 0.3511,
1569 0.3543, 0.3575, 0.3605, 0.3634, 0.3663, 0.3691, 0.3720, 0.3748,
1570 0.3775, 0.3803, 0.3830, 0.3857, 0.3884, 0.3911, 0.3937, 0.3962,
1571 0.3987, 0.4012, 0.4036, 0.4060, 0.4084, 0.4107, 0.4129, 0.4152,
1572 0.4174, 0.4196, 0.4218, 0.4239, 0.4261, 0.4282, 0.4303, 0.4324,
1573 0.4344, 0.4365, 0.4385, 0.4405, 0.4425, 0.4445, 0.4464, 0.4483,
1574 0.4502, 0.4521, 0.4539, 0.4558, 0.4576, 0.4593, 0.4611, 0.4629,
1575 0.4646, 0.4663, 0.4680, 0.4697, 0.4714, 0.4730, 0.4747, 0.4759,
1576 0.4769, 0.4780, 0.4790, 0.4800, 0.4811, 0.4827, 0.4843, 0.4859,
1577 0.4874, 0.4889, 0.4905, 0.4920, 0.4935, 0.4950, 0.4965, 0.4979
1578 }
1579 };
1580
1581 if (fLUT) delete [] fLUT;
1582 fLUT = new Float_t[fLUTbin];
1583
1584 for (Int_t iplan = 0; iplan < kNplan; iplan++) {
1585 for (Int_t ilut = 0; ilut < kNlut; ilut++) {
cfed70d4 1586 fLUT[iplan*kNlut+ilut] = lut[iplan][ilut];
17b26de4 1587 }
1588 }
1589
1590}
1591
1592//_____________________________________________________________________________
1593void AliTRDparameter::SetTiltingAngle(Float_t v)
1594{
1595 //
1596 // Set the tilting angle for the readout pads
1597 //
1598
1599 fTiltingAngle = TMath::Tan(TMath::Pi()/180.0 * v);
1600
1601}
1602
1603//_____________________________________________________________________________
1604Float_t AliTRDparameter::GetTiltingAngle() const
1605{
1606 //
1607 // Get the tilting angle for the readout pads
1608 //
1609
ab35e265 1610 return 180.0 / TMath::Pi() * TMath::ATan(fTiltingAngle);
17b26de4 1611
1612}
1613
1614//_____________________________________________________________________________
1615Float_t AliTRDparameter::GetDiffusionL(Float_t vd, Float_t b)
1616{
1617 //
1618 // Returns the longitudinal diffusion coefficient for a given drift
1619 // velocity <vd> and a B-field <b> for Xe/CO2 (15%).
1620 // The values are according to a GARFIELD simulation.
1621 //
1622
1623 const Int_t kNb = 5;
1624 Float_t p0[kNb] = { 0.007440, 0.007493, 0.007513, 0.007672, 0.007831 };
1625 Float_t p1[kNb] = { 0.019252, 0.018912, 0.018636, 0.018012, 0.017343 };
1626 Float_t p2[kNb] = { -0.005042, -0.004926, -0.004867, -0.004650, -0.004424 };
1627 Float_t p3[kNb] = { 0.000195, 0.000189, 0.000195, 0.000182, 0.000169 };
1628
1629 Int_t ib = ((Int_t) (10 * (b - 0.15)));
1630 ib = TMath::Max( 0,ib);
1631 ib = TMath::Min(kNb,ib);
1632
1633 Float_t diff = p0[ib]
1634 + p1[ib] * vd
1635 + p2[ib] * vd*vd
1636 + p3[ib] * vd*vd*vd;
1637
1638 return diff;
1639
1640}
1641
1642//_____________________________________________________________________________
1643Float_t AliTRDparameter::GetDiffusionT(Float_t vd, Float_t b)
1644{
1645 //
1646 // Returns the transverse diffusion coefficient for a given drift
1647 // velocity <vd> and a B-field <b> for Xe/CO2 (15%).
1648 // The values are according to a GARFIELD simulation.
1649 //
1650
1651 const Int_t kNb = 5;
1652 Float_t p0[kNb] = { 0.009550, 0.009599, 0.009674, 0.009757, 0.009850 };
1653 Float_t p1[kNb] = { 0.006667, 0.006539, 0.006359, 0.006153, 0.005925 };
1654 Float_t p2[kNb] = { -0.000853, -0.000798, -0.000721, -0.000635, -0.000541 };
1655 Float_t p3[kNb] = { 0.000131, 0.000122, 0.000111, 0.000098, 0.000085 };
1656
1657 Int_t ib = ((Int_t) (10 * (b - 0.15)));
1658 ib = TMath::Max( 0,ib);
1659 ib = TMath::Min(kNb,ib);
1660
1661 Float_t diff = p0[ib]
1662 + p1[ib] * vd
1663 + p2[ib] * vd*vd
1664 + p3[ib] * vd*vd*vd;
1665
1666 return diff;
1667
1668}
1669
1670//_____________________________________________________________________________
1671Float_t AliTRDparameter::GetOmegaTau(Float_t vd, Float_t b)
1672{
1673 //
1674 // Returns omega*tau (tan(Lorentz-angle)) for a given drift velocity <vd>
1675 // and a B-field <b> for Xe/CO2 (15%).
1676 // The values are according to a GARFIELD simulation.
1677 //
1678
1679 const Int_t kNb = 5;
1680 Float_t p0[kNb] = { 0.004810, 0.007412, 0.010252, 0.013409, 0.016888 };
1681 Float_t p1[kNb] = { 0.054875, 0.081534, 0.107333, 0.131983, 0.155455 };
1682 Float_t p2[kNb] = { -0.008682, -0.012896, -0.016987, -0.020880, -0.024623 };
1683 Float_t p3[kNb] = { 0.000155, 0.000238, 0.000330, 0.000428, 0.000541 };
1684
1685 Int_t ib = ((Int_t) (10 * (b - 0.15)));
1686 ib = TMath::Max( 0,ib);
1687 ib = TMath::Min(kNb,ib);
1688
1689 Float_t alphaL = p0[ib]
1690 + p1[ib] * vd
1691 + p2[ib] * vd*vd
1692 + p3[ib] * vd*vd*vd;
1693
1694 return TMath::Tan(alphaL);
1695
1696}
1697
1698//_____________________________________________________________________________
0a29d0f1 1699Double_t AliTRDparameter::LUTposition(Int_t iplane, Double_t ampL
1700 , Double_t ampC
1701 , Double_t ampR) const
17b26de4 1702{
1703 //
1704 // Calculates the cluster position using the lookup table.
1705 // Method provided by Bogdan Vulpescu.
1706 //
1707
1708 const Int_t kNplan = AliTRDgeometry::kNplan;
cfed70d4 1709 const Int_t kNlut = 128;
17b26de4 1710
1711 Double_t pos;
1712 Double_t x = 0.0;
1713 Double_t xmin;
1714 Double_t xmax;
1715 Double_t xwid;
1716
1717 Int_t side = 0;
1718 Int_t ix;
1719
1720 Double_t xMin[kNplan] = { 0.006492, 0.006377, 0.006258, 0.006144, 0.006030, 0.005980 };
1721 Double_t xMax[kNplan] = { 0.960351, 0.965870, 0.970445, 0.974352, 0.977667, 0.996101 };
1722
1723 if (ampL > ampR) {
1724 x = (ampL - ampR) / ampC;
1725 side = -1;
1726 }
1727 else if (ampL < ampR) {
1728 x = (ampR - ampL) / ampC;
1729 side = +1;
1730 }
1731
1732 if (ampL != ampR) {
1733
1734 xmin = xMin[iplane] + 0.000005;
1735 xmax = xMax[iplane] - 0.000005;
1736 xwid = (xmax - xmin) / 127.0;
1737
1738 if (x < xmin) {
1739 pos = 0.0000;
1740 }
1741 else if (x > xmax) {
1742 pos = side * 0.5000;
1743 }
1744 else {
1745 ix = (Int_t) ((x - xmin) / xwid);
cfed70d4 1746 pos = side * fLUT[iplane*kNlut+ix];
17b26de4 1747 }
1748
1749 }
1750 else {
1751
1752 pos = 0.0;
1753
1754 }
1755
1756 return pos;
1757
1758}