]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSQADataMakerSim.cxx
Fixing a little overlap. Adding flags to perform material budget studies (Mario)
[u/mrichter/AliRoot.git] / ITS / AliITSQADataMakerSim.cxx
index 61f03f5eea1141a94a49437f7850c984592cfb20..a56d3bc29a697239072549118f661321c51f4ef5 100644 (file)
@@ -25,6 +25,7 @@
 
 // --- ROOT system ---
 #include <TTree.h>
+#include <TMath.h>
 // --- Standard library ---
 
 // --- AliRoot header files ---
@@ -33,7 +34,7 @@
 #include "AliITSQASDDDataMakerSim.h"
 #include "AliITSQASSDDataMakerSim.h"
 #include "AliLog.h"
-#include "AliQA.h"
+#include "AliQAv1.h"
 #include "AliQAChecker.h"
 #include "AliITSQAChecker.h"
 #include "AliRawReader.h"
@@ -42,7 +43,7 @@ ClassImp(AliITSQADataMakerSim)
 
 //____________________________________________________________________________ 
 AliITSQADataMakerSim::AliITSQADataMakerSim(Short_t subDet) :
-AliQADataMakerSim(AliQA::GetDetName(AliQA::kITS), "ITS Quality Assurance Data Maker"),
+AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kITS), "ITS Quality Assurance Data Maker"),
 fSubDetector(subDet),
 fSPDDataMaker(NULL),
 fSDDDataMaker(NULL),
@@ -53,18 +54,18 @@ fSSDDataMaker(NULL)
        AliError("Error: fSubDetector number out of range; return\n");
   }
 
-  // Initialization for RAW data 
+  // Initialization for each subdetector 
   if(fSubDetector == 0 || fSubDetector == 1) {
-    AliDebug(1,"AliITSQADM::Create SPD DataMakerSim\n");
+    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SPD DataMakerSim\n");
        fSPDDataMaker = new AliITSQASPDDataMakerSim(this);
   }
   if(fSubDetector == 0 || fSubDetector == 2) {
-       AliDebug(1,"AliITSQADM::Create SDD DataMakerSim\n");
+    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SDD DataMakerSim\n");
        //printf("AliITSQADM::Create SDD DataMakerSim\n");                  
        fSDDDataMaker = new AliITSQASDDDataMakerSim(this);
   }
   if(fSubDetector == 0 || fSubDetector == 3) {
-       AliDebug(1,"AliITSQADM::Create SSD DataMakerSim\n");
+    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SSD DataMakerSim\n");
        fSSDDataMaker = new AliITSQASSDDataMakerSim(this);
   }
 }
@@ -86,7 +87,6 @@ fSDDDataMaker(NULL),
 fSSDDataMaker(NULL)
 {
   //copy ctor 
-  //printf("AliITSQADataMakerSim::AliITSQADataMakerSim   using copy ctor \n");
   SetName((const char*)qadm.GetName()) ; 
   SetTitle((const char*)qadm.GetTitle());
 }
@@ -104,95 +104,150 @@ AliITSQADataMakerSim& AliITSQADataMakerSim::operator = (const AliITSQADataMakerS
 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();
-  //printf("AliITSQADataMakerSim::StartOfDetectorCycle() now launching fSDDDataMaker->StartOfDetectorCycle... \n");
   if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->StartOfDetectorCycle();
   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->StartOfDetectorCycle();
 }
 
 //____________________________________________________________________________ 
