ATO-98 make macro to calculate correction derivative again working

[u/mrichter/AliRoot.git] / HMPID / AliHMPIDQADataMakerSim.cxx

diff --git a/HMPID/AliHMPIDQADataMakerSim.cxx b/HMPID/AliHMPIDQADataMakerSim.cxx
index 84a6d1221b5d2979c7a8cda73e4edf71558a638b..677f7fbe0ee62ad055d3ea897b56273b671e26c6 100644 (file)
--- a/HMPID/AliHMPIDQADataMakerSim.cxx
+++ b/HMPID/AliHMPIDQADataMakerSim.cxx
@@ -47,14 +47,14 @@ ClassImp(AliHMPIDQADataMakerSim)
            
 //____________________________________________________________________________ 
   AliHMPIDQADataMakerSim::AliHMPIDQADataMakerSim() : 
-  AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kHMPID), "HMPID Quality Assurance Data Maker")
+  AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kHMPID), "HMPID Quality Assurance Data Maker"), fChannel(0)
 {
   // ctor
 }
 
 //____________________________________________________________________________ 
 AliHMPIDQADataMakerSim::AliHMPIDQADataMakerSim(const AliHMPIDQADataMakerSim& qadm) :
-  AliQADataMakerSim() 
+  AliQADataMakerSim(),fChannel(0)
 {
   //copy ctor 
   SetName((const char*)qadm.GetName()) ; 
@@ -84,6 +84,8 @@ void AliHMPIDQADataMakerSim::InitHits()
     hHitMap[iCh]=new TH2F(Form("HMPID HitMap%i",iCh),Form("Ch%i;x_{Hit};y_{Hit};Entries",iCh),162,-1,161,146,-1,145);   
     Add2HitsList(hHitMap[iCh],iCh+1,expert,!image);
   }
+  //
+  ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
 }
 
 //____________________________________________________________________________ 
@@ -108,6 +110,8 @@ void AliHMPIDQADataMakerSim::InitDigits()
   Add2DigitsList(hDigPcEvt,1,expert, !image);
   for(Int_t iMap=0; iMap < 7; iMap++) Add2DigitsList(hDigMap[iMap],2+iMap,expert, !image);
   for(Int_t iH =0; iH < 42 ; iH++) Add2DigitsList(hDigQ[iH]    ,9+iH,expert,!image);
+  //
+  ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
 }
 
 //____________________________________________________________________________ 
@@ -119,103 +123,93 @@ void AliHMPIDQADataMakerSim::InitSDigits()
 
   TH1F   *hSDigits     = new TH1F("hHmpidSDigits",    "SDigits Q  distribution in HMPID;QDC;Entries",  500, 0., 5000.) ; 
   Add2SDigitsList(hSDigits,0, !expert, image);
+  //
+  ClonePerTrigClass(AliQAv1::kSDIGITS); // this should be the last line}
 }
 
 //____________________________________________________________________________ 
 
-void AliHMPIDQADataMakerSim::MakeHits(TClonesArray * data)
+void AliHMPIDQADataMakerSim::MakeHits()
 {
  //
  //filling QA histos for Hits
  //
-  TClonesArray * hits = dynamic_cast<TClonesArray *>(data) ; 
-  if (!hits){
-    AliError("Wrong type of hits container") ; 
-  } else {
-    TIter next(hits); 
-    AliHMPIDHit * hit ; 
-    while ( (hit = dynamic_cast<AliHMPIDHit *>(next())) ) {
-      if(hit->Pid()<500000) GetHitsData(0)->Fill(hit->Q()) ;
-      if(hit->Pid()<500000) GetHitsData(hit->Ch()+1)->Fill(hit->LorsX(),hit->LorsY());
-    }
-  } 
+  
+  TIter next(fHitsArray); 
+  AliHMPIDHit * hit ; 
+  while ( (hit = dynamic_cast<AliHMPIDHit *>(next())) ) {
+    if(hit->Pid()<500000) FillHitsData(0,hit->Q()) ;
+    if(hit->Pid()<500000) FillHitsData(hit->Ch()+1,hit->LorsX(),hit->LorsY());
+  }
+} 
 
-}
 //___________________________________________________________________________
 void AliHMPIDQADataMakerSim::MakeHits(TTree * data)
 {
 //
 //Opening of the Hit TTree 
 //
- TClonesArray *pHits=new TClonesArray("AliHMPIDHit");  data->SetBranchAddress("HMPID",&pHits);
+ if (fHitsArray) 
+   fHitsArray->Clear() ; 
+  else 
+    fHitsArray=new TClonesArray("AliHMPIDHit");  
+  data->SetBranchAddress("HMPID",&fHitsArray);
   for(Int_t iEnt=0;iEnt<data->GetEntriesFast();iEnt++){//entries loop
     data->GetEntry(iEnt);
-    MakeHits(pHits);
+    MakeHits();
   }//entries loop
 }
