]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
Bug corrected in QA checker
authorcoppedis <coppedis@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 29 Jun 2010 20:43:23 +0000 (20:43 +0000)
committercoppedis <coppedis@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 29 Jun 2010 20:43:23 +0000 (20:43 +0000)
ZDC/AliZDCQAChecker.cxx
ZDC/AliZDCQADataMakerRec.cxx

index b952601d6d5218c44129f9bbc74269aa3cb427a6..505b256cad565de9f19dea7d753f378bc8e2297c 100644 (file)
@@ -36,44 +36,51 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
   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;
+
   for(Int_t specie = 0; specie<AliRecoParam::kNSpecies; specie++){
     Int_t count = 0; 
-    //printf("\tAliZDCQAChecker -> specie %d\n\n",specie);
+    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)));
     
-    if(!AliQAv1::Instance()->IsEventSpecieSet(specie))  continue; 
     // ====================================================================
     //         Checks for p-p events
     // ====================================================================
-    if(specie == AliRecoParam::kLowMult){
+    if(AliRecoParam::ConvertIndex(specie) == AliRecoParam::kLowMult){
       if(list[specie]->GetEntries()==0){  
         AliWarning("\t The list to be checked is empty!"); // nothing to check
-        return ;
+        return;
       }
       //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.15;
+      
+      Float_t res=0., percentageDiff=0.20;
+      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
-            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();
+               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();
@@ -86,7 +93,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
              //
              // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
              if(ihitHisto==11){
-               if(sumZNC!=0){
+               if(TMath::Abs(meanZNC)>1.e-10){
                   if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
                     res=1.;
                   else 
@@ -94,7 +101,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumZNA!=0){
+               if(TMath::Abs(meanZNA)>1.e-10){
                   if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
                     res=1.;
                   else percentageDiff=
@@ -102,7 +109,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumZPC!=0){
+               if(TMath::Abs(meanZPC)>1.e-10){
                   if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
                     res=1.;
                   else 
@@ -110,7 +117,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumZPA!=0){
+               if(TMath::Abs(meanZPA)>1.e-10){
                   if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
                     res=1.;
                   else 
@@ -140,7 +147,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
              //
              // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
              if(idigHisto==11){
-               if(sumADCZNC!=0){
+               if(TMath::Abs(sumADCZNC)>1.e-10){
                   if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
                     res=1.;
                   else 
@@ -148,7 +155,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumADCZNA!=0){
+               if(TMath::Abs(sumADCZNA)>1.e-10){
                   if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
                     res=1.;
                   else 
@@ -156,7 +163,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumADCZPC!=0){
+               if(TMath::Abs(sumADCZPC)>1.e-10){
                   if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
                     res=1.;
                   else 
@@ -164,7 +171,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumADCZPA!=0){
+               if(TMath::Abs(sumADCZPA)>1.e-10){
                   if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
                     res=1.;
                   else 
@@ -178,10 +185,6 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
           } 
           // -------------------------------------------------------------------
          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){
@@ -200,8 +203,8 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
            }
            //
            // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
-           if(irawHisto==11){
-             if(sumADCZNC!=0){
+           if(irawHisto==13){
+             if(TMath::Abs(sumADCZNC)>1.e-10){
                if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
                  res=1.;
                else 
@@ -209,7 +212,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZNA!=0){
+             if(TMath::Abs(sumADCZNA)>1.e-10){
                if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
                  res=1.;
                else 
@@ -217,7 +220,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPC!=0){
+             if(TMath::Abs(sumADCZPC)>1.e-10){
                if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
                  res=1.;
                else 
@@ -225,7 +228,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPA!=0){
+             if(TMath::Abs(sumADCZPA)>1.e-10){
                if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
                  res=1.;
                else 
@@ -238,17 +241,13 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
           } 
           // -------------------------------------------------------------------
          else if(index == AliQAv1::kREC) {
-            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.;
-           Int_t irecHisto=0;
            //
             // Check REC HIGH GAIN CHAIN histos
             if(hdata->GetEntries()>0){
-             if(irecHisto==0)       sumZNC = hdata->GetMean();
-             else if(irecHisto==1)  sumZNA = hdata->GetMean();
-             else if(irecHisto==2)  sumZPC = hdata->GetMean();
-             else if(irecHisto==3)  sumZPA = hdata->GetMean();
+             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();
@@ -261,7 +260,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
            //
            // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
            if(irecHisto==11){
-             if(sumZNC!=0){
+             if(TMath::Abs(meanZNC)>1.e-10){
                if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
                  res=1.;
                else 
@@ -269,7 +268,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumZNA!=0){
+             if(TMath::Abs(meanZNA)>1.e-10){
                if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
                  res=1.;
                else 
@@ -277,7 +276,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumZPC!=0){
+             if(TMath::Abs(meanZPC)>1.e-10){
                if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
                  res=1.;
                else 
@@ -285,7 +284,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumZPA!=0){
+             if(TMath::Abs(meanZPA)>1.e-10){
                if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
                  res=1.;
                else 
@@ -298,19 +297,15 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
           } 
           // -------------------------------------------------------------------
          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(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==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();
@@ -320,8 +315,8 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
            }
            //
            // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
-           if(esdInd==13){
-             if(sumADCZNC!=0){
+           if(esdInd==15){
+             if(TMath::Abs(sumADCZNC)>1.e-10){
                if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
                  res=1.;
                else 
@@ -329,7 +324,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZNA!=0){
+             if(TMath::Abs(sumADCZNA)>1.e-10){
                if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
                  res=1.;
                else 
@@ -337,7 +332,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPC!=0){
+             if(TMath::Abs(sumADCZPC)>1.e-10){
                if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
                  res=1.;
                else 
@@ -345,7 +340,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPA!=0){
+             if(TMath::Abs(sumADCZPA)>1.e-10){
                if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
                  res=1.;
                else 
@@ -367,7 +362,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
     // ====================================================================
     //         Checks for A-A events
     // ====================================================================
-    else if (specie == AliRecoParam::kHighMult) {
+    else if (AliRecoParam::ConvertIndex(specie) == AliRecoParam::kHighMult) {
       if(list[specie]->GetEntries()==0){  
         AliWarning("\t The list to be checked is empty!");
         return ;
@@ -377,40 +372,41 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
       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){
-            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();
+               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(hitInd==3){
-               if(meanZNC!=0 && meanZNA!=0){
+             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 
@@ -418,7 +414,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(meanZPC!=0 && meanZPA!=0){
+               if(TMath::Abs(meanZPC)>1.e-10 && TMath::Abs(meanZPA)>1.e-10){
                   if((TMath::Abs(meanZPC-meanZPA)/(meanZPA+meanZPC))<percentageDiff) 
                     res=1.;
                   else 
@@ -428,8 +424,8 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                }
               }
              // --- Check whether (mean PMQi - PMC)/PMC < percentageDiff
-             if(hitInd==11){
-               if(meanZNC!=0){
+             if(ihitHisto==11){
+               if(TMath::Abs(meanZNC)>1.e-10){
                   if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
                     res=1.;
                   else 
@@ -437,7 +433,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(meanZNA!=0){
+               if(TMath::Abs(meanZNA)>1.e-10){
                   if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
                     res=1.;
                   else 
@@ -445,7 +441,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(meanZPC!=0){
+               if(TMath::Abs(meanZPC)>1.e-10){
                   if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
                     res=1.;
                   else 
@@ -453,7 +449,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(meanZPA!=0){
+               if(TMath::Abs(meanZPA)>1.e-10){
                   if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
                     res=1.;
                   else 
@@ -462,29 +458,29 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   count++;
                }
              }
-             hitInd++;
+             ihitHisto++;
             }
             // 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();
+               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(digInd==3){
-               if(sumADCZNC!=0 && sumADCZNA!=0){
+             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 
@@ -492,7 +488,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumADCZPC!=0 && sumADCZPA!=0){
+               if(TMath::Abs(sumADCZPC)>1.e-10 && TMath::Abs(sumADCZPA)>1.e-10){
                   if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) 
                     res=1.;
                   else 
@@ -502,8 +498,8 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                }
               }
              // --- Check whether (sumADC PMQi - PMC)/PMC < percentageDiff
-             if(digInd==11){
-               if(sumADCZNC!=0){
+             if(idigHisto==11){
+               if(TMath::Abs(sumADCZNC)>1.e-10){
                   if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
                     res=1.;
                   else 
@@ -511,7 +507,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumADCZNA!=0){
+               if(TMath::Abs(sumADCZNA)>1.e-10){
                   if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
                     res=1.;
                   else 
@@ -519,7 +515,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumADCZPC!=0){
+               if(TMath::Abs(sumADCZPC)>1.e-10){
                   if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
                     res=1.;
                   else 
@@ -527,7 +523,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   test[specie] += res;
                   count++;
                }
-               if(sumADCZPA!=0){
+               if(TMath::Abs(sumADCZPA)>1.e-10){
                   if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
                     res=1.;
                   else 
@@ -536,36 +532,32 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                   count++;
                }
              }
-              digInd++;
+              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 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();
+             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(rawInd==3){
-             if(sumADCZNC!=0 && sumADCZNA!=0){
+           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 
@@ -573,7 +565,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPC!=0 && sumADCZPA!=0){
+             if(TMath::Abs(sumADCZPC)>1.e-10 && TMath::Abs(sumADCZPA)>1.e-10){
                if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) 
                  res=1.;
                else 
@@ -583,8 +575,8 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
              }
             }
            // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
-           if(rawInd==11){
-             if(sumADCZNC!=0){
+           if(irawHisto==13){
+             if(TMath::Abs(sumADCZNC)>1.e-10){
                if((TMath::Abs(adcQZNC-adcCZNC)/adcCZNC)<percentageDiff) 
                  res=1.;
                else 
@@ -592,7 +584,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZNA!=0){
+             if(TMath::Abs(sumADCZNA)>1.e-10){
                if((TMath::Abs(adcQZNA-adcCZNA)/adcCZNA)<percentageDiff) 
                  res=1.;
                else 
@@ -600,7 +592,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPC!=0){
+             if(TMath::Abs(sumADCZPC)>1.e-10){
                if((TMath::Abs(adcQZPC-adcCZPC)/adcCZPC)<percentageDiff) 
                  res=1.;
                else 
@@ -608,7 +600,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPA!=0){
+             if(TMath::Abs(sumADCZPA)>1.e-10){
                if((TMath::Abs(adcQZPA-adcCZPA)/adcCZPA)<percentageDiff) 
                  res=1.;
                else 
@@ -617,44 +609,40 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                count++;
              }
            }
-           rawInd++;    
+           irawHisto++;         
          }   
           // -------------------------------------------------------------------
           else if(index == AliQAv1::kREC){
-            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.;
-           Int_t recInd=0;
            //
             // Check RAW HIGH GAIN CHAIN histos
             if(hdata->GetEntries()>0){
-             if(recInd==0)       sumZNC = hdata->GetMean();
-             else if(recInd==1)  sumZNA = hdata->GetMean();
-             else if(recInd==2)  sumZPC = hdata->GetMean();
-             else if(recInd==3)  sumZPA = hdata->GetMean();
-             else if(recInd==4)  pmQZNC = hdata->GetMean();
-             else if(recInd==5)  pmQZNA = hdata->GetMean();
-             else if(recInd==6)  pmQZPC = hdata->GetMean();
-             else if(recInd==7)  pmQZPA = hdata->GetMean();
-             else if(recInd==8)  pmCZNC = hdata->GetMean();
-             else if(recInd==9)  pmCZNA = hdata->GetMean();
-             else if(recInd==10) pmCZPC = hdata->GetMean();
-             else if(recInd==11) pmCZPA = hdata->GetMean();
+             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(recInd==3){
-             if(sumZNC!=0 && sumZNA!=0){
-               if((2*TMath::Abs(sumZNC-sumZNA)/(sumZNA+sumZNC))<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(sumZPC!=0 && sumZPA!=0){
-               if((TMath::Abs(sumZPC-sumZPA)/(sumZPA+sumZPC))<percentageDiff) 
+             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;
@@ -663,8 +651,8 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
              }
             }
            // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
-           if(recInd==11){
-             if(sumZNC!=0){
+           if(irecHisto==11){
+             if(TMath::Abs(meanZNC)>1.e-10){
                if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
                  res=1.;
                else 
@@ -672,7 +660,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumZNA!=0){
+             if(TMath::Abs(meanZNA)>1.e-10){
                if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
                  res=1.;
                else 
@@ -680,7 +668,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumZPC!=0){
+             if(TMath::Abs(meanZPC)>1.e-10){
                if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
                  res=1.;
                else 
@@ -688,7 +676,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumZPA!=0){
+             if(TMath::Abs(meanZPA)>1.e-10){
                if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
                  res=1.;
                else 
@@ -697,35 +685,31 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                count++;
              }
            }
-           recInd++;    
+           irecHisto++;         
          }   
           // -------------------------------------------------------------------
           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(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();
+             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(sumADCZNC!=0 && sumADCZNA!=0){
+             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 
@@ -733,7 +717,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPC!=0 && sumADCZPA!=0){
+             if(TMath::Abs(sumADCZPC)>1.e-10 && TMath::Abs(sumADCZPA)>1.e-10){
                if((TMath::Abs(sumADCZPC-sumADCZPA)/(sumADCZPA+sumADCZPC))<percentageDiff) 
                  res=1.;
                else 
@@ -743,8 +727,8 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
              }
             }
            // --- Check whether (sum PMQi - PMC)/PMC < percentageDiff
-           if(esdInd==13){
-             if(sumADCZNC!=0){
+           if(esdInd==15){
+             if(TMath::Abs(sumADCZNC)>1.e-10){
                if((TMath::Abs(pmQZNC-pmCZNC)/pmCZNC)<percentageDiff) 
                  res=1.;
                else 
@@ -752,7 +736,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZNA!=0){
+             if(TMath::Abs(sumADCZNA)>1.e-10){
                if((TMath::Abs(pmQZNA-pmCZNA)/pmCZNA)<percentageDiff) 
                  res=1.;
                else 
@@ -760,7 +744,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPC!=0){
+             if(TMath::Abs(sumADCZPC)>1.e-10){
                if((TMath::Abs(pmQZPC-pmCZPC)/pmCZPC)<percentageDiff) 
                  res=1.;
                else 
@@ -768,7 +752,7 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
                test[specie] += res;
                count++;
              }
-             if(sumADCZPA!=0){
+             if(TMath::Abs(sumADCZPA)>1.e-10){
                if((TMath::Abs(pmQZPA-pmCZPA)/pmCZPA)<percentageDiff) 
                  res=1.;
                else 
@@ -790,18 +774,18 @@ void AliZDCQAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArra
     // ====================================================================
     //         Checks for Calibration events
     // ====================================================================
-    else if (specie == AliRecoParam::kCalib) {
+    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 (specie == AliRecoParam::kCosmic) {
+    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(count!=0) test[specie] = test[specie]/count;
+    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
 }  
index 6f04ac00934ca5e78454d71bfca34be8cef4ef4a..311637874d01b1d2ac19ddf0b2d76e7c2ed63a3e 100644 (file)
@@ -628,7 +628,7 @@ void AliZDCQADataMakerRec::MakeESDs(AliESDEvent * esd)
     Float_t beamEne = esd->GetBeamEnergy();
     zdcESD->GetZNCentroidInPbPb(beamEne, centr_ZNC, centr_ZNA);
   }
-  else printf(" AliZDCQADataMakerRec::MakeESDs: can't calculate centroids for beam type: %s\n\n",beamType.Data());
+  else printf("\n WARNING!!! AliZDCQADataMakerRec::MakeESDs: can't calculate centroids for beam type: %s\n\n",beamType.Data());
   GetESDsData(0)->Fill(centr_ZNC[0], centr_ZNC[1]);
   GetESDsData(1)->Fill(centr_ZNA[0], centr_ZNA[1]);
 
@@ -705,43 +705,42 @@ void AliZDCQADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArr
 {
   //Detector specific actions at end of cycle
   // do the QA checking
+  printf("  AliZDCQADataMakerRec::EndOfDetectorCycle for task %d \n", task);
   if( task == AliQAv1::kRAWS){
-     
-    if (!GetRawsData(4) || !GetRawsData(5) || !GetRawsData(6) || !GetRawsData(7) || 
+     if (!GetRawsData(4) || !GetRawsData(5) || !GetRawsData(6) || !GetRawsData(7) || 
         !GetRawsData(8) || !GetRawsData(9) || !GetRawsData(10) || !GetRawsData(11) || 
        !GetRawsData(12) || !GetRawsData(13)) {
         printf("  WARNING!!! AliZDCQADataMaker Rec -> No histogram for DQM found!\n"); 
-        return;
      }
+     else{
+       TLine* diag = new TLine(7., 7., 1407., 1407.);
+       diag->SetLineColor(kRed);
+       diag->SetLineWidth(2);
      
-     TLine* diag = new TLine(7., 7., 1407., 1407.);
-     diag->SetLineColor(kRed);
-     diag->SetLineWidth(2);
-     
-     ((TH2F*)GetRawsData(6))->GetListOfFunctions()->Add(diag);
-     ((TH2F*)GetRawsData(7))->GetListOfFunctions()->Add(diag);
-     ((TH2F*)GetRawsData(8))->GetListOfFunctions()->Add(diag);
-     ((TH2F*)GetRawsData(9))->GetListOfFunctions()->Add(diag);
-
-     GetRawsData(6)->SetOption("colz");
-     GetRawsData(7)->SetOption("colz");
-     GetRawsData(8)->SetOption("colz");
-     GetRawsData(9)->SetOption("colz");  
+       ((TH2F*)GetRawsData(6))->GetListOfFunctions()->Add(diag);
+       ((TH2F*)GetRawsData(7))->GetListOfFunctions()->Add(diag);
+       ((TH2F*)GetRawsData(8))->GetListOfFunctions()->Add(diag);
+       ((TH2F*)GetRawsData(9))->GetListOfFunctions()->Add(diag);
+      
+       GetRawsData(6)->SetOption("colz");
+       GetRawsData(7)->SetOption("colz");
+       GetRawsData(8)->SetOption("colz");
+       GetRawsData(9)->SetOption("colz");  
      
-     GetRawsData(4)->SetLineColor(kBlue+1);  GetRawsData(4)->SetLineWidth(2);
-     GetRawsData(5)->SetLineColor(kBlue+2);  GetRawsData(5)->SetLineWidth(2);
-     GetRawsData(10)->SetLineColor(kBlue+3); GetRawsData(10)->SetLineWidth(2);
-     GetRawsData(11)->SetLineColor(kBlue+4); GetRawsData(11)->SetLineWidth(2);
-     GetRawsData(12)->SetLineColor(kBlue+5); GetRawsData(12)->SetLineWidth(2);
-     GetRawsData(13)->SetLineColor(kBlue+6); GetRawsData(13)->SetLineWidth(2);
+       GetRawsData(4)->SetLineColor(kBlue+1);  GetRawsData(4)->SetLineWidth(2);
+       GetRawsData(5)->SetLineColor(kBlue+2);  GetRawsData(5)->SetLineWidth(2);
+       GetRawsData(10)->SetLineColor(kBlue+3); GetRawsData(10)->SetLineWidth(2);
+       GetRawsData(11)->SetLineColor(kBlue+4); GetRawsData(11)->SetLineWidth(2);
+       GetRawsData(12)->SetLineColor(kBlue+5); GetRawsData(12)->SetLineWidth(2);
+       GetRawsData(13)->SetLineColor(kBlue+6); GetRawsData(13)->SetLineWidth(2);
      
-     /*GetRawsData(4)->SetDrawOption("LOGY");
-     GetRawsData(5)->SetDrawOption("LOGY");
-     GetRawsData(10)->SetDrawOption("LOGY");
-     GetRawsData(11)->SetDrawOption("LOGY");
-     GetRawsData(12)->SetDrawOption("LOGY");
-     GetRawsData(13)->SetDrawOption("LOGY");*/
-
+       /*GetRawsData(4)->SetDrawOption("LOGY");
+       GetRawsData(5)->SetDrawOption("LOGY");
+       GetRawsData(10)->SetDrawOption("LOGY");
+       GetRawsData(11)->SetDrawOption("LOGY");
+       GetRawsData(12)->SetDrawOption("LOGY");
+       GetRawsData(13)->SetDrawOption("LOGY");*/
+     }
   }
        
   AliQAChecker::Instance()->Run(AliQAv1::kZDC, task, list) ;