New TRF parametrization as a workaround solution

diff --git a/TRD/AliTRDSimParam.cxx b/TRD/AliTRDSimParam.cxx
index 2b30b366324b39bb058f0d947e5a59186d4d0309..123208f15bddf711925dab1d9d3117e91a9c33e7 100644 (file)
--- a/TRD/AliTRDSimParam.cxx
+++ b/TRD/AliTRDSimParam.cxx
@@ -358,7 +358,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*(pow((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++) {