// --- ROOT system ---
#include "TH1.h"
#include <TCanvas.h>
+#include <TH1D.h>
+#include <TH2.h>
+#include <TF1.h>
+
+
// --- AliRoot header files ---
#include "AliITSQASDDChecker.h"
#include "AliLog.h"
#include "AliCDBEntry.h"
-#include "AliCDBManager.h"
+#include "AliQAManager.h"
#include "AliQACheckerBase.h"
#include "TSystem.h"
-
+#include "AliITSCalibrationSDD.h"
+#include "AliITSgeomTGeo.h"
ClassImp(AliITSQASDDChecker)
//__________________________________________________________________
return *this;
}
+AliITSQASDDChecker::~AliITSQASDDChecker()
+{
+ if(fStepBitSDD)
+ {
+ delete[] fStepBitSDD ;
+ fStepBitSDD = NULL;
+ }
+ if(fLowSDDValue)
+ {
+ delete[]fLowSDDValue;
+ fLowSDDValue=NULL;
+ }
+ if(fHighSDDValue)
+ {
+ delete[]fHighSDDValue;
+ fHighSDDValue=NULL;
+ }
+ if(fCalibration)
+ {
+ delete fCalibration;
+ fCalibration=NULL;
+ }
+} // dtor
+
//__________________________________________________________________
-Double_t AliITSQASDDChecker::Check(AliQA::ALITASK_t index, TObjArray * list)
+Double_t AliITSQASDDChecker::Check(AliQAv1::ALITASK_t index, TObjArray * list, const AliDetectorRecoParam * /*recoparam*/)
{
+ AliInfo(Form("AliITSQASDDChecker called with offset: %d\n", fSubDetOffset) );
+
AliDebug(1,Form("AliITSQASDDChecker called with offset: %d\n", fSubDetOffset));
- AliCDBEntry *QARefObj = AliCDBManager::Instance()->Get("ITS/QARef/SDD");
- if( !QARefObj){
- AliError("Calibration object retrieval failed! SDD will not be processed");
- return 1.;
- }
-
- Double_t test = 0.0;
- Int_t offset = 0;
-
- if(index==AliQA::kRAW){ //analizing RAWS
- TH1F *ModPattern = (TH1F*)QARefObj->GetObject();
- if (list->GetEntries() == 0){
- test = 1. ; // nothing to check
- }
- else {
- TIter next(list) ;
- TH1 * hdata ;
- for(offset =0;offset < fSubDetOffset; offset++){
- hdata = dynamic_cast<TH1 *>(next());
- }
- while ( (hdata = dynamic_cast<TH1 *>(next())) && offset >= fSubDetOffset){
- if (hdata) {
- if(offset == fSubDetOffset){ //ModPattern check
- if ( hdata->Integral() == 0 ) {
- AliWarning(Form("Spectrum %s is empty", hdata->GetName())) ;
- return 0.5 ;
- }
- test = hdata->Chi2Test(ModPattern,"UU,p");
- } // ModPattern check
- }
- else{
- AliError("Data type cannot be processed") ;
- }
- offset++;
- } //SDD histo
+ Double_t test = 0.;
+ TH1 *hdata=NULL;
+ Double_t entries=0.;
+ Double_t entries2[2];
+ for(Int_t i=0;i<2;i++)entries2[i]=0.;
+ switch(index)
+ {
- while ( (hdata = dynamic_cast<TH1 *>(next()))) {
- offset++;
- }
- } //else entries !=0
- AliInfo(Form("Test Result for RAWS = %f", test)) ;
- } // if(index==0)
-
-
- if( index==AliQA::kREC){ //analizing RECP
- //printf("analizing recp, offset %d \n",fSubDetOffset);
- if (list->GetEntries() == 0){
- test = 1. ; // nothing to check
- }
- else {
- TIter next(list) ;
- TH1 * hdata ;
- for(offset =0;offset < fSubDetOffset; offset++){
- hdata = dynamic_cast<TH1 *>(next()); // magari TIter++ ??
- //printf("Skipping histo %s, offset %d \n",hdata->GetName(),fSubDetOffset);
+ case AliQAv1::kRAW:
+ AliInfo(Form("Check on %s\n",AliQAv1::GetAliTaskName(index)));
+ if(!fCalibration){
+ AliCDBEntry *calibSDD = AliCDBManager::Instance()->Get("ITS/Calib/CalibSDD");
+ Bool_t cacheStatus = AliCDBManager::Instance()->GetCacheFlag();
+ if(!calibSDD)
+ {
+ AliError("Calibration object retrieval failed! SDD will not be processed");
+ fCalibration = NULL;
+ test= fHighSDDValue[AliQAv1::kWARNING];
+ }
+ fCalibration = (TObjArray *)calibSDD->GetObject();
+
+ if(!cacheStatus)calibSDD->SetObject(NULL);
+ calibSDD->SetOwner(kTRUE);
+ if(!cacheStatus)
+ {
+ delete calibSDD;
+ }
}
-
- while ( (hdata = dynamic_cast<TH1 *>(next())) && offset >= fSubDetOffset ){
- if (hdata) { // offset=9 ModPatternRP
- //printf("Treating histo %s, offset %d \n",hdata->GetName(),fSubDetOffset);
- if( offset == 9 && hdata->GetEntries()>0)test = 0.1;
+ AliInfo("Calib SDD Created\n ");
+
+ if (list->GetEntries() == 0.){ //check if the list is empty
+ //printf("test = %f \t value %f\n",test,fHighSDDValue[AliQAv1::kFATAL]);
+ test=test+fHighSDDValue[AliQAv1::kFATAL];
+ break;
+ }//end if getentries
+ else{
+ TIter next(list);
+ Int_t offset = 0;
+ for(offset =0;offset < fSubDetOffset; offset++){hdata = dynamic_cast<TH1*>(next());}//end for
+ Int_t emptymodules[2];
+ Int_t filledmodules[2];
+ Int_t emptyladders[2];
+ Int_t filledladders[2];
+ for(Int_t i=0;i<2;i++){
+ emptymodules[i]=0;
+ filledmodules[i]=0;
+ emptyladders[i]=0;
+ filledladders[i]=0;
}
+ TH1 *hmodule=NULL;
+ TH2 *hlayer[2];
+ for(Int_t i=0;i<2;i++)hlayer[i]=NULL;
+ while( hdata = dynamic_cast<TH1* >(next()) ){
+ if (hdata){
+ TString hname=hdata->GetName();
+ if(hname.Contains("SDDchargeMap"))continue;
+ if(hname.Contains("SDDModPattern")){
+ hmodule=hdata;
+ entries= hdata->GetEntries();
+ if(entries==0){
+ AliWarning(Form("===================>>>>>> No entries in %s \n",hname.Data()));
+ //printf("test = %f \t value %f\n",test,fHighSDDValue[AliQAv1::kFATAL]);
+ test=test+fStepBitSDD[AliQAv1::kFATAL];
+ }//endif entries
+ else{
+ int modmax=hdata->GetNbinsX();
+ Int_t empty=0;
+ Int_t filled=0;
+ Double_t content=0;
+ for(Int_t i=1;i<=modmax;i++){
+ content=hdata->GetBinContent(i);
+ if(content==0.){empty++;}
+ else if(content!=0.){filled++;}
+ }//end for
+ AliInfo(Form(" %s : empty modules %i \t filled modules %i",hname.Data(), empty, filled));
+ }//end else pattern entries !=0
+ }//end if modpattern
+ else if(hname.Contains("SDDphizL3")||hname.Contains("SDDphizL4")){
+ Int_t layer=0;
+ if(hname.Contains("3"))layer=0;
+ else if(hname.Contains("4"))layer=1;
+ entries2[layer]=hdata->GetEntries();
+ if(entries2[layer]==0){
+ AliWarning(Form("===================>>>>>> No entries in %s \n",hname.Data()));
+ //printf("test = %f \t value %f\n",test,fStepBitSDD[AliQAv1::kFATAL]);
+ test=test+fStepBitSDD[AliQAv1::kFATAL];
+ if(entries==0){
+ //return test;
+ //break;
+ }
+ }//end if getentries
+ else{
+ Int_t layer=0;
+ if(hname.Contains("3"))layer=0;
+ else if(hname.Contains("4"))layer=1;
+ hlayer[layer]=dynamic_cast<TH2*>(hdata);
+ hlayer[layer]->RebinX(2);
+ int modmay=hlayer[layer]->GetNbinsY();
+ TH1D* hproj= hlayer[layer]->ProjectionY();
+ Double_t ladcontent=0;
+ for(Int_t i=1;i<=modmay;i++) {//loop on the ladders
+ ladcontent=hproj->GetBinContent(i);
+ if(ladcontent==0){emptyladders[layer]++;}
+ else if(ladcontent!=0){filledladders[layer]++;}
+ }//end for
+ AliInfo(Form(" %s : empty ladders %i \t filled ladders %i\n",hname.Data(), emptyladders[layer], filledladders[layer]));//end else layer 3
+ delete hproj;
+ hproj=NULL;
+ }//end else entries !=0
+ }//end if layer 3
+ }//end if hdata
+ }//end while
+ if(entries==0.&&entries2[0]==0.&&entries2[1]==0.) break;
else{
- AliError("Data type cannot be processed") ;
+ if(hmodule||(hlayer[0]&&hlayer[1])){
+ Int_t excluded=0;
+ Int_t active=0;
+ Int_t exactive=0;//excluded but taking data
+ //AliITSCalibrationSDD *cal;
+ for(Int_t imod=0;imod<fgknSDDmodules;imod++){
+ Int_t lay=0;
+ Int_t lad=0;
+ Int_t det=0;
+ Int_t module=0;
+ module=imod+fgkmodoffset;
+ AliITSCalibrationSDD * cal=(AliITSCalibrationSDD*)fCalibration->At(imod);
+ if(cal==0) { delete cal; continue;}
+ AliITSgeomTGeo::GetModuleId(module,lay,lad,det);
+ if (cal->IsBad()){
+ excluded++;
+ Double_t content=0.;
+ Double_t contentlayer[2];
+ for(Int_t i=0;i<2;i++)contentlayer[i]=0.;
+ if(hmodule)content=hmodule->GetBinContent(imod+1);//if expert bit is active the histogram has been created
+ contentlayer[lay-3]=hlayer[lay-3]->GetBinContent(det,lad);
+ if(content!=0.||contentlayer[lay-3]!=0.)
+ {
+ filledmodules[lay-3]++;
+ AliWarning(Form("The module %d (layer %i, ladder %i det %i ) excluded from the acquisition, took data \n ",module,lay,lad,det));
+ exactive++;
+ }
+ else if(content==0.&&contentlayer[lay-3]==0.)emptymodules[lay-3]++;
+ //AliInfo(Form("The module %d (layer %i, ladder %i det %i ) is bad, content %f content layer %f filled modules position %d ",module,lay,lad,det,contentlayer[lay-3],content,lay-3) );
+ }//end if bad
+ else
+ {
+ Double_t contentgood=0.;
+ active++;
+ //printf("lay: %i\t det %i \t lad %i \n",lay,det,lad );
+ contentgood=hlayer[lay-3]->GetBinContent(det,lad);
+ if(contentgood==0.){emptymodules[lay-3]++;}
+ else if(contentgood!=0.){filledmodules[lay-3]++;}
+ }
+
+ //delete cal;
+ //cal=NULL;
+ }//end for
+ for(Int_t i=0;i<2;i++){AliInfo(Form(" %s :Layer %i \tempty modules %i \t filled modules %i\n",hlayer[i]->GetName(), i+3,emptymodules[i], filledmodules[i]));}//end else layers
+ if(exactive==0){
+ AliInfo(Form("All the active modules (%i) are in acquisition. The number of excluded modules are %i \n",active,excluded));
+ test=fHighSDDValue[AliQAv1::kINFO];}
+ if(exactive!=0){
+ AliWarning(Form("%i modules excluded from the acquisition took data. Active modules%i \n ",exactive,active));
+ test=fHighSDDValue[AliQAv1::kWARNING];
+ }
+ if(excluded==exactive){
+ AliWarning(Form("All the modules exluded from the acquisition (%d) took data! Active modules %i\n",excluded,active));
+ test=fHighSDDValue[AliQAv1::kWARNING];
+ }
+ if(active==0){
+ AliWarning(Form("No modules took data: excluded %i \t exactive %i \n", excluded, exactive));
+ test=fHighSDDValue[AliQAv1::kFATAL];
+ }
+ }
}
- offset++;
- }
- } // GetEntries loop
- AliInfo(Form("Test Result for RECP = %f", test)) ;
- } // if(index==2) loop
-
+ }//end getentries !=0
+ //delete calSDD;
+
+
+
+ //delete calibSDD;
+ //delete calSDD;
+
+ break;
+ case AliQAv1::kNULLTASK:
+ AliInfo(Form("No Check on %s\n",AliQAv1::GetAliTaskName(index)));
+ test=1.;
+ break;
+ case AliQAv1::kSIM:
+ AliInfo(Form("Check on %s\n",AliQAv1::GetAliTaskName(index)));
+ test=1.;
+ break;
+ case AliQAv1::kREC:
+ AliInfo(Form("Check on %s\n",AliQAv1::GetAliTaskName(index)));
+ test=1.;
+ break;
+ case AliQAv1::kANA:
+ AliInfo(Form("Check on %s\n",AliQAv1::GetAliTaskName(index)));
+ test=1.;
+ break;
+ case AliQAv1::kESD:
+ AliInfo(Form("Check on %s\n",AliQAv1::GetAliTaskName(index)));
+ test=1.;
+ break;
+ case AliQAv1::kNTASK:
+ AliInfo(Form("No Check on %s\n",AliQAv1::GetAliTaskName(index)));
+ test=1.;
+ break;
+
+ }//end switch
+ delete hdata;
return test;
}
//__________________________________________________________________
-void AliITSQASDDChecker::SetTaskOffset(Int_t TaskOffset)
+void AliITSQASDDChecker::SetTaskOffset(Int_t taskoffset)
+{
+ fSubDetOffset = taskoffset;
+}
+
+
+//__________________________________________________________________
+void AliITSQASDDChecker::SetStepBit(Double_t *steprange)
+{
+
+ fStepBitSDD = new Double_t[AliQAv1::kNBIT];
+ for(Int_t bit=0;bit<AliQAv1::kNBIT;bit++)
+ {
+ fStepBitSDD[bit]=steprange[bit];
+ }
+}
+
+//__________________________________________________________________
+void AliITSQASDDChecker::SetSDDLimits(Float_t *lowvalue, Float_t * highvalue)
{
- fSubDetOffset = TaskOffset;
+
+ fLowSDDValue = new Float_t[AliQAv1::kNBIT];
+ fHighSDDValue= new Float_t[AliQAv1::kNBIT];
+
+ for(Int_t bit=0;bit<AliQAv1::kNBIT;bit++)
+ {
+ fLowSDDValue[bit]=lowvalue[bit];
+ fHighSDDValue[bit]= highvalue[bit];
+ }
+
}
+
+