Introducing event specie in QA (Yves)

[u/mrichter/AliRoot.git] / FMD / AliFMDQADataMakerSim.cxx

diff --git a/FMD/AliFMDQADataMakerSim.cxx b/FMD/AliFMDQADataMakerSim.cxx
index eb876bd79f9976c6418bcbb7402abd62ca146988..df3f81cd4c04866c9f61648298ebd932e053ebc4 100644 (file)
--- a/FMD/AliFMDQADataMakerSim.cxx
+++ b/FMD/AliFMDQADataMakerSim.cxx
@@ -27,6 +27,7 @@
 #include "AliFMDHit.h"
 #include "AliQAChecker.h"
 #include "AliFMDParameters.h"
+#include "AliFMDSDigit.h"
 
 //_____________________________________________________________________
 // This is the class that collects the QA data for the FMD during simulation.
@@ -47,33 +48,45 @@ ClassImp(AliFMDQADataMakerSim)
 AliFMDQADataMakerSim::AliFMDQADataMakerSim() 
   :  AliQADataMakerSim(AliQA::GetDetName(AliQA::kFMD),
                       "FMD Quality Assurance Data Maker"),
-     fDigitsArray(0),
-     fHitsArray(0)
+     fSDigitsArray("AliFMDSDigit", 1000),
+     fDigitsArray("AliFMDDigit", 1000),
+     fHitsArray("AliFMDHit", 10)
 {
   // ctor
-  fDigitsArray = new TClonesArray("AliFMDDigit", 1000) ; 
-  fHitsArray   = new TClonesArray("AliFMDHit", 10);
+
 }
 
 //_____________________________________________________________________
-AliFMDQADataMakerSim::AliFMDQADataMakerSim(const AliFMDQADataMakerSim& /*qadm*/) 
-  : AliQADataMakerSim()
+AliFMDQADataMakerSim::AliFMDQADataMakerSim(const AliFMDQADataMakerSim& qadm) 
+  : AliQADataMakerSim(),
+    fSDigitsArray(qadm.fSDigitsArray),
+    fDigitsArray(qadm.fDigitsArray),
+    fHitsArray(qadm.fHitsArray)
 {
   //copy ctor 
+  
   // Parameters: 
   //    qadm    Object to copy from
   
 }
 //_____________________________________________________________________
+AliFMDQADataMakerSim& AliFMDQADataMakerSim::operator = (const AliFMDQADataMakerSim& qadm ) 
+{
+  fSDigitsArray = qadm.fSDigitsArray;
+  fDigitsArray = qadm.fDigitsArray;
+  fHitsArray = qadm.fHitsArray;
+  
+  return *this;
+}
+//_____________________________________________________________________
 AliFMDQADataMakerSim::~AliFMDQADataMakerSim()
 {
-  delete fDigitsArray;
-  delete fHitsArray;
+
 }
 
 //_____________________________________________________________________
 void AliFMDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, 
-                                             TObjArray * list)
+                                             TObjArray ** list)
 {
   //Detector specific actions at end of cycle
   // do the QA checking
@@ -84,6 +97,14 @@ void AliFMDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task,
   AliQAChecker::Instance()->Run(AliQA::kFMD, task, list) ;  
   
 }
+//_____________________________________________________________________
+void AliFMDQADataMakerSim::InitSDigits()
+{
+  // create SDigits histograms in SDigits subdir
+  TH1I* hADCCounts = new TH1I("hADCCounts","Dist of ADC counts",1024,0,1024);
+  hADCCounts->SetXTitle("ADC counts");
+  Add2SDigitsList(hADCCounts, 0);
+}
 
 //____________________________________________________________________ 
 void AliFMDQADataMakerSim::InitHits()
@@ -107,7 +128,7 @@ void AliFMDQADataMakerSim::InitDigits()
 //_____________________________________________________________________
 void AliFMDQADataMakerSim::MakeHits(TClonesArray * hits)
 {
-  TIter next(fHitsArray);
+  TIter next(hits);
   AliFMDHit * hit;
   while ((hit = static_cast<AliFMDHit *>(next()))) 
     GetHitsData(0)->Fill(hit->Edep()/hit->Length()*0.032);
@@ -118,17 +139,21 @@ void AliFMDQADataMakerSim::MakeHits(TTree * hitTree)
 {
   // make QA data from Hit Tree
   
+  fHitsArray.Clear();
+  
   TBranch * branch = hitTree->GetBranch("FMD") ;
   if (!branch) {
     AliWarning("FMD branch in Hit Tree not found") ; 
     return;
   }
-  fHitsArray->Clear();
-  branch->SetAddress(&fHitsArray) ;
+  
+  TClonesArray* hitsAddress = &fHitsArray;
+  
+  branch->SetAddress(&hitsAddress) ;
   
   for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
     branch->GetEntry(ientry);
-    MakeHits(fHitsArray);   //tmp); 
+    MakeHits(hitsAddress);   //tmp); 
   }    
   
 }
@@ -139,7 +164,7 @@ void AliFMDQADataMakerSim::MakeDigits(TClonesArray * digits)
   // makes data from Digits
   if(!digits) return;
 
-  for(Int_t i = 0 ; i < fDigitsArray->GetEntriesFast() ; i++) {
+  for(Int_t i = 0 ; i < digits->GetEntriesFast() ; i++) {
     //Raw ADC counts
     AliFMDDigit* digit = static_cast<AliFMDDigit*>(digits->At(i));
     GetDigitsData(0)->Fill(digit->Counts());
@@ -150,15 +175,45 @@ void AliFMDQADataMakerSim::MakeDigits(TClonesArray * digits)
 void AliFMDQADataMakerSim::MakeDigits(TTree * digitTree)
 {
   
-  fDigitsArray->Clear();
+  fDigitsArray.Clear();
   TBranch * branch = digitTree->GetBranch("FMD") ;
   if (!branch)    {
       AliWarning("FMD branch in Digit Tree not found") ; 
       return;
   } 
-  branch->SetAddress(&fDigitsArray) ;
+  TClonesArray* digitAddress = &fDigitsArray;
+  branch->SetAddress(&digitAddress) ;
+  branch->GetEntry(0) ; 
+  MakeDigits(digitAddress) ; 
+}
+
+//_____________________________________________________________________
+void AliFMDQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
+{
+  // makes data from Digits
+  if(!sdigits) return;
+  
+  for(Int_t i = 0 ; i < sdigits->GetEntriesFast() ; i++) {
+    //Raw ADC counts
+    AliFMDSDigit* sdigit = static_cast<AliFMDSDigit*>(sdigits->At(i));
+    GetSDigitsData(0)->Fill(sdigit->Counts());
+  }
+}
+
+//_____________________________________________________________________
+void AliFMDQADataMakerSim::MakeSDigits(TTree * sdigitTree)
+{
+  
+  fSDigitsArray.Clear();
+  TBranch * branch = sdigitTree->GetBranch("FMD") ;
+  if (!branch)    {
+    AliWarning("FMD branch in SDigit Tree not found") ; 
+    return;
+  } 
+  TClonesArray* sdigitAddress = &fSDigitsArray;
+  branch->SetAddress(&sdigitAddress) ;
   branch->GetEntry(0) ; 
-  MakeDigits(fDigitsArray) ; 
+  MakeSDigits(sdigitAddress) ; 
 }
 
 //_____________________________________________________________________