From 0929ca8ea83ae66b21741589083a11bd70640d32 Mon Sep 17 00:00:00 2001
From: cblume <Christoph.Blume@cern.ch>
Date: Sun, 10 Jul 2016 13:01:45 +0200
Subject: [PATCH] bugfix by Raphaelle

See ALIROOT-6776
---
 TRD/TRDbase/AliTRDCalDet.cxx     | 36 ++++++++++++++++++++++++++++++++
 TRD/TRDbase/AliTRDCalDet.h       |  1 +
 TRD/TRDbase/AliTRDCalibraFit.cxx |  3 ++-
 3 files changed, 39 insertions(+), 1 deletion(-)

diff --git a/TRD/TRDbase/AliTRDCalDet.cxx b/TRD/TRDbase/AliTRDCalDet.cxx
index 537779a274e..375f7582321 100644
--- a/TRD/TRDbase/AliTRDCalDet.cxx
+++ b/TRD/TRDbase/AliTRDCalDet.cxx
@@ -207,6 +207,42 @@ Double_t AliTRDCalDet::GetRMSRobust(Double_t robust) const
   delete [] index;
   return mean;
 }
+//____________________________________________________________________________________________
+Double_t AliTRDCalDet::GetMeanRobust(Double_t robust) const
+{
+  //
+  // Calculate the robust mean
+  //
+
+  // sorted
+  Int_t *index = new Int_t[kNdet];
+  TMath::Sort((Int_t)kNdet,fData,index);
+ 
+  // reject
+  Double_t reject = (Int_t) (kNdet*(1-robust)/2.0);
+  if(reject <= 0.0) reject = 0.0;
+  if(reject >= kNdet) reject = 0.0;
+  //printf("Rejecter %f\n",reject);
+
+  Double_t *ddata = new Double_t[kNdet];
+  Int_t nPoints = 0;
+  for (Int_t i=0;i<kNdet;i++) {
+    Bool_t rej = kFALSE;
+    for(Int_t k = 0; k < reject; k++){
+      if(i==index[k]) rej = kTRUE;
+      if(i==index[kNdet-(k+1)]) rej  = kTRUE;
+    }
+    if(!rej){
+      ddata[nPoints]= fData[i];
+      nPoints++;
+    }
+  }
+  //printf("Number of points %d\n",nPoints);
+  Double_t mean = TMath::Mean(nPoints,ddata);
+  delete [] ddata;
+  delete [] index;
+  return mean;
+}
 
 //______________________________________________________________________________________________
 Double_t AliTRDCalDet::GetLTM(Double_t *sigma
diff --git a/TRD/TRDbase/AliTRDCalDet.h b/TRD/TRDbase/AliTRDCalDet.h
index 520e7ffec60..c68bdf61f00 100644
--- a/TRD/TRDbase/AliTRDCalDet.h
+++ b/TRD/TRDbase/AliTRDCalDet.h
@@ -42,6 +42,7 @@ class AliTRDCalDet : public TNamed {
 
   // statistic
   Double_t GetMean(AliTRDCalDet * const outlierDet=0) const;
+  Double_t GetMeanRobust(Double_t robust=0.92) const;
   Double_t GetRMS(AliTRDCalDet * const outlierDet=0) const;
   Double_t GetRMSRobust(Double_t robust=0.92) const;
   Double_t GetMedian(AliTRDCalDet * const outlierDet=0) const;
diff --git a/TRD/TRDbase/AliTRDCalibraFit.cxx b/TRD/TRDbase/AliTRDCalibraFit.cxx
index 8faef785409..77bce2cea3f 100644
--- a/TRD/TRDbase/AliTRDCalibraFit.cxx
+++ b/TRD/TRDbase/AliTRDCalibraFit.cxx
@@ -1239,7 +1239,7 @@ Bool_t AliTRDCalibraFit::AnalyseLinearFitters(AliTRDCalibraVdriftLinearFit *cali
     fStatisticMean += entriesCurrent;     
 
     // Check the fit
-    if((-(param[1])) <= 0.000001) {
+    if((-(param[1])) <= 0.01) {
       NotEnoughStatisticLinearFitter();
       continue;
     }
@@ -1255,6 +1255,7 @@ Bool_t AliTRDCalibraFit::AnalyseLinearFitters(AliTRDCalibraVdriftLinearFit *cali
     fCurrentCoef[0]  = -param[1];
     // here the database must be the one of the reconstruction for the lorentz angle....
     fCurrentCoef2[0] = (param[0]+fCurrentCoef[1]*fCurrentCoef2[1])/fCurrentCoef[0];
+    //printf("fCurrentCoef2[0] %f, param[0] %f, fCurrentCoef[1]*fCurrentCoef2[1] %f,fCurrentCoef[0] %f\n",fCurrentCoef2[0],param[0],fCurrentCoef[1]*fCurrentCoef2[1],fCurrentCoef[0]);
     fCurrentCoefE    = error[1];
     fCurrentCoefE2   = error[0];
     if((TMath::Abs(fCurrentCoef2[0]) > 0.0000001) && (TMath::Abs(param[0]) > 0.0000001)){