First implementation of EMCAL trigger QA from Nicola Arbor
[u/mrichter/AliRoot.git] / ZDC / AliZDCQAChecker.cxx
index ca19f47..c6a4a9c 100644 (file)
@@ -18,6 +18,9 @@
 #include <TH1F.h> 
 #include <TIterator.h> 
 #include <TString.h> 
+#include <TPaveText.h> 
+#include <TObjArray.h>
+#include "TList.h"
 
 // --- Standard library ---
 
 ClassImp(AliZDCQAChecker)
 
 //____________________________________________________________________________
-Double_t * AliZDCQAChecker::Check(AliQA::ALITASK_t index, TObjArray ** list) 
+void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArray ** list,
+      const AliDetectorRecoParam * /*recoParam*/) 
 {
   // Checks the QA histograms on the input list: 
   //
-  Double_t * test   = new Double_t[AliRecoParam::kNSpecies] ;
-  Int_t *    ntests = new Int_t[AliRecoParam::kNSpecies]  ; 
-  const char* taskName = AliQA::GetAliTaskName(index);
+  const char* taskName = AliQAv1::GetAliTaskName(index);
+  //printf("\n\tAliZDCQAChecker -> checking QA histos for task %s\n",taskName);
   //
-  
-  //YS Int_t beamType=0; // 0 -> protons, 1 -> ions
-  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
-    test[specie] = 1.0 ; 
-    ntests[specie] = 0 ; 
-    if ( !AliQA::Instance()->IsEventSpecieSet(specie) ) 
-      continue ; 
+  Int_t ihitHisto=0, idigHisto=0;
+  Int_t irecHisto=0, irawHisto=0, esdInd=0;
+
+  TObjArray messages;
+  messages.SetOwner(kTRUE);
+
+  for(Int_t specie = 0; specie<AliRecoParam::kNSpecies; specie++){
+    Int_t count = 0; 
+    if(!AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::ConvertIndex(specie))) 
+      continue ;
+    //printf("\tAliZDCQAChecker -> specie %d, AliRecoParam::ConvertIndex(specie) %d, AliRecoParam::kLowMult %d, IsEventSpecieSet(specie) %d\n",
+    //  specie, AliRecoParam::ConvertIndex(specie) ,AliRecoParam::kLowMult,
+    //  AliQAv1::Instance(AliQAv1::GetDetIndex(GetName()))->IsEventSpecieSet(AliRecoParam::ConvertIndex(specie)));
+    
     // ====================================================================
     //         Checks for p-p events
     // ====================================================================
-    //if (beamType==0){
-    if ( specie == AliRecoParam::kLowMult) {
+    if(AliRecoParam::ConvertIndex(specie) == AliRecoParam::kLowMult){
       if(list[specie]->GetEntries()==0){  
-        AliWarning("\tAliZDCQAChecker->The list to be checked is empty!"); // nothing to check
-        return test;
+        AliWarning("\t The list to be checked is empty!"); // nothing to check
+        return;
       }
-      //AliInfo(Form("\n\tAliZDCQAChecker-> checking QA histograms for task %s\n\n",taskName));
+      //AliDebug(AliQAv1::GetQADebugLevel(), Form("\n\tAliZDCQAChecker-> checking QA histograms for task %s\n\n",taskName));
       TIter next(list[specie]); 
-      ntests[specie] = 0; 
+      count = 0; 
       TH1 * hdata;       
+      
+      Float_t res=0., percentageDiff=0.30;
+      Float_t meanZNA=0., meanZNC=0., meanZPA=0., meanZPC=0.;
+      Float_t pmCZNA=0., pmCZNC=0., pmCZPA=0., pmCZPC=0.;
+      Float_t pmQZNA=0., pmQZNC=0., pmQZPA=0., pmQZPC=0.;
+      Float_t sumADCZNA=0., sumADCZNC=0., sumADCZPA=0., sumADCZPC=0.;
+      Float_t adcCZNA=0., adcCZNC=0., adcCZPA=0., adcCZPC=0.;
+      Float_t adcQZNA=0., adcQZNC=0., adcQZPA=0., adcQZPC=0.;
+      Int_t nentries = -99;
+      
       while((hdata = dynamic_cast<TH1 *>(next()))){
         if(hdata){ 
           // -------------------------------------------------------------------
-          if(index == AliQA::kSIM){
-            //AliInfo(Form("\tAliZDCQAChecker-> checking histo %s",hdata->GetName()));
-            // Check DIGITS histos
-            if(!(strncmp(hdata->GetName(),"hDig",4))){
-              if(hdata->GetEntries()>0){
-                test[specie] += 1.; 
-                ntests[specie]++;
-              }
-            }
+          if(index == AliQAv1::kSIM){
+            //AliDebug(AliQAv1::GetQADebugLevel(), Form("\tAliZDCQAChecker-> checking histo %s",hdata->GetName()));
             // Check HITS histos
-            else{ 
+            //
+           if(!(strncmp(hdata->GetName(),"hHits",5))){
               if(hdata->GetEntries()>0){
-                test[specie] += 1.; 
-                ntests[specie]++;
-              }
+               if(ihitHisto==0)      meanZNC = hdata->GetMean();
+               else if(ihitHisto==1) meanZNA = hdata->GetMean();
+               else if(ihitHisto==2) meanZPC = hdata->GetMean();
+               else if(ihitHisto==3) meanZPA = hdata->GetMean();
+               else if(ihitHisto==4) pmQZNC = hdata->GetMean();
+               else if(ihitHisto==5) pmQZNA = hdata->GetMean();
+               else if(ihitHisto==6) pmQZPC = hdata->GetMean();
+               else if(ihitHisto==7) pmQZPA = hdata->GetMean();
+               else if(ihitHisto==8)  pmCZNC = hdata->GetMean();
+               else if(ihitHisto==9)  pmCZNA = hdata->GetMean();
+               else if(ihitHisto==10) pmCZPC = hdata->GetMean();
+               else if(ihitHisto==11) pmCZPA = hdata->GetMean();
+             }
+             //
+             // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
+             if(ihitHisto==11){
+               if(TMath::Abs(meanZNC)>1.e-10){
+                  if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
+                    res=1.;
+                  else 
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
+               if(TMath::Abs(meanZNA)>1.e-10){
+                  if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
+                    res=1.;
+                  else percentageDiff=
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
+               if(TMath::Abs(meanZPC)>1.e-10){
+                  if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
+                    res=1.;
+                  else 
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
+               if(TMath::Abs(meanZPA)>1.e-10){
+                  if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
+                    res=1.;
+                  else 
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
+             }
+             ihitHisto++;
             }
-            // -------------------------------------------------------------------
-          } else if(index == AliQA::kRAW) {
-            if(hdata->GetEntries()!=0){
-              if(hdata->GetMean()>10.) 
-                test[specie] += 1.; 
-              else 
-                test[specie] = 0.5; 
-              ntests[specie]++;
+            // Check DIGIT HIGH GAIN CHAIN histos
+            else if(!(strncmp(hdata->GetName(),"hDig",4))){ 
+              if(hdata->GetEntries()>0){
+               if(idigHisto==0)      sumADCZNC = hdata->GetMean();
+               else if(idigHisto==1) sumADCZNA = hdata->GetMean();
+               else if(idigHisto==2) sumADCZPC = hdata->GetMean();
+               else if(idigHisto==3) sumADCZPA = hdata->GetMean();
+               else if(idigHisto==4) pmQZNC = hdata->GetMean();
+               else if(idigHisto==5) pmQZNA = hdata->GetMean();
+               else if(idigHisto==6) pmQZPC = hdata->GetMean();
+               else if(idigHisto==7) pmQZPA = hdata->GetMean();
+               else if(idigHisto==8)  pmCZNC = hdata->GetMean();
+               else if(idigHisto==9)  pmCZNA = hdata->GetMean();
+               else if(idigHisto==10) pmCZPC = hdata->GetMean();
+               else if(idigHisto==11) pmCZPA = hdata->GetMean();
+             }
+             //
+             // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
+             if(idigHisto==11){
+               if(TMath::Abs(sumADCZNC)>1.e-10){
+                  if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
+                    res=1.;
+                  else 
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
+               if(TMath::Abs(sumADCZNA)>1.e-10){
+                  if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
+                    res=1.;
+                  else 
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
+               if(TMath::Abs(sumADCZPC)>1.e-10){
+                  if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
+                    res=1.;
+                  else 
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
+               if(TMath::Abs(sumADCZPA)>1.e-10){
+                  if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
+                    res=1.;
+                  else 
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
+             }
+             idigHisto++;            
             }
-            // -------------------------------------------------------------------
-          } else if(index == AliQA::kESD) {
-            Int_t    esdInd=0;
-            if(hdata->GetEntries()!=0){
-              if(esdInd>1){
-                if(hdata->GetMean()>10.) 
-                  test[specie] += 1.; 
-                else 
-                  test[specie] = 0.5; 
-                ntests[specie]++;
-              }
+          } 
+          // -------------------------------------------------------------------
+         else if(index == AliQAv1::kRAW) {
+           //
+            // Check RAW HIGH GAIN CHAIN histos
+            
+           messages.Clear();
+           
+           Bool_t iDetPM = kTRUE;
+           // --- Checks
+           if(irawHisto==16) nentries = Int_t (hdata->GetEntries());
+           if(irawHisto==18){ 
+             Float_t resADC=0.;
+             for(int ibin=1; ibin<=hdata->GetNbinsX(); ibin++){
+                if((hdata->GetBinContent(ibin))>10.){
+                  res=1.;
+                }
+                else if((hdata->GetBinContent(ibin))<10.){
+                  res=0.5;
+                  if(ibin==1 || ibin==6 || ibin==11 || ibin==12 || ibin==13 || ibin==18){
+                    iDetPM = kFALSE;
+                  }
+                }
+                //
+                resADC += res;
+                test[specie] += res;
+                count++;
+             }
+             if(nentries != -99) messages.Add(new TObjString(Form("#entries %d",nentries)));
+             else messages.Add(new TObjString("#entries not known"));
+             //
+             Float_t rv=1.;
+             if(hdata->GetNbinsX() != 0) rv = resADC/hdata->GetNbinsX();
+             if(rv == 1.) messages.Add(new TObjString("ADCs are OK!")); 
+             else if(iDetPM==kFALSE){
+               messages.Add(new TObjString("Problem with some ADC!"));
+                messages.Add(new TObjString("IF THIS IS A PHYSICS RUN"));
+             }
+             else if(iDetPM==kTRUE) messages.Add(new TObjString("Minor problem with ADCs"));
+             SetupHisto(messages, *hdata, rv);
+           }
+           else if(irawHisto==19){
+             // Reference values from RUN 145456
+             Double_t refTDCs[6] = {-322.5,-319.1,-320.9,-319.2,-319.7,-319.2};
+             Float_t resTDC=0.;
+             for(int ibin=1; ibin<=hdata->GetNbinsX(); ibin++){
+                if(TMath::Abs((hdata->GetBinContent(ibin))-refTDCs[ibin-1])<3.){
+                  res=1.;
+                }
+                else if(TMath::Abs((hdata->GetBinContent(ibin))-refTDCs[ibin-1])<5.){
+                  res=0.8;
+                }
+                else res=0.5;
+                //
+                resTDC += res;
+                test[specie] += res;
+                count++;
+             }
+             Float_t rv=1.;
+             if(hdata->GetNbinsX() != 0) rv = resTDC/hdata->GetNbinsX();
+             if(rv == 1.) messages.Add(new TObjString("TDCs are OK!")); 
+             else if(rv<1. && rv>0.75) messages.Add(new TObjString("Minor problem with TDCs"));
+             else{
+               messages.Add(new TObjString("Serious problem in ZDC timing"));
+                messages.Add(new TObjString("IF THIS IS A PHYSICS RUN"));
+             }
+             SetupHisto(messages, *hdata, rv);
+           }
+           irawHisto++;
+           
+          } 
+          // -------------------------------------------------------------------
+         else if(index == AliQAv1::kREC) {
+           //
+            // Check REC HIGH GAIN CHAIN histos
+            if(hdata->GetEntries()>0){
+             if(irecHisto==0)       meanZNC = hdata->GetMean();
+             else if(irecHisto==1)  meanZNA = hdata->GetMean();
+             else if(irecHisto==2)  meanZPC = hdata->GetMean();
+             else if(irecHisto==3)  meanZPA = hdata->GetMean();
+             else if(irecHisto==4)  pmQZNC = hdata->GetMean();
+             else if(irecHisto==5)  pmQZNA = hdata->GetMean();
+             else if(irecHisto==6)  pmQZPC = hdata->GetMean();
+             else if(irecHisto==7)  pmQZPA = hdata->GetMean();
+             else if(irecHisto==8)  pmCZNC = hdata->GetMean();
+             else if(irecHisto==9)  pmCZNA = hdata->GetMean();
+             else if(irecHisto==10) pmCZPC = hdata->GetMean();
+             else if(irecHisto==11) pmCZPA = hdata->GetMean();
+           }
+           //
+           // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
+           if(irecHisto==11){
+             if(TMath::Abs(meanZNC)>1.e-10){
+               if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(meanZNA)>1.e-10){
+               if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(meanZPC)>1.e-10){
+               if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(meanZPA)>1.e-10){
+               if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+           }
+           irecHisto++;            
+          } 
+          // -------------------------------------------------------------------
+         else if(index == AliQAv1::kESD) {
+           //
+            // Check ESD HIGH GAIN CHAIN histos
+            if(hdata->GetEntries()>0){
+             if(esdInd==0)      sumADCZNC = hdata->GetMean();
+             else if(esdInd==1) sumADCZNA = hdata->GetMean();
+             else if(esdInd==2) sumADCZPC = hdata->GetMean();
+             else if(esdInd==3) sumADCZPA = hdata->GetMean();
+             else if(esdInd==6)  pmQZNC = hdata->GetMean();
+             else if(esdInd==7)  pmQZNA = hdata->GetMean();
+             else if(esdInd==8)  pmQZPC = hdata->GetMean();
+             else if(esdInd==9) pmQZPA = hdata->GetMean();
+             else if(esdInd==10) pmCZNC = hdata->GetMean();
+             else if(esdInd==11) pmCZNA = hdata->GetMean();
+             else if(esdInd==12) pmCZPC = hdata->GetMean();
+             else if(esdInd==13) pmCZPA = hdata->GetMean();
+           }
+           //
+           // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
+           if(esdInd==15){
+             if(TMath::Abs(sumADCZNC)>1.e-10){
+               if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(sumADCZNA)>1.e-10){
+               if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(sumADCZPC)>1.e-10){
+               if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(sumADCZPA)>1.e-10){
+               if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
             }
             esdInd++;
-            // -------------------------------------------------------------------
-         } else {
+         }  
+        else {
            AliWarning(Form("\n\t No ZDC QA for %s task\n",taskName)); 
-           return NULL;
+           return ;
          }
-        } else
-            AliError("AliZDCQAChecker-> No histos!!!\n");
+        }//if(hdata) 
+       else AliError("AliZDCQAChecker-> No histos!!!\n");
       }
-    }
+    } // LowMult (p-p)
     // ====================================================================
     //         Checks for A-A events
     // ====================================================================
-    //else if (beamType==1){ 
-    else if ( specie == AliRecoParam::kHighMult) {
+    if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kHighMult) {
       if(list[specie]->GetEntries()==0){  
-        AliWarning("\tAliZDCQAChecker->The list to be checked is empty!");
-        return test;
+        AliWarning("\t The list to be checked is empty!");
+        return ;
       }
-      //AliInfo(Form("\n\tAliZDCQAChecker-> checking QA histograms for task %s\n\n",taskName));
-      
+      //AliDebug(AliQAv1::GetQADebugLevel(), Form("\n\tAliZDCQAChecker-> checking QA histograms for task %s\n\n",taskName));
+      //
       TIter next(list[specie]); 
+      count = 0; 
       TH1 * hdata;       
+      
+      Float_t res=0., percentageDiff=0.10;
+      Float_t meanZNA=0., meanZNC=0., meanZPA=0., meanZPC=0.;
+      Float_t pmCZNA=0., pmCZNC=0., pmCZPA=0., pmCZPC=0.;
+      Float_t pmQZNA=0., pmQZNC=0., pmQZPA=0., pmQZPC=0.;
+      Float_t sumADCZNA=0., sumADCZNC=0., sumADCZPA=0., sumADCZPC=0.;
+      Float_t adcCZNA=0., adcCZNC=0., adcCZPA=0., adcCZPC=0.;
+      Float_t adcQZNA=0., adcQZNC=0., adcQZPA=0., adcQZPC=0.;
+      Int_t nentries=-99;
+      
       while((hdata = dynamic_cast<TH1 *>(next()))){
         if(hdata){ 
-          //AliInfo(Form("\tAliZDCQAChecker-> checking histo %s",hdata->GetName()));
-          ntests[specie] = 0; 
-          Double_t meanX=0., meanY=0.;
-          Double_t meanZNA=0., rmsZNA=0., meanZNC=0.;
-          Double_t meanZPA=0., rmsZPA=0., meanZPC=0.;
-          Double_t pmCZNA=0., pmCZNC=0., pmCZPA=0., pmCZPC=0.;
-          Double_t pmQZNA=0., pmQZNC=0., pmQZPA=0., pmQZPC=0.;
-          Float_t  res=0.;
-          Int_t    testgood=0;
+          //AliDebug(AliQAv1::GetQADebugLevel(), Form("\tAliZDCQAChecker-> checking histo %s",hdata->GetName()));
           // -------------------------------------------------------------------
-          if(index == AliQA::kSIM){
-            Int_t    digInd=0;
-            // Check DIGITS histos
-            if (!(strncmp(hdata->GetName(),"hDig",4))){
-              // [1] check response of ZNC vs. ZNA
-              if(digInd==0 || digInd==1){
-                if(digInd==0){
-                  if(hdata->GetEntries() != 0.){
-                    testgood=1;
-                    meanZNA = hdata->GetMean();
-                    rmsZNA = hdata->GetRMS();
-                  }
-                }
-                else{
-                  if(hdata->GetEntries() != 0.){
-                    testgood=1;
-                    meanZNC = hdata->GetMean();
-                  }
-                  else testgood=0;
-                  // check if the response m.v. of ZNA and ZNC are equal (@ 1sigma level)
-                  if(testgood==1){
-                    if(TMath::Abs(meanZNA-meanZNC)<rmsZNA) 
-                      res=1.;
-                    else 
-                      res=.5;
-                    testgood=0;
-                    test[specie] += res;
-                    ntests[specie]++;
-                    //AliInfo(Form("\t %d performed tests, results %1.2f\n",ntests[specie],test[specie]/ntests[specie]));
-                  }
-                  else res=0.;
-                }
-              }
-              // [2] check response of ZPC vs. ZPA
-              else if(digInd==2 || digInd==3){
-                if(digInd==2){
-                  if(hdata->GetEntries() != 0.){
-                    testgood=1;
-                    meanZPA = hdata->GetMean();
-                    rmsZPA = hdata->GetRMS();
-                  }
-                }
-                else{
-                  if(hdata->GetEntries() != 0.){
-                    testgood=1;
-                    meanZPC = hdata->GetMean();
-                  }
-               // check if the response m.v. of ZPA and ZPC are equal (@ 3sigma level)
-               if(testgood==1){
-                 if(TMath::Abs(meanZPA-meanZPC)<(3.*rmsZPA)) 
-                   res=1.;
-                 else 
-                   res=.5;
-                 test[specie] += res;
-                 ntests[specie]++;
-                 testgood=0;
-                 //AliInfo(Form("\t %d performed tests, results %1.2f\n",ntests[specie],test[specie]/ntests[specie]);
-               }
-               else res=0.;
-                }
+          if(index == AliQAv1::kSIM){
+            // Check HITS histos
+            if (!(strncmp(hdata->GetName(),"hHits",5))){
+              if(hdata->GetEntries()>0){
+               if(ihitHisto==0)      meanZNC = hdata->GetMean();
+               else if(ihitHisto==1) meanZNA = hdata->GetMean();
+               else if(ihitHisto==2) meanZPC = hdata->GetMean();
+               else if(ihitHisto==3) meanZPA = hdata->GetMean();
+               else if(ihitHisto==4) pmQZNC = hdata->GetMean();
+               else if(ihitHisto==5) pmQZNA = hdata->GetMean();
+               else if(ihitHisto==6) pmQZPC = hdata->GetMean();
+               else if(ihitHisto==7) pmQZPA = hdata->GetMean();
+               else if(ihitHisto==8)  pmCZNC = hdata->GetMean();
+               else if(ihitHisto==9)  pmCZNA = hdata->GetMean();
+               else if(ihitHisto==10) pmCZPC = hdata->GetMean();
+               else if(ihitHisto==11) pmCZPA = hdata->GetMean();
+             }
+             //
+             // --- Check whether 2*|Mean ZNA - Mean ZNC|/(Mean ZNA + Mean ZNC) < percentageDiff
+             // --- and 2*|Mean ZPA - Mean ZPC|/(Mean ZPA + Mean ZPC) < 2*percentageDiff
+             if(ihitHisto==3){
+               if(TMath::Abs(meanZNC)>1.e-10 && TMath::Abs(meanZNA)>1.e-10){
+                  if((2*TMath::Abs(meanZNC-meanZNA)/(meanZNA+meanZNC))<percentageDiff) 
+                    res=1.;
+                  else 
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
+               if(TMath::Abs(meanZPC)>1.e-10 && TMath::Abs(meanZPA)>1.e-10){
+                  if((TMath::Abs(meanZPC-meanZPA)/(meanZPA+meanZPC))<percentageDiff) 
+                    res=1.;
+                  else 
+                    res=.5;
+                  test[specie] += res;
+                  count++;
+               }
               }
-              // [2] check PMC responses vs. summed PMQ responses
-              else if(digInd>3 && digInd<12){
-                if(digInd==4) pmQZNC = hdata->GetMean();
-                else if(digInd==5) pmQZNA = hdata->GetMean();
-                else if(digInd==6) pmQZPC = hdata->GetMean();
-                else if(digInd==7) pmQZPA = hdata->GetMean();
-                else if(digInd==8){
-                  pmCZNC = hdata->GetMean();
-                  if(TMath::Abs(pmQZNC-pmCZNC)<(0.1*(pmQZNC+pmCZNC)/2)) 
+             // --- Check whether (mean PMQi - PMC)/PMC < percentageDiff
+             if(ihitHisto==11){
+               if(TMath::Abs(meanZNC)>1.e-10){
+                  if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
                     res=1.;
                   else 
-                    res=0.5;
+                    res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                  //AliInfo(Form("\t %d performed tests, results %1.2f\n",ntests[specie],test[specie]/ntests[specie]));
-                }
-                else if(digInd==9){
-                  pmCZNA = hdata->GetMean();
-                  if(TMath::Abs(pmQZNA-pmCZNA)<(0.1*(pmQZNA+pmCZNA)/2)) 
+                  count++;
+               }
+               if(TMath::Abs(meanZNA)>1.e-10){
+                  if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
                     res=1.;
                   else 
-                    res=0.5;
+                    res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                  //AliInfo(Form("\t %d performed tests, results %1.2f\n",ntests[specie],test[specie]/ntests[specie]));
-                }
-                else if(digInd==10){
-                  pmCZPC = hdata->GetMean();
-                  if(TMath::Abs(pmQZPC-pmCZPC)<(0.1*(pmQZPC+pmCZPC)/2)) 
+                  count++;
+               }
+               if(TMath::Abs(meanZPC)>1.e-10){
+                  if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
                     res=1.;
                   else 
-                    res=0.5;
+                    res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                  //AliInfor(Form("\t %d performed tests, results %1.2f\n",ntests[specie],test[specie]/ntests[specie]);
-                }
-                else if(digInd==11){
-                  pmCZPA = hdata->GetMean();
-                  if(TMath::Abs(pmQZPA-pmCZPA)<(0.1*(pmQZPA+pmCZPA)/2)) 
+                  count++;
+               }
+               if(TMath::Abs(meanZPA)>1.e-10){
+                  if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
                     res=1.;
                   else 
-                    res=0.5;
+                    res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                  //AliInfor(Form("\t %d performed tests, results %1.2f\n",ntests[specie],test[specie]/ntests[specie]));
-                }
-              }
-              digInd++;
-            }
-            // Check HITS histos
-            else{ 
-              // hits histos
-              meanX = hdata->GetMean(1);
-              meanY = hdata->GetMean(2);
-              // check if the spot is centered
-              if((TMath::Abs(meanX)<0.2) && (TMath::Abs(meanY)<0.2)) 
-                res=1.;
-              else 
-                res=0.5;
-              test[specie] += res;
-              ntests[specie]++;
-              //AliInfo(Form("\t %d performed tests, results %1.2f\n",ntests[specie],test[specie]/ntests[specie]);
+                  count++;
+               }
+             }
+             ihitHisto++;
             }
-          }
-          // -------------------------------------------------------------------
-          else if(index == AliQA::kRAW){
-            Int_t    rawInd=0;
-            // [1] check response of ZNC vs. ZNA
-            if(rawInd==0 || rawInd==1){
-              if(rawInd==0){
-                if(hdata->GetEntries() != 0.){
-                  testgood=1;
-                  meanZNA = hdata->GetMean();
-                  rmsZNA = hdata->GetRMS();
-                }
-              }
-              else{
-                if(hdata->GetEntries() != 0.){
-                  testgood=1;
-                  meanZNC = hdata->GetMean();
-                }
-                else testgood=0;
-                // check if the response m.v. of ZNA and ZNC are equal (@ 1sigma level)
-                if(testgood==1){
-                  if(TMath::Abs(meanZNA-meanZNC)<rmsZNA) 
+            // Check DIGITS histos
+            else if (!(strncmp(hdata->GetName(),"hDig",4))){
+              if(hdata->GetEntries()>0){
+               if(idigHisto==0)      sumADCZNC = hdata->GetMean();
+               else if(idigHisto==1) sumADCZNA = hdata->GetMean();
+               else if(idigHisto==2) sumADCZPC = hdata->GetMean();
+               else if(idigHisto==3) sumADCZPA = hdata->GetMean();
+               else if(idigHisto==4) adcQZNC = hdata->GetMean();
+               else if(idigHisto==5) adcQZNA = hdata->GetMean();
+               else if(idigHisto==6) adcQZPC = hdata->GetMean();
+               else if(idigHisto==7) adcQZPA = hdata->GetMean();
+               else if(idigHisto==8)  adcCZNC = hdata->GetMean();
+               else if(idigHisto==9)  adcCZNA = hdata->GetMean();
+               else if(idigHisto==10) adcCZPC = hdata->GetMean();
+               else if(idigHisto==11) adcCZPA = hdata->GetMean();
+             }
+             //
+             // --- Check whether 2*|Mean ZNA - Mean ZNC|/(Mean ZNA + Mean ZNC) < percentageDiff
+             // --- and 2*|Mean ZPA - Mean ZPC|/(Mean ZPA + Mean ZPC) < 2*percentageDiff
+             if(idigHisto==3){
+               if(TMath::Abs(sumADCZNC)>1.e-10 && TMath::Abs(sumADCZNA)>1.e-10){
+                  if((2*TMath::Abs(sumADCZNC-sumADCZNA)/(sumADCZNA+sumADCZNC))<percentageDiff) 
                     res=1.;
                   else 
                     res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                  testgood=0;
-                }
-                else 
-                  res=0.;
-              }
-            }
-            // [2] check response of ZPC vs. ZPA
-            else if(rawInd==2 || rawInd==3){
-              if(rawInd==2){
-                if(hdata->GetEntries() != 0.){
-                  testgood=1;
-                  meanZPA = hdata->GetMean();
-                  rmsZPA = hdata->GetRMS();
-                }
-              }
-              else{
-                if(hdata->GetEntries() != 0.){
-                  testgood=1;
-                  meanZPC = hdata->GetMean();
-                }
-                // check if the response m.v. of ZPA and ZPC are equal (@ 3sigma level)
-                if(testgood==1){
-                  if(TMath::Abs(meanZPA-meanZPC)<(3.*rmsZPA)) 
+                  count++;
+               }
+               if(TMath::Abs(sumADCZPC)>1.e-10 && TMath::Abs(sumADCZPA)>1.e-10){
+                  if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) 
                     res=1.;
                   else 
                     res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                  testgood=0;
-                }
-                else 
-                  res=0.;
-              }
-            }
-            // [2] check PMC responses vs. summed PMQ responses
-            else if(rawInd>3 && rawInd<12){
-              if(rawInd==4) pmQZNC = hdata->GetMean();
-              else if(rawInd==5) pmQZNA = hdata->GetMean();
-              else if(rawInd==6) pmQZPC = hdata->GetMean();
-              else if(rawInd==7) pmQZPA = hdata->GetMean();
-              else if(rawInd==8){
-                pmCZNC = hdata->GetMean();
-                if(TMath::Abs(pmQZNC-pmCZNC)<(0.1*(pmQZNC+pmCZNC)/2)) 
-                  res=1.;
-                else 
-                  res=0.5;
-                test[specie] += res;
-                ntests[specie]++;
-              }
-              else if(rawInd==9){
-                pmCZNA = hdata->GetMean();
-                if(TMath::Abs(pmQZNA-pmCZNA)<(0.1*(pmQZNA+pmCZNA)/2)) 
-                  res=1.;
-                else 
-                  res=0.5;
-                test[specie] += res;
-                ntests[specie]++;
+                  count++;
+               }
               }
-              else if(rawInd==10){
-                pmCZPC = hdata->GetMean();
-                if(TMath::Abs(pmQZPC-pmCZPC)<(0.1*(pmQZPC+pmCZPC)/2)) 
-                  res=1.;
-                else 
-                  res=0.5;
-                test[specie] += res;
-                ntests[specie]++;
-              }
-              else if(rawInd==11){
-                pmCZPA = hdata->GetMean();
-                if(TMath::Abs(pmQZPA-pmCZPA)<(0.1*(pmQZPA+pmCZPA)/2)) 
-                  res=1.;
-                else 
-                  res=0.5;
-                test[specie] += res;
-                ntests[specie]++;
-              }
-            }
-            rawInd++;
-          }
-          // -------------------------------------------------------------------
-          else if(index == AliQA::kESD){
-            Double_t eneQZNC, eneQZNA ;  
-            Double_t eneQZPC, eneQZPA ;  
-            Double_t eneCZNC, eneCZNA ;  
-            Double_t eneCZPC, eneCZPA ;            
-            Int_t esdInd=0;
-            if(esdInd<2){
-              // hits histos
-              meanX = hdata->GetMean(1);
-              meanY = hdata->GetMean(2);
-              // check if the spot is centered
-              if((TMath::Abs(meanX)<0.2) && (TMath::Abs(meanY)<0.2)) 
-                res=1.;
-              else 
-                res=0.5;
-              test[specie] += res;
-              ntests[specie]++;
-            }
-            //
-            else{
-              // [1] check response of ZNC vs. ZNA
-              if(esdInd==0 || esdInd==1){
-                if(esdInd==0){
-                  if(hdata->GetEntries() != 0.){
-                    testgood=1;
-                    meanZNA = hdata->GetMean();
-                    rmsZNA = hdata->GetRMS();
-                  }
-                }
-                else{
-                  if(hdata->GetEntries() != 0.){
-                    testgood=1;
-                    meanZNC = hdata->GetMean();
-                  }
-                  else testgood=0;
-                  // check if the response m.v. of ZNA and ZNC are equal (@ 1sigma level)
-                  if(testgood==1){
-                    if(TMath::Abs(meanZNA-meanZNC)<rmsZNA) 
-                      res=1.;
-                    else 
-                      res=.5;
-                    testgood=0;
-                    test[specie] += res;
-                    ntests[specie]++;
-                  }
-                  else res=0.;
-                }
-              }
-              // [2] check response of ZPC vs. ZPA
-              else if(esdInd==2 || esdInd==3){
-                if(esdInd==2){
-                  if(hdata->GetEntries() != 0.){
-                    testgood=1;
-                    meanZPA = hdata->GetMean();
-                    rmsZPA = hdata->GetRMS();
-                  }
-                }
-                else{
-                  if(hdata->GetEntries() != 0.){
-                    testgood=1;
-                    meanZPC = hdata->GetMean();
-                  }
-                  // check if the response m.v. of ZPA and ZPC are equal (@ 3sigma level)
-                  if(testgood==1){
-                    if(TMath::Abs(meanZPA-meanZPC)<(3.*rmsZPA)) 
-                      res=1.;
-                    else 
-                      res=.5;
-                    test[specie] += res;
-                    ntests[specie]++;
-                    testgood=0;
-                  }
-                  else res=0.;
-                }
-              }
-              // [2] check eneC responses vs. summed eneQ responses
-              else if(esdInd>3 && esdInd<12){
-                if(esdInd==4) eneQZNC = hdata->GetMean();
-                else if(esdInd==5) eneQZNA = hdata->GetMean();
-                else if(esdInd==6) eneQZPC = hdata->GetMean();
-                else if(esdInd==7) eneQZPA = hdata->GetMean();
-                else if(esdInd==8){
-                  eneCZNC = hdata->GetMean();
-                  if(TMath::Abs(eneQZNC-eneCZNC)<(0.1*(eneQZNC+eneCZNC)/2)) 
+             // --- Check whether (sumADC PMQi - PMC)/PMC < percentageDiff
+             if(idigHisto==11){
+               if(TMath::Abs(sumADCZNC)>1.e-10){
+                  if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
                     res=1.;
                   else 
-                    res=0.5;
+                    res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                }
-                else if(esdInd==9){
-                  eneCZNA = hdata->GetMean();
-                  if(TMath::Abs(eneQZNA-eneCZNA)<(0.1*(eneQZNA+eneCZNA)/2)) 
+                  count++;
+               }
+               if(TMath::Abs(sumADCZNA)>1.e-10){
+                  if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
                     res=1.;
                   else 
-                    res=0.5;
+                    res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                }
-                else if(esdInd==10){
-                  eneCZPC = hdata->GetMean();
-                  if(TMath::Abs(eneQZPC-eneCZPC)<(0.1*(eneQZPC+eneCZPC)/2)) 
+                  count++;
+               }
+               if(TMath::Abs(sumADCZPC)>1.e-10){
+                  if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
                     res=1.;
                   else 
-                    res=0.5;
+                    res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                }
-                else if(esdInd==11){
-                  eneCZPA = hdata->GetMean();
-                  if(TMath::Abs(eneQZPA-eneCZPA)<(0.1*(eneQZPA+eneCZPA)/2)) 
+                  count++;
+               }
+               if(TMath::Abs(sumADCZPA)>1.e-10){
+                  if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
                     res=1.;
                   else 
-                    res=0.5;
+                    res=.5;
                   test[specie] += res;
-                  ntests[specie]++;
-                }
-              }
-              esdInd++;
+                  count++;
+               }
+             }
+              idigHisto++;
             }
-            //AliInfor(Form("\t %d performed tests, results %1.2f\n",ntests[specie],test[specie]/ntests[specie]));
-          }
-          else {
-          AliWarning(Form("\n\t No ZDC QA for %s task\n",taskName)); 
-          return NULL;
           }
-        } else
-            AliError("\t AliZDCQAChecker->No histos!!!\n");
-        }
-    } else {
-      AliError(Form("Checking not implemented for %s, %s", 
-                    AliRecoParam::GetEventSpecieName(AliRecoParam::kCosmic), 
-                    AliRecoParam::GetEventSpecieName(AliRecoParam::kCalib))) ; 
-    }
+          // -------------------------------------------------------------------
+          else if(index == AliQAv1::kRAW){
+           //
+            // Check RAW HIGH GAIN CHAIN histos
+            
+           messages.Clear();
+           
+           Bool_t iDetPM = kTRUE;
+           // --- Checks
+           if(irawHisto==16) nentries = Int_t (hdata->GetEntries());
+           if(irawHisto==18){ 
+             Float_t resADC=0.;
+             for(int ibin=1; ibin<=hdata->GetNbinsX(); ibin++){
+                if((hdata->GetBinContent(ibin))>10.){
+                  res=1.;
+                }
+                else if((hdata->GetBinContent(ibin))<10. && 
+                  (ibin==1 || ibin==6 || ibin==11 || ibin==12 || ibin==13 || ibin==18)){
+                  res=0.5;
+                  iDetPM = kFALSE;
+                }
+                //
+                resADC += res;
+                test[specie] += res;
+                count++;
+             }
+             if(nentries != -99) messages.Add(new TObjString(Form("#entries %d",nentries)));
+             else messages.Add(new TObjString("#entries not known"));
+             //
+             Float_t rv=1.;
+             if(hdata->GetNbinsX() != 0) rv = resADC/hdata->GetNbinsX();
+             if(rv > 0.98) messages.Add(new TObjString("ADCs are OK!")); 
+             else if(iDetPM==kFALSE){
+               messages.Add(new TObjString("Problem with ADCs!"));
+                messages.Add(new TObjString("IF THIS IS A PHYSICS RUN"));
+             }
+             else if(iDetPM==kTRUE) messages.Add(new TObjString("Minor problem with ADCs"));
+             SetupHisto(messages, *hdata, rv);
+           }
+           else if(irawHisto==19){
+             // Reference values from RUN 137161
+             Double_t refTDCs[6] = {-320.7,-319.0,-318.6,-319.9,-321.3,-320.8};
+             Float_t resTDC=0.;
+             for(int ibin=1; ibin<=hdata->GetNbinsX(); ibin++){
+                if(TMath::Abs((hdata->GetBinContent(ibin))-refTDCs[ibin-1])<3.){
+                  res=1.;
+                }
+                else if(TMath::Abs((hdata->GetBinContent(ibin))-refTDCs[ibin-1])<4.){
+                  res=0.8;
+                }
+                else res=0.5;
+                //
+                resTDC += res;
+                test[specie] += res;
+                count++;
+             }
+             Float_t rv=1.;
+             if(hdata->GetNbinsX() != 0) rv = resTDC/hdata->GetNbinsX();
+             if(rv == 1.) messages.Add(new TObjString("TDCs are OK!")); 
+             else if(rv<1 && rv>0.75) messages.Add(new TObjString("Minor problem with TDCs"));
+             else{
+               messages.Add(new TObjString("Serious problem in ZDC timing"));
+                messages.Add(new TObjString("IF THIS IS A PHYSICS RUN"));
+             }
+             SetupHisto(messages, *hdata, rv);
+           }
+           irawHisto++;
+           
+         }   
+          // -------------------------------------------------------------------
+          else if(index == AliQAv1::kREC){
+           //
+            // Check RAW HIGH GAIN CHAIN histos
+            if(hdata->GetEntries()>0){
+             if(irecHisto==0)       meanZNC = hdata->GetMean();
+             else if(irecHisto==1)  meanZNA = hdata->GetMean();
+             else if(irecHisto==2)  meanZPC = hdata->GetMean();
+             else if(irecHisto==3)  meanZPA = hdata->GetMean();
+             else if(irecHisto==4)  pmQZNC = hdata->GetMean();
+             else if(irecHisto==5)  pmQZNA = hdata->GetMean();
+             else if(irecHisto==6)  pmQZPC = hdata->GetMean();
+             else if(irecHisto==7)  pmQZPA = hdata->GetMean();
+             else if(irecHisto==8)  pmCZNC = hdata->GetMean();
+             else if(irecHisto==9)  pmCZNA = hdata->GetMean();
+             else if(irecHisto==10) pmCZPC = hdata->GetMean();
+             else if(irecHisto==11) pmCZPA = hdata->GetMean();
+           }
+            //
+           // --- Check whether 2*|Mean ZNA - Mean ZNC|/(Mean ZNA + Mean ZNC) < percentageDiff
+           // --- and 2*|Mean ZPA - Mean ZPC|/(Mean ZPA + Mean ZPC) < 2*percentageDiff
+           if(irecHisto==3){
+             if(TMath::Abs(meanZNC)>1.e-10 && TMath::Abs(meanZNA)>1.e-10){
+               if((2*TMath::Abs(meanZNC-meanZNA)/(meanZNA+meanZNC))<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(meanZPC)>1.e-10 && TMath::Abs(meanZPA)>1.e-10){
+               if((TMath::Abs(meanZPC-meanZPA)/(meanZPA+meanZPC))<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+            }
+           // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
+           if(irecHisto==11){
+             if((TMath::Abs(meanZNC)>1.e-10) && (pmCZNC>1.e-10)){
+               if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if((TMath::Abs(meanZNA)>1.e-10) && (pmCZNA>1.e-10)){
+               if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if((TMath::Abs(meanZPC)>1.e-10) && (pmCZPC>1.e-10)){
+               if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if((TMath::Abs(meanZPA)>1.e-10) && (pmCZPA>1.e-10)){
+               if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+           }
+           irecHisto++;         
+         }   
+          // -------------------------------------------------------------------
+          else if(index == AliQAv1::kESD){
+           //
+            // Check ESD HIGH GAIN CHAIN histos
+            if(hdata->GetEntries()>0){
+             if(esdInd==0)      sumADCZNC = hdata->GetMean();
+             else if(esdInd==1) sumADCZNA = hdata->GetMean();
+             else if(esdInd==2) sumADCZPC = hdata->GetMean();
+             else if(esdInd==3) sumADCZPA = hdata->GetMean();
+             else if(esdInd==6)  pmQZNC = hdata->GetMean();
+             else if(esdInd==7)  pmQZNA = hdata->GetMean();
+             else if(esdInd==8) pmQZPC = hdata->GetMean();
+             else if(esdInd==9) pmQZPA = hdata->GetMean();
+             else if(esdInd==10) pmCZNC = hdata->GetMean();
+             else if(esdInd==11) pmCZNA = hdata->GetMean();
+             else if(esdInd==12) pmCZPC = hdata->GetMean();
+             else if(esdInd==13) pmCZPA = hdata->GetMean();
+           }
+           //
+           // --- Check whether 2*|Mean ZNA - Mean ZNC|/(Mean ZNA + Mean ZNC) < percentageDiff
+           // --- and 2*|Mean ZPA - Mean ZPC|/(Mean ZPA + Mean ZPC) < 2*percentageDiff
+           if(esdInd==5){
+             if(TMath::Abs(sumADCZNC)>1.e-10 && TMath::Abs(sumADCZNA)>1.e-10){
+               if((2*TMath::Abs(sumADCZNC-sumADCZNA)/(sumADCZNA+sumADCZNC))<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(sumADCZPC)>1.e-10 && TMath::Abs(sumADCZPA)>1.e-10){
+               if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+            }
+           // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
+           if(esdInd==15){
+             if(TMath::Abs(sumADCZNC)>1.e-10){
+               if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(sumADCZNA)>1.e-10){
+               if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+/*           if(TMath::Abs(sumADCZPC)>1.e-10){
+               if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+             if(TMath::Abs(sumADCZPA)>1.e-10){
+               if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
+                 res=1.;
+               else 
+                 res=.5;
+               test[specie] += res;
+               count++;
+             }
+*/
+            }
+            esdInd++;
+          }  
+         else {
+            AliWarning(Form("\n\t No ZDC QA for %s task\n",taskName)); 
+            return ;
+          } 
+        }//if(hdata) 
+       else AliError("\t  No histos found for ZDC!!!\n");
+      }
+    } // HighMult (Pb-Pb) 
+    // ====================================================================
+    //         Checks for Calibration events
+    // ====================================================================
+    if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCalib) {
+      AliWarning(Form("\n\t No check implemented in ZDC QA for %s task in CALIBRATION events\n",taskName)); 
+      return ;
+    } // Calibration
+    // ====================================================================
+    //         Checks for cosmic events
+    // ====================================================================
+    else if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCosmic) {
+      AliWarning(Form("\n\t No check needed in ZDC QA for %s task in COSMIC events\n",taskName)); 
+      return ; 
+    } // Cosmic
+    if(TMath::Abs(count)>1.e-10) test[specie] = test[specie]/count;
+    AliDebug(AliQAv1::GetQADebugLevel(), Form("\n\t ZDC QA check result = %1.2f\n",test[specie]));
+  } // Loop on species
+}  
+
+//___________________________________________________________________ 
+void AliZDCQAChecker::SetupHisto(const TObjArray& messages, TH1& histo, Float_t& code)
+{
+  //
+  /// Add text to histos
+  //
+
+  TPaveText* text = new TPaveText(0.70,0.70,0.99,0.99,"NDC");
+    
+  TIter next(&messages);
+  TObjString* str;
+    
+  while ( ( str = static_cast<TObjString*>(next()) ) ){
+    text->AddText(str->String());
   }
-   
-  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
-    if(ntests[specie]!=0) test[specie] = test[specie]/ntests[specie];
-        printf("\n\tAliZDCQAChecker-> QA check result = %1.2f\n",test[specie]);
+
+  TString defaultText = "";
+
+  Int_t color = 0;
+  if(code==1.){
+    color = kGreen;
+    defaultText = "Everything is fine!";
+  }  
+  else if(code<1. && code>=0.85){  
+    color = kYellow;
+    defaultText = "To be monitored in next runs";
   }
-  delete [] ntests ; 
-  return test; 
-}  
+  else if(code<0.85 && code>=0.6){
+    color = kOrange;
+    defaultText = "notify the expert DURING THE DAY!";
+  }
+  else if(code<0.6){
+    color = kRed;
+    defaultText = "PLEASE CALL THE EXPERT!!!!";
+  }
+
+
+  text->AddText(defaultText.Data());
+  text->SetFillColor(color);
+                      
+  //histo.SetFillStyle(1001);
+  //histo.SetFillColor(color);
 
+  histo.SetStats(kFALSE);
+    
+  histo.GetListOfFunctions()->Add(text);
+}