Check PHOS and EMCAL clusters

diff --git a/STEER/CheckESD.C b/STEER/CheckESD.C
index aac7a8e3dee34d96ac7b0487282bd943ce966193..a918f27dcadffdb184f1d9ef64ba23cd93b4f28f 100644 (file)
--- a/STEER/CheckESD.C
+++ b/STEER/CheckESD.C
@@ -221,7 +221,7 @@ Bool_t CheckESD(const char* gAliceFileName = "galice.root",
 
   // calorimeters
   TH1F* hEPHOS = CreateHisto("hEPHOS", "PHOS", 100, 0, 5, "E [GeV]", "N");
-  TH1F* hEEMCAL = CreateHisto("hEEMCAL", "EMCAL", 100, 0, 2, "E [GeV]", "N");
+  TH1F* hEEMCAL = CreateHisto("hEEMCAL", "EMCAL", 100, 0, 50, "E [GeV]", "N");
 
   // muons
   TH1F* hPtMUON = CreateHisto("hPtMUON", "MUON", 100, 0, 20, 
@@ -435,6 +435,22 @@ Bool_t CheckESD(const char* gAliceFileName = "galice.root",
       selCascades.Remove(particle);
       nRecCascades++;
     }
+
+    // loop over the PHOS clusters
+    {
+    Int_t firstPHOSCluster = esd->GetFirstPHOSCluster();
+    Int_t lastPHOSCluster  = firstPHOSCluster + esd->GetNumberOfPHOSClusters();
+    for (Int_t iCluster=firstPHOSCluster; iCluster<lastPHOSCluster; iCluster++)
+      hEPHOS->Fill(esd->GetCaloCluster(iCluster)->GetClusterEnergy());
+    }
+
+    // loop over the EMCAL clusters
+    {
+    Int_t firstEMCALCluster = esd->GetFirstEMCALCluster();
+    Int_t lastEMCALCluster  = firstEMCALCluster + esd->GetNumberOfEMCALClusters();
+    for (Int_t iCluster=firstEMCALCluster; iCluster<lastEMCALCluster; iCluster++)
+      hEEMCAL->Fill(esd->GetCaloCluster(iCluster)->GetClusterEnergy());
+    }
   }
 
   // perform checks