]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSQADataMakerSim.cxx
Added fit macro from M. Putis
[u/mrichter/AliRoot.git] / ITS / AliITSQADataMakerSim.cxx
index c6b3784386b1ac5d8cfd70c9b025829a8155bd8e..faef0a88f9ed471b2ca50017cdf43d02a1639a14 100644 (file)
@@ -54,7 +54,7 @@ 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(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SPD DataMakerSim\n");
        fSPDDataMaker = new AliITSQASPDDataMakerSim(this);
@@ -114,83 +114,99 @@ void AliITSQADataMakerSim::StartOfDetectorCycle()
 //____________________________________________________________________________ 
 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
   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
-    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),fSDDDataMaker->GetOffset(task),fSSDDataMaker->GetOffset(task)); //Setting the offset for the QAChecker list           
+    if (! IsValidEventSpecie(specie, list) ){
+      continue ; 
     }
-    else
-      if(subdet!=0){
-       Int_t offset=GetDetTaskOffset(subdet, task);
-       qacb->SetDetTaskOffset(subdet,offset);
-      }
-    qac->Run( AliQAv1::kITS , task, list); 
-  }
+    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()
 {  
-  // Initialization for RAW data 
+
+  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(TClonesArray * digits)
+void AliITSQADataMakerSim::MakeDigits()
 { 
   // Fill QA for digits   
-  if(fSubDetector == 0 || fSubDetector == 1) {
-    Int_t rv = fSPDDataMaker->MakeDigits(digits) ; 
-    if ( rv != 0 )
-      fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
-  }
+  if(fSubDetector == 0 || fSubDetector == 1) 
+    fSPDDataMaker->MakeDigits() ; 
+
   
-  if(fSubDetector == 0 || fSubDetector == 2) {
-    Int_t rv = fSDDDataMaker->MakeDigits(digits) ; 
-    if ( rv != 0 )
-      fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
-  }
+  if(fSubDetector == 0 || fSubDetector == 2) 
+    fSDDDataMaker->MakeDigits() ; 
 
-  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
+  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits();
 }
 
 //____________________________________________________________________________
 void AliITSQADataMakerSim::MakeDigits(TTree * digits)
 { 
   // Fill QA for digits   
-  if(fSubDetector == 0 || fSubDetector == 1) {
-    Int_t rv = fSPDDataMaker->MakeDigits(digits) ; 
-    if ( rv != 0 )
-      fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
-  }
-  
-  if(fSubDetector == 0 || fSubDetector == 2){
-    Int_t rv = fSDDDataMaker->MakeDigits(digits) ; 
-    if ( rv != 0 )
-      fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
-  }
+  if(fSubDetector == 0 || fSubDetector == 1) 
+    fSPDDataMaker->MakeDigits(digits) ; 
+
+  if(fSubDetector == 0 || fSubDetector == 2)
+    fSDDDataMaker->MakeDigits(digits) ; 
 
   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
 }
@@ -198,38 +214,40 @@ void AliITSQADataMakerSim::MakeDigits(TTree * 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) {
+  if(fSubDetector == 0 || fSubDetector == 2)
     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitSDigits\n");
-       fSDDDataMaker->InitSDigits();
+
+    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->InitSDigits();
+
+    fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
+    fSSDDataMaker->InitSDigits();
   }
 }
 
 //____________________________________________________________________________ 
-void AliITSQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
+void AliITSQADataMakerSim::MakeSDigits()
 {
   // Fill QA for recpoints
-  if(fSubDetector == 0 || fSubDetector == 1){
-    Int_t rv = fSPDDataMaker->MakeSDigits(sdigits) ; 
-    if ( rv != 0 )
-      fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
-  }
+  if(fSubDetector == 0 || fSubDetector == 1)
+    fSPDDataMaker->MakeSDigits() ; 
+
   
-  if(fSubDetector == 0 || fSubDetector == 2) {
-    Int_t rv = fSDDDataMaker->MakeSDigits(sdigits) ; 
-    if ( rv != 0 )
-      fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
-  }
+  if(fSubDetector == 0 || fSubDetector == 2) 
+    fSDDDataMaker->MakeSDigits() ; 
 
-  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
+
+  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits();
 }
 
 //____________________________________________________________________________ 
@@ -237,15 +255,11 @@ void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
 {
   // Fill QA for recpoints
   if(fSubDetector == 0 || fSubDetector == 1){
-    Int_t rv = fSPDDataMaker->MakeSDigits(sdigits) ; 
-    if ( rv != 0 )
-      fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
-  }
+    fSPDDataMaker->MakeSDigits(sdigits) ; 
+ }
   
   if(fSubDetector == 0 || fSubDetector == 2){
-    Int_t rv = fSDDDataMaker->MakeSDigits(sdigits) ; 
-    if ( rv != 0 )
-      fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
+    fSDDDataMaker->MakeSDigits(sdigits) ; 
   }
 
   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
@@ -254,6 +268,7 @@ void AliITSQADataMakerSim::MakeSDigits(TTree * 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");
@@ -261,31 +276,31 @@ void AliITSQADataMakerSim::InitHits()
   }
   if(fSubDetector == 0 || fSubDetector == 2) {
     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitHits\n");
-       fSDDDataMaker->InitHits();
+
+    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->InitHits();
+
+    fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
+    fSSDDataMaker->InitHits();
   }
 }
 
 //____________________________________________________________________________ 
-void AliITSQADataMakerSim::MakeHits(TClonesArray * hits)
+void AliITSQADataMakerSim::MakeHits()
 {
   // Fill QA for hits
   if(fSubDetector == 0 || fSubDetector == 1) {
-    Int_t rv = fSPDDataMaker->MakeHits(hits) ; 
-    if ( rv != 0 )
-      fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
-  }
+    fSPDDataMaker->MakeHits() ; 
+    }
   
   if(fSubDetector == 0 || fSubDetector == 2) {
-    Int_t rv = fSDDDataMaker->MakeHits(hits) ; 
-    if ( rv != 0 )
-      fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
+    fSDDDataMaker->MakeHits() ; 
   }
 
-  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
+  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits();
 }
 
 //____________________________________________________________________________ 
@@ -293,14 +308,10 @@ void AliITSQADataMakerSim::MakeHits(TTree * hits)
 {
   // Fill QA for hits
   if(fSubDetector == 0 || fSubDetector == 1) {
-    Int_t rv = fSPDDataMaker->MakeHits(hits) ; 
-    if ( rv != 0 )
-      fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
-  }
+    fSPDDataMaker->MakeHits(hits) ; 
+   }
   if(fSubDetector == 0 || fSubDetector == 2) {
-    Int_t rv = fSDDDataMaker->MakeHits(hits) ; 
-    if ( rv != 0 )
-      fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
+    fSDDDataMaker->MakeHits(hits) ; 
   }
 
   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
@@ -333,3 +344,31 @@ Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t t
     }
   //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;
+}