2 /**************************************************************************
3 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
5 * Author: The ALICE Off-line Project. *
6 * Contributors are mentioned in the code where appropriate. *
8 * Permission to use, copy, modify and distribute this software and its *
9 * documentation strictly for non-commercial purposes is hereby granted *
10 * without fee, provided that the above copyright notice appears in all *
11 * copies and that both the copyright notice and this permission notice *
12 * appear in the supporting documentation. The authors make no claims *
13 * about the suitability of this software for any purpose. It is *
14 * provided "as is" without express or implied warranty. *
15 **************************************************************************/
19 // *****************************************
20 // Checks the quality assurance
21 // by comparing with reference data
22 // P. Cerello Apr 2008
25 // --- ROOT system ---
29 // --- AliRoot header files ---
30 #include "AliITSQASDDChecker.h"
32 #include "AliCDBEntry.h"
33 #include "AliCDBManager.h"
34 #include "AliQACheckerBase.h"
38 ClassImp(AliITSQASDDChecker)
39 //__________________________________________________________________
40 AliITSQASDDChecker& AliITSQASDDChecker::operator = (const AliITSQASDDChecker& qac )
43 this->~AliITSQASDDChecker();
44 new(this) AliITSQASDDChecker(qac);
48 //__________________________________________________________________
49 const Double_t AliITSQASDDChecker::Check(AliQA::ALITASK_t index, TObjArray * list)
51 AliDebug(1,Form("AliITSQASDDChecker called with offset: %d\n", fSubDetOffset));
52 AliCDBEntry *QARefObj = AliCDBManager::Instance()->Get("ITS/QARef/SDD");
54 AliError("Calibration object retrieval failed! SDD will not be processed");
61 if(index==0){ //analizing RAWS
62 TH1F *ModPattern = (TH1F*)QARefObj->GetObject();
63 if (list->GetEntries() == 0){
64 test = 1. ; // nothing to check
69 for(offset =0;offset < fSubDetOffset; offset++){
70 hdata = dynamic_cast<TH1 *>(next());
73 while ( (hdata = dynamic_cast<TH1 *>(next())) && offset >= fSubDetOffset){
75 if(offset == fSubDetOffset){ //ModPattern check
76 if ( hdata->Integral() == 0 ) {
77 AliWarning(Form("Spectrum %s is empty", hdata->GetName())) ;
80 test = hdata->Chi2Test(ModPattern,"UU,p");
84 AliError("Data type cannot be processed") ;
89 while ( (hdata = dynamic_cast<TH1 *>(next()))) {
93 AliInfo(Form("Test Result for RAWS = %f", test)) ;
97 if( index==2){ //analizing RECP
98 //printf("analizing recp, offset %d \n",fSubDetOffset);
99 if (list->GetEntries() == 0){
100 test = 1. ; // nothing to check
105 for(offset =0;offset < fSubDetOffset; offset++){
106 hdata = dynamic_cast<TH1 *>(next()); // magari TIter++ ??
107 //printf("Skipping histo %s, offset %d \n",hdata->GetName(),fSubDetOffset);
110 while ( (hdata = dynamic_cast<TH1 *>(next())) && offset >= fSubDetOffset ){
111 if (hdata) { // offset=9 ModPatternRP
112 //printf("Treating histo %s, offset %d \n",hdata->GetName(),fSubDetOffset);
113 if( offset == 9 && hdata->GetEntries()>0)test = 0.1;
116 AliError("Data type cannot be processed") ;
121 AliInfo(Form("Test Result for RECP = %f", test)) ;
122 } // if(index==2) loop
127 //__________________________________________________________________
128 void AliITSQASDDChecker::SetTaskOffset(Int_t TaskOffset)
130 fSubDetOffset = TaskOffset;