some warnings fixed according Federico's request
[u/mrichter/AliRoot.git] / T0 / AliT0QAChecker.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16
17
18
19 //...
20 //  Checks the quality assurance. 
21 //  By comparing with reference data
22 //  Skeleton for T0
23 //...
24
25 // --- ROOT system ---
26 #include <Riostream.h>
27 #include <TClass.h>
28 #include <TH1F.h> 
29 #include <TH1I.h> 
30 #include <TIterator.h> 
31 #include <TKey.h> 
32 #include <TFile.h> 
33 #include <TMath.h>
34
35 // --- Standard library ---
36
37 // --- AliRoot header files ---
38 #include "AliLog.h"
39 #include "AliQA.h"
40 #include "AliQAChecker.h"
41 #include "AliT0QAChecker.h"
42
43 ClassImp(AliT0QAChecker)
44
45 //__________________________________________________________________
46 AliT0QAChecker& AliT0QAChecker::operator = (const AliT0QAChecker& qac )
47 {
48   // Equal operator.
49   this->~AliT0QAChecker();
50   new(this) AliT0QAChecker(qac);
51   return *this;
52 }
53 //__________________________________________________________________
54
55 const Double_t AliT0QAChecker::Check(AliQA::ALITASK_t index,TObjArray * list)
56 {
57
58   // Super-basic check on the QA histograms on the input list:
59   // look whether they are empty!
60
61   cout<<" GetAliTaskName "<<AliQA::GetAliTaskName(index)<<" index "<<index<<endl;
62
63
64   Double_t test = 0.0  ;
65   Int_t count = 0 ;
66   Double_t nent[100];
67   memset(nent,0,100*sizeof(Double_t));
68   Double_t w[100];
69   memset(w,1,100*sizeof(Double_t));
70   TH1 *fhRecLEDAmp[24];  TH1 * fhRecQTC[24];
71   TH1 *fhOnlineMean = 0x0;  
72   TH1 * fhRecMean = 0x0;
73   TString dataType = AliQA::GetAliTaskName(index);
74
75   if (list->GetEntries() == 0){
76     test = 1. ; // nothing to check
77   }
78   else {
79     
80     TIter next(list) ;
81     TH1 * hdata ;
82     
83     count = 0 ;
84     while ( (hdata = dynamic_cast<TH1 *>(next())) ) {
85       if (hdata) {
86         nent[count] = hdata->GetEntries();
87         AliDebug(10,Form("count %i %s -> %f",count, hdata->GetName(),nent[count])) ;
88
89         if(index==2){
90           if(count>23 && count<48)fhRecLEDAmp[count-24] = hdata;
91           if(count>47 && count<72)fhRecQTC[count-48] = hdata;
92           if(count == 72)  fhOnlineMean = hdata; 
93           if(count == 73)  fhRecMean = hdata; 
94         }
95         count++ ;
96         Double_t rv = 0.;
97         if(hdata->GetEntries()>0) rv = 1;
98         //   AliInfo(Form("%s -> %f", hdata->GetName(), rv)) ;
99         test += rv ;
100      }
101       else{
102         AliError("Data type cannot be processed") ;
103       }
104
105     }
106
107     if (count != 0) {
108       if (test==0) {
109         AliWarning("Histograms are there, but they are all empty: setting flag to kWARNING");
110         test = 0.5;  //upper limit value to set kWARNING flag for a task
111       }
112       else {
113
114         if(index == 2){
115           if ( TMath::Abs(fhRecMean->GetMean() - fhOnlineMean->GetMean()) > 5) 
116             AliWarning(Form("rec mean %f -> online mean %f",fhRecMean->GetMean(), fhOnlineMean->GetMean())) ;
117           Double_t meanLED, meanQTC;
118           for (Int_t idet=0; idet<24; idet++) {
119             meanLED = fhRecLEDAmp[idet]->GetMean();
120             meanQTC = fhRecQTC[idet]->GetMean();
121             if (TMath::Abs(meanLED-meanQTC)> 1.) 
122             AliWarning(Form("Amplitude measurements are different in channel %i : Amp LED %f -> Amp QTC %f",idet,meanLED, meanQTC)) ;
123           }
124         }       
125         else
126           {
127             AliDebug(10,Form(" MaxElement %f ", TMath::MaxElement(count,nent)));        
128             if(TMath::MaxElement(count,nent) > 1000) {
129               Double_t mean = TMath::Mean(count,nent,w);
130               AliDebug(10,Form(" Mean %f ", mean));     
131               for (Int_t i=0; i<count; i++) 
132                 {
133                   Double_t diff = TMath::Abs(nent[i]-mean);
134                   if (diff > 0.1*mean )
135                     AliWarning(Form("Problem in Number of entried in hist %i  is %f\n", i, nent[i])) ; 
136           }
137         }
138           }
139       }
140     }
141   } //  if (list->GetEntries() != 0
142   AliInfo(Form("Test Result = %f", test)) ;
143   return test ;
144 }