-void AliITSQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t 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"); 
-  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->EndOfDetectorCycle(task, list);  
-  //printf("AliITSQADataMakerSim::StartOfDetectorCycle() now launching fSDDDataMaker->EndOfDetectorCycle... \n");
-  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->EndOfDetectorCycle(task, list);
-  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->EndOfDetectorCycle(task, list);
-  
-  AliQAChecker *qac = AliQAChecker::Instance();
-  AliITSQAChecker *qacb = (AliITSQAChecker *) qac->GetDetQAChecker(0);
-  qac->Run( AliQA::kITS , task, list);  //temporary skipping the checking
+  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 RAW data 
-       if(fSubDetector == 0 || fSubDetector == 1) {
-         AliDebug(1,"AliITSQADM:: SPD InitDigits\n");
-         fSPDDataMaker->InitDigits();
-       }
-       if(fSubDetector == 0 || fSubDetector == 2) {
-         AliDebug(1,"AliITSQADM:: SDD InitDigits\n");
-//       printf("AliITSQADataMakerSim::InitDigits()   launching AliITSQADM:: SDD InitDigits\n");
-         fSDDDataMaker->InitDigits();
-       }
-       if(fSubDetector == 0 || fSubDetector == 3) {
-         AliDebug(1,"AliITSQADM:: SSD InitDigits\n");
-         fSSDDataMaker->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(TClonesArray * digits)
+void AliITSQADataMakerSim::MakeDigits()
 { 
-  // Fill QA for RAW   
-  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeDigits(digits);
-  //printf("AliITSQADataMakerSim::MakeDigits()   launching fSDDDataMaker->MakeDigits with TClonesArray\n");
-  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeDigits(digits);
-  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
+  // 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 RAW   
-  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeDigits(digits);
-  //printf("AliITSQADataMakerSim::MakeDigits()   launching fSDDDataMaker->MakeDigits with TTree\n");
-  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeDigits(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 RECPOINTS
+  // Initialization for SDigits
+  fSDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(70);
   if(fSubDetector == 0 || fSubDetector == 1) {
-       AliDebug(1,"AliITSQADM:: SPD InitSDigits\n");
+    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitSDigits\n");
+
     fSPDDataMaker->InitSDigits();
   }
-  if(fSubDetector == 0 || fSubDetector == 2) {
-       AliDebug(1,"AliITSQADM:: SDD InitSDigits\n");
-//     printf("AliITSQADataMakerSim::InitSDigits()   launching AliITSQADM:: SDD InitSDigits\n");
-       fSDDDataMaker->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(1,"AliITSQADM:: SSD InitSDigits\n");
-       fSSDDataMaker->InitSDigits();
+    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(TClonesArray * sdigits)
+void AliITSQADataMakerSim::MakeSDigits()
 {
-  // Fill QA for recpoints
-  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeSDigits(sdigits); 
-  //printf("AliITSQADataMakerSim::MakeSDigits()   launching fSDDDataMaker->MakeSDigits with TClonesArray\n");
-  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeSDigits(sdigits);
-  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
+  // 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();
+  //
 }
 
 //____________________________________________________________________________ 
@@ -200,46 +255,119 @@ void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
 {
   // Fill QA for recpoints
   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeSDigits(sdigits); 
-  //printf("AliITSQADataMakerSim::MakeSDigits()   launching fSDDDataMaker->MakeSDigits with TTree\n");
-  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->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 RECPOINTS
+  // Initialization for hits
+  fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(50);
   if(fSubDetector == 0 || fSubDetector == 1) {
-       AliDebug(1,"AliITSQADM:: SPD InitHits\n");
+    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitHits\n");
     fSPDDataMaker->InitHits();
   }
   if(fSubDetector == 0 || fSubDetector == 2) {
-       AliDebug(1,"AliITSQADM:: SDD InitHits\n");
-//     printf("AliITSQADataMakerSim::InitHits()   launching AliITSQADM:: SDD InitHits\n");
-       fSDDDataMaker->InitHits();
+    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(1,"AliITSQADM:: SSD InitHits\n");
-       fSSDDataMaker->InitHits();
+    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(TClonesArray * hits)
+void AliITSQADataMakerSim::MakeHits()
 {
-  // Fill QA for recpoints
-  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeHits(hits);
-  //printf("AliITSQADataMakerSim::MakeHits()   launching fSDDDataMaker->MakeHits with TClonesArray\n");
-  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeHits(hits);
-  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
+  // 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 recpoints
+  // Fill QA for hits
   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeHits(hits);
-  //printf("AliITSQADataMakerSim::MakeHits()   launching fSDDDataMaker->MakeHits with TTree\n");
   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;
 }