TRD/AliTRDparameter.cxx

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