// 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() const
+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");
+ ResetEventTrigClasses(); // reset triggers list to select all histos
+ ResetEvCountCycle();
+ //
+ 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() ));
+ ResetEventTrigClasses(); // reset triggers list to select all histos
+ //
// launch the QA checking
- AliDebug(1,"AliITSDM instantiates checker with Run(AliQA::kITS, task, list)\n");
-
- AliQAChecker::Instance()->Run( AliQA::kITS , task, list);
+ for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
+ if (! IsValidEventSpecie(specie, list) ) continue;
+ SetEventSpecie(AliRecoParam::ConvertIndex(specie));
+ //
+ 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 for
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::InitDigits()
+{
+
+ // Initialization for Digits data
+ fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(60);
+ 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();
+ }
+ //
+ ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::MakeDigits()
+{
+ // Fill QA for digits
+ if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeDigits();
+ if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeDigits();
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits();
+ //
+ IncEvCountCycleDigits();
+ IncEvCountTotalDigits();
+ //
+}
+
+//____________________________________________________________________________
+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);
+ //
+ IncEvCountCycleDigits();
+ IncEvCountTotalDigits();
+ //
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::InitSDigits()
+{
+ // Initialization for SDigits
+ fSDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(70);
+ 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();
+ }
+ //
+ ClonePerTrigClass(AliQAv1::kSDIGITS); // this should be the last line
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::MakeSDigits()
+{
+ // Fill QA for sdigits
+ if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeSDigits();
+ if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeSDigits();
+ if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits();
+ //
+ IncEvCountCycleSDigits();
+ IncEvCountTotalSDigits();
+ //
+}
+
+//____________________________________________________________________________
+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);
+ //
+ IncEvCountCycleSDigits();
+ IncEvCountTotalSDigits();
+ //
+}
+
+//____________________________________________________________________________
+void AliITSQADataMakerSim::InitHits()
+{
+ // Initialization for hits
+ fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(50);
+ 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();
+ }
+ //
+ ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
}
+//____________________________________________________________________________
+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();
+ //
+ IncEvCountCycleHits();
+ IncEvCountTotalHits();
+ //
+}
+
+//____________________________________________________________________________
+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);
+ //
+ IncEvCountCycleHits();
+ IncEvCountTotalHits();
+ //
+}
+
+//_________________________________________________________________
+Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task)
+{
+
+ //return the offset for each subdetector
+ 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)
+{
+ //return of the number of histograms for each task and for each sub detector
+ 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;
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff