Further trigger QA development

diff --git a/PWGGA/PHOSTasks/PHOS_TriggerQA/AliAnalysisTaskPHOSTriggerQA.cxx b/PWGGA/PHOSTasks/PHOS_TriggerQA/AliAnalysisTaskPHOSTriggerQA.cxx
index c515914c068b311073336cc1d8bc5955bcf9c07f..3ca686f415e94daaecf8f6bfece52da322ef6b34 100644 (file)
--- a/PWGGA/PHOSTasks/PHOS_TriggerQA/AliAnalysisTaskPHOSTriggerQA.cxx
+++ b/PWGGA/PHOSTasks/PHOS_TriggerQA/AliAnalysisTaskPHOSTriggerQA.cxx
@@ -84,7 +84,7 @@ void AliAnalysisTaskPHOSTriggerQA::UserCreateOutputObjects()
   Int_t nTrMax = 1200;
   Float_t trMax = 600.;
 
-  fOutputContainer->Add(new TH1F("hNtr","Number of fired 4x4 regions",nTrMax,0.,trMax));
+  fOutputContainer->Add(new TH1F("hNtr","Number of fired 4x4 regions per event",nTrMax,0.,trMax));
 
   for(Int_t sm=1; sm<4; sm++) {
 
@@ -96,6 +96,10 @@ void AliAnalysisTaskPHOSTriggerQA::UserCreateOutputObjects()
     snprintf(titl,55,"SM%d 4x4 occupancy",sm);
     fOutputContainer->Add(new TH2F(key,titl,nRows,0.,nRows,nCols,0.,nCols));
 
+    snprintf(key,55,"h4x4CluSM%d",sm);
+    snprintf(titl,55,"SM%d 4x4 occupancy associated with clusters (E>2GeV)",sm);
+    fOutputContainer->Add(new TH2F(key,titl,nRows,0.,nRows,nCols,0.,nCols));
+
     snprintf(key,55,"hCluSM%d",sm);
     snprintf(titl,55,"SM%d cluster occupancy",sm);
     fOutputContainer->Add(new TH2F(key,titl,nRows,0.,nRows,nCols,0.,nCols));
@@ -108,9 +112,22 @@ void AliAnalysisTaskPHOSTriggerQA::UserCreateOutputObjects()
     snprintf(titl,55,"SM%d cluster energy",sm);
     fOutputContainer->Add(new TH1F(key,titl,nPtPhot,0.,ptPhotMax));
 
+    for(Int_t iTRU=1; iTRU<=8; iTRU++) {
+      snprintf(key,55,"hPhotAllSM%dTRU%d",sm,iTRU);
+      snprintf(titl,55,"SM%d: clusters energy in TRU%d",sm,iTRU);
+      fOutputContainer->Add(new TH1F(key,titl,nPtPhot,0.,ptPhotMax));
+    }
+    
     snprintf(key,55,"hPhotTrigSM%d",sm);
     snprintf(titl,55,"SM%d triggered cluster energy",sm);
     fOutputContainer->Add(new TH1F(key,titl,nPtPhot,0.,ptPhotMax));
+    
+    for(Int_t iTRU=1; iTRU<=8; iTRU++) {
+      snprintf(key,55,"hPhotTrigSM%dTRU%d",sm,iTRU);
+      snprintf(titl,55,"SM%d: triggered clusters energy in TRU%d",sm,iTRU);
+      fOutputContainer->Add(new TH1F(key,titl,nPtPhot,0.,ptPhotMax));
+    }
+    
   }
   
   PostData(1, fOutputContainer);
@@ -137,14 +154,9 @@ void AliAnalysisTaskPHOSTriggerQA::UserExec(Option_t *)
   AliESDCaloTrigger* trgESD = event->GetCaloTrigger("PHOS");
   trgESD->Reset();
   
-  if (!trgESD->GetEntries()) {
-    PostData(1, fOutputContainer);
-    return;
-  }
-  
-  FillHistogram("hNev",1.); // triggered events
+  if(trgESD->GetEntries()) FillHistogram("hNev",1.); // triggered events
   FillHistogram("hNtr",trgESD->GetEntries());
-
+  
   TString trigClasses = event->GetFiredTriggerClasses();
   printf("\nEvent %d: %d non-zero trigger digits %s\n",
         fEventCounter,trgESD->GetEntries(),trigClasses.Data());
