1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
20 // Checks the quality assurance.
21 // By comparing with reference data
25 // --- ROOT system ---
26 #include <Riostream.h>
30 #include <TIterator.h>
35 // --- Standard library ---
37 // --- AliRoot header files ---
40 #include "AliQAChecker.h"
41 #include "AliT0QAChecker.h"
43 ClassImp(AliT0QAChecker)
45 //__________________________________________________________________
46 AliT0QAChecker& AliT0QAChecker::operator = (const AliT0QAChecker& qac )
49 this->~AliT0QAChecker();
50 new(this) AliT0QAChecker(qac);
53 //__________________________________________________________________
55 const Double_t AliT0QAChecker::Check(AliQA::ALITASK_t index,TObjArray * list)
58 // Super-basic check on the QA histograms on the input list:
59 // look whether they are empty!
61 cout<<" GetAliTaskName "<<AliQA::GetAliTaskName(index)<<" index "<<index<<endl;
67 memset(nent,0,100*sizeof(Double_t));
69 memset(w,1.,100*sizeof(Double_t));
70 TH1 *fhRecLEDAmp[24]; TH1 * fhRecQTC[24];
71 TH1 *fhOnlineMean; TH1 * fhRecMean;
72 TString dataType = AliQA::GetAliTaskName(index);
74 if (list->GetEntries() == 0){
75 test = 1. ; // nothing to check
83 while ( (hdata = dynamic_cast<TH1 *>(next())) ) {
85 nent[count] = hdata->GetEntries();
86 AliDebug(10,Form("count %i %s -> %f",count, hdata->GetName(),nent[count])) ;
89 if(count>23 && count<48)fhRecLEDAmp[count-24] = hdata;
90 if(count>47 && count<72)fhRecQTC[count-48] = hdata;
91 if(count == 72) fhOnlineMean = hdata;
92 if(count == 73) fhRecMean = hdata;
96 if(hdata->GetEntries()>0) rv = 1;
97 // AliInfo(Form("%s -> %f", hdata->GetName(), rv)) ;
101 AliError("Data type cannot be processed") ;
108 AliWarning("Histograms are there, but they are all empty: setting flag to kWARNING");
109 test = 0.5; //upper limit value to set kWARNING flag for a task
114 if ( TMath::Abs(fhRecMean->GetMean() - fhOnlineMean->GetMean()) > 5)
115 AliWarning(Form("rec mean %f -> online mean %f",fhRecMean->GetMean(), fhOnlineMean->GetMean())) ;
116 Double_t meanLED, meanQTC;
117 for (Int_t idet=0; idet<24; idet++) {
118 meanLED = fhRecLEDAmp[idet]->GetMean();
119 meanQTC = fhRecQTC[idet]->GetMean();
120 if (TMath::Abs(meanLED-meanQTC)> 1.)
121 AliWarning(Form("Amplitude measurements are different in channel %i : Amp LED %f -> Amp QTC %f",idet,meanLED, meanQTC)) ;
126 AliDebug(10,Form(" MaxElement %f ", TMath::MaxElement(count,nent)));
127 if(TMath::MaxElement(count,nent) > 1000) {
128 Double_t mean = TMath::Mean(count,nent,w);
129 AliDebug(10,Form(" Mean %f ", mean));
130 for (Int_t i=0; i<count; i++)
132 Double_t diff = TMath::Abs(nent[i]-mean);
133 if (diff > 0.1*mean )
134 AliWarning(Form("Problem in Number of entried in hist %i is %f\n", i, nent[i])) ;
140 } // if (list->GetEntries() != 0
141 AliInfo(Form("Test Result = %f", test)) ;