// by comparing with reference data
// contained in a DB
// -------------------------------------------------------------
-// W. Ferrarese Nov 2007
+// W. Ferrarese + P. Cerello Feb 2008
// INFN Torino
// --- ROOT system ---
-#include <TH2D.h>
#include <TTree.h>
-
+#include <TMath.h>
// --- Standard library ---
// --- AliRoot header files ---
#include "AliITSQADataMakerSim.h"
+#include "AliITSQASPDDataMakerSim.h"
+#include "AliITSQASDDDataMakerSim.h"
+#include "AliITSQASSDDataMakerSim.h"
#include "AliLog.h"
+#include "AliQAv1.h"
#include "AliQAChecker.h"
-
+#include "AliITSQAChecker.h"
+#include "AliRawReader.h"
ClassImp(AliITSQADataMakerSim)
//____________________________________________________________________________
-AliITSQADataMakerSim::AliITSQADataMakerSim() :
- AliQADataMakerSim(AliQA::GetDetName(AliQA::kITS), "SDD Quality Assurance Data Maker")
-{
- // ctor
+AliITSQADataMakerSim::AliITSQADataMakerSim(Short_t subDet) :
+AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kITS), "ITS Quality Assurance Data Maker"),
+fSubDetector(subDet),
+fSPDDataMaker(NULL),
+fSDDDataMaker(NULL),
+fSSDDataMaker(NULL)
+{
+ //ctor used to discriminate OnLine-Offline analysis
+ if(fSubDetector < 0 || fSubDetector > 3) {
+ AliError("Error: fSubDetector number out of range; return\n");
+ }
+
+ // Initialization for each subdetector
+ if(fSubDetector == 0 || fSubDetector == 1) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SPD DataMakerSim\n");
+ fSPDDataMaker = new AliITSQASPDDataMakerSim(this);
+ }
+ if(fSubDetector == 0 || fSubDetector == 2) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SDD DataMakerSim\n");
+ //printf("AliITSQADM::Create SDD DataMakerSim\n");
+ fSDDDataMaker = new AliITSQASDDDataMakerSim(this);
+ }
+ if(fSubDetector == 0 || fSubDetector == 3) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SSD DataMakerSim\n");
+ fSSDDataMaker = new AliITSQASSDDataMakerSim(this);
+ }
}
//____________________________________________________________________________
-AliITSQADataMakerSim::AliITSQADataMakerSim(Int_t /*ldc */, Bool_t /*kMode */) :
- AliQADataMakerSim(AliQA::GetDetName(AliQA::kITS), "SDD Quality Assurance Data Maker")
-{
- //ctor used to discriminate OnLine-Offline analysis
+AliITSQADataMakerSim::~AliITSQADataMakerSim(){
+ // destructor
+ if(fSPDDataMaker)delete fSPDDataMaker;
+ if(fSDDDataMaker)delete fSDDDataMaker;
+ if(fSSDDataMaker)delete fSSDDataMaker;
}
//____________________________________________________________________________
AliITSQADataMakerSim::AliITSQADataMakerSim(const AliITSQADataMakerSim& qadm) :
- AliQADataMakerSim()
+AliQADataMakerSim(),
+fSubDetector(qadm.fSubDetector),
+fSPDDataMaker(NULL),
+fSDDDataMaker(NULL),
+fSSDDataMaker(NULL)
{
//copy ctor
SetName((const char*)qadm.GetName()) ;
- SetTitle((const char*)qadm.GetTitle());
+ SetTitle((const char*)qadm.GetTitle());
}
//__________________________________________________________________
void AliITSQADataMakerSim::StartOfDetectorCycle()
{
//Detector specific actions at start of cycle
- AliDebug(1,"AliITSQADM::Start of ITS Cycle\n");
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Start of ITS Cycle\n");
+
+ if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->StartOfDetectorCycle();
+ if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->StartOfDetectorCycle();
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->StartOfDetectorCycle();
}
//____________________________________________________________________________
-void AliITSQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX task, TObjArray *list)
+void AliITSQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
{
+
+ AliInfo(Form("End of Dedetctor Cycle called for %s\n",AliQAv1::GetTaskName(task).Data() ));
+
// launch the QA checking
- AliDebug(1,"AliITSDM instantiates checker with Run(AliQA::kITS, task, list)\n");
+ for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
+ if (! IsValidEventSpecie(specie, list) ){
+ continue ;
+ }
+ else{
+ Int_t idnumber=list[specie]->GetUniqueID();
+ //printf("specie %s \t id number == %d\n",AliRecoParam::GetEventSpecieName(specie),idnumber);
+ if(idnumber==0)
+ {
+ //AliInfo(Form("No check for %s\n",AliQAv1::GetTaskName(task).Data() ))
+ continue;
+ } //skip kDigitsR and not filled TobjArray specie
+ else{
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
+ if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->EndOfDetectorCycle(task, list[specie]);
+ if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->EndOfDetectorCycle(task, list[specie]);
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->EndOfDetectorCycle(task, list[specie]);
+
+
+ AliQAChecker *qac = AliQAChecker::Instance();
+ AliITSQAChecker *qacb = (AliITSQAChecker *) qac->GetDetQAChecker(0);
+ Int_t subdet=GetSubDet();
+ qacb->SetSubDet(subdet);
+
+ if(subdet== 0 ){
+ qacb->SetTaskOffset(fSPDDataMaker->GetOffset(task,specie),fSDDDataMaker->GetOffset(task,specie),fSSDDataMaker->GetOffset(task,specie)); //Setting the offset for the QAChecker list
+ qacb->SetHisto(fSPDDataMaker->GetTaskHisto(task), fSDDDataMaker->GetTaskHisto(task), fSSDDataMaker->GetTaskHisto(task));
+ }
+ else
+ if(subdet!=0){
+ Int_t offset=GetDetTaskOffset(subdet, task);
+ qacb->SetDetTaskOffset(subdet,offset);
+ Int_t histo=GetDetTaskHisto(subdet, task);
+ qacb->SetDetHisto(subdet,histo);
+ }
+ qac->Run( AliQAv1::kITS , task, list);
+ }//end else unique id
+ }//end else event specie
+ }//end for
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::InitDigits()
+{
+
+ fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(60);
+ // Initialization for Digits data
+ if(fSubDetector == 0 || fSubDetector == 1) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitDigits\n");
+
+ fSPDDataMaker->InitDigits();
+ }
+ if(fSubDetector == 0 || fSubDetector == 2) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitDigits\n");
+
+ fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
+ fSDDDataMaker->InitDigits();
+ }
+ if(fSubDetector == 0 || fSubDetector == 3) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitDigits\n");
+
+ fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
+ fSSDDataMaker->InitDigits();
+ }
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::MakeDigits()
+{
+ // Fill QA for digits
+ if(fSubDetector == 0 || fSubDetector == 1)
+ fSPDDataMaker->MakeDigits() ;
+
- AliQAChecker::Instance()->Run( AliQA::kITS , task, list);
+ if(fSubDetector == 0 || fSubDetector == 2)
+ fSDDDataMaker->MakeDigits() ;
+
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits();
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::MakeDigits(TTree * digits)
+{
+ // Fill QA for digits
+ if(fSubDetector == 0 || fSubDetector == 1)
+ fSPDDataMaker->MakeDigits(digits) ;
+
+ if(fSubDetector == 0 || fSubDetector == 2)
+ fSDDDataMaker->MakeDigits(digits) ;
+
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
}
+//____________________________________________________________________________
+void AliITSQADataMakerSim::InitSDigits()
+{
+ fSDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(70);
+ // Initialization for RECPOINTS
+ if(fSubDetector == 0 || fSubDetector == 1) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitSDigits\n");
+
+ fSPDDataMaker->InitSDigits();
+ }
+ if(fSubDetector == 0 || fSubDetector == 2){
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitSDigits\n");
+
+ fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
+ fSDDDataMaker->InitSDigits();
+ }
+ if(fSubDetector == 0 || fSubDetector == 3) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitSDigits\n");
+
+ fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
+ fSSDDataMaker->InitSDigits();
+ }
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::MakeSDigits()
+{
+ // Fill QA for recpoints
+ if(fSubDetector == 0 || fSubDetector == 1)
+ fSPDDataMaker->MakeSDigits() ;
+
+
+ if(fSubDetector == 0 || fSubDetector == 2)
+ fSDDDataMaker->MakeSDigits() ;
+
+
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits();
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
+{
+ // Fill QA for recpoints
+ if(fSubDetector == 0 || fSubDetector == 1){
+ fSPDDataMaker->MakeSDigits(sdigits) ;
+ }
+
+ if(fSubDetector == 0 || fSubDetector == 2){
+ fSDDDataMaker->MakeSDigits(sdigits) ;
+ }
+
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::InitHits()
+{
+ fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(50);
+ // Initialization for hits
+ if(fSubDetector == 0 || fSubDetector == 1) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitHits\n");
+ fSPDDataMaker->InitHits();
+ }
+ if(fSubDetector == 0 || fSubDetector == 2) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitHits\n");
+
+ fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
+ fSDDDataMaker->InitHits();
+ }
+ if(fSubDetector == 0 || fSubDetector == 3) {
+ AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitHits\n");
+
+ fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
+ fSSDDataMaker->InitHits();
+ }
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::MakeHits()
+{
+ // Fill QA for hits
+ if(fSubDetector == 0 || fSubDetector == 1) {
+ fSPDDataMaker->MakeHits() ;
+ }
+
+ if(fSubDetector == 0 || fSubDetector == 2) {
+ fSDDDataMaker->MakeHits() ;
+ }
+
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits();
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::MakeHits(TTree * hits)
+{
+ // Fill QA for hits
+ if(fSubDetector == 0 || fSubDetector == 1) {
+ fSPDDataMaker->MakeHits(hits) ;
+ }
+ if(fSubDetector == 0 || fSubDetector == 2) {
+ fSDDDataMaker->MakeHits(hits) ;
+ }
+
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
+}
+
+//_________________________________________________________________
+Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task)
+{
+ switch(subdet)
+ {
+
+ Int_t offset;
+ case 1:
+ offset=fSPDDataMaker->GetOffset(task);
+ return offset;
+ break;
+ case 2:
+ offset=fSDDDataMaker->GetOffset(task);
+ return offset;
+ break;
+ case 3:
+ offset=fSSDDataMaker->GetOffset(task);
+ return offset;
+ break;
+ default:
+ AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
+ offset=0;
+ return offset;
+ break;
+ }
+ //return offset;
+}
+
+//_________________________________________________________________
+Int_t AliITSQADataMakerSim::GetDetTaskHisto(Int_t subdet,AliQAv1::TASKINDEX_t task)
+{
+ switch(subdet)
+ {
+
+ Int_t histo;
+ case 1:
+ histo=fSPDDataMaker->GetOffset(task);
+ return histo;
+ break;
+ case 2:
+ histo=fSDDDataMaker->GetOffset(task);
+ return histo;
+ break;
+ case 3:
+ histo=fSSDDataMaker->GetOffset(task);
+ return histo;
+ break;
+ default:
+ AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
+ histo=0;
+ return histo;
+ break;
+ }
+ //return offset;
+}