-//____________________________________________________________________________
-void AliHMPIDQADataMakerSim::MakeDigits(TClonesArray * data)
+//___________________________________________________________________________
+void AliHMPIDQADataMakerSim::MakeDigits()
 {
- //
- //filling QA histos for Digits
- //
-
-  TObjArray *chamber = dynamic_cast<TObjArray*>(data);
-  if ( !chamber) {
-    AliError("Wrong type of digits container") ; 
-  } else {
-    for(Int_t i =0; i< chamber->GetEntries(); i++)
-      {
-       TClonesArray * digits = dynamic_cast<TClonesArray*>(chamber->At(i)); 
-       GetDigitsData(0)->Fill(i,digits->GetEntriesFast()/(48.*80.*6.));
-       TIter next(digits); 
-       AliHMPIDDigit * digit; 
-       while ( (digit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
-         GetDigitsData(1)->Fill(10.*i+digit->Pc(),1./(48.*80.));
-          GetDigitsData(2+i)->Fill(digit->PadChX(),digit->PadChY());
-          GetDigitsData(9+i*6+digit->Pc())->Fill(digit->Q());
-       }  
-      }
-  }
-}
+  //
+  //filling QA histos for Digits
+  //
+   
+  Int_t i = fChannel ; 
+  FillDigitsData(0,i,fDigitsArray->GetEntriesFast()/(48.*80.*6.));
+  TIter next(fDigitsArray); 
+  AliHMPIDDigit * digit; 
+  while ( (digit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
+    FillDigitsData(1,10.*i+digit->Pc(),1./(48.*80.));
+    FillDigitsData(2+i,digit->PadChX(),digit->PadChY());
+    FillDigitsData(9+i*6+digit->Pc(),digit->Q());
+  }  
+}  
 //___________________________________________________________________________
 void AliHMPIDQADataMakerSim::MakeDigits(TTree * data)
 {
-//
-//Opening the Digit Tree
-//
- TObjArray *pObjDig=new TObjArray(AliHMPIDParam::kMaxCh+1);
-  for(Int_t iCh=AliHMPIDParam::kMinCh;iCh<=AliHMPIDParam::kMaxCh;iCh++){
-    TClonesArray *pCA=new TClonesArray("AliHMPIDDigit");
-    pObjDig->AddAt(pCA,iCh);
-  }
-
-  pObjDig->SetOwner(kTRUE);
-
+  //
+  //Opening the Digit Tree
+  //
+  
+  if(fDigitsArray) 
+    fDigitsArray->Clear() ; 
+  else
+    fDigitsArray=new TClonesArray("AliHMPIDDigit");
+  
   for(Int_t iCh=AliHMPIDParam::kMinCh;iCh<=AliHMPIDParam::kMaxCh;iCh++){
-    data->SetBranchAddress(Form("HMPID%i",iCh),&(*pObjDig)[iCh]);
+    fChannel = iCh ; 
+    data->SetBranchAddress(Form("HMPID%i",iCh),&fDigitsArray);
+    data->GetEntry(0);
+    MakeDigits();
+    fDigitsArray->Clear() ; 
   }
-  data->GetEntry(0);
-
-   MakeDigits((TClonesArray *)pObjDig);
 }
+
 //____________________________________________________________________________
 
-void AliHMPIDQADataMakerSim::MakeSDigits(TClonesArray * data)
+void AliHMPIDQADataMakerSim::MakeSDigits()
 {
  //
  //filling QA histos for SDigits
  //
-  TClonesArray * sdigits = dynamic_cast<TClonesArray *>(data) ; 
-  if (!sdigits) {
-    AliError("Wrong type of sdigits container") ; 
-  } else {
-    TIter next(sdigits) ; 
-    AliHMPIDDigit * sdigit ; 
-    while ( (sdigit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
-           GetSDigitsData(0)->Fill(sdigit->Q());
-    } 
-  }
+ 
+  TIter next(fSDigitsArray) ; 
+  AliHMPIDDigit * sdigit ; 
+  while ( (sdigit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
+    FillSDigitsData(0,sdigit->Q());
+  } 
 }
 //___________________________________________________________________________
 void AliHMPIDQADataMakerSim::MakeSDigits(TTree * data)
@@ -223,16 +217,19 @@ void AliHMPIDQADataMakerSim::MakeSDigits(TTree * data)
  //
  // Opening the SDigit Tree
  //
- TClonesArray * sdigits = new TClonesArray("AliHMPIDDigit", 1000) ;
+ if (fSDigitsArray)
+   fSDigitsArray->Clear() ; 
+  else 
+    fSDigitsArray = new TClonesArray("AliHMPIDDigit", 1000) ;
 
   TBranch * branch = data->GetBranch("HMPID") ;
   if ( ! branch ) {
     AliError("HMPID SDigit Tree not found") ;
     return;
   }
-  branch->SetAddress(&sdigits) ;
+  branch->SetAddress(&fSDigitsArray) ;
   branch->GetEntry(0) ;
-  MakeSDigits(sdigits) ;
+  MakeSDigits() ;
 }
 //____________________________________________________________________________
 void AliHMPIDQADataMakerSim::StartOfDetectorCycle()
@@ -245,6 +242,7 @@ void AliHMPIDQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjA
 {
   //Detector specific actions at end of cycle
   // do the QA checking
+  ResetEventTrigClasses(); // reset triggers list to select all histos
   AliQAChecker::Instance()->Run(AliQAv1::kHMPID, task, obj) ;  
 }