/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/


// --- ROOT system ---
#include <TH1F.h> 
#include <TIterator.h> 
#include <TString.h> 

// --- Standard library ---

// --- AliRoot header files ---
#include "AliLog.h"
#include "AliZDCQAChecker.h"

ClassImp(AliZDCQAChecker)

//____________________________________________________________________________
Double_t * AliZDCQAChecker::Check(AliQAv1::ALITASK_t index, TObjArray ** list) 
{
  // 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 = AliQAv1::GetAliTaskName(index);
  //
  for(Int_t specie = 0; specie < AliRecoParam::kNSpecies; specie++) {
    test[specie] = 1.0; 
    ntests[specie] = 0; 
    if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) )  continue ; 
    // ====================================================================
    // 	Checks for p-p events
    // ====================================================================
    if (specie == AliRecoParam::kLowMult) {
      if(list[specie]->GetEntries()==0){  
        AliWarning("\t The list to be checked is empty!"); // nothing to check
        return test;
      }
      //AliDebug(AliQAv1::GetQADebugLevel(), Form("\n\tAliZDCQAChecker-> checking QA histograms for task %s\n\n",taskName));
      TIter next(list[specie]); 
      ntests[specie] = 0; 
      TH1 * hdata;	  
      Float_t res=0., percentageDiff=0.15;
      while((hdata = dynamic_cast<TH1 *>(next()))){
        if(hdata){ 
          // -------------------------------------------------------------------
          if(index == AliQAv1::kSIM){
            //AliDebug(AliQAv1::GetQADebugLevel(), Form("\tAliZDCQAChecker-> checking histo %s",hdata->GetName()));
            // Check HITS histos
            Float_t sumZNA=0., sumZNC=0., sumZPA=0., sumZPC=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 ihitHisto=0, idigHisto=0;
            //
	    if(!(strncmp(hdata->GetName(),"hHits",5))){
              if(hdata->GetEntries()>0){
	        if(ihitHisto==0)      sumZNC = hdata->GetMean();
		else if(ihitHisto==1) sumZNA = hdata->GetMean();
		else if(ihitHisto==2) sumZPC = hdata->GetMean();
		else if(ihitHisto==3) sumZPA = 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(sumZNC!=0){
                  if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumZNA!=0){
                  if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
                    res=1.;
                  else percentageDiff=
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumZPC!=0){
                  if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumZPA!=0){
                  if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	      }
	      ihitHisto++;
            }
            // 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(sumADCZNC!=0){
                  if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumADCZNA!=0){
                  if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumADCZPC!=0){
                  if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumADCZPA!=0){
                  if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	      }
	      idigHisto++;	      
            }
          } 
          // -------------------------------------------------------------------
	  else if(index == AliQAv1::kRAW) {
            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 irawHisto=0;
	    //
            // Check RAW HIGH GAIN CHAIN histos
            if(hdata->GetEntries()>0){
	      if(irawHisto==0)      sumADCZNC = hdata->GetMean();
	      else if(irawHisto==1) sumADCZNA = hdata->GetMean();
	      else if(irawHisto==2) sumADCZPC = hdata->GetMean();
	      else if(irawHisto==3) sumADCZPA = hdata->GetMean();
	      else if(irawHisto==4) adcQZNC = hdata->GetMean();
	      else if(irawHisto==5) adcQZNA = hdata->GetMean();
	      else if(irawHisto==6) adcQZPC = hdata->GetMean();
	      else if(irawHisto==7) adcQZPA = hdata->GetMean();
	      else if(irawHisto==8)  adcCZNC = hdata->GetMean();
	      else if(irawHisto==9)  adcCZNA = hdata->GetMean();
	      else if(irawHisto==10) adcCZPC = hdata->GetMean();
	      else if(irawHisto==11) adcCZPA = hdata->GetMean();
	    }
	    //
	    // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
	    if(irawHisto==11){
	      if(sumADCZNC!=0){
            	if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZNA!=0){
            	if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPC!=0){
            	if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPA!=0){
            	if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	    }
	    irawHisto++;	    
          } 
          // -------------------------------------------------------------------
	  else if(index == AliQAv1::kESD) {
            Float_t sumADCZNA=0., sumADCZNC=0., sumADCZPA=0., sumADCZPC=0.;
            Float_t pmCZNA=0., pmCZNC=0., pmCZPA=0., pmCZPC=0.;
            Float_t pmQZNA=0., pmQZNC=0., pmQZPA=0., pmQZPC=0.;
            Int_t esdInd=0;
	    //
            // Check ESD HIGH GAIN CHAIN histos
            if(hdata->GetEntries()!=0){
	      if(esdInd==2)      sumADCZNC = hdata->GetMean();
	      else if(esdInd==3) sumADCZNA = hdata->GetMean();
	      else if(esdInd==4) sumADCZPC = hdata->GetMean();
	      else if(esdInd==5) 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==13){
	      if(sumADCZNC!=0){
            	if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZNA!=0){
            	if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPC!=0){
            	if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPA!=0){
            	if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
            }
            esdInd++;
         }  
	 else {
           AliWarning(Form("\n\t No ZDC QA for %s task\n",taskName)); 
           return NULL;
         }
        }//if(hdata) 
	else AliError("AliZDCQAChecker-> No histos!!!\n");
      }
    } // LowMult (p-p)
    // ====================================================================
    // 	Checks for A-A events
    // ====================================================================
    else if (specie == AliRecoParam::kHighMult) {
      if(list[specie]->GetEntries()==0){  
        AliWarning("\t The list to be checked is empty!");
        return test;
      }
      //AliDebug(AliQAv1::GetQADebugLevel(), Form("\n\tAliZDCQAChecker-> checking QA histograms for task %s\n\n",taskName));
      //
      TIter next(list[specie]); 
      ntests[specie] = 0; 
      TH1 * hdata;	  
      Float_t res=0., percentageDiff=0.10;
      while((hdata = dynamic_cast<TH1 *>(next()))){
        if(hdata){ 
          //AliDebug(AliQAv1::GetQADebugLevel(), Form("\tAliZDCQAChecker-> checking histo %s",hdata->GetName()));
          // -------------------------------------------------------------------
          if(index == AliQAv1::kSIM){
            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 hitInd=0, digInd=0;
            // Check HITS histos
            if (!(strncmp(hdata->GetName(),"hHits",5))){
              if(hdata->GetEntries()>0){
	        if(hitInd==0)      meanZNC = hdata->GetMean();
		else if(hitInd==1) meanZNA = hdata->GetMean();
		else if(hitInd==2) meanZPC = hdata->GetMean();
		else if(hitInd==3) meanZPA = hdata->GetMean();
	        else if(hitInd==4) pmQZNC = hdata->GetMean();
	        else if(hitInd==5) pmQZNA = hdata->GetMean();
	        else if(hitInd==6) pmQZPC = hdata->GetMean();
	        else if(hitInd==7) pmQZPA = hdata->GetMean();
	        else if(hitInd==8)  pmCZNC = hdata->GetMean();
	        else if(hitInd==9)  pmCZNA = hdata->GetMean();
	        else if(hitInd==10) pmCZPC = hdata->GetMean();
	        else if(hitInd==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(hitInd==3){
	        if(meanZNC!=0 && meanZNA!=0){
                  if((2*TMath::Abs(meanZNC-meanZNA)/(meanZNA+meanZNC))<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(meanZPC!=0 && meanZPA!=0){
                  if((TMath::Abs(meanZPC-meanZPA)/(meanZPA+meanZPC))<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
              }
	      // --- Check whether (mean PMQi - PMC)/PMC < percentageDiff
	      if(hitInd==11){
	        if(meanZNC!=0){
                  if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(meanZNA!=0){
                  if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(meanZPC!=0){
                  if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(meanZPA!=0){
                  if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	      }
	      hitInd++;
            }
            // Check DIGITS histos
            else if (!(strncmp(hdata->GetName(),"hDig",4))){
              if(hdata->GetEntries()>0){
	        if(digInd==0)      sumADCZNC = hdata->GetMean();
		else if(digInd==1) sumADCZNA = hdata->GetMean();
		else if(digInd==2) sumADCZPC = hdata->GetMean();
		else if(digInd==3) sumADCZPA = hdata->GetMean();
	        else if(digInd==4) adcQZNC = hdata->GetMean();
	        else if(digInd==5) adcQZNA = hdata->GetMean();
	        else if(digInd==6) adcQZPC = hdata->GetMean();
	        else if(digInd==7) adcQZPA = hdata->GetMean();
	        else if(digInd==8)  adcCZNC = hdata->GetMean();
	        else if(digInd==9)  adcCZNA = hdata->GetMean();
	        else if(digInd==10) adcCZPC = hdata->GetMean();
	        else if(digInd==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(digInd==3){
	        if(sumADCZNC!=0 && sumADCZNA!=0){
                  if((2*TMath::Abs(sumADCZNC-sumADCZNA)/(sumADCZNA+sumADCZNC))<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumADCZPC!=0 && sumADCZPA!=0){
                  if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
              }
	      // --- Check whether (sumADC PMQi - PMC)/PMC < percentageDiff
	      if(digInd==11){
	        if(sumADCZNC!=0){
                  if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumADCZNA!=0){
                  if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumADCZPC!=0){
                  if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	        if(sumADCZPA!=0){
                  if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
                    res=1.;
                  else 
                    res=.5;
                  test[specie] += res;
                  ntests[specie]++;
		}
	      }
              digInd++;
            }
          }
          // -------------------------------------------------------------------
          else if(index == AliQAv1::kRAW){
            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 rawInd=0;
	    //
            // Check RAW HIGH GAIN CHAIN histos
            if(hdata->GetEntries()>0){
	      if(rawInd==0)      sumADCZNC = hdata->GetMean();
	      else if(rawInd==1) sumADCZNA = hdata->GetMean();
	      else if(rawInd==2) sumADCZPC = hdata->GetMean();
	      else if(rawInd==3) sumADCZPA = hdata->GetMean();
	      else if(rawInd==4) adcQZNC = hdata->GetMean();
	      else if(rawInd==5) adcQZNA = hdata->GetMean();
	      else if(rawInd==6) adcQZPC = hdata->GetMean();
	      else if(rawInd==7) adcQZPA = hdata->GetMean();
	      else if(rawInd==8)  adcCZNC = hdata->GetMean();
	      else if(rawInd==9)  adcCZNA = hdata->GetMean();
	      else if(rawInd==10) adcCZPC = hdata->GetMean();
	      else if(rawInd==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(rawInd==3){
	      if(sumADCZNC!=0 && sumADCZNA!=0){
            	if((2*TMath::Abs(sumADCZNC-sumADCZNA)/(sumADCZNA+sumADCZNC))<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPC!=0 && sumADCZPA!=0){
            	if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
            }
	    // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
	    if(rawInd==11){
	      if(sumADCZNC!=0){
            	if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZNA!=0){
            	if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPC!=0){
            	if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPA!=0){
            	if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	    }
	    rawInd++;	 
	  }   
          // -------------------------------------------------------------------
          else if(index == AliQAv1::kESD){
            Float_t sumADCZNA=0., sumADCZNC=0., sumADCZPA=0., sumADCZPC=0.;
            Float_t pmCZNA=0., pmCZNC=0., pmCZPA=0., pmCZPC=0.;
            Float_t pmQZNA=0., pmQZNC=0., pmQZPA=0., pmQZPC=0.;
            Int_t esdInd=0;
	    //
            // Check ESD HIGH GAIN CHAIN histos
            if(hdata->GetEntries()!=0){
	      if(esdInd==2)      sumADCZNC = hdata->GetMean();
	      else if(esdInd==3) sumADCZNA = hdata->GetMean();
	      else if(esdInd==4) sumADCZPC = hdata->GetMean();
	      else if(esdInd==5) 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(sumADCZNC!=0 && sumADCZNA!=0){
            	if((2*TMath::Abs(sumADCZNC-sumADCZNA)/(sumADCZNA+sumADCZNC))<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPC!=0 && sumADCZPA!=0){
            	if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
            }
	    // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
	    if(esdInd==13){
	      if(sumADCZNC!=0){
            	if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZNA!=0){
            	if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPC!=0){
            	if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
	      if(sumADCZPA!=0){
            	if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
            	  res=1.;
            	else 
            	  res=.5;
            	test[specie] += res;
            	ntests[specie]++;
	      }
            }
            esdInd++;
          }  
	  else {
            AliWarning(Form("\n\t No ZDC QA for %s task\n",taskName)); 
            return NULL;
          } 
        }//if(hdata) 
	else AliError("\t  No histos found for ZDC!!!\n");
      }
    } // HighMult (Pb-Pb) 
    else {
       AliError(Form("Checking not implemented for %s, %s", 
                    AliRecoParam::GetEventSpecieName(AliRecoParam::kCosmic), 
                    AliRecoParam::GetEventSpecieName(AliRecoParam::kCalib))) ; 
    }
  } // Loop on species
   
  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
    if(ntests[specie]!=0) test[specie] = test[specie]/ntests[specie];
    AliDebug(AliQAv1::GetQADebugLevel(), Form("\n\t ZDC QA check result = %1.2f\n",test[specie]));
  }
  delete [] ntests ; 
  return test; 
}