One header file for PYQUEN only.

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

diff --git a/TRD/AliTRDSimParam.cxx b/TRD/AliTRDSimParam.cxx
index 123208f15bddf711925dab1d9d3117e91a9c33e7..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,25 +354,25 @@ 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 };
 
-       // 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*(pow((x-t0)/t1, 2.5)*(exp(-(x-t0)/t1))+k2*exp(-(x-t0)/t2));
-       }
+//     // 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];