// --- ROOT system ---
-#include <TClass.h>
#include <TH1F.h>
-#include <TH1I.h>
#include <TIterator.h>
-#include <TKey.h>
-#include <TFile.h>
+#include <TString.h>
// --- Standard library ---
// --- AliRoot header files ---
#include "AliLog.h"
-#include "AliQA.h"
-#include "AliQAChecker.h"
#include "AliZDCQAChecker.h"
ClassImp(AliZDCQAChecker)
-//__________________________________________________________________
-AliZDCQAChecker& AliZDCQAChecker::operator = (const AliZDCQAChecker& qac )
-{
- // Equal operator.
- this->~AliZDCQAChecker();
- new(this) AliZDCQAChecker(qac);
- return *this;
-}
-
-
//____________________________________________________________________________
-const 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:
+ //
+ const char* taskName = AliQAv1::GetAliTaskName(index);
+ //printf("\n\tAliZDCQAChecker -> checking QA histos for task %s\n",taskName);
+ //
+ Int_t ihitHisto=0, idigHisto=0;
+ Int_t irecHisto=0, irawHisto=0, esdInd=0;
- // Super-basic check on the QA histograms on the input list:
- // look whether they are empty!
-
- Double_t test = 0.0 ;
- Int_t count = 0 ;
-
- if(list->GetEntries() == 0){
- test = 1. ; // nothing to check
- }
- else {
- TIter next(list) ;
- TH1 * hdata ;
- count = 0 ;
- while((hdata = dynamic_cast<TH1 *>(next()))){
- if(hdata){
- Double_t rv = 0.;
- if(hdata->GetEntries()>0)rv=1;
- AliInfo(Form("%s -> %f", hdata->GetName(), rv)) ;
- count++ ;
- test += rv ;
+ 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(AliRecoParam::ConvertIndex(specie) == AliRecoParam::kLowMult){
+ if(list[specie]->GetEntries()==0){
+ AliWarning("\t The list to be checked is empty!"); // nothing to check
+ return;
}
- else{
- AliError("Data type cannot be processed") ;
+ //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.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.;
+
+ 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
+ //
+ 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 (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++;
+ }
+ // 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 == AliQAv1::kRAW) {
+ //
+ // 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==6) adcQZNC = hdata->GetMean();
+ else if(irawHisto==7) adcQZNA = hdata->GetMean();
+ else if(irawHisto==8) adcQZPC = hdata->GetMean();
+ else if(irawHisto==9) adcQZPA = hdata->GetMean();
+ else if(irawHisto==10) adcCZNC = hdata->GetMean();
+ else if(irawHisto==11) adcCZNA = hdata->GetMean();
+ else if(irawHisto==12) adcCZPC = hdata->GetMean();
+ else if(irawHisto==13) adcCZPA = hdata->GetMean();
+ }
+ //
+ // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
+ if(irawHisto==13){
+ 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++;
+ }
+ }
+ 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==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==8) pmQZNC = hdata->GetMean();
+ else if(esdInd==9) pmQZNA = hdata->GetMean();
+ else if(esdInd==10) pmQZPC = hdata->GetMean();
+ else if(esdInd==11) pmQZPA = hdata->GetMean();
+ else if(esdInd==12) pmCZNC = hdata->GetMean();
+ else if(esdInd==13) pmCZNA = hdata->GetMean();
+ else if(esdInd==14) pmCZPC = hdata->GetMean();
+ else if(esdInd==15) 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 {
+ AliWarning(Form("\n\t No ZDC QA for %s task\n",taskName));
+ return ;
+ }
+ }//if(hdata)
+ else AliError("AliZDCQAChecker-> No histos!!!\n");
}
- }
- if (count != 0) {
- if (test==0) {
- AliWarning("Histograms are there, but they are all empty: setting flag to kWARNING");
- test = 0.5; //upper limit value to set kWARNING flag for a task
+ } // LowMult (p-p)
+ // ====================================================================
+ // Checks for A-A events
+ // ====================================================================
+ else if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kHighMult) {
+ if(list[specie]->GetEntries()==0){
+ AliWarning("\t The list to be checked is empty!");
+ return ;
}
- else {
- test /= count ;
+ //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.;
+
+ while((hdata = dynamic_cast<TH1 *>(next()))){
+ if(hdata){
+ //AliDebug(AliQAv1::GetQADebugLevel(), Form("\tAliZDCQAChecker-> checking histo %s",hdata->GetName()));
+ // -------------------------------------------------------------------
+ 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++;
+ }
+ }
+ // --- 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=.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++;
+ }
+ }
+ ihitHisto++;
+ }
+ // 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;
+ 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 (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=.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 == AliQAv1::kRAW){
+ //
+ // 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==6) adcQZNC = hdata->GetMean();
+ else if(irawHisto==7) adcQZNA = hdata->GetMean();
+ else if(irawHisto==8) adcQZPC = hdata->GetMean();
+ else if(irawHisto==9) adcQZPA = hdata->GetMean();
+ else if(irawHisto==10) adcCZNC = hdata->GetMean();
+ else if(irawHisto==11) adcCZNA = hdata->GetMean();
+ else if(irawHisto==12) adcCZPC = hdata->GetMean();
+ else if(irawHisto==13) 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(irawHisto==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;
+ 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(irawHisto==13){
+ 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++;
+ }
+ }
+ 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){
+ 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==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==8) pmQZNC = hdata->GetMean();
+ else if(esdInd==9) pmQZNA = hdata->GetMean();
+ else if(esdInd==10) pmQZPC = hdata->GetMean();
+ else if(esdInd==11) pmQZPA = hdata->GetMean();
+ else if(esdInd==12) pmCZNC = hdata->GetMean();
+ else if(esdInd==13) pmCZNA = hdata->GetMean();
+ else if(esdInd==14) pmCZPC = hdata->GetMean();
+ else if(esdInd==15) 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");
}
- }
- }
-
- AliInfo(Form("Test Result = %f", test)) ;
- return test ;
+ } // HighMult (Pb-Pb)
+ // ====================================================================
+ // Checks for Calibration events
+ // ====================================================================
+ else if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCalib) {
+ AliWarning(Form("\n\t No check implemented in ZDC QA for %s task\n",taskName));
+ return ;
+ } // Calibration
+ // ====================================================================
+ // Checks for cosmic events
+ // ====================================================================
+ else if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCosmic) {
+ AliWarning(Form("\n\t No check implemented in ZDC QA for %s task\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
}