From 102f1568aabbe92a659d5a0267f73f18d60cd34d Mon Sep 17 00:00:00 2001
From: cblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Thu, 4 Jun 2009 09:47:51 +0000
Subject: [PATCH] New QA checking procedure by Sylwester

---
 TRD/AliTRDQAChecker.cxx | 43 ++++++++++++++++++++++++++++++++++++-----
 TRD/AliTRDQAChecker.h   | 17 ++++++++--------
 2 files changed, 47 insertions(+), 13 deletions(-)

diff --git a/TRD/AliTRDQAChecker.cxx b/TRD/AliTRDQAChecker.cxx
index caf703bbc5f..14ae0b68c4b 100644
--- a/TRD/AliTRDQAChecker.cxx
+++ b/TRD/AliTRDQAChecker.cxx
@@ -13,11 +13,15 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
-/*
-  Checks the quality assurance. 
-  By comparing with reference data
-  S.Radomski Uni-Heidelberg October 2007
-*/
+/* $Id$ */
+
+////////////////////////////////////////////////////////////////////////////
+//                                                                        //
+//  Checks the quality assurance.                                         //
+//  By comparing with reference data                                      //
+//  S.Radomski Uni-Heidelberg October 2007                                //
+//                                                                        //
+////////////////////////////////////////////////////////////////////////////
 
 // --- ROOT system ---
 #include <TClass.h>
@@ -38,3 +42,32 @@
 ClassImp(AliTRDQAChecker)
 
 //__________________________________________________________________
+
+Double_t * AliTRDQAChecker::Check(AliQAv1::ALITASK_t index, TObjArray ** list) 
+{
+
+  // Super-basic check on the QA histograms on the input list: 
+
+  Double_t * test  = new Double_t[AliRecoParam::kNSpecies];
+  for(Int_t i=0; i<AliRecoParam::kNSpecies; i++) test[i] = 0.5; 
+
+  //Int_t count[AliRecoParam::kNSpecies] = { 0 }; 
+
+  if (index != AliQAv1::kRECPOINTS) return test;
+
+  const Double_t lowAmp = 30;
+  const Double_t highAmp = 50;
+
+  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
+    
+    TH1D *hist = (TH1D*) list[specie]->At(12);
+    if (!hist) continue;
+    
+    Double_t value = hist->Integral(hist->FindBin(lowAmp), hist->FindBin(highAmp));
+    test[specie] = value / hist->GetSum();
+
+  }
+  return test ; 
+}  
+
+//____________________________________________________________________________
diff --git a/TRD/AliTRDQAChecker.h b/TRD/AliTRDQAChecker.h
index d7fa324b1c1..84be75aff3d 100644
--- a/TRD/AliTRDQAChecker.h
+++ b/TRD/AliTRDQAChecker.h
@@ -3,14 +3,15 @@
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 
-//////////////////////////////////////////////////////
-//
-//  Checks the quality assurance. 
-//  By comparing with reference data
-//  S. Radomski Uni-Heidelberg October 2007
-//
-///////////////////////////////////////////////////////
+/* $Id$ */
 
+////////////////////////////////////////////////////////////////////////////
+//                                                                        //
+//  Checks the quality assurance.                                         //
+//  By comparing with reference data                                      //
+//  S.Radomski Uni-Heidelberg October 2007                                //
+//                                                                        //
+////////////////////////////////////////////////////////////////////////////
 
 // --- ROOT system ---
 class TFile ; 
@@ -35,7 +36,7 @@ public:
 
   virtual Double_t * Check(AliQAv1::ALITASK_t /*index*/) {return NULL;}
   virtual Double_t * Check(TList * /*list*/) {return NULL;}
-  virtual Double_t * Check(AliQAv1::ALITASK_t /*index*/, TObjArray ** /*list*/) {return NULL;}
+  virtual Double_t * Check(AliQAv1::ALITASK_t /*index*/, TObjArray ** /*list*/);
   virtual Double_t * Check(AliQAv1::ALITASK_t /*index*/, TNtupleD** /*nt*/)     {return NULL;}
 
 private:
-- 
2.39.3