]> 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
Possibility to change sign of Pz event by event, needed for pA, Ap simulation in
[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 fDiffusionT = 0.0;
53 fDiffusionL = 0.0;
17b26de4 54 fElAttachProp = 0.0;
17b26de4 55 fOmegaTau = 0.0;
17b26de4 56 fDriftVelocity = 0.0;
ccb4315c 57 fSamplingFrequency = 0.0;
17b26de4 58 fPadCoupling = 0.0;
59 fTimeCoupling = 0.0;
17b26de4 60 fField = 0.0;
61 fTiltingAngle = 0.0;
62 fPRFbin = 0;
63 fPRFlo = 0.0;
64 fPRFhi = 0.0;
65 fPRFwid = 0.0;
66 fPRFpad = 0;
67 fTRFbin = 0;
68 fTRFlo = 0.0;
69 fTRFhi = 0.0;
70 fTRFwid = 0.0;
71 fTCnexp = 0;
72
17b26de4 73 fLUT = 0;
74 fClusMaxThresh = 0;
75 fClusSigThresh = 0;
76
ccb4315c 77 fTimeStruct1 = 0;
78 fTimeStruct2 = 0;
a328fff9 79 fAnodeWireOffset = 0.0;
80
ccb4315c 81 fDiffusionOn = kFALSE;
82 fElAttachOn = kFALSE;
83 fExBOn = kFALSE;
84 fPRFOn = kFALSE;
85 fTRFOn = kFALSE;
86 fCTOn = kFALSE;
87 fTCOn = kFALSE;
88 fLUTOn = kFALSE;
89 fTimeStructOn = kFALSE;
90
17b26de4 91}
92
93//_____________________________________________________________________________
94AliTRDparameter::AliTRDparameter(const Text_t *name, const Text_t *title)
95 :TNamed(name,title)
96{
97 //
98 // AliTRDparameter constructor
99 //
100
101 fGeo = new AliTRDgeometryFull();
67689223 102 fPRFsmp = 0;
103 fTRFsmp = 0;
104 fCTsmp = 0;
105 fGasGain = 0.0;
106 fNoise = 0.0;
107 fChipGain = 0.0;
108 fADCoutRange = 0.0;
109 fADCinRange = 0.0;
110 fADCthreshold = 0;
7e582e9f 111 fADCbaseline = 0;
67689223 112 fDiffusionT = 0.0;
113 fDiffusionL = 0.0;
67689223 114 fElAttachProp = 0.0;
115 fExBOn = 0;
116 fOmegaTau = 0.0;
67689223 117 fDriftVelocity = 0.0;
ccb4315c 118 fSamplingFrequency = 0.0;
67689223 119 fPadCoupling = 0.0;
120 fTimeCoupling = 0.0;
67689223 121 fField = 0.0;
122 fTiltingAngle = 0.0;
123 fPRFbin = 0;
124 fPRFlo = 0.0;
125 fPRFhi = 0.0;
126 fPRFwid = 0.0;
127 fPRFpad = 0;
128 fTRFbin = 0;
129 fTRFlo = 0.0;
130 fTRFhi = 0.0;
131 fTRFwid = 0.0;
132 fTCnexp = 0;
133
67689223 134 fLUT = 0;
135 fClusMaxThresh = 0;
136 fClusSigThresh = 0;
17b26de4 137
ccb4315c 138 fTimeStruct1 = 0;
139 fTimeStruct2 = 0;
a328fff9 140 fAnodeWireOffset = 0.0;
141
ccb4315c 142 fDiffusionOn = kFALSE;
143 fElAttachOn = kFALSE;
144 fExBOn = kFALSE;
145 fPRFOn = kFALSE;
146 fTRFOn = kFALSE;
147 fCTOn = kFALSE;
148 fTCOn = kFALSE;
149 fLUTOn = kFALSE;
150 fTimeStructOn = kFALSE;
151
17b26de4 152 Init();
153
154}
155
156
157//_____________________________________________________________________________
73ae7b59 158AliTRDparameter::AliTRDparameter(const AliTRDparameter &p):TNamed(p)
17b26de4 159{
160 //
161 // AliTRDparameter copy constructor
162 //
163
164 ((AliTRDparameter &) p).Copy(*this);
165
166}
167
168///_____________________________________________________________________________
169AliTRDparameter::~AliTRDparameter()
170{
171 //
172 // AliTRDparameter destructor
173 //
174
175 if (fTRFsmp) {
176 delete [] fTRFsmp;
177 fTRFsmp = 0;
178 }
179
180 if (fPRFsmp) {
181 delete [] fPRFsmp;
182 fPRFsmp = 0;
183 }
184
185 if (fCTsmp) {
186 delete [] fCTsmp;
187 fCTsmp = 0;
188 }
189
190 if (fLUT) {
191 delete [] fLUT;
192 fLUT = 0;
193 }
194
195 if (fGeo) {
196 delete fGeo;
197 fGeo = 0;
198 }
199
ccb4315c 200 if (fTimeStruct1) {
201 delete [] fTimeStruct1;
202 fTimeStruct1 = 0;
203 }
204
205 if (fTimeStruct2) {
206 delete [] fTimeStruct2;
207 fTimeStruct2 = 0;
a328fff9 208 }
209
17b26de4 210}
211
212//_____________________________________________________________________________
213AliTRDparameter &AliTRDparameter::operator=(const AliTRDparameter &p)
214{
215 //
216 // Assignment operator
217 //
218
219 if (this != &p) ((AliTRDparameter &) p).Copy(*this);
220 return *this;
221
222}
223
224//_____________________________________________________________________________
e0d47c25 225void AliTRDparameter::Copy(TObject &p) const
17b26de4 226{
227 //
228 // Copy function
229 //
230
231 Int_t iBin;
232
233 ((AliTRDparameter &) p).fGasGain = fGasGain;
234 ((AliTRDparameter &) p).fNoise = fNoise;
235 ((AliTRDparameter &) p).fChipGain = fChipGain;
236 ((AliTRDparameter &) p).fADCoutRange = fADCoutRange;
237 ((AliTRDparameter &) p).fADCinRange = fADCinRange;
238 ((AliTRDparameter &) p).fADCthreshold = fADCthreshold;
7e582e9f 239 ((AliTRDparameter &) p).fADCbaseline = fADCbaseline;
17b26de4 240 ((AliTRDparameter &) p).fDiffusionOn = fDiffusionOn;
241 ((AliTRDparameter &) p).fDiffusionT = fDiffusionT;
242 ((AliTRDparameter &) p).fDiffusionL = fDiffusionL;
243 ((AliTRDparameter &) p).fElAttachOn = fElAttachOn;
244 ((AliTRDparameter &) p).fElAttachProp = fElAttachProp;
245 ((AliTRDparameter &) p).fExBOn = fExBOn;
246 ((AliTRDparameter &) p).fOmegaTau = fOmegaTau;
247 ((AliTRDparameter &) p).fLorentzFactor = fLorentzFactor;
248 ((AliTRDparameter &) p).fDriftVelocity = fDriftVelocity;
ccb4315c 249 ((AliTRDparameter &) p).fSamplingFrequency = fSamplingFrequency;
17b26de4 250 ((AliTRDparameter &) p).fPadCoupling = fPadCoupling;
251 ((AliTRDparameter &) p).fTimeCoupling = fTimeCoupling;
17b26de4 252 ((AliTRDparameter &) p).fField = fField;
253 ((AliTRDparameter &) p).fPRFOn = fPRFOn;
254 ((AliTRDparameter &) p).fTRFOn = fTRFOn;
255 ((AliTRDparameter &) p).fCTOn = fCTOn;
256 ((AliTRDparameter &) p).fTCOn = fTCOn;
257 ((AliTRDparameter &) p).fTiltingAngle = fTiltingAngle;
258 ((AliTRDparameter &) p).fPRFbin = fPRFbin;
259 ((AliTRDparameter &) p).fPRFlo = fPRFlo;
260 ((AliTRDparameter &) p).fPRFhi = fPRFhi;
261 ((AliTRDparameter &) p).fPRFwid = fPRFwid;
262 ((AliTRDparameter &) p).fPRFpad = fPRFpad;
263 if (((AliTRDparameter &) p).fPRFsmp) delete [] ((AliTRDparameter &) p).fPRFsmp;
264 ((AliTRDparameter &) p).fPRFsmp = new Float_t[fPRFbin];
265 for (iBin = 0; iBin < fPRFbin; iBin++) {
266 ((AliTRDparameter &) p).fPRFsmp[iBin] = fPRFsmp[iBin];
267 }
268 ((AliTRDparameter &) p).fTRFbin = fTRFbin;
269 ((AliTRDparameter &) p).fTRFlo = fTRFlo;
270 ((AliTRDparameter &) p).fTRFhi = fTRFhi;
271 ((AliTRDparameter &) p).fTRFwid = fTRFwid;
272 if (((AliTRDparameter &) p).fTRFsmp) delete [] ((AliTRDparameter &) p).fTRFsmp;
273 ((AliTRDparameter &) p).fTRFsmp = new Float_t[fTRFbin];
274 for (iBin = 0; iBin < fTRFbin; iBin++) {
275 ((AliTRDparameter &) p).fTRFsmp[iBin] = fTRFsmp[iBin];
276 }
277 if (((AliTRDparameter &) p).fCTsmp) delete [] ((AliTRDparameter &) p).fCTsmp;
278 ((AliTRDparameter &) p).fCTsmp = new Float_t[fTRFbin];
279 for (iBin = 0; iBin < fTRFbin; iBin++) {
280 ((AliTRDparameter &) p).fCTsmp[iBin] = fCTsmp[iBin];
281 }
282 ((AliTRDparameter &) p).fTCnexp = fTCnexp;
283
284 ((AliTRDparameter &) p).fLUTOn = fLUTOn;
285 ((AliTRDparameter &) p).fLUTbin = fLUTbin;
286 if (((AliTRDparameter &) p).fLUT) delete [] ((AliTRDparameter &) p).fLUT;
287 ((AliTRDparameter &) p).fLUT = new Float_t[fLUTbin];
288 for (iBin = 0; iBin < fLUTbin; iBin++) {
289 ((AliTRDparameter &) p).fLUT[iBin] = fLUT[iBin];
290 }
291 ((AliTRDparameter &) p).fClusMaxThresh = fClusMaxThresh;
292 ((AliTRDparameter &) p).fClusSigThresh = fClusSigThresh;
293
a328fff9 294 ((AliTRDparameter &) p).fAnodeWireOffset = fAnodeWireOffset;
a328fff9 295 ((AliTRDparameter &) p).fTimeStructOn = fTimeStructOn;
ccb4315c 296 if (((AliTRDparameter &) p).fTimeStruct1)
297 delete [] ((AliTRDparameter &) p).fTimeStruct1;
298 ((AliTRDparameter &) p).fTimeStruct1 = new Float_t[38*11];
682619b9 299 for (Int_t i = 0; i < 38*11; i++) {
ccb4315c 300 ((AliTRDparameter &) p).fTimeStruct1[i] = fTimeStruct1[i];
301 }
302 if (((AliTRDparameter &) p).fTimeStruct2)
303 delete [] ((AliTRDparameter &) p).fTimeStruct2;
304 ((AliTRDparameter &) p).fTimeStruct2 = new Float_t[38*11];
305 for (Int_t i = 0; i < 38*11; i++) {
306 ((AliTRDparameter &) p).fTimeStruct2[i] = fTimeStruct2[i];
307 }
a328fff9 308
17b26de4 309}
310
311//_____________________________________________________________________________
312void AliTRDparameter::Init()
313{
314 //
315 // Initializes the parameter
316 //
317 // The maximum number of pads
318 // and the position of pad 0,0,0
319 //
320 // chambers seen from the top:
321 // +----------------------------+
322 // | |
323 // | | ^
324 // | | rphi|
325 // | | |
326 // |0 | |
327 // +----------------------------+ +------>
328 // z
329 // chambers seen from the side: ^
330 // +----------------------------+ drift|
331 // |0 | |
332 // | | |
333 // +----------------------------+ +------>
334 // z
335 //
336 // IMPORTANT: time bin 0 is now the first one in the drift region
337 // closest to the readout !!!
338 //
339
340 //
341 // ----------------------------------------------------------------------------
342 // The pad definition
343 // ----------------------------------------------------------------------------
344 //
345
346 // The pad size in column direction (rphi-direction)
ccb4315c 347
348 //SetColPadSize(0,0.65);
349 //SetColPadSize(1,0.68);
350 //SetColPadSize(2,0.71);
351 //SetColPadSize(3,0.74);
352 //SetColPadSize(4,0.77);
353 //SetColPadSize(5,0.80);
354
a328fff9 355 SetColPadSize(0,0.664);
356 SetColPadSize(1,0.695);
357 SetColPadSize(2,0.726);
358 SetColPadSize(3,0.756);
359 SetColPadSize(4,0.788);
360 SetColPadSize(5,0.818);
73ae7b59 361
17b26de4 362 // The pad row (z-direction)
363 SetNRowPad();
364
ccb4315c 365 // !CHANGED!
366 // Time position of anode wire plane
367 for (Int_t iplan = 0; iplan < AliTRDgeometry::Nplan(); iplan++) {
368 fTime0[iplan] = AliTRDgeometry::Rmin() + AliTRDgeometry::CraHght() + AliTRDgeometry::CdrHght()
369 + AliTRDgeometry::CamHght()/2. + iplan * (AliTRDgeometry::Cheight() + AliTRDgeometry::Cspace());
370 }
17b26de4 371
372 //
373 // ----------------------------------------------------------------------------
ccb4315c 374 // The digitization parameters
17b26de4 375 // ----------------------------------------------------------------------------
376 //
377
ccb4315c 378 // Use drift time maps
379 fTimeStructOn = kTRUE;
380
381 // The Sampling Frequency. Default is 10MHz
382 fSamplingFrequency = 10.0;
383
384 // The drift velocity (in drift region) and the time structure of the
385 // drift cells. Default is 1.5 cm/mus
386 SetDriftVelocity(1.5);
387
388 // Additional time bins before and after the drift region.
389 // Default is to only sample the drift region
390 SetExpandTimeBin(0,0);
391
17b26de4 392 // The default parameter for the digitization
a328fff9 393 fGasGain = 4000.;
17b26de4 394 fChipGain = 12.4;
395 fNoise = 1000.;
396 fADCoutRange = 1023.; // 10-bit ADC
7e582e9f 397 fADCinRange = 2000.; // 2V input range
17b26de4 398 fADCthreshold = 1;
7e582e9f 399 fADCbaseline = 0;
17b26de4 400
17b26de4 401 // Diffusion on
ccb4315c 402 fDiffusionOn = kTRUE;
17b26de4 403
404 // E x B effects
ccb4315c 405 fExBOn = kTRUE;
17b26de4 406
407 // Propability for electron attachment
ccb4315c 408 fElAttachOn = kFALSE;
17b26de4 409 fElAttachProp = 0.0;
410
411 // The pad response function
ccb4315c 412 fPRFOn = kTRUE;
17b26de4 413
414 // The time response function
ccb4315c 415 fTRFOn = kTRUE;
17b26de4 416
417 // The cross talk
ccb4315c 418 fCTOn = kTRUE;
17b26de4 419
420 // The tail cancelation
ccb4315c 421 fTCOn = kTRUE;
17b26de4 422
423 // The number of exponentials
81a6930e 424 fTCnexp = 1;
17b26de4 425
b258a7c7 426 // The pad coupling factor
427 //fPadCoupling = 0.3;
682619b9 428 // Use 0.46 instead which reproduces better the test beam
b258a7c7 429 // data, even tough it is not understood why.
682619b9 430 fPadCoupling = 0.46;
17b26de4 431
432 // The time coupling factor (same number as for the TPC)
433 fTimeCoupling = 0.4;
434
a328fff9 435 // Distance of first Anode wire from first pad edge
436 fAnodeWireOffset = 0.25;
437
17b26de4 438 // The tilting angle for the readout pads
73ae7b59 439 SetTiltingAngle(2.0);
17b26de4 440
441 // The magnetic field strength in Tesla
81a6930e 442 Double_t x[3] = { 0.0, 0.0, 0.0 };
443 Double_t b[3];
444 gAlice->Field(x,b); // b[] is in kilo Gauss
445 fField = b[2] * 0.1; // Tesla
17b26de4 446
447 //
448 // ----------------------------------------------------------------------------
449 // The clusterization parameter
450 // ----------------------------------------------------------------------------
451 //
452
453 // The default parameter for the clustering
454 fClusMaxThresh = 3;
455 fClusSigThresh = 1;
456
457 // Use the LUT
ccb4315c 458 fLUTOn = kTRUE;
17b26de4 459
460 ReInit();
461
462}
463
464//_____________________________________________________________________________
465void AliTRDparameter::ReInit()
466{
467 //
468 // Reinitializes the parameter class after a change
469 //
470
ccb4315c 471 Float_t maxdrift = AliTRDgeometry::DrThick() + AliTRDgeometry::AmThick()/2.0;
a328fff9 472
ccb4315c 473 // Number of time bins:
474 fTimeMax = (Int_t) ( maxdrift / fDriftVelocity * fSamplingFrequency );
17b26de4 475
476 // The range and the binwidth for the sampled TRF
477 fTRFbin = 100;
478 // Start 0.2 mus before the signal
ccb4315c 479 fTRFlo = -0.2;
480 // End the maximum drift time after the signal
481 fTRFhi = 2.2;
482 //
17b26de4 483 fTRFwid = (fTRFhi - fTRFlo) / ((Float_t) fTRFbin);
484
485 // Transverse and longitudinal diffusion coefficients (Xe/CO2)
486 fDiffusionT = GetDiffusionT(fDriftVelocity,fField);
487 fDiffusionL = GetDiffusionL(fDriftVelocity,fField);
488
ccb4315c 489 // omega * tau = tan(Lorentz-angle)
17b26de4 490 fOmegaTau = GetOmegaTau(fDriftVelocity,fField);
491
492 // The Lorentz factor
493 if (fExBOn) {
494 fLorentzFactor = 1.0 / (1.0 + fOmegaTau*fOmegaTau);
495 }
496 else {
497 fLorentzFactor = 1.0;
498 }
499
500 // Create the sampled PRF
501 SamplePRF();
502
503 // Create the sampled TRF
504 SampleTRF();
505
506 // Create the LUT
507 FillLUT();
508
509}
510
511//_____________________________________________________________________________
afc51ac2 512void AliTRDparameter::SetNRowPad(Int_t p, Int_t c, Int_t npad)
17b26de4 513{
514 //
515 // Redefines the number of pads in raw direction for
516 // a given plane and chamber number
517 //
518
519 for (Int_t isect = 0; isect < AliTRDgeometry::Nsect(); isect++) {
520
521 fRowMax[p][c][isect] = npad;
522
523 fRowPadSize[p][c][isect] = (fGeo->GetChamberLength(p,c)
524 - 2.* AliTRDgeometry::RpadW())
525 / ((Float_t) npad);
526
527 }
528
529}
530
531//_____________________________________________________________________________
532void AliTRDparameter::SetNRowPad()
533{
534 //
535 // Defines the number of pads in row direction
536 //
537
538 Int_t isect;
539 Int_t icham;
540 Int_t iplan;
541
542 Int_t rowMax[kNplan][kNcham] = { { 16, 16, 12, 16, 16 }
543 , { 16, 16, 12, 16, 16 }
544 , { 16, 16, 12, 16, 16 }
545 , { 16, 16, 12, 16, 16 }
ad973dd4 546 , { 16, 16, 12, 16, 16 }
547 , { 16, 16, 12, 16, 16 } };
17b26de4 548
549 Float_t rpadW = AliTRDgeometry::RpadW();
550
551 for (isect = 0; isect < kNsect; isect++) {
552 for (icham = 0; icham < kNcham; icham++) {
553 for (iplan = 0; iplan < kNplan; iplan++) {
554
555 fRowMax[iplan][icham][isect] = rowMax[iplan][icham];
556
557 fRowPadSize[iplan][icham][isect] = (fGeo->GetChamberLength(iplan,icham)
558 - 2.*rpadW)
559 / ((Float_t) rowMax[iplan][icham]);
560
561 Float_t row0 = rpadW - fGeo->GetChamberLength(iplan,0)
562 - fGeo->GetChamberLength(iplan,1)
563 - fGeo->GetChamberLength(iplan,2) / 2.;
564 for (Int_t ic = 0; ic < icham; ic++) {
4bd34b78 565 row0 += fGeo->GetChamberLength(iplan,ic);
17b26de4 566 }
81a6930e 567
568 fRow0[iplan][icham][isect] = row0;
569 // For new chamber ordering
570 //fRow0[iplan][kNcham-icham-1][isect] = row0;
17b26de4 571
572 }
573 }
574 }
575
576}
577
578//_____________________________________________________________________________
afc51ac2 579void AliTRDparameter::SetColPadSize(Int_t p, Float_t s)
17b26de4 580{
581 //
582 // Redefines the pad size in column direction
583 //
584
585 Float_t cpadW = AliTRDgeometry::CpadW();
586
587 fColPadSize[p] = s;
588 fCol0[p] = - fGeo->GetChamberWidth(p)/2. + cpadW;
589 fColMax[p] = ((Int_t) ((fGeo->GetChamberWidth(p) - 2.*cpadW) / s));
590
591}
592
17b26de4 593//_____________________________________________________________________________
0a29d0f1 594Float_t AliTRDparameter::CrossTalk(Float_t time) const
17b26de4 595{
596 //
597 // Applies the pad-pad capacitive cross talk
598 //
599
600 Int_t iBin = ((Int_t) ((time - fTRFlo) / fTRFwid));
601 if ((iBin >= 0) && (iBin < fTRFbin)) {
602 return fCTsmp[iBin];
603 }
604 else {
605 return 0.0;
606 }
607
608}
609
a328fff9 610//_____________________________________________________________________________
ccb4315c 611void AliTRDparameter::PrintDriftVelocity()
612{
613 printf("<AliTRDparameter::PrintDriftVelocity> Driftvelocity = %.3f\n", fDriftVelocity);
614}
615
616//_____________________________________________________________________________
617Float_t AliTRDparameter::TimeStruct(Float_t dist, Float_t z) const
a328fff9 618{
619 //
620 // Applies the time structure of the drift cells (by C.Lippmann).
621 // The drift time of electrons to the anode wires depends on the
622 // distance to the wire (z) and on the position in the drift region.
682619b9 623 //
624 // input :
ccb4315c 625 // dist = radial distance from (cathode) pad plane [cm]
682619b9 626 // z = distance from anode wire (parallel to cathode planes) [cm]
a328fff9 627 //
ccb4315c 628 // output :
629 // tdrift = the drift time of an electron at the given position
630 //
631 // We interpolate between the drift time values at the two drift
632 // velocities fVDlo and fVDhi, being smaller and larger than
633 // fDriftVelocity. We use the two stored drift time maps fTimeStruct1
634 // and fTimeStruct2, calculated for the two mentioned drift velocities.
635 //
a328fff9 636
682619b9 637 // indices:
638 Int_t r1 = (Int_t)(10*dist);
a328fff9 639 Int_t r2 = r1+1;
ccb4315c 640 if (r1<0) r1 = 0;
641 if (r1>37) r1 = 37;
642 const Int_t z1 = (Int_t)(100*z/2.5);
643 const Int_t z2 = z1+1;
a328fff9 644
682619b9 645 if (r1<0 || r1>37 || z1<0 || z1>10) {
a328fff9 646 printf("<AliTRDparameter::TimeStruct> Warning. Indices out of range: ");
682619b9 647 printf("dist=%.2f, z=%.2f, r1=%d, z1=%d\n",dist,z,r1,z1);
a328fff9 648 }
649
ccb4315c 650 const Float_t y111 = fTimeStruct1[r1+38*z1];
651 const Float_t y221 = (r2 <= 37 && z2 <= 10) ? fTimeStruct1[r2+38*z2] : fTimeStruct1[37+38*10];
652 const Float_t y121 = (z2 <= 10) ? fTimeStruct1[r1+38*z2] : fTimeStruct1[r1+38*10];
653 const Float_t y211 = (r2 <= 37) ? fTimeStruct1[r2+38*z1] : fTimeStruct1[37+38*z1];
654
655 // 2D Interpolation, lower drift time map
656 const Float_t y11 = (y211-y111)*10*dist + y111 - (y211-y111)*r1;
657 const Float_t y21 = (y221-y121)*10*dist + y121 - (y221-y121)*r1;
682619b9 658
ccb4315c 659 const Float_t y112 = fTimeStruct2[r1+38*z1];
660 const Float_t y222 = (r2 <= 37 && z2 <= 10) ? fTimeStruct2[r2+38*z2] : fTimeStruct2[37+38*10];
661 const Float_t y122 = (z2 <= 10) ? fTimeStruct2[r1+38*z2] : fTimeStruct2[r1+38*10];
662 const Float_t y212 = (r2 <= 37) ? fTimeStruct2[r2+38*z1] : fTimeStruct2[37+38*z1];
682619b9 663
ccb4315c 664 // 2D Interpolation, larger drift time map
665 const Float_t y12 = (y212-y112)*10*dist + y112 - (y212-y112)*r1;
666 const Float_t y22 = (y222-y122)*10*dist + y122 - (y222-y122)*r1;
682619b9 667
ccb4315c 668 // dist now is the drift distance to the anode wires (negative if electrons are
682619b9 669 // between anode wire plane and cathode pad plane)
ccb4315c 670 dist -= AliTRDgeometry::AmThick()/2.0;
682619b9 671
ccb4315c 672 // Get the drift times for the drift velocities fVDlo and fVDhi
673 const Float_t tdrift1 =
674 ( TMath::Abs(dist)>0.005 || z>0.005 ) ? (y21-y11)*100*z/2.5+y11-(y21-y11)*z1 : 0.0;
675 const Float_t tdrift2 =
676 ( TMath::Abs(dist)>0.005 || z>0.005 ) ? (y22-y12)*100*z/2.5+y12-(y22-y12)*z1 : 0.0;
682619b9 677
ccb4315c 678 // 1D Interpolation between the values at fVDlo and fVDhi
679 Float_t a = (tdrift2 - tdrift1) / (fVDhi - fVDlo);
680 Float_t b = tdrift2 - a * fVDhi;
682619b9 681
682619b9 682
ccb4315c 683 //printf("(%.2f, %.2f): %f, %f -> %f\n",
684 // dist+AliTRDgeometry::AmThick()/2.0, z, tdrift1, tdrift2, a*fDriftVelocity+b);
a328fff9 685
ccb4315c 686 return a * fDriftVelocity + b;
a328fff9 687
688}
689
17b26de4 690//_____________________________________________________________________________
691Int_t AliTRDparameter::Diffusion(Float_t driftlength, Float_t *xyz)
692{
693 //
694 // Applies the diffusion smearing to the position of a single electron
695 //
696
697 Float_t driftSqrt = TMath::Sqrt(driftlength);
698 Float_t sigmaT = driftSqrt * fDiffusionT;
699 Float_t sigmaL = driftSqrt * fDiffusionL;
700 xyz[0] = gRandom->Gaus(xyz[0], sigmaL * fLorentzFactor);
701 xyz[1] = gRandom->Gaus(xyz[1], sigmaT * fLorentzFactor);
702 xyz[2] = gRandom->Gaus(xyz[2], sigmaT);
703
704 return 1;
705
706}
707
708//_____________________________________________________________________________
0a29d0f1 709Int_t AliTRDparameter::ExB(Float_t driftlength, Float_t *xyz) const
17b26de4 710{
711 //
712 // Applies E x B effects to the position of a single electron
713 //
714
715 xyz[0] = xyz[0];
716 xyz[1] = xyz[1] + fOmegaTau * driftlength;
717 xyz[2] = xyz[2];
718
719 return 1;
720
721}
722
723//_____________________________________________________________________________
724Int_t AliTRDparameter::PadResponse(Float_t signal, Float_t dist
0a29d0f1 725 , Int_t plane, Float_t *pad) const
17b26de4 726{
727 //
728 // Applies the pad response
729 //
730
731 const Int_t kNplan = AliTRDgeometry::kNplan;
732
733 Int_t iBin = ((Int_t) (( - dist - fPRFlo) / fPRFwid));
734 Int_t iOff = plane * fPRFbin;
735
736 Int_t iBin0 = iBin - fPRFpad + iOff;
737 Int_t iBin1 = iBin + iOff;
738 Int_t iBin2 = iBin + fPRFpad + iOff;
739
740 pad[0] = 0.0;
741 pad[1] = 0.0;
742 pad[2] = 0.0;
743 if ((iBin1 >= 0) && (iBin1 < (fPRFbin*kNplan))) {
744
745 if (iBin0 >= 0) {
746 pad[0] = signal * fPRFsmp[iBin0];
747 }
748 pad[1] = signal * fPRFsmp[iBin1];
749 if (iBin2 < (fPRFbin*kNplan)) {
750 pad[2] = signal * fPRFsmp[iBin2];
751 }
752
753 return 1;
754
755 }
756 else {
757
758 return 0;
759
760 }
761
762}
763
764//_____________________________________________________________________________
0a29d0f1 765Float_t AliTRDparameter::TimeResponse(Float_t time) const
17b26de4 766{
767 //
768 // Applies the preamp shaper time response
769 //
770
81a6930e 771 Int_t iBin = ((Int_t) ((time) / fTRFwid));
772 //Int_t iBin = ((Int_t) ((time - fTRFlo) / fTRFwid));
17b26de4 773 if ((iBin >= 0) && (iBin < fTRFbin)) {
774 return fTRFsmp[iBin];
775 }
776 else {
777 return 0.0;
778 }
779
780}
781
782//_____________________________________________________________________________
783Float_t AliTRDparameter::Col0Tilted(Float_t col0, Float_t rowOffset
784 , Int_t plane)
785{
786 //
787 // Calculates col0 for tilted pads
788 //
789
790 Float_t diff = fTiltingAngle * rowOffset;
a328fff9 791 return (col0 + TMath::Power(-1.0,(plane+1)) * diff);
17b26de4 792
793}
794
795//_____________________________________________________________________________
796void AliTRDparameter::SampleTRF()
797{
798 //
799 // Samples the time response function
800 //
801 // New TRF from Venelin Angelov, simulated with CADENCE
802 // Pad-ground capacitance = 25 pF
803 // Pad-pad cross talk capacitance = 6 pF
804 //
805
806 Int_t ipos1;
807 Int_t ipos2;
808 Float_t diff;
809
810 const Int_t kNpasa = 252;
811
812 Float_t time[kNpasa] = { -0.220000, -0.210000, -0.200000, -0.190000
813 , -0.180000, -0.170000, -0.160000, -0.150000
814 , -0.140000, -0.130000, -0.120000, -0.110000
815 , -0.100000, -0.090000, -0.080000, -0.070000
816 , -0.060000, -0.050000, -0.040000, -0.030000
817 , -0.020000, -0.010000, -0.000000, 0.010000
818 , 0.020000, 0.030000, 0.040000, 0.050000
819 , 0.060000, 0.070000, 0.080000, 0.090000
820 , 0.100000, 0.110000, 0.120000, 0.130000
821 , 0.140000, 0.150000, 0.160000, 0.170000
822 , 0.180000, 0.190000, 0.200000, 0.210000
823 , 0.220000, 0.230000, 0.240000, 0.250000
824 , 0.260000, 0.270000, 0.280000, 0.290000
825 , 0.300000, 0.310000, 0.320000, 0.330000
826 , 0.340000, 0.350000, 0.360000, 0.370000
827 , 0.380000, 0.390000, 0.400000, 0.410000
828 , 0.420000, 0.430000, 0.440000, 0.450000
829 , 0.460000, 0.470000, 0.480000, 0.490000
830 , 0.500000, 0.510000, 0.520000, 0.530000
831 , 0.540000, 0.550000, 0.560000, 0.570000
832 , 0.580000, 0.590000, 0.600000, 0.610000
833 , 0.620000, 0.630000, 0.640000, 0.650000
834 , 0.660000, 0.670000, 0.680000, 0.690000
835 , 0.700000, 0.710000, 0.720000, 0.730000
836 , 0.740000, 0.750000, 0.760000, 0.770000
837 , 0.780000, 0.790000, 0.800000, 0.810000
838 , 0.820000, 0.830000, 0.840000, 0.850000
839 , 0.860000, 0.870000, 0.880000, 0.890000
840 , 0.900000, 0.910000, 0.920000, 0.930000
841 , 0.940000, 0.950000, 0.960000, 0.970000
842 , 0.980000, 0.990000, 1.000000, 1.010000
843 , 1.020000, 1.030000, 1.040000, 1.050000
844 , 1.060000, 1.070000, 1.080000, 1.090000
845 , 1.100000, 1.110000, 1.120000, 1.130000
846 , 1.140000, 1.150000, 1.160000, 1.170000
847 , 1.180000, 1.190000, 1.200000, 1.210000
848 , 1.220000, 1.230000, 1.240000, 1.250000
849 , 1.260000, 1.270000, 1.280000, 1.290000
850 , 1.300000, 1.310000, 1.320000, 1.330000
851 , 1.340000, 1.350000, 1.360000, 1.370000
852 , 1.380000, 1.390000, 1.400000, 1.410000
853 , 1.420000, 1.430000, 1.440000, 1.450000
854 , 1.460000, 1.470000, 1.480000, 1.490000
855 , 1.500000, 1.510000, 1.520000, 1.530000
856 , 1.540000, 1.550000, 1.560000, 1.570000
857 , 1.580000, 1.590000, 1.600000, 1.610000
858 , 1.620000, 1.630000, 1.640000, 1.650000
859 , 1.660000, 1.670000, 1.680000, 1.690000
860 , 1.700000, 1.710000, 1.720000, 1.730000
861 , 1.740000, 1.750000, 1.760000, 1.770000
862 , 1.780000, 1.790000, 1.800000, 1.810000
863 , 1.820000, 1.830000, 1.840000, 1.850000
864 , 1.860000, 1.870000, 1.880000, 1.890000
865 , 1.900000, 1.910000, 1.920000, 1.930000
866 , 1.940000, 1.950000, 1.960000, 1.970000
867 , 1.980000, 1.990000, 2.000000, 2.010000
868 , 2.020000, 2.030000, 2.040000, 2.050000
869 , 2.060000, 2.070000, 2.080000, 2.090000
870 , 2.100000, 2.110000, 2.120000, 2.130000
871 , 2.140000, 2.150000, 2.160000, 2.170000
872 , 2.180000, 2.190000, 2.200000, 2.210000
873 , 2.220000, 2.230000, 2.240000, 2.250000
874 , 2.260000, 2.270000, 2.280000, 2.290000 };
875
876 Float_t signal[kNpasa] = { 0.000000, 0.000000, 0.000000, 0.000000
877 , 0.000000, 0.000000, 0.000000, 0.000396
878 , 0.005096, 0.022877, 0.061891, 0.126614
879 , 0.215798, 0.324406, 0.444507, 0.566817
880 , 0.683465, 0.787089, 0.873159, 0.937146
881 , 0.979049, 0.999434, 1.000000, 0.983579
882 , 0.954134, 0.913364, 0.866365, 0.813703
883 , 0.759910, 0.706116, 0.653454, 0.603624
884 , 0.556625, 0.514156, 0.475085, 0.439977
885 , 0.408834, 0.380578, 0.355549, 0.333352
886 , 0.313647, 0.296093, 0.280351, 0.266195
887 , 0.253397, 0.241789, 0.231257, 0.221574
888 , 0.212627, 0.204417, 0.196772, 0.189581
889 , 0.182956, 0.176784, 0.171008, 0.165515
890 , 0.160419, 0.155606, 0.151076, 0.146716
891 , 0.142639, 0.138845, 0.135221, 0.131767
892 , 0.128482, 0.125368, 0.122424, 0.119592
893 , 0.116931, 0.114326, 0.111891, 0.109513
894 , 0.107248, 0.105096, 0.103058, 0.101019
895 , 0.099151, 0.097282, 0.095527, 0.093715
896 , 0.092129, 0.090544, 0.088958, 0.087429
897 , 0.086014, 0.084598, 0.083239, 0.081880
898 , 0.080634, 0.079388, 0.078143, 0.077010
899 , 0.075878, 0.074745, 0.073669, 0.072593
900 , 0.071574, 0.070612, 0.069649, 0.068686
901 , 0.067780, 0.066874, 0.066025, 0.065176
902 , 0.064326, 0.063533, 0.062684, 0.061948
903 , 0.061212, 0.060419, 0.059740, 0.059003
904 , 0.058324, 0.057644, 0.057022, 0.056342
905 , 0.055663, 0.055096, 0.054473, 0.053851
906 , 0.053284, 0.052718, 0.052152, 0.051585
907 , 0.051019, 0.050566, 0.050000, 0.049490
908 , 0.048981, 0.048528, 0.048018, 0.047508
909 , 0.047055, 0.046602, 0.046149, 0.045696
910 , 0.045300, 0.044904, 0.044451, 0.044054
911 , 0.043658, 0.043205, 0.042865, 0.042469
912 , 0.042072, 0.041733, 0.041336, 0.040997
913 , 0.040657, 0.040260, 0.039921, 0.039581
914 , 0.039241, 0.038958, 0.038618, 0.038335
915 , 0.037995, 0.037656, 0.037373, 0.037089
916 , 0.036806, 0.036467, 0.036183, 0.035900
917 , 0.035617, 0.035334, 0.035108, 0.034824
918 , 0.034541, 0.034315, 0.034032, 0.033805
919 , 0.033522, 0.033296, 0.033069, 0.032786
920 , 0.032559, 0.032333, 0.032106, 0.031880
921 , 0.031653, 0.031427, 0.031200, 0.030974
922 , 0.030804, 0.030578, 0.030351, 0.030125
923 , 0.029955, 0.029785, 0.029558, 0.029332
924 , 0.029162, 0.028992, 0.028766, 0.028596
925 , 0.028426, 0.028199, 0.028086, 0.027860
926 , 0.027746, 0.027633, 0.027463, 0.027293
927 , 0.027180, 0.027067, 0.026954, 0.026954
928 , 0.026840, 0.026727, 0.026727, 0.026614
929 , 0.026614, 0.026614, 0.026557, 0.026501
930 , 0.026501, 0.026501, 0.026501, 0.026501
931 , 0.026501, 0.026501, 0.026501, 0.026387
932 , 0.026387, 0.026387, 0.026387, 0.026387
933 , 0.026387, 0.026387, 0.026387, 0.026387
934 , 0.026387, 0.026387, 0.026387, 0.026387
935 , 0.026387, 0.026274, 0.026274, 0.026274
936 , 0.026274, 0.026274, 0.026274, 0.026274
937 , 0.026274, 0.026274, 0.026274, 0.026274
938 , 0.026274, 0.026274, 0.026274, 0.026161 };
939
940 Float_t xtalk[kNpasa] = { 0.000000, 0.000000, 0.000000, 0.000000
941 , 0.000000, 0.000000, 0.000000, 0.000113
942 , 0.000793, 0.003058, 0.007305, 0.013194
943 , 0.019706, 0.025821, 0.030634, 0.033465
944 , 0.034145, 0.032729, 0.029615, 0.025198
945 , 0.019989, 0.014496, 0.009003, 0.003964
946 , -0.000510, -0.004190, -0.007191, -0.009400
947 , -0.010872, -0.011835, -0.012288, -0.012288
948 , -0.012005, -0.011495, -0.010872, -0.010136
949 , -0.009343, -0.008607, -0.007871, -0.007191
950 , -0.006512, -0.005946, -0.005379, -0.004926
951 , -0.004473, -0.004077, -0.003737, -0.003398
952 , -0.003114, -0.002831, -0.002605, -0.002378
953 , -0.002208, -0.002039, -0.001869, -0.001699
954 , -0.001585, -0.001472, -0.001359, -0.001246
955 , -0.001132, -0.001019, -0.001019, -0.000906
956 , -0.000906, -0.000793, -0.000793, -0.000680
957 , -0.000680, -0.000680, -0.000566, -0.000566
958 , -0.000566, -0.000566, -0.000453, -0.000453
959 , -0.000453, -0.000453, -0.000453, -0.000453
960 , -0.000340, -0.000340, -0.000340, -0.000340
961 , -0.000340, -0.000340, -0.000340, -0.000340
962 , -0.000340, -0.000340, -0.000340, -0.000340
963 , -0.000340, -0.000227, -0.000227, -0.000227
964 , -0.000227, -0.000227, -0.000227, -0.000227
965 , -0.000227, -0.000227, -0.000227, -0.000227
966 , -0.000227, -0.000227, -0.000227, -0.000227
967 , -0.000227, -0.000227, -0.000227, -0.000227
968 , -0.000227, -0.000227, -0.000227, -0.000227
969 , -0.000227, -0.000227, -0.000227, -0.000227
970 , -0.000227, -0.000227, -0.000227, -0.000227
971 , -0.000227, -0.000227, -0.000227, -0.000227
972 , -0.000227, -0.000227, -0.000227, -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.000113
977 , -0.000113, -0.000113, -0.000113, -0.000113
978 , -0.000113, -0.000113, -0.000113, -0.000113
979 , -0.000113, -0.000113, -0.000113, -0.000113
980 , -0.000113, -0.000113, -0.000113, -0.000113
981 , -0.000113, -0.000113, -0.000113, -0.000113
982 , -0.000113, -0.000113, -0.000113, -0.000113
983 , -0.000113, -0.000113, -0.000113, -0.000113
984 , -0.000113, -0.000113, -0.000113, -0.000113
985 , -0.000113, -0.000113, -0.000113, -0.000113
986 , -0.000113, -0.000113, -0.000113, -0.000113
987 , -0.000113, -0.000113, -0.000113, -0.000113
988 , -0.000113, -0.000113, -0.000113, -0.000113
989 , -0.000113, -0.000113, -0.000113, -0.000113
990 , -0.000113, -0.000113, -0.000113, -0.000113
991 , -0.000113, -0.000113, -0.000113, -0.000113
992 , -0.000113, -0.000113, -0.000113, 0.000000
993 , 0.000000, 0.000000, 0.000000, 0.000000
994 , 0.000000, 0.000000, 0.000000, 0.000000
995 , 0.000000, 0.000000, 0.000000, 0.000000
996 , 0.000000, 0.000000, 0.000000, 0.000000
997 , 0.000000, 0.000000, 0.000000, 0.000000
998 , 0.000000, 0.000000, 0.000000, 0.000000
999 , 0.000000, 0.000000, 0.000000, 0.000000
1000 , 0.000000, 0.000000, 0.000000, 0.000000
1001 , 0.000000, 0.000000, 0.000000, 0.000000
1002 , 0.000000, 0.000000, 0.000000, 0.000000 };
1003
1004 // increase CrossTalk to measurements
1005 for (Int_t ipasa = 0; ipasa < kNpasa; ipasa++) {
1006 xtalk[ipasa] *= 1.75;
1007 }
1008
1009 if (fTRFsmp) delete [] fTRFsmp;
1010 fTRFsmp = new Float_t[fTRFbin];
1011 if (fCTsmp) delete [] fCTsmp;
1012 fCTsmp = new Float_t[fTRFbin];
1013
ccb4315c 1014 Float_t loTRF = TMath::Max(fTRFlo, time[0]);
1015 Float_t hiTRF = TMath::Min(fTRFhi, time[kNpasa-1]);
17b26de4 1016 Float_t binWidth = (hiTRF - loTRF) / ((Float_t) fTRFbin);
1017
1018 // Take the linear interpolation
1019 for (Int_t iBin = 0; iBin < fTRFbin; iBin++) {
1020
1021 Float_t bin = (((Float_t) iBin) + 0.5) * binWidth + loTRF;
1022 ipos1 = ipos2 = 0;
1023 diff = 0;
1024 do {
1025 diff = bin - time[ipos2++];
1026 } while (diff > 0);
1027 ipos2--;
1028 if (ipos2 >= kNpasa) ipos2 = kNpasa - 1;
1029 ipos1 = ipos2 - 1;
1030
1031 fTRFsmp[iBin] = signal[ipos2]
1032 + diff * (signal[ipos2] - signal[ipos1])
1033 / ( time[ipos2] - time[ipos1]);
1034
1035 fCTsmp[iBin] = xtalk[ipos2]
1036 + diff * (xtalk[ipos2] - xtalk[ipos1])
1037 / ( time[ipos2] - time[ipos1]);
1038
1039 }
1040
1041}
1042
a328fff9 1043//_____________________________________________________________________________
1044void AliTRDparameter::SampleTimeStruct()
1045{
ccb4315c 1046
a328fff9 1047 //
1048 // Samples the timing structure of a drift cell
1049 // Drift Time data calculated with Garfield (by C.Lippmann)
1050 //
1051
ccb4315c 1052 // drift time maps are saved for some drift velocity values (in drift region):
1053 Float_t fVDsmp[8];
1054 fVDsmp[0] = 1.032;
1055 fVDsmp[1] = 1.158;
1056 fVDsmp[2] = 1.299;
1057 fVDsmp[3] = 1.450;
1058 fVDsmp[4] = 1.610;
1059 fVDsmp[5] = 1.783;
1060 fVDsmp[6] = 1.959;
1061 fVDsmp[7] = 2.134;
1062
1063 if ( fDriftVelocity < fVDsmp[0] ) {
1064 printf("<AliTRDparameter::SampleTimeStruct> !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
1065 printf("<AliTRDparameter::SampleTimeStruct> Drift Velocity too small (%.3f<%.3f)\n", fDriftVelocity, fVDsmp[0]);
1066 printf("<AliTRDparameter::SampleTimeStruct> !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
1067 fDriftVelocity = fVDsmp[0];
1068 } else if ( fDriftVelocity > fVDsmp[7] ) {
1069 printf("<AliTRDparameter::SampleTimeStruct> !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
1070 printf("<AliTRDparameter::SampleTimeStruct> Drift Velocity too large (%.3f>%.3f)\n", fDriftVelocity,fVDsmp[6]);
1071 printf("<AliTRDparameter::SampleTimeStruct> !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
1072 fDriftVelocity = fVDsmp[7];
1073 }
1074
a328fff9 1075 const Int_t ktimebin = 38;
682619b9 1076 const Int_t kZbin = 11;
ccb4315c 1077
1078 // Garfield simulation at UD = -1500V; vd = 1.032cm/microsec, <driftfield> = 525V/cm
1079 Float_t time1500[ktimebin][kZbin] =
1080 {{0.09170, 0.09205, 0.09306, 0.09475, 0.09716, 0.10035,
1081 0.10445, 0.10993, 0.11838, 0.13986, 0.37858},
1082 {0.06588, 0.06626, 0.06739, 0.06926, 0.07186, 0.07524,
1083 0.07951, 0.08515, 0.09381, 0.11601, 0.35673},
1084 {0.03946, 0.04003, 0.04171, 0.04435, 0.04780, 0.05193,
1085 0.05680, 0.06306, 0.07290, 0.09873, 0.34748},
1086 {0.01151, 0.01283, 0.01718, 0.02282, 0.02880, 0.03479,
1087 0.04098, 0.04910, 0.06413, 0.10567, 0.36897},
1088 {0.01116, 0.01290, 0.01721, 0.02299, 0.02863, 0.03447,
1089 0.04074, 0.04984, 0.06839, 0.11625, 0.37277},
1090 {0.03919, 0.03974, 0.04131, 0.04380, 0.04703, 0.05102,
1091 0.05602, 0.06309, 0.07651, 0.10938, 0.36838},
1092 {0.06493, 0.06560, 0.06640, 0.06802, 0.07051, 0.07392,
1093 0.07853, 0.08510, 0.09690, 0.12621, 0.38058},
1094 {0.09174, 0.09186, 0.09225, 0.09303, 0.09477, 0.00000,
1095 0.11205, 0.11952, 0.13461, 0.16984, 0.43017},
1096 {0.14356, 0.14494, 0.14959, 0.16002, 0.18328, 0.27981,
1097 0.22785, 0.21240, 0.21948, 0.25965, 0.52392},
1098 {0.23120, 0.23366, 0.24046, 0.25422, 0.28071, 0.36914,
1099 0.32999, 0.31208, 0.31772, 0.35804, 0.62249},
1100 {0.32686, 0.32916, 0.33646, 0.35053, 0.37710, 0.46292,
1101 0.42773, 0.40948, 0.41497, 0.45527, 0.71955},
1102 {0.42353, 0.42583, 0.43317, 0.44727, 0.47380, 0.55884,
1103 0.52479, 0.50650, 0.51194, 0.55225, 0.81658},
1104 {0.52038, 0.52271, 0.53000, 0.54415, 0.57064, 0.65545,
1105 0.62172, 0.60341, 0.60885, 0.64915, 0.91339},
1106 {0.61724, 0.61953, 0.62694, 0.64098, 0.66756, 0.75226,
1107 0.71862, 0.70030, 0.70575, 0.74604, 1.01035},
1108 {0.71460, 0.71678, 0.72376, 0.73786, 0.76447, 0.84913,
1109 0.81551, 0.79720, 0.80264, 0.84292, 1.10723},
1110 {0.81101, 0.81334, 0.82066, 0.83475, 0.86127, 0.94599,
1111 0.91240, 0.89408, 0.89952, 0.93981, 1.20409},
1112 {0.90788, 0.91023, 0.91752, 0.93163, 0.95815, 1.04293,
1113 1.00929, 0.99096, 0.99640, 1.03669, 1.30106},
1114 {1.00477, 1.00707, 1.01449, 1.02852, 1.05504, 1.13976,
1115 1.10617, 1.08784, 1.09329, 1.13358, 1.39796},
1116 {1.10166, 1.10397, 1.11130, 1.12541, 1.15257, 1.23672,
1117 1.20307, 1.18472, 1.19018, 1.23046, 1.49486},
1118 {1.19854, 1.20084, 1.20818, 1.22235, 1.24885, 1.33355,
1119 1.29992, 1.28156, 1.28707, 1.32735, 1.59177},
1120 {1.29544, 1.29780, 1.30507, 1.31917, 1.34575, 1.43073,
1121 1.39681, 1.37851, 1.38396, 1.42377, 1.68854},
1122 {1.39236, 1.39462, 1.40205, 1.41607, 1.44259, 1.52745,
1123 1.49368, 1.47541, 1.48083, 1.52112, 1.78546},
1124 {1.49314, 1.49149, 1.49885, 1.51297, 1.53949, 1.62420,
1125 1.59016, 1.57231, 1.57772, 1.61800, 1.88048},
1126 {1.58610, 1.58839, 1.59572, 1.60983, 1.63635, 1.72109,
1127 1.68651, 1.66921, 1.67463, 1.71489, 1.97918},
1128 {1.68400, 1.68529, 1.69261, 1.70671, 1.73331, 1.81830,
1129 1.78341, 1.76605, 1.77150, 1.81179, 2.07608},
1130 {1.77991, 1.78215, 1.78952, 1.80385, 1.83014, 1.91486,
1131 1.88128, 1.86215, 1.86837, 1.90865, 2.17304},
1132 {1.87674, 1.87904, 1.88647, 1.90052, 1.92712, 2.01173,
1133 1.97812, 1.95905, 1.96527, 2.00710, 2.26979},
1134 {1.97369, 1.97594, 1.98326, 1.99869, 2.02442, 2.10865,
1135 2.07501, 2.05666, 2.06214, 2.10243, 2.36669},
1136 {2.07052, 2.07281, 2.08016, 2.09425, 2.12132, 2.20555,
1137 2.17182, 2.15341, 2.15904, 2.19933, 2.46363},
1138 {2.16742, 2.16971, 2.17707, 2.19114, 2.21766, 2.30240,
1139 2.26877, 2.25015, 2.25573, 2.29586, 2.56060},
1140 {2.26423, 2.26659, 2.27396, 2.28803, 2.31456, 2.40828,
1141 2.36567, 2.34705, 2.35282, 2.39765, 2.65744},
1142 {2.36153, 2.36349, 2.37330, 2.38501, 2.41159, 2.49940,
1143 2.46257, 2.44420, 2.44843, 2.48987, 2.75431},
1144 {2.46558, 2.46035, 2.46822, 2.48181, 2.50849, 2.59630,
1145 2.55947, 2.54112, 2.54513, 2.58677, 2.85094},
1146 {2.56248, 2.55723, 2.56486, 2.57871, 2.60520, 2.68998,
1147 2.65626, 2.63790, 2.64316, 2.68360, 2.94813},
1148 {2.65178, 2.65441, 2.66153, 2.67556, 2.70210, 2.78687,
1149 2.75319, 2.73463, 2.74032, 2.78060, 3.04503},
1150 {2.74868, 2.75131, 2.75870, 2.77245, 2.79385, 2.88700,
1151 2.85009, 2.83177, 2.83723, 2.87750, 3.14193},
1152 {2.84574, 2.84789, 2.85560, 2.86935, 2.89075, 2.98060,
1153 2.94576, 2.92868, 2.93356, 2.97436, 3.23868},
1154 {2.94239, 2.94469, 2.95221, 2.96625, 2.99345, 3.07747,
1155 3.04266, 3.02545, 3.03051, 3.07118, 3.33555}};
1156
1157 // Garfield simulation at UD = -1600V; vd = 1.158cm/microsec, <driftfield> = 558V/cm
1158 Float_t time1600[ktimebin][kZbin] =
1159 {{0.09169, 0.09203, 0.09305, 0.09473, 0.09714, 0.10032,
1160 0.10441, 0.10990, 0.11835, 0.13986, 0.37845},
1161 {0.06589, 0.06626, 0.06738, 0.06924, 0.07184, 0.07521,
1162 0.07947, 0.08512, 0.09379, 0.11603, 0.35648},
1163 {0.03947, 0.04003, 0.04171, 0.04434, 0.04778, 0.05190,
1164 0.05678, 0.06304, 0.07292, 0.09876, 0.34759},
1165 {0.01111, 0.01281, 0.01718, 0.02281, 0.02879, 0.03477,
1166 0.04097, 0.04910, 0.06415, 0.10573, 0.36896},
1167 {0.01120, 0.01311, 0.01721, 0.02279, 0.02862, 0.03446,
1168 0.04074, 0.04981, 0.06825, 0.11595, 0.37255},
1169 {0.03919, 0.03980, 0.04132, 0.04380, 0.04704, 0.05102,
1170 0.05602, 0.06302, 0.07633, 0.10896, 0.36743},
1171 {0.06531, 0.06528, 0.06631, 0.06805, 0.07053, 0.07392,
1172 0.07853, 0.08505, 0.09669, 0.12578, 0.37967},
1173 {0.09157, 0.09171, 0.09216, 0.09301, 0.09475, 0.00000,
1174 0.11152, 0.11879, 0.13352, 0.16802, 0.42750},
1175 {0.13977, 0.14095, 0.14509, 0.15433, 0.17534, 0.26406,
1176 0.21660, 0.20345, 0.21113, 0.25067, 0.51434},
1177 {0.21816, 0.22041, 0.22631, 0.23850, 0.26208, 0.34340,
1178 0.30755, 0.29237, 0.29878, 0.33863, 0.60258},
1179 {0.30344, 0.30547, 0.31241, 0.32444, 0.34809, 0.42696,
1180 0.39464, 0.37919, 0.38546, 0.42530, 0.68926},
1181 {0.38969, 0.39164, 0.39810, 0.41059, 0.43441, 0.51246,
1182 0.48112, 0.46562, 0.47191, 0.51172, 0.77558},
1183 {0.47592, 0.47799, 0.48442, 0.49689, 0.52061, 0.59855,
1184 0.56752, 0.55201, 0.55826, 0.59808, 0.86202},
1185 {0.56226, 0.56428, 0.57074, 0.58324, 0.60696, 0.68483,
1186 0.65388, 0.63837, 0.64461, 0.68445, 0.94830},
1187 {0.64881, 0.65063, 0.65709, 0.66958, 0.69331, 0.77117,
1188 0.74023, 0.72472, 0.73098, 0.77079, 1.03486},
1189 {0.73506, 0.73698, 0.74344, 0.75596, 0.77964, 0.85751,
1190 0.82658, 0.81107, 0.81731, 0.85712, 1.12106},
1191 {0.82132, 0.82333, 0.82979, 0.84228, 0.86608, 0.94386,
1192 0.91293, 0.89742, 0.90367, 0.94335, 1.20737},
1193 {0.90767, 0.90968, 0.91614, 0.92863, 0.95236, 1.03021,
1194 0.99928, 0.98377, 0.99001, 1.02984, 1.29371},
1195 {0.99410, 0.99602, 1.00257, 1.01498, 1.03869, 1.11720,
1196 1.08563, 1.07011, 1.07637, 1.11621, 1.37873},
1197 {1.08036, 1.08240, 1.08884, 1.10138, 1.12504, 1.20301,
1198 1.17198, 1.15647, 1.16272, 1.20255, 1.46651},
1199 {1.16670, 1.16872, 1.17525, 1.18783, 1.21139, 1.28934,
1200 1.25833, 1.24281, 1.24909, 1.28889, 1.55275},
1201 {1.25307, 1.25510, 1.26153, 1.27404, 1.29773, 1.37584,
1202 1.34469, 1.32916, 1.33536, 1.37524, 1.63915},
1203 {1.33942, 1.34146, 1.34788, 1.36040, 1.38410, 1.46438,
1204 1.43105, 1.41537, 1.42176, 1.46158, 1.72538},
1205 {1.42579, 1.42782, 1.43458, 1.44674, 1.47043, 1.55085,
1206 1.51675, 1.50168, 1.50810, 1.54793, 1.81174},
1207 {1.51207, 1.51454, 1.52060, 1.53307, 1.55684, 1.63478,
1208 1.60336, 1.58820, 1.59446, 1.63414, 1.89814},
1209 {1.59856, 1.60047, 1.60693, 1.61942, 1.64317, 1.72257,
1210 1.69008, 1.67454, 1.68080, 1.72063, 1.98433},
1211 {1.68481, 1.68682, 1.69330, 1.70584, 1.72949, 1.80752,
1212 1.77643, 1.76089, 1.76716, 1.80692, 2.07069},
1213 {1.77117, 1.77319, 1.77969, 1.79260, 1.81583, 1.89376,
1214 1.86226, 1.84720, 1.85355, 1.89256, 2.15343},
1215 {1.85748, 1.85967, 1.86605, 1.87848, 1.90222, 1.98010,
1216 1.94913, 1.93271, 1.93981, 1.97968, 2.24355},
1217 {1.94386, 1.94587, 1.95233, 1.96484, 1.98854, 2.06646,
1218 2.03542, 2.01755, 2.02617, 2.06604, 2.32993},
1219 {2.03022, 2.03230, 2.03868, 2.05134, 2.07488, 2.15367,
1220 2.12178, 2.10391, 2.11252, 2.15432, 2.41623},
1221 {2.11656, 2.11857, 2.12505, 2.13772, 2.16147, 2.23919,
1222 2.20817, 2.19265, 2.20744, 2.23872, 2.49996},
1223 {2.20291, 2.20611, 2.21137, 2.22387, 2.24758, 2.32563,
1224 2.29450, 2.27901, 2.28525, 2.32507, 2.58897},
1225 {2.28922, 2.29172, 2.29774, 2.31345, 2.33400, 2.41287,
1226 2.38086, 2.36535, 2.37160, 2.40869, 2.67113},
1227 {2.37572, 2.37764, 2.38410, 2.39803, 2.42046, 2.49817,
1228 2.46721, 2.45171, 2.45794, 2.49505, 2.76061},
1229 {2.46190, 2.46396, 2.47043, 2.48340, 2.50665, 2.58453,
1230 2.55357, 2.53728, 2.54430, 2.58407, 2.84816},
1231 {2.54833, 2.55032, 2.55679, 2.56976, 2.59312, 2.67103,
1232 2.63993, 2.62364, 2.63062, 2.67040, 2.93444},
1233 {2.63456, 2.63660, 2.64304, 2.65555, 2.67938, 2.75739,
1234 2.72629, 2.71064, 2.71688, 2.75671, 3.01886}};
682619b9 1235
ccb4315c 1236 // Garfield simulation at UD = -1700V; vd = 1.299cm/microsec, <driftfield> = 590V/cm
1237 Float_t time1700[ktimebin][kZbin] =
1238 {{0.09167, 0.09201, 0.09302, 0.09471, 0.09712, 0.10029,
1239 0.10438, 0.10986, 0.11832, 0.13986, 0.37824},
1240 {0.06591, 0.06626, 0.06736, 0.06923, 0.07183, 0.07519,
1241 0.07944, 0.08511, 0.09378, 0.11603, 0.35625},
1242 {0.03946, 0.04003, 0.04170, 0.04433, 0.04777, 0.05189,
1243 0.05676, 0.06301, 0.07291, 0.09880, 0.34724},
1244 {0.01110, 0.01281, 0.01718, 0.02280, 0.02879, 0.03476,
1245 0.04096, 0.04910, 0.06417, 0.10582, 0.36861},
1246 {0.01115, 0.01294, 0.01721, 0.02276, 0.02862, 0.03447,
1247 0.04074, 0.04980, 0.06812, 0.11565, 0.37231},
1248 {0.03920, 0.03974, 0.04133, 0.04381, 0.04706, 0.05105,
1249 0.05603, 0.06299, 0.07618, 0.10860, 0.36646},
1250 {0.06498, 0.06529, 0.06634, 0.06808, 0.07055, 0.07395,
1251 0.07852, 0.08500, 0.09650, 0.12532, 0.37850},
1252 {0.09143, 0.09159, 0.09207, 0.09297, 0.09473, 0.00000,
1253 0.11102, 0.11809, 0.13245, 0.16627, 0.42496},
1254 {0.13646, 0.13750, 0.14112, 0.14926, 0.16806, 0.24960,
1255 0.20627, 0.19536, 0.20366, 0.24256, 0.50557},
1256 {0.20678, 0.20848, 0.21384, 0.22450, 0.24552, 0.32018,
1257 0.28740, 0.27477, 0.28196, 0.32128, 0.58475},
1258 {0.28287, 0.28461, 0.29020, 0.30108, 0.32224, 0.39467,
1259 0.36500, 0.35217, 0.35926, 0.39860, 0.66194},
1260 {0.35972, 0.36145, 0.36713, 0.37797, 0.39912, 0.47091,
1261 0.44212, 0.42925, 0.43632, 0.47563, 0.73892},
1262 {0.43667, 0.43841, 0.44413, 0.45494, 0.47607, 0.54780,
1263 0.51912, 0.50627, 0.51334, 0.55254, 0.81595},
1264 {0.51365, 0.51540, 0.52101, 0.53193, 0.55305, 0.62463,
1265 0.59617, 0.58328, 0.59035, 0.62965, 0.89303},
1266 {0.59064, 0.59240, 0.59801, 0.60893, 0.63009, 0.70169,
1267 0.67317, 0.66028, 0.66735, 0.70666, 0.96995},
1268 {0.66765, 0.66939, 0.67501, 0.68592, 0.70724, 0.77863,
1269 0.75016, 0.73728, 0.74435, 0.78366, 1.04696},
1270 {0.74464, 0.74636, 0.75200, 0.76293, 0.78405, 0.85561,
1271 0.82716, 0.81427, 0.82133, 0.86064, 1.12396},
1272 {0.82165, 0.82340, 0.82902, 0.83991, 0.86104, 0.93266,
1273 0.90414, 0.89128, 0.89834, 0.93763, 1.20100},
1274 {0.89863, 0.90042, 0.90659, 0.91705, 0.93805, 1.00960,
1275 0.98115, 0.96825, 0.97533, 1.01462, 1.27801},
1276 {0.97563, 0.97740, 0.98310, 0.99391, 1.01504, 1.08659,
1277 1.05814, 1.04526, 1.05233, 1.09163, 1.35503},
1278 {1.05276, 1.05451, 1.06002, 1.07090, 1.09099, 1.16357,
1279 1.13516, 1.12225, 1.12933, 1.16863, 1.43195},
1280 {1.12977, 1.13138, 1.13700, 1.14792, 1.16797, 1.24061,
1281 1.21212, 1.19926, 1.20626, 1.24554, 1.50900},
1282 {1.20664, 1.20839, 1.21400, 1.22490, 1.24606, 1.31772,
1283 1.28914, 1.27382, 1.28329, 1.32262, 1.58550},
1284 {1.28367, 1.28541, 1.29099, 1.30189, 1.32312, 1.39460,
1285 1.36612, 1.34924, 1.36030, 1.39961, 1.66310},
1286 {1.36064, 1.36249, 1.36799, 1.37896, 1.40004, 1.48030,
1287 1.44314, 1.43032, 1.43731, 1.47659, 1.73442},
1288 {1.43762, 1.43937, 1.44497, 1.45618, 1.47704, 1.54932,
1289 1.52012, 1.50725, 1.51430, 1.55357, 1.81708},
1290 {1.51462, 1.51937, 1.52203, 1.53316, 1.55403, 1.62572,
1291 1.59713, 1.58424, 1.59128, 1.63061, 1.89406},
1292 {1.59162, 1.59338, 1.59947, 1.60989, 1.63103, 1.70270,
1293 1.67411, 1.66124, 1.66799, 1.70759, 1.97103},
1294 {1.66874, 1.67037, 1.67597, 1.68687, 1.70814, 1.77969,
1295 1.75112, 1.73806, 1.74530, 1.78457, 2.04794},
1296 {1.74693, 1.74749, 1.75297, 1.76476, 1.78500, 1.85667,
1297 1.82811, 1.81504, 1.82101, 1.86161, 2.12492},
1298 {1.82260, 1.82437, 1.82995, 1.84174, 1.86202, 1.93372,
1299 1.90509, 1.89202, 1.89930, 1.93859, 2.20189},
1300 {1.89964, 1.90135, 1.90693, 1.91789, 1.93952, 2.01080,
1301 1.98207, 1.96921, 1.97628, 2.01384, 2.27887},
1302 {1.97662, 1.97917, 1.98611, 1.99487, 2.01601, 2.08778,
1303 2.05846, 2.04623, 2.05330, 2.09244, 2.35585},
1304 {2.05359, 2.05615, 2.06309, 2.07187, 2.09867, 2.16459,
1305 2.13610, 2.12322, 2.13029, 2.16942, 2.43199},
1306 {2.13063, 2.13233, 2.13795, 2.14886, 2.17008, 2.24199,
1307 2.21310, 2.20020, 2.20727, 2.24659, 2.50983},
1308 {2.20761, 2.20931, 2.21955, 2.22624, 2.24708, 2.32147,
1309 2.29009, 2.27725, 2.28276, 2.32359, 2.58680},
1310 {2.28459, 2.29108, 2.29202, 2.30286, 2.32007, 2.39559,
1311 2.36683, 2.35422, 2.36119, 2.40058, 2.66081},
1312 {2.36153, 2.36806, 2.36889, 2.37985, 2.40092, 2.47828,
1313 2.44381, 2.43099, 2.43819, 2.47750, 2.73779}};
1314
1315 // Garfield simulation at UD = -1800V; vd = 1.450cm/microsec, <driftfield> = 623V/cm
1316 Float_t time1800[ktimebin][kZbin] =
1317 {{0.09166, 0.09199, 0.09300, 0.09470, 0.09709, 0.10026,
1318 0.10434, 0.10983, 0.11831, 0.13987, 0.37802},
1319 {0.06585, 0.06623, 0.06735, 0.06921, 0.07180, 0.07520,
1320 0.07941, 0.08507, 0.09376, 0.11604, 0.35624},
1321 {0.03945, 0.04004, 0.04169, 0.04432, 0.04776, 0.05187,
1322 0.05674, 0.06300, 0.07290, 0.09884, 0.34704},
1323 {0.01108, 0.01287, 0.01717, 0.02280, 0.02880, 0.03476,
1324 0.04095, 0.04909, 0.06419, 0.10589, 0.36843},
1325 {0.01115, 0.01287, 0.01720, 0.02276, 0.02862, 0.03448,
1326 0.04073, 0.04973, 0.06799, 0.11535, 0.37224},
1327 {0.03918, 0.03975, 0.04134, 0.04382, 0.04707, 0.05105,
1328 0.05603, 0.06296, 0.07605, 0.10822, 0.36560},
1329 {0.06498, 0.06532, 0.06635, 0.06809, 0.07058, 0.07395,
1330 0.07855, 0.08495, 0.09632, 0.12488, 0.37730},
1331 {0.09130, 0.09160, 0.09199, 0.09300, 0.09472, 0.00000,
1332 0.11059, 0.11747, 0.13146, 0.16462, 0.42233},
1333 {0.13364, 0.13449, 0.13767, 0.14481, 0.16147, 0.23635,
1334 0.19706, 0.18812, 0.19704, 0.23520, 0.49749},
1335 {0.19698, 0.19844, 0.20311, 0.21236, 0.23082, 0.29920,
1336 0.26936, 0.25927, 0.26732, 0.30601, 0.56871},
1337 {0.26518, 0.26670, 0.27160, 0.28099, 0.29955, 0.36597,
1338 0.33885, 0.32858, 0.33653, 0.37524, 0.63801},
1339 {0.33441, 0.33553, 0.34040, 0.34987, 0.36841, 0.43428,
1340 0.40797, 0.39763, 0.40556, 0.44425, 0.70698},
1341 {0.40296, 0.40447, 0.40934, 0.41881, 0.43737, 0.50306,
1342 0.47695, 0.46662, 0.47455, 0.51329, 0.77600},
1343 {0.47296, 0.47344, 0.47830, 0.48779, 0.50632, 0.57200,
1344 0.54593, 0.53559, 0.54351, 0.58222, 0.84489},
1345 {0.54089, 0.54264, 0.54727, 0.55673, 0.57529, 0.64094,
1346 0.61490, 0.60457, 0.61249, 0.65118, 0.91394},
1347 {0.60987, 0.61138, 0.61624, 0.62573, 0.64428, 0.70989,
1348 0.68397, 0.67354, 0.68147, 0.72015, 0.98291},
1349 {0.67883, 0.68035, 0.68521, 0.69469, 0.71324, 0.77896,
1350 0.75287, 0.74251, 0.75043, 0.78912, 1.04458},
1351 {0.74780, 0.74932, 0.75421, 0.76367, 0.78221, 0.84785,
1352 0.82185, 0.81148, 0.81941, 0.85811, 1.12085},
1353 {0.81690, 0.81830, 0.82316, 0.83263, 0.85120, 0.91683,
1354 0.89077, 0.88045, 0.88837, 0.92707, 1.18976},
1355 {0.88574, 0.88726, 0.89228, 0.90198, 0.92017, 0.98578,
1356 0.95974, 0.94947, 0.95734, 0.99604, 1.25873},
1357 {0.95493, 0.95624, 0.96110, 0.97058, 0.98913, 1.05481,
1358 1.02873, 1.01839, 1.02631, 1.06503, 1.32772},
1359 {1.02392, 1.02524, 1.03008, 1.03955, 1.05810, 1.12378,
1360 1.09757, 1.08605, 1.09530, 1.13399, 1.39669},
1361 {1.09270, 1.09418, 1.09911, 1.10854, 1.12714, 1.19281,
1362 1.16502, 1.15633, 1.16427, 1.20271, 1.46574},
1363 {1.16168, 1.16323, 1.16801, 1.17772, 1.19604, 1.26190,
1364 1.23399, 1.22531, 1.23323, 1.27194, 1.53475},
1365 {1.23061, 1.23214, 1.23698, 1.24669, 1.26503, 1.33073,
1366 1.30461, 1.29428, 1.30220, 1.34091, 1.60372},
1367 {1.29960, 1.30110, 1.30596, 1.31544, 1.33398, 1.39962,
1368 1.37228, 1.36323, 1.37121, 1.40988, 1.67273},
1369 {1.36851, 1.37007, 1.37512, 1.38441, 1.40297, 1.46865,
1370 1.44256, 1.43222, 1.44017, 1.47878, 1.74155},
1371 {1.43752, 1.43904, 1.44773, 1.45338, 1.47220, 1.53759,
1372 1.51136, 1.50119, 1.50914, 1.54775, 1.81050},
1373 {1.50646, 1.50802, 1.51288, 1.52237, 1.54097, 1.60697,
1374 1.58049, 1.57018, 1.57811, 1.61678, 1.87947},
1375 {1.57545, 1.57720, 1.58185, 1.59134, 1.60996, 1.67787,
1376 1.64929, 1.63914, 1.64707, 1.68570, 1.94851},
1377 {1.64442, 1.64617, 1.65081, 1.66035, 1.67893, 1.74684,
1378 1.71826, 1.70745, 1.71604, 1.75310, 2.01748},
1379 {1.71337, 1.71513, 1.71978, 1.72932, 1.74645, 1.81346,
1380 1.78739, 1.77642, 1.78501, 1.82151, 2.08644},
1381 {1.78238, 1.78410, 1.78876, 1.79824, 1.81678, 1.88332,
1382 1.85639, 1.84262, 1.85397, 1.89270, 2.15533},
1383 {1.85135, 1.85306, 1.85778, 1.86728, 1.88580, 1.95615,
1384 1.92536, 1.91171, 1.92283, 1.96165, 2.22428},
1385 {1.92774, 1.92184, 1.92672, 1.93618, 1.95477, 2.02048,
1386 1.99427, 1.98068, 1.99192, 2.03062, 2.29338},
1387 {1.98929, 1.99081, 1.99567, 2.00515, 2.02373, 2.08987,
1388 2.06332, 2.05249, 2.05939, 2.09928, 2.36227},
1389 {2.05829, 2.05978, 2.06464, 2.07414, 2.09272, 2.15850,
1390 2.12928, 2.12194, 2.12987, 2.16825, 2.43083},
1391 {2.12726, 2.12869, 2.13360, 2.14425, 2.16160, 2.22872,
1392 2.20118, 2.19078, 2.19876, 2.23416, 2.50007}};
1393
1394 // Garfield simulation at UD = -1900V; vd = 1.610cm/microsec, <driftfield> = 655V/cm
1395 Float_t time1900[ktimebin][kZbin] =
1396 {{0.09166, 0.09198, 0.09298, 0.09467, 0.09707, 0.10023,
1397 0.10431, 0.10980, 0.11828, 0.13988, 0.37789},
1398 {0.06584, 0.06622, 0.06735, 0.06920, 0.07179, 0.07514,
1399 0.07938, 0.08505, 0.09374, 0.11606, 0.35599},
1400 {0.03945, 0.04002, 0.04169, 0.04432, 0.04775, 0.05185,
1401 0.05672, 0.06298, 0.07290, 0.09888, 0.34695},
1402 {0.01109, 0.01281, 0.01717, 0.02279, 0.02878, 0.03476,
1403 0.04094, 0.04909, 0.06421, 0.10597, 0.36823},
1404 {0.01115, 0.01287, 0.01720, 0.02294, 0.02862, 0.03448,
1405 0.04074, 0.04973, 0.06783, 0.11506, 0.37206},
1406 {0.03940, 0.03975, 0.04135, 0.04386, 0.04708, 0.05106,
1407 0.05604, 0.06293, 0.07592, 0.10787, 0.36484},
1408 {0.06500, 0.06534, 0.06638, 0.06811, 0.07060, 0.07413,
1409 0.07852, 0.08491, 0.09614, 0.12446, 0.37626},
1410 {0.09119, 0.09140, 0.09194, 0.09293, 0.09471, 0.00000,
1411 0.11013, 0.11685, 0.13050, 0.16302, 0.41991},
1412 {0.13111, 0.13190, 0.13466, 0.14091, 0.15554, 0.22409,
1413 0.18846, 0.18167, 0.19113, 0.22854, 0.48995},
1414 {0.18849, 0.18975, 0.19380, 0.20185, 0.21797, 0.28050,
1415 0.25368, 0.24575, 0.25446, 0.29249, 0.55442},
1416 {0.24995, 0.25125, 0.25563, 0.26366, 0.27986, 0.34065,
1417 0.31605, 0.30815, 0.31680, 0.35482, 0.61697},
1418 {0.31187, 0.31324, 0.31745, 0.32580, 0.34205, 0.40217,
1419 0.37825, 0.37031, 0.37897, 0.41696, 0.67890},
1420 {0.37401, 0.37531, 0.37955, 0.38777, 0.40395, 0.46408,
1421 0.44037, 0.43242, 0.44108, 0.47906, 0.74122},
1422 {0.43610, 0.43741, 0.44161, 0.44986, 0.46604, 0.52614,
1423 0.50248, 0.49452, 0.50316, 0.54116, 0.80326},
1424 {0.49820, 0.49988, 0.50372, 0.51196, 0.52814, 0.58822,
1425 0.56459, 0.55661, 0.56527, 0.60326, 0.86526},
1426 {0.56032, 0.56161, 0.56582, 0.57408, 0.59024, 0.65032,
1427 0.62670, 0.61872, 0.62737, 0.66537, 0.92749},
1428 {0.62240, 0.62371, 0.62792, 0.63624, 0.65236, 0.71241,
1429 0.68881, 0.68081, 0.68947, 0.72750, 0.98941},
1430 {0.68449, 0.68581, 0.69002, 0.69828, 0.71444, 0.77452,
1431 0.75089, 0.74295, 0.75157, 0.78957, 1.05157},
1432 {0.74660, 0.74790, 0.75212, 0.76036, 0.77654, 0.83748,
1433 0.81299, 0.80501, 0.81193, 0.85168, 1.11375},
1434 {0.80870, 0.81017, 0.81423, 0.82246, 0.83867, 0.89908,
1435 0.87509, 0.86660, 0.87577, 0.91376, 1.17586},
1436 {0.87080, 0.87233, 0.87632, 0.88458, 0.90074, 0.96083,
1437 0.93718, 0.92922, 0.93787, 0.97588, 1.23794},
1438 {0.93291, 0.93422, 0.93844, 0.94667, 0.96293, 1.02295,
1439 0.99929, 0.99127, 0.99997, 1.03797, 1.29995},
1440 {0.99500, 0.99645, 1.00308, 1.00877, 1.02497, 1.08504,
1441 1.06140, 1.05343, 1.06203, 1.10006, 1.36216},
1442 {1.05712, 1.05926, 1.06262, 1.07092, 1.08706, 1.14754,
1443 1.12350, 1.11550, 1.12417, 1.16218, 1.42427},
1444 {1.11921, 1.12059, 1.12473, 1.13297, 1.14916, 1.21140,
1445 1.18560, 1.17284, 1.18625, 1.22414, 1.48629},
1446 {1.18140, 1.18262, 1.18690, 1.19508, 1.21125, 1.27139,
1447 1.24164, 1.23495, 1.24838, 1.28634, 1.54852},
1448 {1.24340, 1.24473, 1.24901, 1.25732, 1.27336, 1.33358,
1449 1.30793, 1.30179, 1.31047, 1.34848, 1.61066},
1450 {1.30551, 1.30684, 1.31104, 1.32056, 1.33553, 1.39609,
1451 1.37004, 1.36392, 1.37045, 1.41057, 1.67259},
1452 {1.36755, 1.36892, 1.37315, 1.39148, 1.39755, 1.45820,
1453 1.43215, 1.42602, 1.43467, 1.47268, 1.73477},
1454 {1.42966, 1.43101, 1.43549, 1.45359, 1.45976, 1.52031,
1455 1.49601, 1.48811, 1.49677, 1.53477, 1.79691},
1456 {1.49180, 1.49321, 1.49760, 1.51570, 1.52175, 1.58185,
1457 1.55771, 1.55023, 1.55888, 1.59672, 1.85501},
1458 {1.55391, 1.55527, 1.55943, 1.57782, 1.58391, 1.64395,
1459 1.62008, 1.61233, 1.62085, 1.65883, 1.92091},
1460 {1.61599, 1.61732, 1.62154, 1.63993, 1.64612, 1.70608,
1461 1.68237, 1.67108, 1.68301, 1.72110, 1.97931},
1462 {1.67808, 1.67948, 1.68364, 1.70204, 1.70823, 1.76858,
1463 1.74404, 1.73539, 1.74512, 1.78321, 2.04522},
1464 {1.74019, 1.74152, 1.74573, 1.76415, 1.77015, 1.83040,
1465 1.80615, 1.79366, 1.80723, 1.84509, 2.10742},
1466 {1.80235, 1.80362, 1.80783, 1.82626, 1.83227, 1.89246,
1467 1.86795, 1.85405, 1.86938, 1.90720, 2.16953},
1468 {1.86442, 1.86572, 1.86994, 1.88837, 1.89438, 1.95445,
1469 1.93006, 1.92283, 1.93148, 1.96931, 2.23147},
1470 {1.92700, 1.92783, 1.93197, 1.95049, 1.95649, 2.01668,
1471 1.99217, 1.98486, 1.99352, 2.03143, 2.29358}};
1472
1473 // Garfield simulation at UD = -2000V; vd = 1.783cm/microsec, <driftfield> = 688V/cm
1474 Float_t time2000[ktimebin][kZbin] =
1475 {{0.09176, 0.09196, 0.09296, 0.09465, 0.09704, 0.10020,
1476 0.10427, 0.10977, 0.11825, 0.13988, 0.37774},
1477 {0.06583, 0.06620, 0.06735, 0.06918, 0.07177, 0.07513,
1478 0.07936, 0.08503, 0.09372, 0.11606, 0.35586},
1479 {0.03944, 0.04001, 0.04170, 0.04431, 0.04774, 0.05184,
1480 0.05670, 0.06296, 0.07291, 0.09893, 0.34680},
1481 {0.01108, 0.01281, 0.01719, 0.02279, 0.02879, 0.03474,
1482 0.04093, 0.04908, 0.06422, 0.10605, 0.36800},
1483 {0.01114, 0.01287, 0.01720, 0.02276, 0.02863, 0.03449,
1484 0.04073, 0.04970, 0.06774, 0.11478, 0.37179},
1485 {0.03925, 0.03977, 0.04135, 0.04386, 0.04711, 0.05108,
1486 0.05604, 0.06290, 0.07580, 0.10748, 0.36386},
1487 {0.06501, 0.06536, 0.06640, 0.06814, 0.07062, 0.07398,
1488 0.07852, 0.08487, 0.09598, 0.12405, 0.37519},
1489 {0.09109, 0.09127, 0.09188, 0.09292, 0.09472, 0.00000,
1490 0.10964, 0.11630, 0.12960, 0.16150, 0.41765},
1491 {0.12898, 0.12968, 0.13209, 0.13749, 0.15034, 0.21286,
1492 0.18088, 0.17590, 0.18591, 0.22254, 0.48315},
1493 {0.18122, 0.18227, 0.18574, 0.19263, 0.20674, 0.26376,
1494 0.23960, 0.23375, 0.24316, 0.28047, 0.54179},
1495 {0.23674, 0.23784, 0.24142, 0.24847, 0.26264, 0.31810,
1496 0.29602, 0.29008, 0.29944, 0.33675, 0.59795},
1497 {0.29279, 0.29382, 0.29742, 0.30448, 0.31865, 0.37364,
1498 0.35215, 0.34629, 0.35555, 0.39286, 0.65411},
1499 {0.34875, 0.34987, 0.35346, 0.36054, 0.37472, 0.42956,
1500 0.40825, 0.40229, 0.41167, 0.44894, 0.71033},
1501 {0.40484, 0.40594, 0.40954, 0.41660, 0.43077, 0.48560,
1502 0.46433, 0.45840, 0.46772, 0.50500, 0.76632},
1503 {0.46090, 0.46201, 0.46560, 0.47267, 0.48684, 0.54167,
1504 0.52041, 0.51449, 0.52382, 0.56108, 0.82227},
1505 {0.51698, 0.51809, 0.52173, 0.52874, 0.54291, 0.59776,
1506 0.57646, 0.57052, 0.57986, 0.61717, 0.87836},
1507 {0.57306, 0.57418, 0.57782, 0.58513, 0.59899, 0.65380,
1508 0.63255, 0.62661, 0.63594, 0.67325, 0.93460},
1509 {0.62912, 0.63024, 0.63383, 0.64103, 0.65506, 0.70988,
1510 0.68484, 0.68267, 0.69202, 0.72878, 0.99046},
1511 {0.68521, 0.68633, 0.68990, 0.69699, 0.71115, 0.76595,
1512 0.74468, 0.73872, 0.74814, 0.78538, 1.04674},
1513 {0.74127, 0.74239, 0.74605, 0.75303, 0.77022, 0.82204,
1514 0.80078, 0.79484, 0.80416, 0.84147, 1.10261},
1515 {0.79736, 0.79846, 0.80206, 0.80947, 0.82330, 0.87813,
1516 0.85688, 0.85087, 0.86023, 0.89752, 1.15874},
1517 {0.85342, 0.85454, 0.85815, 0.86519, 0.87936, 0.93417,
1518 0.91293, 0.90428, 0.91631, 0.95360, 1.20760},
1519 {0.90949, 0.91061, 0.91423, 0.92128, 0.93544, 0.99026,
1520 0.96807, 0.96305, 0.97239, 1.00967, 1.27078},
1521 {0.96556, 0.96669, 0.97111, 0.97734, 0.99151, 1.04664,
1522 1.02508, 1.01879, 1.02846, 1.06167, 1.32695},
1523 {1.02167, 1.02279, 1.02656, 1.03341, 1.04759, 1.10242,
1524 1.08115, 1.07003, 1.08453, 1.12184, 1.38304},
1525 {1.07776, 1.07883, 1.08242, 1.08950, 1.10384, 1.16422,
1526 1.13725, 1.13133, 1.14061, 1.17793, 1.43910},
1527 {1.13379, 1.13492, 1.13864, 1.14567, 1.15973, 1.21455,
1528 1.19323, 1.18734, 1.19668, 1.23401, 1.49528},
1529 {1.18988, 1.19098, 1.19457, 1.20164, 1.21582, 1.27064,
1530 1.24937, 1.24044, 1.25275, 1.29004, 1.55137},
1531 {1.24592, 1.24706, 1.25087, 1.25773, 1.27188, 1.32670,
1532 1.30544, 1.29953, 1.30883, 1.34613, 1.60743},
1533 {1.30202, 1.30313, 1.30673, 1.31381, 1.32797, 1.38278,
1534 1.36151, 1.35167, 1.36490, 1.40221, 1.66306},
1535 {1.35809, 1.35921, 1.36282, 1.36986, 1.38403, 1.43888,
1536 1.41760, 1.41174, 1.42083, 1.45830, 1.71915},
1537 {1.41419, 1.41528, 1.41890, 1.42595, 1.44011, 1.49496,
1538 1.47368, 1.46769, 1.47706, 1.51436, 1.77523},
1539 {1.47131, 1.47141, 1.47494, 1.48850, 1.49620, 1.55137,
1540 1.52977, 1.51820, 1.53315, 1.57042, 1.83158},
1541 {1.52635, 1.52750, 1.53103, 1.53814, 1.55228, 1.60736,
1542 1.58503, 1.57986, 1.58920, 1.62649, 1.88767},
1543 {1.58418, 1.58355, 1.58711, 1.59526, 1.60833, 1.66316,
1544 1.63345, 1.63261, 1.64556, 1.68204, 1.94359},
1545 {1.64027, 1.63958, 1.64489, 1.65024, 1.66443, 1.71925,
1546 1.69794, 1.69201, 1.70143, 1.73865, 1.99968},
1547 {1.69450, 1.69566, 1.69940, 1.70697, 1.71841, 1.77819,
1548 1.75396, 1.74814, 1.75743, 1.79083, 2.05427},
1549 {1.75054, 1.75221, 1.75527, 1.76306, 1.77662, 1.83428,
1550 1.81006, 1.81173, 1.81345, 1.85076, 2.10289}};
1551
1552 // Garfield simulation at UD = -2100V; vd = 1.959cm/microsec, <driftfield> = 720V/cm
1553 Float_t time2100[ktimebin][kZbin] =
1554 {{0.09160, 0.09194, 0.09294, 0.09462, 0.09701, 0.10017,
1555 0.10424, 0.10974, 0.11823, 0.13988, 0.37762},
1556 {0.06585, 0.06619, 0.06731, 0.06916, 0.07174, 0.07509,
1557 0.07933, 0.08500, 0.09370, 0.11609, 0.35565},
1558 {0.03960, 0.04001, 0.04171, 0.04430, 0.04774, 0.05182,
1559 0.05668, 0.06294, 0.07291, 0.09896, 0.34676},
1560 {0.01109, 0.01280, 0.01716, 0.02279, 0.02876, 0.03474,
1561 0.04096, 0.04908, 0.06424, 0.10612, 0.36790},
1562 {0.01114, 0.01285, 0.01719, 0.02287, 0.02863, 0.03449,
1563 0.04073, 0.04964, 0.06759, 0.11446, 0.37162},
1564 {0.03922, 0.03977, 0.04146, 0.04386, 0.04711, 0.05109,
1565 0.05605, 0.06287, 0.07575, 0.10713, 0.36298},
1566 {0.06504, 0.06538, 0.06641, 0.06818, 0.07064, 0.07426,
1567 0.07852, 0.08483, 0.09581, 0.12363, 0.37424},
1568 {0.09103, 0.09129, 0.09186, 0.09291, 0.09476, 0.00000,
1569 0.10923, 0.11578, 0.12873, 0.16005, 0.41525},
1570 {0.12723, 0.12777, 0.12988, 0.13458, 0.14579, 0.20264,
1571 0.17421, 0.17078, 0.18132, 0.21708, 0.47699},
1572 {0.17508, 0.17601, 0.17897, 0.18487, 0.19698, 0.24881,
1573 0.22737, 0.22337, 0.23348, 0.27000, 0.53032},
1574 {0.22571, 0.22663, 0.22969, 0.23570, 0.24787, 0.29826,
1575 0.27871, 0.27462, 0.28471, 0.32122, 0.58166},
1576 {0.27664, 0.27759, 0.28067, 0.28669, 0.29891, 0.34898,
1577 0.32982, 0.32570, 0.33576, 0.37229, 0.63268},
1578 {0.32766, 0.32862, 0.33170, 0.33778, 0.34988, 0.39973,
1579 0.38088, 0.37675, 0.38680, 0.42333, 0.68159},
1580 {0.37872, 0.37966, 0.38275, 0.38875, 0.40093, 0.45073,
1581 0.43192, 0.42780, 0.43786, 0.47438, 0.73480},
1582 {0.42974, 0.43070, 0.43378, 0.43982, 0.45196, 0.50177,
1583 0.48297, 0.47884, 0.48889, 0.52544, 0.78581},
1584 {0.48081, 0.48175, 0.48482, 0.49084, 0.50302, 0.55290,
1585 0.53398, 0.52988, 0.53994, 0.57647, 0.83687},
1586 {0.53645, 0.53295, 0.53586, 0.54188, 0.55408, 0.60398,
1587 0.58504, 0.58092, 0.59100, 0.62768, 0.88773},
1588 {0.58345, 0.58409, 0.58690, 0.59292, 0.60510, 0.65562,
1589 0.63609, 0.63197, 0.64203, 0.67856, 0.93907},
1590 {0.63397, 0.63490, 0.63795, 0.64403, 0.65613, 0.70612,
1591 0.68714, 0.68301, 0.69294, 0.72955, 0.99000},
1592 {0.68496, 0.68592, 0.68899, 0.69504, 0.70733, 0.75708,
1593 0.73818, 0.73405, 0.74412, 0.78064, 1.04100},
1594 {0.73600, 0.73696, 0.74003, 0.74624, 0.75828, 0.80805,
1595 0.78904, 0.78512, 0.79517, 0.83152, 1.09205},
1596 {0.78709, 0.78801, 0.79108, 0.79709, 0.80931, 0.85906,
1597 0.84027, 0.83614, 0.84621, 0.88269, 1.14058},
1598 {0.83808, 0.83905, 0.84215, 0.84816, 0.86031, 0.91011,
1599 0.89139, 0.88718, 0.89725, 0.93377, 1.19413},
1600 {0.88916, 0.89010, 0.89320, 0.89920, 0.91136, 0.96117,
1601 0.94235, 0.93822, 0.94828, 0.98480, 1.24538},
1602 {0.94036, 0.94113, 0.94422, 0.95023, 0.96241, 1.01220,
1603 0.99310, 0.98927, 0.99933, 1.03585, 1.29629},
1604 {0.99139, 0.99220, 0.99525, 1.00127, 1.01344, 1.06324,
1605 1.04451, 1.04033, 1.04836, 1.08690, 1.34727},
1606 {1.04261, 1.04325, 1.04628, 1.05232, 1.06448, 1.12090,
1607 1.09546, 1.09136, 1.10142, 1.13795, 1.39831},
1608 {1.09331, 1.09429, 1.09742, 1.10336, 1.11557, 1.16547,
1609 1.14658, 1.13642, 1.15247, 1.18898, 1.44936},
1610 {1.14436, 1.14539, 1.14847, 1.15443, 1.16662, 1.21794,
1611 1.19763, 1.19329, 1.20349, 1.23956, 1.50043},
1612 {1.19533, 1.19651, 1.19943, 1.20548, 1.21666, 1.26753,
1613 1.24862, 1.24434, 1.25455, 1.29106, 1.55142},
1614 {1.24638, 1.24756, 1.25046, 1.25648, 1.26764, 1.31858,
1615 1.29967, 1.29538, 1.30499, 1.34211, 1.60250},
1616 {1.29747, 1.29847, 1.30175, 1.30753, 1.31869, 1.36969,
1617 1.35069, 1.34656, 1.35663, 1.39316, 1.64644},
1618 {1.35537, 1.34952, 1.35255, 1.35869, 1.36973, 1.41387,
1619 1.40173, 1.39761, 1.40768, 1.44396, 1.70238},
1620 {1.39956, 1.40056, 1.40380, 1.40961, 1.42178, 1.46492,
1621 1.45278, 1.45423, 1.45872, 1.49522, 1.75557},
1622 {1.45080, 1.45159, 1.45463, 1.46109, 1.47287, 1.52263,
1623 1.50382, 1.50050, 1.50977, 1.54502, 1.80670},
1624 {1.50165, 1.50264, 1.50570, 1.51214, 1.52233, 1.57370,
1625 1.55484, 1.55155, 1.56080, 1.59731, 1.85778},
1626 {1.55269, 1.55364, 1.55675, 1.56274, 1.57491, 1.62598,
1627 1.60590, 1.60259, 1.61185, 1.64836, 1.90883},
1628 {1.60368, 1.60469, 1.60779, 1.61373, 1.62596, 1.67738,
1629 1.65651, 1.65249, 1.66290, 1.69936, 1.95959}};
1630
1631 // Garfield simulation at UD = -2200V; vd = 2.134cm/microsec, <driftfield> = 753V/cm
1632 Float_t time2200[ktimebin][kZbin] =
1633 {{0.09162, 0.09194, 0.09292, 0.09460, 0.09702, 0.10014,
1634 0.10421, 0.10971, 0.11820, 0.13990, 0.37745},
1635 {0.06581, 0.06618, 0.06730, 0.06915, 0.07173, 0.07507,
1636 0.07931, 0.08497, 0.09368, 0.11609, 0.35560},
1637 {0.03947, 0.04001, 0.04167, 0.04429, 0.04772, 0.05183,
1638 0.05667, 0.06293, 0.07292, 0.09900, 0.34673},
1639 {0.01111, 0.01280, 0.01716, 0.02279, 0.02876, 0.03473,
1640 0.04091, 0.04907, 0.06426, 0.10620, 0.36766},
1641 {0.01113, 0.01285, 0.01719, 0.02276, 0.02863, 0.03452,
1642 0.04076, 0.04960, 0.06745, 0.11419, 0.37139},
1643 {0.03923, 0.03978, 0.04137, 0.04387, 0.04713, 0.05110,
1644 0.05605, 0.06284, 0.07551, 0.10677, 0.36210},
1645 {0.06505, 0.06540, 0.06644, 0.06820, 0.07069, 0.07401,
1646 0.07852, 0.08479, 0.09565, 0.12325, 0.37313},
1647 {0.09107, 0.09127, 0.09181, 0.09291, 0.09474, 0.00000,
1648 0.10883, 0.11528, 0.12789, 0.15865, 0.41313},
1649 {0.12559, 0.12622, 0.12800, 0.13206, 0.14166, 0.19331,
1650 0.16832, 0.16632, 0.17724, 0.21218, 0.47098},
1651 {0.16992, 0.17070, 0.17325, 0.17831, 0.18871, 0.23557,
1652 0.21690, 0.21451, 0.22514, 0.26082, 0.52034},
1653 {0.21640, 0.21722, 0.21987, 0.22500, 0.23540, 0.28097,
1654 0.26399, 0.26154, 0.27214, 0.30784, 0.56734},
1655 {0.26318, 0.26400, 0.26679, 0.27181, 0.28220, 0.32739,
1656 0.31090, 0.30842, 0.31902, 0.35474, 0.61415},
1657 {0.31001, 0.31085, 0.31348, 0.31866, 0.32903, 0.37412,
1658 0.35777, 0.35546, 0.36588, 0.40159, 0.66103},
1659 {0.35687, 0.35769, 0.36033, 0.36556, 0.37588, 0.42094,
1660 0.40523, 0.40214, 0.41273, 0.44841, 0.70785},
1661 {0.40372, 0.40457, 0.40723, 0.41234, 0.42273, 0.46778,
1662 0.45148, 0.44903, 0.45961, 0.49526, 0.75486},
1663 {0.45062, 0.45139, 0.45404, 0.45966, 0.46958, 0.51470,
1664 0.49833, 0.49584, 0.50644, 0.54211, 0.80160},
1665 {0.49742, 0.49825, 0.50088, 0.50605, 0.51644, 0.56148,
1666 0.54518, 0.54270, 0.55330, 0.58897, 0.84854},
1667 {0.54427, 0.54510, 0.54774, 0.55290, 0.56329, 0.60846,
1668 0.59203, 0.58955, 0.60014, 0.63578, 0.89528},
1669 {0.59119, 0.59199, 0.59471, 0.59975, 0.61014, 0.65533,
1670 0.63889, 0.63636, 0.64699, 0.68269, 0.94197},
1671 {0.63866, 0.63880, 0.64145, 0.64664, 0.65701, 0.70639,
1672 0.68574, 0.68325, 0.69385, 0.72949, 0.98900},
1673 {0.68483, 0.68566, 0.68831, 0.69347, 0.70386, 0.74890,
1674 0.73260, 0.73010, 0.74069, 0.77638, 1.03320},
1675 {0.73168, 0.73255, 0.73515, 0.74031, 0.75072, 0.79576,
1676 0.77117, 0.77501, 0.78755, 0.82318, 1.08006},
1677 {0.77854, 0.78310, 0.78200, 0.79525, 0.79756, 0.84281,
1678 0.81803, 0.82393, 0.83441, 0.87008, 1.12692},
1679 {0.82541, 0.82642, 0.82916, 0.83528, 0.84442, 0.89086,
1680 0.87315, 0.87079, 0.88125, 0.91694, 1.17648},
1681 {0.87226, 0.87308, 0.87602, 0.88086, 0.89128, 0.93772,
1682 0.92001, 0.91751, 0.92811, 0.95587, 1.22328},
1683 {0.91921, 0.91994, 0.92256, 0.92772, 0.94713, 0.98566,
1684 0.96690, 0.96436, 0.97495, 1.01064, 1.26882},
1685 {0.96790, 0.96679, 0.96941, 0.97463, 0.99399, 1.03001,
1686 1.01376, 1.01112, 1.02181, 1.05749, 1.31568},
1687 {1.01278, 1.01390, 1.01674, 1.02147, 1.03182, 1.07820,
1688 1.06056, 1.05798, 1.06867, 1.10433, 1.36390},
1689 {1.05964, 1.06076, 1.06331, 1.06833, 1.07870, 1.13297,
1690 1.10742, 1.10520, 1.11543, 1.15120, 1.41069},
1691 {1.10664, 1.10762, 1.10997, 1.11519, 1.12556, 1.17531,
1692 1.15427, 1.14620, 1.16229, 1.19805, 1.45758},
1693 {1.15352, 1.15421, 1.15683, 1.16218, 1.17242, 1.21910,
1694 1.20035, 1.19863, 1.20579, 1.24473, 1.50412},
1695 {1.20019, 1.20115, 1.20369, 1.20892, 1.21928, 1.26913,
1696 1.24721, 1.24549, 1.25605, 1.29159, 1.54920},
1697 {1.24707, 1.24846, 1.25052, 1.25602, 1.26608, 1.31558,
1698 1.29448, 1.29232, 1.30293, 1.33675, 1.59798},
1699 {1.29391, 1.29475, 1.29738, 1.30255, 1.31294, 1.36244,
1700 1.34167, 1.33918, 1.34979, 1.38229, 1.64496},
1701 {1.34078, 1.34304, 1.34424, 1.35565, 1.35980, 1.40930,
1702 1.38853, 1.38229, 1.39664, 1.42863, 1.69162},
1703 {1.38762, 1.38847, 1.39110, 1.39627, 1.40666, 1.45183,
1704 1.43539, 1.43289, 1.44348, 1.47549, 1.73876},
1705 {1.43524, 1.43533, 1.43796, 1.44310, 1.45371, 1.49305,
1706 1.48224, 1.47941, 1.49034, 1.52601, 1.78552},
1707 {1.48122, 1.48219, 1.48482, 1.48991, 1.50030, 1.53991,
1708 1.52898, 1.52653, 1.53653, 1.57282, 1.82386}};
1709
1710 if (fTimeStruct1) delete [] fTimeStruct1;
1711 fTimeStruct1 = new Float_t[ktimebin*kZbin];
1712
1713 if (fTimeStruct2) delete [] fTimeStruct2;
1714 fTimeStruct2 = new Float_t[ktimebin*kZbin];
1715
1716 for (Int_t ctrt = 0; ctrt<ktimebin; ctrt++) {
a328fff9 1717 for (Int_t ctrz = 0; ctrz<kZbin; ctrz++) {
ccb4315c 1718 if ( fDriftVelocity > fVDsmp[6] ) {
1719 fTimeStruct1[ctrt+ctrz*ktimebin] = time2100[ctrt][ctrz];
1720 fTimeStruct2[ctrt+ctrz*ktimebin] = time2200[ctrt][ctrz];
1721 fVDlo = fVDsmp[6];
1722 fVDhi = fVDsmp[7];
1723 } else if ( fDriftVelocity > fVDsmp[5] ) {
1724 fTimeStruct1[ctrt+ctrz*ktimebin] = time2000[ctrt][ctrz];
1725 fTimeStruct2[ctrt+ctrz*ktimebin] = time2100[ctrt][ctrz];
1726 fVDlo = fVDsmp[5];
1727 fVDhi = fVDsmp[6];
1728 } else if ( fDriftVelocity > fVDsmp[4] ) {
1729 fTimeStruct1[ctrt+ctrz*ktimebin] = time1900[ctrt][ctrz];
1730 fTimeStruct2[ctrt+ctrz*ktimebin] = time2000[ctrt][ctrz];
1731 fVDlo = fVDsmp[4];
1732 fVDhi = fVDsmp[5];
1733 } else if ( fDriftVelocity > fVDsmp[3] ) {
1734 fTimeStruct1[ctrt+ctrz*ktimebin] = time1800[ctrt][ctrz];
1735 fTimeStruct2[ctrt+ctrz*ktimebin] = time1900[ctrt][ctrz];
1736 fVDlo = fVDsmp[3];
1737 fVDhi = fVDsmp[4];
1738 } else if ( fDriftVelocity > fVDsmp[2] ) {
1739 fTimeStruct1[ctrt+ctrz*ktimebin] = time1700[ctrt][ctrz];
1740 fTimeStruct2[ctrt+ctrz*ktimebin] = time1800[ctrt][ctrz];
1741 fVDlo = fVDsmp[2];
1742 fVDhi = fVDsmp[3];
1743 } else if ( fDriftVelocity > fVDsmp[1] ) {
1744 fTimeStruct1[ctrt+ctrz*ktimebin] = time1600[ctrt][ctrz];
1745 fTimeStruct2[ctrt+ctrz*ktimebin] = time1700[ctrt][ctrz];
1746 fVDlo = fVDsmp[1];
1747 fVDhi = fVDsmp[2];
1748 } else if ( fDriftVelocity > (fVDsmp[0] - 1.e-5) ) {
1749 fTimeStruct1[ctrt+ctrz*ktimebin] = time1500[ctrt][ctrz];
1750 fTimeStruct2[ctrt+ctrz*ktimebin] = time1600[ctrt][ctrz];
1751 fVDlo = fVDsmp[0];
1752 fVDhi = fVDsmp[1];
1753 }
a328fff9 1754 }
ccb4315c 1755 }
1756
1757 //printf("fVDlo = %f, fVDhi = %f\n", fVDlo, fVDhi);
1758
1759 ReInit();
a328fff9 1760
1761}
1762
17b26de4 1763//_____________________________________________________________________________
1764void AliTRDparameter::SamplePRF()
1765{
1766 //
1767 // Samples the pad response function
1768 //
1769
1770 const Int_t kNplan = AliTRDgeometry::kNplan;
1771 const Int_t kPRFbin = 61;
1772
a328fff9 1773 Float_t prf[kNplan][kPRFbin] = { {0.0267681, 0.031056, 0.036111, 0.0419945, 0.0489578,
1774 0.0571508, 0.0666984, 0.078064, 0.0912508, 0.106579,
1775 0.124429, 0.144803, 0.168187, 0.194645, 0.224143,
1776 0.25717, 0.293146, 0.331943, 0.373128, 0.415798,
1777 0.459149, 0.502028, 0.54349, 0.582109, 0.617335,
1778 0.648306, 0.674189, 0.695082, 0.710222, 0.719711,
1779 0.723751, 0.719711, 0.710222, 0.695082, 0.674189,
1780 0.648306, 0.617335, 0.582109, 0.54349, 0.502028,
1781 0.459149, 0.415798, 0.373128, 0.331943, 0.293146,
1782 0.25717, 0.224143, 0.194645, 0.168187, 0.144803,
1783 0.124429, 0.106579, 0.0912508, 0.078064, 0.0666984,
1784 0.0571508, 0.0489578, 0.0419945, 0.036111, 0.031056,
1785 0.0267681},
1786 {0.0240831, 0.0281427, 0.0328987, 0.0384801, 0.0451249,
1787 0.0530119, 0.0623139, 0.0733268, 0.0864095, 0.101704,
1788 0.119473, 0.140308, 0.164191, 0.191333, 0.222141,
1789 0.256551, 0.29433, 0.335294, 0.37892, 0.424206,
1790 0.470202, 0.515621, 0.559255, 0.600068, 0.636821,
1791 0.668836, 0.695767, 0.717076, 0.73233, 0.741818,
1792 0.745491, 0.741818, 0.73233, 0.717076, 0.695767,
1793 0.668836, 0.636821, 0.600068, 0.559255, 0.515621,
1794 0.470202, 0.424206, 0.37892, 0.335294, 0.29433,
1795 0.256551, 0.222141, 0.191333, 0.164191, 0.140308,
1796 0.119473, 0.101704, 0.0864095, 0.0733268, 0.0623139,
1797 0.0530119, 0.0451249, 0.0384801, 0.0328987, 0.0281427,
1798 0.0240831},
1799 {0.0206855, 0.0243011, 0.0285603, 0.0335902, 0.0396015,
1800 0.0467843, 0.0553424, 0.0655316, 0.0776688, 0.0921782,
1801 0.109269, 0.129259, 0.152459, 0.179369, 0.210091,
1802 0.244681, 0.283099, 0.325176, 0.370396, 0.417732,
1803 0.466086, 0.514192, 0.560518, 0.603805, 0.643028,
1804 0.67738, 0.706095, 0.728728, 0.745307, 0.755731,
1805 0.759954, 0.755731, 0.745307, 0.728728, 0.706095,
1806 0.67738, 0.643028, 0.603805, 0.560518, 0.514192,
1807 0.466086, 0.417732, 0.370396, 0.325176, 0.283099,
1808 0.244681, 0.210091, 0.179369, 0.152459, 0.129259,
1809 0.109269, 0.0921782, 0.0776688, 0.0655316, 0.0553424,
1810 0.0467843, 0.0396015, 0.0335902, 0.0285603, 0.0243011,
1811 0.0206855},
1812 {0.0186168, 0.0219999, 0.0260102, 0.0307769, 0.0364947,
1813 0.0433655, 0.0516213, 0.0615466, 0.0734611, 0.0877121,
1814 0.104666, 0.124855, 0.14853, 0.176033, 0.207639,
1815 0.243511, 0.283633, 0.327786, 0.37537, 0.425281,
1816 0.476227, 0.526727, 0.575268, 0.620462, 0.66101,
1817 0.69611, 0.725313, 0.748317, 0.764969, 0.775206,
1818 0.779006, 0.775206, 0.764969, 0.748317, 0.725313,
1819 0.69611, 0.66101, 0.620462, 0.575268, 0.526727,
1820 0.476227, 0.425281, 0.37537, 0.327786, 0.283633,
1821 0.243511, 0.207639, 0.176033, 0.14853, 0.124855,
1822 0.104666, 0.0877121, 0.0734611, 0.0615466, 0.0516213,
1823 0.0433655, 0.0364947, 0.0307769, 0.0260102, 0.0219999,
1824 0.0186168, },
1825 {0.0159737, 0.0189921, 0.0225916, 0.0268927, 0.0320634,
1826 0.0382995, 0.0458393, 0.0549765, 0.0660512, 0.0794439,
1827 0.095565, 0.114844, 0.137713, 0.164586, 0.195824,
1828 0.231681, 0.272223, 0.31727, 0.366248, 0.418078,
1829 0.471358, 0.524425, 0.575561, 0.623193, 0.666055,
1830 0.703243, 0.734192, 0.7586, 0.776331, 0.787347,
1831 0.791646, 0.787347, 0.776331, 0.7586, 0.734192,
1832 0.703243, 0.666055, 0.623193, 0.575561, 0.524425,
1833 0.471358, 0.418078, 0.366248, 0.31727, 0.272223,
1834 0.231681, 0.195824, 0.164586, 0.137713, 0.114844,
1835 0.095565, 0.0794439, 0.0660512, 0.0549765, 0.0458393,
1836 0.0382995, 0.0320634, 0.0268927, 0.0225916, 0.0189921,
1837 0.0159737},
1838 {0.0143532, 0.0171745, 0.0205622, 0.024635, 0.0295461,
1839 0.0354915, 0.0427198, 0.0515368, 0.0623058, 0.0754416,
1840 0.0913994, 0.110662, 0.133721, 0.161057, 0.193093,
1841 0.230138, 0.27229, 0.319327, 0.370597, 0.424961,
1842 0.480832, 0.536333, 0.589555, 0.638805, 0.682781,
1843 0.720621, 0.751845, 0.776243, 0.793764, 0.804426,
1844 0.808259, 0.804426, 0.793764, 0.776243, 0.751845,
1845 0.720621, 0.682781, 0.638805, 0.589555, 0.536333,
1846 0.480832, 0.424961, 0.370597, 0.319327, 0.27229,
1847 0.230138, 0.193093, 0.161057, 0.133721, 0.110662,
1848 0.0913994, 0.0754416, 0.0623058, 0.0515368, 0.0427198,
1849 0.0354915, 0.0295461, 0.024635, 0.0205622, 0.0171745,
1850 0.0143532} };
17b26de4 1851
1852 // More sampling precision with linear interpolation
1853 fPRFlo = -1.5;
1854 fPRFhi = 1.5;
1855 Float_t pad[kPRFbin];
1856 Int_t sPRFbin = kPRFbin;
1857 Float_t sPRFwid = (fPRFhi - fPRFlo) / ((Float_t) sPRFbin);
1858 for (Int_t iPad = 0; iPad < sPRFbin; iPad++) {
1859 pad[iPad] = ((Float_t) iPad + 0.5) * sPRFwid + fPRFlo;
1860 }
1861 fPRFbin = 500;
1862 fPRFwid = (fPRFhi - fPRFlo) / ((Float_t) fPRFbin);
1863 fPRFpad = ((Int_t) (1.0 / fPRFwid));
1864
1865 if (fPRFsmp) delete [] fPRFsmp;
1866 fPRFsmp = new Float_t[kNplan*fPRFbin];
1867
1868 Int_t ipos1;
1869 Int_t ipos2;
1870 Float_t diff;
1871
1872 for (Int_t iPla = 0; iPla < kNplan; iPla++) {
1873
1874 for (Int_t iBin = 0; iBin < fPRFbin; iBin++) {
1875
1876 Float_t bin = (((Float_t) iBin) + 0.5) * fPRFwid + fPRFlo;
1877 ipos1 = ipos2 = 0;
1878 diff = 0;
1879 do {
1880 diff = bin - pad[ipos2++];
8e64dd77 1881 } while ((diff > 0) && (ipos2 < kPRFbin));
1882 if (ipos2 == kPRFbin) {
1883 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2-1];
1884 }
1885 else if (ipos2 == 1) {
1886 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2-1];
1887 }
1888 else {
1889 ipos2--;
1890 if (ipos2 >= kPRFbin) ipos2 = kPRFbin - 1;
1891 ipos1 = ipos2 - 1;
1892 fPRFsmp[iPla*fPRFbin+iBin] = prf[iPla][ipos2]
1893 + diff * (prf[iPla][ipos2] - prf[iPla][ipos1])
1894 / sPRFwid;
1895 }
17b26de4 1896
1897 }
1898 }
1899
1900}
1901
1902//_____________________________________________________________________________
1903void AliTRDparameter::FillLUT()
1904{
1905 //
1906 // Create the LUT
1907 //
1908
1909 const Int_t kNplan = AliTRDgeometry::kNplan;
1910 const Int_t kNlut = 128;
1911
1912 fLUTbin = kNplan * kNlut;
1913
1914 // The lookup table from Bogdan
1915 Float_t lut[kNplan][kNlut] = {
1916 {
1917 0.0070, 0.0150, 0.0224, 0.0298, 0.0374, 0.0454, 0.0533, 0.0611,
1918 0.0684, 0.0755, 0.0827, 0.0900, 0.0975, 0.1049, 0.1120, 0.1187,
1919 0.1253, 0.1318, 0.1385, 0.1453, 0.1519, 0.1584, 0.1646, 0.1704,
1920 0.1762, 0.1821, 0.1879, 0.1938, 0.1996, 0.2053, 0.2108, 0.2160,
1921 0.2210, 0.2260, 0.2310, 0.2361, 0.2411, 0.2461, 0.2509, 0.2557,
1922 0.2602, 0.2646, 0.2689, 0.2732, 0.2774, 0.2816, 0.2859, 0.2901,
1923 0.2942, 0.2983, 0.3022, 0.3061, 0.3099, 0.3136, 0.3172, 0.3207,
1924 0.3242, 0.3278, 0.3312, 0.3347, 0.3382, 0.3416, 0.3450, 0.3483,
1925 0.3515, 0.3547, 0.3579, 0.3609, 0.3639, 0.3669, 0.3698, 0.3727,
1926 0.3756, 0.3785, 0.3813, 0.3842, 0.3870, 0.3898, 0.3926, 0.3952,
1927 0.3979, 0.4005, 0.4032, 0.4057, 0.4082, 0.4108, 0.4132, 0.4157,
1928 0.4181, 0.4205, 0.4228, 0.4252, 0.4275, 0.4299, 0.4322, 0.4345,
1929 0.4367, 0.4390, 0.4412, 0.4434, 0.4456, 0.4478, 0.4499, 0.4520,
1930 0.4541, 0.4562, 0.4583, 0.4603, 0.4623, 0.4643, 0.4663, 0.4683,
1931 0.4702, 0.4722, 0.4741, 0.4758, 0.4774, 0.4790, 0.4805, 0.4824,
1932 0.4844, 0.4863, 0.4883, 0.4902, 0.4921, 0.4940, 0.4959, 0.4978
1933 },
1934 {
1935 0.0072, 0.0156, 0.0235, 0.0313, 0.0394, 0.0478, 0.0561, 0.0642,
1936 0.0718, 0.0792, 0.0868, 0.0947, 0.1025, 0.1101, 0.1172, 0.1241,
1937 0.1309, 0.1378, 0.1449, 0.1518, 0.1586, 0.1650, 0.1710, 0.1770,
1938 0.1830, 0.1891, 0.1952, 0.2011, 0.2070, 0.2125, 0.2177, 0.2229,
1939 0.2280, 0.2332, 0.2383, 0.2435, 0.2484, 0.2533, 0.2581, 0.2627,
1940 0.2670, 0.2714, 0.2757, 0.2799, 0.2842, 0.2884, 0.2927, 0.2968,
1941 0.3008, 0.3048, 0.3086, 0.3123, 0.3159, 0.3195, 0.3231, 0.3266,
1942 0.3301, 0.3335, 0.3370, 0.3404, 0.3438, 0.3471, 0.3504, 0.3536,
1943 0.3567, 0.3598, 0.3628, 0.3657, 0.3686, 0.3715, 0.3744, 0.3772,
1944 0.3800, 0.3828, 0.3856, 0.3884, 0.3911, 0.3938, 0.3965, 0.3991,
1945 0.4016, 0.4042, 0.4067, 0.4092, 0.4116, 0.4140, 0.4164, 0.4187,
1946 0.4211, 0.4234, 0.4257, 0.4280, 0.4302, 0.4325, 0.4347, 0.4369,
1947 0.4391, 0.4413, 0.4434, 0.4456, 0.4477, 0.4497, 0.4518, 0.4538,
1948 0.4558, 0.4578, 0.4598, 0.4618, 0.4637, 0.4656, 0.4675, 0.4694,
1949 0.4713, 0.4732, 0.4750, 0.4766, 0.4781, 0.4797, 0.4813, 0.4832,
1950 0.4851, 0.4870, 0.4888, 0.4906, 0.4925, 0.4942, 0.4960, 0.4978
1951 },
1952 {
1953 0.0075, 0.0163, 0.0246, 0.0328, 0.0415, 0.0504, 0.0592, 0.0674,
1954 0.0753, 0.0832, 0.0914, 0.0996, 0.1077, 0.1154, 0.1225, 0.1296,
1955 0.1369, 0.1442, 0.1515, 0.1585, 0.1652, 0.1714, 0.1776, 0.1839,
1956 0.1902, 0.1965, 0.2025, 0.2085, 0.2141, 0.2194, 0.2247, 0.2299,
1957 0.2352, 0.2405, 0.2457, 0.2507, 0.2557, 0.2604, 0.2649, 0.2693,
1958 0.2737, 0.2780, 0.2823, 0.2867, 0.2909, 0.2951, 0.2992, 0.3033,
1959 0.3072, 0.3110, 0.3146, 0.3182, 0.3218, 0.3253, 0.3288, 0.3323,
1960 0.3357, 0.3392, 0.3426, 0.3459, 0.3492, 0.3524, 0.3555, 0.3586,
1961 0.3616, 0.3645, 0.3674, 0.3703, 0.3731, 0.3759, 0.3787, 0.3815,
1962 0.3843, 0.3870, 0.3897, 0.3925, 0.3950, 0.3976, 0.4002, 0.4027,
1963 0.4052, 0.4076, 0.4101, 0.4124, 0.4148, 0.4171, 0.4194, 0.4217,
1964 0.4239, 0.4262, 0.4284, 0.4306, 0.4328, 0.4350, 0.4371, 0.4393,
1965 0.4414, 0.4435, 0.4455, 0.4476, 0.4496, 0.4516, 0.4536, 0.4555,
1966 0.4575, 0.4594, 0.4613, 0.4632, 0.4650, 0.4669, 0.4687, 0.4705,
1967 0.4723, 0.4741, 0.4758, 0.4773, 0.4789, 0.4804, 0.4821, 0.4839,
1968 0.4857, 0.4875, 0.4893, 0.4910, 0.4928, 0.4945, 0.4961, 0.4978
1969 },
1970 {
1971 0.0078, 0.0171, 0.0258, 0.0345, 0.0438, 0.0532, 0.0624, 0.0708,
1972 0.0791, 0.0875, 0.0962, 0.1048, 0.1130, 0.1206, 0.1281, 0.1356,
1973 0.1432, 0.1508, 0.1582, 0.1651, 0.1716, 0.1780, 0.1845, 0.1910,
1974 0.1975, 0.2038, 0.2099, 0.2155, 0.2210, 0.2263, 0.2317, 0.2371,
1975 0.2425, 0.2477, 0.2528, 0.2578, 0.2626, 0.2671, 0.2715, 0.2759,
1976 0.2803, 0.2846, 0.2890, 0.2933, 0.2975, 0.3016, 0.3056, 0.3095,
1977 0.3132, 0.3168, 0.3204, 0.3239, 0.3274, 0.3309, 0.3344, 0.3378,
1978 0.3412, 0.3446, 0.3479, 0.3511, 0.3543, 0.3574, 0.3603, 0.3633,
1979 0.3662, 0.3690, 0.3718, 0.3747, 0.3774, 0.3802, 0.3829, 0.3857,
1980 0.3883, 0.3910, 0.3936, 0.3962, 0.3987, 0.4012, 0.4037, 0.4061,
1981 0.4085, 0.4109, 0.4132, 0.4155, 0.4177, 0.4200, 0.4222, 0.4244,
1982 0.4266, 0.4288, 0.4309, 0.4331, 0.4352, 0.4373, 0.4394, 0.4414,
1983 0.4435, 0.4455, 0.4475, 0.4494, 0.4514, 0.4533, 0.4552, 0.4571,
1984 0.4590, 0.4608, 0.4626, 0.4645, 0.4662, 0.4680, 0.4698, 0.4715,
1985 0.4733, 0.4750, 0.4766, 0.4781, 0.4796, 0.4812, 0.4829, 0.4846,
1986 0.4863, 0.4880, 0.4897, 0.4914, 0.4930, 0.4946, 0.4963, 0.4979
1987 },
1988 {
1989 0.0081, 0.0178, 0.0270, 0.0364, 0.0463, 0.0562, 0.0656, 0.0744,
1990 0.0831, 0.0921, 0.1013, 0.1102, 0.1183, 0.1261, 0.1339, 0.1419,
1991 0.1499, 0.1576, 0.1648, 0.1715, 0.1782, 0.1849, 0.1917, 0.1984,
1992 0.2048, 0.2110, 0.2167, 0.2223, 0.2278, 0.2333, 0.2389, 0.2444,
1993 0.2497, 0.2548, 0.2598, 0.2645, 0.2691, 0.2735, 0.2780, 0.2824,
1994 0.2868, 0.2912, 0.2955, 0.2997, 0.3038, 0.3078, 0.3116, 0.3152,
1995 0.3188, 0.3224, 0.3259, 0.3294, 0.3329, 0.3364, 0.3398, 0.3432,
1996 0.3465, 0.3497, 0.3529, 0.3561, 0.3591, 0.3620, 0.3649, 0.3677,
1997 0.3705, 0.3733, 0.3761, 0.3788, 0.3816, 0.3843, 0.3869, 0.3896,
1998 0.3922, 0.3948, 0.3973, 0.3998, 0.4022, 0.4047, 0.4070, 0.4094,
1999 0.4117, 0.4139, 0.4162, 0.4184, 0.4206, 0.4227, 0.4249, 0.4270,
2000 0.4291, 0.4313, 0.4334, 0.4354, 0.4375, 0.4395, 0.4415, 0.4435,
2001 0.4455, 0.4474, 0.4493, 0.4512, 0.4531, 0.4550, 0.4568, 0.4586,
2002 0.4604, 0.4622, 0.4639, 0.4657, 0.4674, 0.4691, 0.4708, 0.4725,
2003 0.4742, 0.4758, 0.4773, 0.4788, 0.4803, 0.4819, 0.4836, 0.4852,
2004 0.4869, 0.4885, 0.4901, 0.4917, 0.4933, 0.4948, 0.4964, 0.4979
2005 },
2006 {
2007 0.0085, 0.0189, 0.0288, 0.0389, 0.0497, 0.0603, 0.0699, 0.0792,
2008 0.0887, 0.0985, 0.1082, 0.1170, 0.1253, 0.1336, 0.1421, 0.1505,
2009 0.1587, 0.1662, 0.1733, 0.1803, 0.1874, 0.1945, 0.2014, 0.2081,
2010 0.2143, 0.2201, 0.2259, 0.2316, 0.2374, 0.2431, 0.2487, 0.2541,
2011 0.2593, 0.2642, 0.2689, 0.2735, 0.2781, 0.2826, 0.2872, 0.2917,
2012 0.2961, 0.3003, 0.3045, 0.3086, 0.3125, 0.3162, 0.3198, 0.3235,
2013 0.3270, 0.3306, 0.3342, 0.3377, 0.3411, 0.3446, 0.3479, 0.3511,
2014 0.3543, 0.3575, 0.3605, 0.3634, 0.3663, 0.3691, 0.3720, 0.3748,
2015 0.3775, 0.3803, 0.3830, 0.3857, 0.3884, 0.3911, 0.3937, 0.3962,
2016 0.3987, 0.4012, 0.4036, 0.4060, 0.4084, 0.4107, 0.4129, 0.4152,
2017 0.4174, 0.4196, 0.4218, 0.4239, 0.4261, 0.4282, 0.4303, 0.4324,
2018 0.4344, 0.4365, 0.4385, 0.4405, 0.4425, 0.4445, 0.4464, 0.4483,
2019 0.4502, 0.4521, 0.4539, 0.4558, 0.4576, 0.4593, 0.4611, 0.4629,
2020 0.4646, 0.4663, 0.4680, 0.4697, 0.4714, 0.4730, 0.4747, 0.4759,
2021 0.4769, 0.4780, 0.4790, 0.4800, 0.4811, 0.4827, 0.4843, 0.4859,
2022 0.4874, 0.4889, 0.4905, 0.4920, 0.4935, 0.4950, 0.4965, 0.4979
2023 }
2024 };
2025
2026 if (fLUT) delete [] fLUT;
2027 fLUT = new Float_t[fLUTbin];
2028
2029 for (Int_t iplan = 0; iplan < kNplan; iplan++) {
2030 for (Int_t ilut = 0; ilut < kNlut; ilut++) {
cfed70d4 2031 fLUT[iplan*kNlut+ilut] = lut[iplan][ilut];
17b26de4 2032 }
2033 }
2034
2035}
2036
2037//_____________________________________________________________________________
2038void AliTRDparameter::SetTiltingAngle(Float_t v)
2039{
2040 //
2041 // Set the tilting angle for the readout pads
2042 //
2043
2044 fTiltingAngle = TMath::Tan(TMath::Pi()/180.0 * v);
2045
2046}
2047
2048//_____________________________________________________________________________
2049Float_t AliTRDparameter::GetTiltingAngle() const
2050{
2051 //
2052 // Get the tilting angle for the readout pads
2053 //
2054
ab35e265 2055 return 180.0 / TMath::Pi() * TMath::ATan(fTiltingAngle);
17b26de4 2056
2057}
2058
2059//_____________________________________________________________________________
2060Float_t AliTRDparameter::GetDiffusionL(Float_t vd, Float_t b)
2061{
2062 //
2063 // Returns the longitudinal diffusion coefficient for a given drift
2064 // velocity <vd> and a B-field <b> for Xe/CO2 (15%).
2065 // The values are according to a GARFIELD simulation.
2066 //
2067
2068 const Int_t kNb = 5;
2069 Float_t p0[kNb] = { 0.007440, 0.007493, 0.007513, 0.007672, 0.007831 };
2070 Float_t p1[kNb] = { 0.019252, 0.018912, 0.018636, 0.018012, 0.017343 };
2071 Float_t p2[kNb] = { -0.005042, -0.004926, -0.004867, -0.004650, -0.004424 };
2072 Float_t p3[kNb] = { 0.000195, 0.000189, 0.000195, 0.000182, 0.000169 };
2073
2074 Int_t ib = ((Int_t) (10 * (b - 0.15)));
2075 ib = TMath::Max( 0,ib);
2076 ib = TMath::Min(kNb,ib);
2077
2078 Float_t diff = p0[ib]
2079 + p1[ib] * vd
2080 + p2[ib] * vd*vd
2081 + p3[ib] * vd*vd*vd;
2082
2083 return diff;
2084
2085}
2086
2087//_____________________________________________________________________________
2088Float_t AliTRDparameter::GetDiffusionT(Float_t vd, Float_t b)
2089{
2090 //
2091 // Returns the transverse diffusion coefficient for a given drift
2092 // velocity <vd> and a B-field <b> for Xe/CO2 (15%).
2093 // The values are according to a GARFIELD simulation.
2094 //
2095
2096 const Int_t kNb = 5;
2097 Float_t p0[kNb] = { 0.009550, 0.009599, 0.009674, 0.009757, 0.009850 };
2098 Float_t p1[kNb] = { 0.006667, 0.006539, 0.006359, 0.006153, 0.005925 };
2099 Float_t p2[kNb] = { -0.000853, -0.000798, -0.000721, -0.000635, -0.000541 };
2100 Float_t p3[kNb] = { 0.000131, 0.000122, 0.000111, 0.000098, 0.000085 };
2101
2102 Int_t ib = ((Int_t) (10 * (b - 0.15)));
2103 ib = TMath::Max( 0,ib);
2104 ib = TMath::Min(kNb,ib);
2105
2106 Float_t diff = p0[ib]
2107 + p1[ib] * vd
2108 + p2[ib] * vd*vd
2109 + p3[ib] * vd*vd*vd;
2110
2111 return diff;
2112
2113}
2114
2115//_____________________________________________________________________________
2116Float_t AliTRDparameter::GetOmegaTau(Float_t vd, Float_t b)
2117{
2118 //
2119 // Returns omega*tau (tan(Lorentz-angle)) for a given drift velocity <vd>
2120 // and a B-field <b> for Xe/CO2 (15%).
2121 // The values are according to a GARFIELD simulation.
2122 //
2123
2124 const Int_t kNb = 5;
2125 Float_t p0[kNb] = { 0.004810, 0.007412, 0.010252, 0.013409, 0.016888 };
2126 Float_t p1[kNb] = { 0.054875, 0.081534, 0.107333, 0.131983, 0.155455 };
2127 Float_t p2[kNb] = { -0.008682, -0.012896, -0.016987, -0.020880, -0.024623 };
2128 Float_t p3[kNb] = { 0.000155, 0.000238, 0.000330, 0.000428, 0.000541 };
2129
2130 Int_t ib = ((Int_t) (10 * (b - 0.15)));
2131 ib = TMath::Max( 0,ib);
2132 ib = TMath::Min(kNb,ib);
2133
2134 Float_t alphaL = p0[ib]
2135 + p1[ib] * vd
2136 + p2[ib] * vd*vd
2137 + p3[ib] * vd*vd*vd;
2138
2139 return TMath::Tan(alphaL);
2140
2141}
2142
2143//_____________________________________________________________________________
0a29d0f1 2144Double_t AliTRDparameter::LUTposition(Int_t iplane, Double_t ampL
2145 , Double_t ampC
2146 , Double_t ampR) const
17b26de4 2147{
2148 //
2149 // Calculates the cluster position using the lookup table.
2150 // Method provided by Bogdan Vulpescu.
2151 //
2152
2153 const Int_t kNplan = AliTRDgeometry::kNplan;
cfed70d4 2154 const Int_t kNlut = 128;
17b26de4 2155
2156 Double_t pos;
2157 Double_t x = 0.0;
2158 Double_t xmin;
2159 Double_t xmax;
2160 Double_t xwid;
2161
2162 Int_t side = 0;
2163 Int_t ix;
2164
2165 Double_t xMin[kNplan] = { 0.006492, 0.006377, 0.006258, 0.006144, 0.006030, 0.005980 };
2166 Double_t xMax[kNplan] = { 0.960351, 0.965870, 0.970445, 0.974352, 0.977667, 0.996101 };
2167
2168 if (ampL > ampR) {
2169 x = (ampL - ampR) / ampC;
2170 side = -1;
2171 }
2172 else if (ampL < ampR) {
2173 x = (ampR - ampL) / ampC;
2174 side = +1;
2175 }
2176
2177 if (ampL != ampR) {
2178
2179 xmin = xMin[iplane] + 0.000005;
2180 xmax = xMax[iplane] - 0.000005;
2181 xwid = (xmax - xmin) / 127.0;
2182
2183 if (x < xmin) {
2184 pos = 0.0000;
2185 }
2186 else if (x > xmax) {
2187 pos = side * 0.5000;
2188 }
2189 else {
2190 ix = (Int_t) ((x - xmin) / xwid);
cfed70d4 2191 pos = side * fLUT[iplane*kNlut+ix];
17b26de4 2192 }
2193
2194 }
2195 else {
2196
2197 pos = 0.0;
2198
2199 }
2200
2201 return pos;
2202
2203}