TRD/AliTRDSimParam.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 // Class containing constant simulation parameters                           //
  21 //                                                                           //
  22 // Request an instance with AliTRDSimParam::Instance()                 //
  23 // Then request the needed values                                            //
  24 //                                                                           //
  25 ///////////////////////////////////////////////////////////////////////////////
  26
  27 #include <TMath.h>
  28
  29 #include "AliRun.h"
  30
  31 #include "AliTRDSimParam.h"
  32
  33 ClassImp(AliTRDSimParam)
  34
  35 AliTRDSimParam* AliTRDSimParam::fgInstance = 0;
  36 Bool_t AliTRDSimParam::fgTerminated = kFALSE;
  37
  38 //_ singleton implementation __________________________________________________
  39 AliTRDSimParam* AliTRDSimParam::Instance()
  40 {
  41   //
  42   // Singleton implementation
  43   // Returns an instance of this class, it is created if neccessary
  44   //
  45
  46   if (fgTerminated != kFALSE)
  47     return 0;
  48
  49   if (fgInstance == 0)
  50     fgInstance = new AliTRDSimParam();
  51
  52   return fgInstance;
  53 }
  54
  55 void AliTRDSimParam::Terminate()
  56 {
  57   //
  58   // Singleton implementation
  59   // Deletes the instance of this class and sets the terminated flag, instances cannot be requested anymore
  60   // This function can be called several times.
  61   //
  62
  63   fgTerminated = kTRUE;
  64
  65   if (fgInstance != 0)
  66   {
  67     delete fgInstance;
  68     fgInstance = 0;
  69   }
  70 }
  71
  72 //_____________________________________________________________________________
  73 AliTRDSimParam::AliTRDSimParam()
  74 {
  75   //
  76   // default constructor
  77   //
  78
  79   fGasGain            = 0.0;
  80   fNoise              = 0.0;
  81   fChipGain           = 0.0;
  82
  83   fADCoutRange        = 0.0;
  84   fADCinRange         = 0.0;
  85   fADCthreshold       = 0;
  86   fADCbaseline        = 0;
  87
  88   fDiffusionOn        = kFALSE;
  89
  90   fElAttachOn         = kFALSE;
  91   fElAttachProp       = 0.0;
  92
  93   fTRFOn              = kFALSE;
  94   fTRFsmp             = 0;
  95   fTRFbin             = 0;
  96   fTRFlo              = 0.0;
  97   fTRFhi              = 0.0;
  98   fTRFwid             = 0.0;
  99
 100   fCTOn               = kFALSE;
 101   fCTsmp              = 0;
 102
 103   fAnodeWireOffset    = 0.0;
 104   fPadCoupling        = 0.0;
 105   fTimeCoupling       = 0.0;
 106   fTimeStructOn       = kFALSE;
 107
 108   fPRFOn              = kFALSE;
 109
 110   Init();
 111 }
 112
 113 //_____________________________________________________________________________
 114 void AliTRDSimParam::Init()
 115 {
 116   //
 117   // default initializiation
 118   //
 119
 120   // The default parameter for the digitization
 121   fGasGain        = 4000.;
 122   fChipGain       = 12.4;
 123   fNoise          = 1000.;
 124   fADCoutRange    = 1023.;          // 10-bit ADC
 125   fADCinRange     = 2000.;          // 2V input range
 126   fADCthreshold   = 1;
 127   fADCbaseline    = 0;
 128
 129   // Diffusion on
 130   fDiffusionOn    = kTRUE;
 131
 132   // Propability for electron attachment
 133   fElAttachOn     = kFALSE;
 134   fElAttachProp   = 0.0;
 135
 136   // The time response function
 137   fTRFOn          = kTRUE;
 138
 139   // The cross talk
 140   fCTOn           = kTRUE;
 141
 142   // The pad coupling factor
 143   //fPadCoupling    = 0.3;
 144   // Use 0.46 instead which reproduces better the test beam
 145   // data, even tough it is not understood why.
 146   fPadCoupling    = 0.46;
 147
 148   // The time coupling factor (same number as for the TPC)
 149   fTimeCoupling   = 0.4;
 150
 151   // Distance of first Anode wire from first pad edge
 152   fAnodeWireOffset = 0.25;
 153
 154   // Use drift time maps
 155   fTimeStructOn = kTRUE;
 156
 157   // The pad response function
 158   fPRFOn          = kTRUE;
 159
 160   ReInit();
 161 }
 162
 163 //_____________________________________________________________________________
 164 AliTRDSimParam::~AliTRDSimParam()
 165 {
 166   //
 167   // destructor
 168   //
 169
 170   if (fTRFsmp) {
 171     delete [] fTRFsmp;
 172     fTRFsmp = 0;
 173   }
 174
 175   if (fCTsmp) {
 176     delete [] fCTsmp;
 177     fCTsmp  = 0;
 178   }
 179 }
 180
 181 //_____________________________________________________________________________
 182 AliTRDSimParam::AliTRDSimParam(const AliTRDSimParam &p):TObject(p)
 183 {
 184   //
 185   // copy constructor
 186   //
 187
 188   ((AliTRDSimParam &) p).Copy(*this);
 189 }
 190
 191
 192 //_____________________________________________________________________________
 193 AliTRDSimParam &AliTRDSimParam::operator=(const AliTRDSimParam &p)
 194 {
 195   //
 196   // Assignment operator
 197   //
 198
 199   if (this != &p) ((AliTRDSimParam &) p).Copy(*this);
 200   return *this;
 201 }
 202
 203 //_____________________________________________________________________________
 204 void AliTRDSimParam::Copy(TObject &p) const
 205 {
 206   //
 207   // Copy function
 208   //
 209
 210   AliTRDSimParam* target = dynamic_cast<AliTRDSimParam*> (&p);
 211   if (!target)
 212     return;
 213
 214   target->fGasGain            = fGasGain;
 215   //target->fField              = fField;
 216   target->fNoise              = fNoise;
 217   target->fChipGain           = fChipGain;
 218
 219   target->fADCoutRange        = fADCoutRange;
 220   target->fADCinRange         = fADCinRange;
 221   target->fADCthreshold       = fADCthreshold;
 222   target->fADCbaseline        = fADCbaseline;
 223
 224   target->fDiffusionOn        = fDiffusionOn;
 225
 226   target->fElAttachOn         = fElAttachOn;
 227   target->fElAttachProp       = fElAttachProp;
 228
 229   target->fTRFOn              = fTRFOn;
 230   if (target->fTRFsmp)
 231     delete[] target->fTRFsmp;
 232   target->fTRFsmp = new Float_t[fTRFbin];
 233   for (Int_t iBin = 0; iBin < fTRFbin; iBin++) {
 234     target->fTRFsmp[iBin] = fTRFsmp[iBin];
 235   }
 236   target->fTRFbin             = fTRFbin;
 237   target->fTRFlo              = fTRFlo;
 238   target->fTRFhi              = fTRFhi;
 239   target->fTRFwid             = fTRFwid;
 240
 241   target->fCTOn               = fCTOn;
 242   if (target->fCTsmp)
 243     delete[] target->fCTsmp;
 244   target->fCTsmp  = new Float_t[fTRFbin];
 245   for (Int_t iBin = 0; iBin < fTRFbin; iBin++) {
 246     target->fCTsmp[iBin]  = fCTsmp[iBin];
 247   }
 248
 249   target->fAnodeWireOffset    = fAnodeWireOffset;
 250   target->fPadCoupling        = fPadCoupling;
 251   target->fTimeCoupling       = fTimeCoupling;
 252
 253   target->fPRFOn              = fPRFOn;
 254 }
 255
 256 //_____________________________________________________________________________
 257 void AliTRDSimParam::ReInit()
 258 {
 259   //
 260   // Reinitializes the parameter class after a change
 261   //
 262
 263   // The range and the binwidth for the sampled TRF
 264   fTRFbin = 200;
 265   // Start 0.2 mus before the signal
 266   fTRFlo  = -0.4;
 267   // End the maximum drift time after the signal
 268   fTRFhi  = 3.58;
 269   //
 270   fTRFwid = (fTRFhi - fTRFlo) / ((Float_t) fTRFbin);
 271
 272   // Create the sampled TRF
 273   SampleTRF();
 274 }
 275
 276
 277 //_____________________________________________________________________________
 278 void AliTRDSimParam::SampleTRF()
 279 {
 280   //
 281   // Samples the new time response function.
 282   // From Antons measurements with Fe55 source, adjusted by C. Lippmann.
 283   // time bins are -0.4, -0.38, -0.36, ...., 3.54, 3.56, 3.58 microseconds
 284   //
 285
 286   const Int_t kNpasa     = 200;  // kNpasa should be equal to fTRFbin!
 287
 288   Float_t signal[kNpasa]={0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000,
 289                           0.0002, 0.0007, 0.0026, 0.0089, 0.0253, 0.0612, 0.1319,
 290                           0.2416, 0.3913, 0.5609, 0.7295, 0.8662, 0.9581, 1.0000,
 291                           0.9990, 0.9611, 0.8995, 0.8269, 0.7495, 0.6714, 0.5987,
 292                           0.5334, 0.4756, 0.4249, 0.3811, 0.3433, 0.3110, 0.2837,
 293                           0.2607, 0.2409, 0.2243, 0.2099, 0.1974, 0.1868, 0.1776,
 294                           0.1695, 0.1627, 0.1566, 0.1509, 0.1457, 0.1407, 0.1362,
 295                           0.1317, 0.1274, 0.1233, 0.1196, 0.1162, 0.1131, 0.1102,
 296                           0.1075, 0.1051, 0.1026, 0.1004, 0.0979, 0.0956, 0.0934,
 297                           0.0912, 0.0892, 0.0875, 0.0858, 0.0843, 0.0829, 0.0815,
 298                           0.0799, 0.0786, 0.0772, 0.0757, 0.0741, 0.0729, 0.0718,
 299                           0.0706, 0.0692, 0.0680, 0.0669, 0.0655, 0.0643, 0.0630,
 300                           0.0618, 0.0607, 0.0596, 0.0587, 0.0576, 0.0568, 0.0558,
 301                           0.0550, 0.0541, 0.0531, 0.0522, 0.0513, 0.0505, 0.0497,
 302                           0.0490, 0.0484, 0.0474, 0.0465, 0.0457, 0.0449, 0.0441,
 303                           0.0433, 0.0425, 0.0417, 0.0410, 0.0402, 0.0395, 0.0388,
 304                           0.0381, 0.0374, 0.0368, 0.0361, 0.0354, 0.0348, 0.0342,
 305                           0.0336, 0.0330, 0.0324, 0.0318, 0.0312, 0.0306, 0.0301,
 306                           0.0296, 0.0290, 0.0285, 0.0280, 0.0275, 0.0270, 0.0265,
 307                           0.0260, 0.0256, 0.0251, 0.0246, 0.0242, 0.0238, 0.0233,
 308                           0.0229, 0.0225, 0.0221, 0.0217, 0.0213, 0.0209, 0.0206,
 309                           0.0202, 0.0198, 0.0195, 0.0191, 0.0188, 0.0184, 0.0181,
 310                           0.0178, 0.0175, 0.0171, 0.0168, 0.0165, 0.0162, 0.0159,
 311                           0.0157, 0.0154, 0.0151, 0.0148, 0.0146, 0.0143, 0.0140,
 312                           0.0138, 0.0135, 0.0133, 0.0131, 0.0128, 0.0126, 0.0124,
 313                           0.0121, 0.0119, 0.0120, 0.0115, 0.0113, 0.0111, 0.0109,
 314                           0.0107, 0.0105, 0.0103, 0.0101, 0.0100, 0.0098, 0.0096,
 315                           0.0094, 0.0092, 0.0091, 0.0089, 0.0088, 0.0086, 0.0084,
 316                           0.0083, 0.0081, 0.0080, 0.0078};
 317
 318   Float_t xtalk[kNpasa];
 319
 320   for (Int_t ipasa = 3; ipasa < kNpasa; ipasa++) {
 321     //signal[ipasa] /= 0.44; // normalise area to 1
 322     // With undershoot, positive peak corresponds to ~3% of the main signal:
 323     xtalk[ipasa] = 0.2*(signal[ipasa-2]-signal[ipasa-3]);
 324   }
 325
 326   xtalk[0] =  0.0;      signal[0] = 0.0;
 327   xtalk[1] =  0.0;      signal[1] = 0.0;
 328   xtalk[2] =  0.0;      signal[2] = 0.0;
 329
 330   if (fTRFsmp) delete [] fTRFsmp;
 331   fTRFsmp = new Float_t[fTRFbin];
 332   if (fCTsmp)  delete [] fCTsmp;
 333   fCTsmp  = new Float_t[fTRFbin];
 334
 335   for (Int_t iBin = 0; iBin < fTRFbin; iBin++) {
 336     fTRFsmp[iBin] = signal[iBin];
 337     fCTsmp[iBin]  = xtalk[iBin];
 338   }
 339
 340 }
 341
 342 //_____________________________________________________________________________
 343 Double_t AliTRDSimParam::TimeResponse(Double_t time) const
 344 {
 345   //
 346   // Applies the preamp shaper time response
 347   // (We assume a signal rise time of 0.2us = fTRFlo/2.
 348   //
 349
 350   Int_t iBin = ((Int_t) ((time - fTRFlo/2.) / fTRFwid));
 351   if ((iBin >= 0) && (iBin < fTRFbin)) {
 352     return fTRFsmp[iBin];
 353   }
 354   else {
 355     return 0.0;
 356   }
 357
 358 }
 359
 360 //_____________________________________________________________________________
 361 Double_t AliTRDSimParam::CrossTalk(Double_t time) const
 362 {
 363   //
 364   // Applies the pad-pad capacitive cross talk
 365   //
 366
 367   Int_t iBin = ((Int_t) ((time - fTRFlo) / fTRFwid));
 368   if ((iBin >= 0) && (iBin < fTRFbin)) {
 369     return fCTsmp[iBin];
 370   }
 371   else {
 372     return 0.0;
 373   }
 374
 375 }
 376