1 /**************************************************************************
2 * Copyright(c) 2007-2009, 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 **************************************************************************/
18 // *****************************************
19 // Checks the quality assurance
20 // by comparing with reference data
21 // W.Ferrarese P.Cerello Mag 2008
24 // --- ROOT system ---
26 #include <Riostream.h>
28 // --- AliRoot header files ---
29 #include "AliITSQAChecker.h"
30 #include "AliITSQASPDChecker.h"
31 #include "AliITSQASDDChecker.h"
32 #include "AliITSQASSDChecker.h"
34 ClassImp(AliITSQAChecker)
36 //____________________________________________________________________________
37 AliITSQAChecker::AliITSQAChecker(Bool_t kMode, Short_t subDet, Short_t ldc) :
38 AliQACheckerBase("ITS","SDD Quality Assurance Checker"),
48 fSPDChecker(0), // SPD Checker
49 fSDDChecker(0), // SDD Checker
50 fSSDChecker(0) // SSD Checker
53 // Standard constructor
54 fkOnline = kMode; fDet = subDet; fLDC = ldc;
55 if(fDet == 0 || fDet == 1) {
56 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQAChecker::Create SPD Checker\n");
57 fSPDChecker = new AliITSQASPDChecker();
59 if(fDet == 0 || fDet == 2) {
60 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQAChecker::Create SDD Checker\n");
61 fSDDChecker = new AliITSQASDDChecker();
63 if(fDet == 0 || fDet == 3) {
64 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQAChecker::Create SSD Checker\n");
65 fSSDChecker = new AliITSQASSDChecker();
67 InitQACheckerLimits();
70 //____________________________________________________________________________
71 AliITSQAChecker::AliITSQAChecker(const AliITSQAChecker& qac):
72 AliQACheckerBase(qac.GetName(), qac.GetTitle()),
73 fkOnline(qac.fkOnline),
76 fSPDOffset(qac.fSPDOffset),
77 fSDDOffset(qac.fSDDOffset),
78 fSSDOffset(qac.fSSDOffset),
79 fSPDHisto(qac.fSPDHisto),
80 fSDDHisto(qac.fSDDHisto),
81 fSSDHisto(qac.fSSDHisto),
82 fSPDChecker(qac.fSPDChecker),
83 fSDDChecker(qac.fSDDChecker),
84 fSSDChecker(qac.fSSDChecker)
87 AliError("Copy should not be used with this class\n");
89 //____________________________________________________________________________
90 AliITSQAChecker& AliITSQAChecker::operator=(const AliITSQAChecker& qac){
91 // assignment operator
92 this->~AliITSQAChecker();
93 new(this)AliITSQAChecker(qac);
98 //____________________________________________________________________________
99 AliITSQAChecker::~AliITSQAChecker(){
101 if(fSPDChecker)delete fSPDChecker;
102 if(fSDDChecker)delete fSDDChecker;
103 if(fSSDChecker)delete fSSDChecker;
106 //____________________________________________________________________________
107 Double_t * AliITSQAChecker::Check(AliQAv1::ALITASK_t index, TObjArray ** list, const AliDetectorRecoParam * recoParam)
111 // basic checks on the QA histograms on the input list
112 //for the ITS subdetectorQA (Raws Digits Hits RecPoints SDigits) return the worst value of the three result
113 if(index == AliQAv1::kESD){
115 Double_t * rv = new Double_t[AliRecoParam::kNSpecies] ;
116 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
118 if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) )
120 AliDebug(AliQAv1::GetQADebugLevel(),"Checker for ESD");
123 // The following flags are set to kTRUE if the corresponding
124 // QA histograms exceed a given quality threshold
125 Bool_t cluMapSA = kFALSE;
126 Bool_t cluMapMI = kFALSE;
127 Bool_t cluMI = kFALSE;
128 Bool_t cluSA = kFALSE;
129 Bool_t verSPDZ = kFALSE;
130 if (list[specie]->GetEntries() == 0) {
131 rv[specie] = 0.; // nothing to check
134 Double_t *stepbit=new Double_t[AliQAv1::kNBIT];
135 Double_t histonumb= list[specie]->GetEntries();
136 CreateStepForBit(histonumb,stepbit);
137 TIter next1(list[specie]);
140 Bool_t skipped[6]={kFALSE,kFALSE,kFALSE,kFALSE,kFALSE,kFALSE};
141 // look for layers that we wanted to skip
142 while ( (hdata = dynamic_cast<TH1 *>(next1())) ) {
144 TString hname = hdata->GetName();
145 if(!hname.Contains("hESDSkippedLayers")) continue;
146 for(Int_t k=1; k<7; k++) {
147 if(hdata->GetBinContent(k)>0) {
153 TIter next(list[specie]);
154 while ( (hdata = dynamic_cast<TH1 *>(next())) ) {
156 TString hname = hdata->GetName();
157 Double_t entries = hdata->GetEntries();
159 if(!(entries>0.))++empty;
160 AliDebug(AliQAv1::GetQADebugLevel(),Form("ESD hist name %s - entries %12.1g",hname.Data(),entries));
161 if(hname.Contains("hESDClusterMapSA") && entries>0.){
163 AliDebug(AliQAv1::GetQADebugLevel(),Form("Processing histogram %s",hname.Data()));
164 // Check if there are layers with anomalously low
165 // contributing points to SA reconstructed tracks
166 for(Int_t k=1;k<7;k++){
167 // check if the layer was skipped
168 if(skipped[k-1]) continue;
169 if(hdata->GetBinContent(k)<0.5*(entries/6.)){
171 AliDebug(AliQAv1::GetQADebugLevel(),Form("SA tracks have few points on layer %d - look at histogram hESDClustersSA",k));
176 else if(hname.Contains("hESDClusterMapMI") && entries>0.){
177 // Check if there are layers with anomalously low
178 // contributing points to MI reconstructed tracks
179 AliDebug(AliQAv1::GetQADebugLevel(),Form("Processing histogram %s",hname.Data()));
181 for(Int_t k=1;k<7;k++){
182 // check if the layer was skipped
183 if(skipped[k-1]) continue;
184 if(hdata->GetBinContent(k)<0.5*(entries/6.)){
186 AliDebug(AliQAv1::GetQADebugLevel(),Form("MI tracks have few points on layer %d - look at histogram hESDClustersMI",k));
191 else if(hname.Contains("hESDClustersMI") && entries>0.){
192 // Check if 6 clusters MI tracks are the majority
193 AliDebug(AliQAv1::GetQADebugLevel(),Form("Processing histogram %s",hname.Data()));
195 Double_t maxlaytracks = hdata->GetBinContent(7-nskipped);
196 for(Int_t k=2; k<7-nskipped; k++){
197 if(hdata->GetBinContent(k)>maxlaytracks){
199 AliDebug(AliQAv1::GetQADebugLevel(),Form("MI Tracks with %d clusters are more than tracks with %d clusters. Look at histogram hESDClustersMI",k-1,6-nskipped));
204 else if(hname.Contains("hESDClustersSA") && entries>0.){
205 // Check if 6 clusters SA tracks are the majority
206 AliDebug(AliQAv1::GetQADebugLevel(),Form("Processing histogram %s",hname.Data()));
208 Double_t maxlaytracks = hdata->GetBinContent(7-nskipped);
209 for(Int_t k=2; k<7-nskipped; k++){
210 if(hdata->GetBinContent(k)>maxlaytracks){
212 AliDebug(AliQAv1::GetQADebugLevel(), Form("SA Tracks with %d clusters are more than tracks with %d clusters. Look at histogram hESDClustersSA",k-1,6-nskipped));
217 else if(hname.Contains("hSPDVertexZ") && entries>0.){
218 // Check if average Z vertex coordinate is -5 < z < 5 cm
219 AliDebug(AliQAv1::GetQADebugLevel(),Form("Processing histogram %s",hname.Data()));
221 if(hdata->GetMean()<-5. && hdata->GetMean()>5.){
223 AliDebug(AliQAv1::GetQADebugLevel(),Form("Average z vertex coordinate is at z= %10.4g cm",hdata->GetMean()));
228 AliError("ESD Checker - invalid data type");
234 rv[specie] = 2500.; // set to error
235 AliWarning(Form("All ESD histograms are empty - specie=%d",specie));
238 rv[specie] = 2500.-1500.*(static_cast<Double_t>(tested-empty)/static_cast<Double_t>(tested)); // INFO if all histos are filled
239 if(cluMapSA)rv[specie]-=200.;
240 if(cluMapMI)rv[specie]-=200.;
241 if(cluMI)rv[specie]-=200.;
242 if(cluSA)rv[specie]-=200.;
243 if(verSPDZ)rv[specie]-=199.; // down to 1 if everything is OK
247 // AliDebug(AliQAv1::GetQADebugLevel(), Form("ESD - Tested %d histograms, Return value %f \n",tested,rv[specie]));
248 AliInfo(Form("ESD - Tested %d histograms, Return value %f \n",tested,rv[specie]));
253 Double_t * retval = new Double_t[AliRecoParam::kNSpecies] ;
254 //____________________________________________________________________________
256 Double_t spdCheck[AliRecoParam::kNSpecies] ;
257 Double_t sddCheck[AliRecoParam::kNSpecies] ;
258 Double_t ssdCheck[AliRecoParam::kNSpecies] ;
262 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
263 if ( AliQAv1::Instance()->IsEventSpecieSet(specie) ) {
264 Double_t histotot=list[specie]->GetEntries();
270 retval[specie] = 0.0 ;//
272 if(fDet == 0 || fDet == 1) {
273 fSPDChecker->SetTaskOffset(fSPDOffset);
274 //printf("spdoffset = %i \n",fSPDOffset );
275 Double_t histoSPD=double(GetSPDHisto());
277 Double_t *stepSPD=new Double_t[AliQAv1::kNBIT];
278 CreateStepForBit(histoSPD,stepSPD);
279 fSPDChecker->SetStepBit(stepSPD);
280 spdCheck[specie] = fSPDChecker->Check(index, list[specie], recoParam);
281 if(spdCheck[specie]>fUpTestValue[AliQAv1::kFATAL]||spdCheck[specie]<0.)
283 AliInfo(Form("SPD check result for %s is out of range (%f)!!! Retval of specie %s is sit to -1\n ",AliQAv1::GetAliTaskName(index),spdCheck[specie],AliRecoParam::GetEventSpecieName(specie)));
284 spdCheck[specie]=fUpTestValue[AliQAv1::kFATAL];
286 //if(spdCheck[specie]<0.5)AliInfo(Form("SPD check result for %s (%s) is < 0.5 .The result is %f ",AliQAv1::GetAliTaskName(index),AliRecoParam::GetEventSpecieName(specie),spdCheck[specie]) );
288 retval[specie]=spdCheck[specie];
291 if(fDet == 0 || fDet == 2) {
292 fSDDChecker->SetTaskOffset(fSDDOffset);
293 Double_t histoSDD=double(GetSDDHisto());
294 Double_t *stepSDD=new Double_t[AliQAv1::kNBIT];
295 CreateStepForBit(histoSDD,stepSDD);
296 fSDDChecker->SetStepBit(stepSDD);
297 sddCheck[specie] = fSDDChecker->Check(index, list[specie], recoParam);
298 if(sddCheck[specie]>fUpTestValue[AliQAv1::kFATAL]||sddCheck[specie]<0.)
300 AliInfo(Form("SDD check result for %s is out of range (%f)!!! Retval of specie %s is sit to -1\n ",AliQAv1::GetAliTaskName(index),sddCheck[specie],AliRecoParam::GetEventSpecieName(specie)));
301 sddCheck[specie]=fUpTestValue[AliQAv1::kFATAL];
303 //if(sddCheck[specie]<0.5)AliInfo(Form("SDD check result for %s (%s) is < 0.5 .The result is %f\f ",AliQAv1::GetAliTaskName(index),AliRecoParam::GetEventSpecieName(specie),sddCheck[specie]) );
305 if(sddCheck[specie]>retval[specie])retval[specie]=sddCheck[specie];
308 if(fDet == 0 || fDet == 3) {
309 fSSDChecker->SetTaskOffset(fSSDOffset);
310 Double_t histoSSD=double(GetSSDHisto());
311 Double_t *stepSSD=new Double_t[AliQAv1::kNBIT];
312 CreateStepForBit(histoSSD,stepSSD);
313 fSSDChecker->SetStepBit(stepSSD);
314 ssdCheck[specie] = fSSDChecker->Check(index, list[specie], recoParam);
315 if(ssdCheck[specie]>fUpTestValue[AliQAv1::kFATAL]||ssdCheck[specie]<0.)
317 AliInfo(Form("SSD check result for %s is out of range (%f)!!! Retval of specie %s is sit to -1\n ",AliQAv1::GetAliTaskName(index),ssdCheck[specie],AliRecoParam::GetEventSpecieName(specie)));
318 ssdCheck[specie]=fUpTestValue[AliQAv1::kFATAL];
320 //if(ssdCheck[specie]<0.5)AliInfo(Form("SSD check result for %s (%s) is < 0.5 . The result is %f ",AliQAv1::GetAliTaskName(index),AliRecoParam::GetEventSpecieName(specie),ssdCheck[specie]) );
322 if(ssdCheck[specie]>retval[specie])retval[specie]=ssdCheck[specie];
325 AliInfo(Form("Check result for %s: \n\t SPD %f \n\t SDD %f \n\t SSD %f \n Check result %f \n ",AliQAv1::GetAliTaskName(index),spdCheck[specie],sddCheck[specie],ssdCheck[specie],retval[specie]));
326 // here merging part for common ITS QA result
329 }//end if event specie
335 //____________________________________________________________________________
336 void AliITSQAChecker::SetTaskOffset(Int_t SPDOffset, Int_t SDDOffset, Int_t SSDOffset)
338 //Setting the 3 offsets for each task called
339 fSPDOffset = SPDOffset;
340 fSDDOffset = SDDOffset;
341 fSSDOffset = SSDOffset;
344 //____________________________________________________________________________
345 void AliITSQAChecker::SetHisto(Int_t SPDhisto, Int_t SDDhisto, Int_t SSDhisto)
347 //Setting the 3 offsets for each task called
348 fSPDHisto = SPDhisto;
349 fSDDHisto = SDDhisto;
350 fSSDHisto = SSDhisto;
353 //____________________________________________________________________________
354 void AliITSQAChecker::SetDetTaskOffset(Int_t subdet,Int_t offset)
358 SetSPDTaskOffset(offset);
361 SetSDDTaskOffset(offset);
364 SetSSDTaskOffset(offset);
367 AliWarning("No specific (SPD,SDD or SSD) subdetector correspond to to this number!!! all offsets set to zero for all the detectors\n");
368 SetTaskOffset(0, 0, 0);
373 //____________________________________________________________________________
374 void AliITSQAChecker::SetDetHisto(Int_t subdet,Int_t histo)
387 AliWarning("No specific (SPD,SDD or SSD) subdetector correspond to to this number!!! all offsets set to zero for all the detectors\n");
393 //_____________________________________________________________________________
395 void AliITSQAChecker::InitQACheckerLimits()
398 AliInfo("Setting of tolerance values\n");
400 Float_t lowtolerancevalue[AliQAv1::kNBIT];
402 Float_t hightolerancevalue[AliQAv1::kNBIT];
403 for(Int_t bit=0;bit<AliQAv1::kNBIT;bit++)
405 lowtolerancevalue[bit]=(bit*1000.);
406 hightolerancevalue[bit]=((bit+1.)*1000.);
408 SetHiLo(hightolerancevalue,lowtolerancevalue);
409 // AliInfo(Form("Range Value \n INFO -> %f < value < %f \n WARNING -> %f < value <= %f \n ERROR -> %f < value <= %f \n FATAL -> %f <= value < %f \n", fLowTestValue[AliQAv1::kINFO], fUpTestValue[AliQAv1::kINFO], fLowTestValue[AliQAv1::kWARNING], fUpTestValue[AliQAv1::kWARNING], fLowTestValue[AliQAv1::kERROR], fUpTestValue[AliQAv1::kERROR], fLowTestValue[AliQAv1::kFATAL], fUpTestValue[AliQAv1::kFATAL] ));
411 if(fDet == 0 || fDet == 1) {
412 fSPDChecker->SetSPDLimits( lowtolerancevalue,hightolerancevalue );
414 if(fDet == 0 || fDet == 2) {
415 fSDDChecker->SetSDDLimits( lowtolerancevalue,hightolerancevalue );
417 if(fDet == 0 || fDet == 3) {
418 fSSDChecker->SetSSDLimits( lowtolerancevalue,hightolerancevalue );
426 //_____________________________________________________________________________
428 void AliITSQAChecker::CreateStepForBit(Double_t histonumb,Double_t *steprange)
430 for(Int_t bit=0;bit < AliQAv1::kNBIT; bit++)
432 //printf("%i\t %f \t %f \t %f \n",bit, fUpTestValue[bit],fLowTestValue[AliQAv1::kINFO],histonumb);
433 steprange[bit]=double((fUpTestValue[bit] - fLowTestValue[AliQAv1::kINFO])/histonumb);
434 //printf("%i\t %f \t %f \t %f \t %f\n",bit, fUpTestValue[bit],fLowTestValue[AliQAv1::kINFO],histonumb,steprange[bit] );
436 //AliInfo(Form("StepBitValue:numner of histo %f\n\t INFO %f \t WARNING %f \t ERROR %f \t FATAL %f \n",histonumb, steprange[AliQAv1::kINFO],steprange[AliQAv1::kWARNING],steprange[AliQAv1::kERROR],steprange[AliQAv1::kFATAL]));
440 //_____________________________________________________________________________
441 void AliITSQAChecker::SetQA(AliQAv1::ALITASK_t index, Double_t * value) const
444 AliQAv1 * qa = AliQAv1::Instance(index) ;
447 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
449 if (! qa->IsEventSpecieSet(AliRecoParam::ConvertIndex(specie)))
451 if ( value == NULL ) { // No checker is implemented, set all QA to Fatal
452 qa->Set(AliQAv1::kFATAL, specie) ;
454 if ( value[specie] > fLowTestValue[AliQAv1::kFATAL] && value[specie] <= fUpTestValue[AliQAv1::kFATAL] )
455 qa->Set(AliQAv1::kFATAL, AliRecoParam::ConvertIndex(specie)) ;
456 else if ( value[specie] > fLowTestValue[AliQAv1::kERROR] && value[specie] <= fUpTestValue[AliQAv1::kERROR] )
457 qa->Set(AliQAv1::kERROR, AliRecoParam::ConvertIndex(specie)) ;
458 else if ( value[specie] > fLowTestValue[AliQAv1::kWARNING] && value[specie] <= fUpTestValue[AliQAv1::kWARNING] )
459 qa->Set(AliQAv1::kWARNING, AliRecoParam::ConvertIndex(specie)) ;
460 else if ( value[specie] > fLowTestValue[AliQAv1::kINFO] && value[specie] <= fUpTestValue[AliQAv1::kINFO] )
461 qa->Set(AliQAv1::kINFO, AliRecoParam::ConvertIndex(specie)) ;
462 //else if(value[specie]==0) qa->Set(AliQAv1::kFATAL, AliRecoParam::ConvertIndex(specie)) ; //no ckeck has been done
464 qa->ShowStatus(AliQAv1::kITS,index,AliRecoParam::ConvertIndex(specie));