MultiAODInputHandler: an input handler for aod event muxing.

[u/mrichter/AliRoot.git] / TRD / AliTRDSimParam.cxx

diff --git a/TRD/AliTRDSimParam.cxx b/TRD/AliTRDSimParam.cxx
index 2b30b366324b39bb058f0d947e5a59186d4d0309..d6f0c3552ceaae98b76f94da61b870c11f00a9a3 100644 (file)
--- a/TRD/AliTRDSimParam.cxx
+++ b/TRD/AliTRDSimParam.cxx
@@ -82,7 +82,6 @@ AliTRDSimParam::AliTRDSimParam()
   ,fChipGain(0.0)
   ,fADCoutRange(0.0)
   ,fADCinRange(0.0)
-  ,fADCthreshold(0)
   ,fADCbaseline(0)
   ,fDiffusionOn(kFALSE)
   ,fElAttachOn(kFALSE)
@@ -122,9 +121,8 @@ void AliTRDSimParam::Init()
   fNoise             = 1250.0;
   fADCoutRange       = 1023.0;          // 10-bit ADC
   fADCinRange        = 2000.0;          // 2V input range
-  fADCthreshold      =    3;
   // Go back to 0 again, just to be consistent with reconstruction
-  fADCbaseline       =    0;
+  fADCbaseline       =      0;
   //fADCbaseline       =   10;
 
   // Diffusion on
@@ -188,7 +186,6 @@ AliTRDSimParam::AliTRDSimParam(const AliTRDSimParam &p)
   ,fChipGain(p.fChipGain)
   ,fADCoutRange(p.fADCoutRange)
   ,fADCinRange(p.fADCinRange)
-  ,fADCthreshold(p.fADCthreshold)
   ,fADCbaseline(p.fADCbaseline)
   ,fDiffusionOn(p.fDiffusionOn)
   ,fElAttachOn(p.fElAttachOn)
@@ -263,7 +260,6 @@ void AliTRDSimParam::Copy(TObject &p) const
   target->fChipGain           = fChipGain;  
   target->fADCoutRange        = fADCoutRange;
   target->fADCinRange         = fADCinRange;
-  target->fADCthreshold       = fADCthreshold;
   target->fADCbaseline        = fADCbaseline; 
   target->fDiffusionOn        = fDiffusionOn; 
   target->fElAttachOn         = fElAttachOn;
@@ -358,7 +354,27 @@ void AliTRDSimParam::SampleTRF()
                         , 0.0094, 0.0092, 0.0091, 0.0089, 0.0088, 0.0086, 0.0084
                         , 0.0083, 0.0081, 0.0080, 0.0078 };
 
-  Float_t xtalk[kNpasa];
+//     // Andronic & Bercuci parametrization
+//     // define new TRF parametrization
+//             // normalizing constant to Fe signal
+//     const Float_t k1=1.055;
+//     // time constants
+//     const Float_t t1=0.04;
+//     const Float_t t2=.9;
+//     // the relative fraction of the long component
+//     const Float_t k2=.15;
+//     // time offset for Fe
+//     const Float_t t0=-.29;
+//     Float_t x = t0; Int_t index;
+//     for(int i=0; i<kNpasa; i++) signal[i] = 0.;
+//   for(int i=0; i<kNpasa; i++){
+//             index = i+6;
+//             if(index >= kNpasa) break;
+//             x += .02;
+//             signal[index]=k1*(TMath::Power((x-t0)/t1, 2.5)*(exp(-(x-t0)/t1))+k2*exp(-(x-t0)/t2));
+//     }
+
+       Float_t xtalk[kNpasa];
 
   // With undershoot, positive peak corresponds to ~3% of the main signal:
   for (Int_t ipasa = 3; ipasa < kNpasa; ipasa++) {