@@ -183,12 +195,18 @@ void AliAnalysisTaskPHOSTriggerQA::UserExec(Option_t *)
       AliESDCaloCluster *c1 = event->GetCaloCluster(i);
       if(!c1->IsPHOS()) continue;
       
+      if(c1->E()<0.3) continue; 
+      if(c1->GetNCells()<3) continue ; 
+      
       Int_t maxId, relid[4];
       MaxEnergyCellPos(phsCells,c1,maxId);
       
       fPHOSGeo->AbsToRelNumbering(maxId, relid);
       snprintf(key,55,"hPhotAllSM%d",relid[0]);
       FillHistogram(key,c1->E());
+
+      snprintf(key,55,"hPhotAllSM%dTRU%d",relid[0],GetTRUNum(relid[2]-1,relid[3]-1));
+      FillHistogram(key,c1->E());
       
       snprintf(key,55,"hCluSM%d",relid[0]);
       FillHistogram(key,relid[2]-1,relid[3]-1);
@@ -197,10 +215,17 @@ void AliAnalysisTaskPHOSTriggerQA::UserExec(Option_t *)
 
        snprintf(key,55,"hPhotTrigSM%d",relid[0]);
        FillHistogram(key,c1->E());
+       
+       snprintf(key,55,"hPhotTrigSM%dTRU%d",relid[0],GetTRUNum(relid[2]-1,relid[3]-1));
+       FillHistogram(key,c1->E());
 
        snprintf(key,55,"hCluTSM%d",relid[0]);
        FillHistogram(key,relid[2]-1,relid[3]-1);
        
+       if(c1->E()>2.) { // Eclu > 2 GeV
+         snprintf(key,55,"h4x4CluSM%d",trelid[0]);
+         FillHistogram(key,trelid[2]-1,trelid[3]-1);
+       }       
        continue;
       }
       
@@ -283,3 +308,42 @@ Bool_t AliAnalysisTaskPHOSTriggerQA::Matched(Int_t *trig_relid, Int_t *cluster_r
 
   return kTRUE;
 }
+
+//_______________________________________________________________________________
+Int_t AliAnalysisTaskPHOSTriggerQA::GetTRUNum(Int_t cellX, Int_t cellZ)
+{
+  //Return TRU region number for given cell.
+  //cellX: [0-63], cellZ: [0-55]
+ 
+  Int_t iTRU=-111;
+
+  //RCU0: TRU 1,2
+  if(0<=cellX&&cellX<16) { 
+
+    if(0<=cellZ&&cellZ<28) iTRU=2;
+    else iTRU=1;
+  }
+
+  //RCU1: TRU 3,4
+  if(16<=cellX&&cellX<32) { 
+    
+    if(0<=cellZ&&cellZ<28) iTRU=4;
+    else iTRU=3;
+  }
+
+  //RCU2: TRU 5,6
+  if(32<=cellX&&cellX<48) { 
+    
+    if(0<=cellZ&&cellZ<28) iTRU=6;
+    else iTRU=5;
+  }
+  
+  //RCU3: TRU 7,8
+  if(48<=cellX&&cellX<64) { 
+    
+    if(0<=cellZ&&cellZ<28) iTRU=8;
+    else iTRU=7;
+  }
+  
+  return iTRU;
+}

diff --git a/PWGGA/PHOSTasks/PHOS_TriggerQA/AliAnalysisTaskPHOSTriggerQA.h b/PWGGA/PHOSTasks/PHOS_TriggerQA/AliAnalysisTaskPHOSTriggerQA.h
index 64ab691f1ef96db34506da27d96e5c7e29f512f4..0a68c15fb7906469c385f8db09b77664558f7925 100644 (file)
--- a/PWGGA/PHOSTasks/PHOS_TriggerQA/AliAnalysisTaskPHOSTriggerQA.h
+++ b/PWGGA/PHOSTasks/PHOS_TriggerQA/AliAnalysisTaskPHOSTriggerQA.h
@@ -35,6 +35,8 @@ private:
   void   MaxEnergyCellPos(AliESDCaloCells *cells, AliESDCaloCluster* clu, Int_t& maxId);
   Bool_t Matched(Int_t *trig_relid,Int_t *cluster_relid); //is cluster position coincides with 4x4 position?
 
+  Int_t GetTRUNum(Int_t cellX, Int_t cellZ);
+
 private:
 
   TList * fOutputContainer;   //final histogram container