#include "AliRunLoader.h"
#include "AliLoader.h"
#include "AliESDEvent.h"
+#include "AliESDv0.h"
+#include "AliESDcascade.h"
+#include "AliESDMuonTrack.h"
+#include "AliESDCaloCluster.h"
#include "AliRun.h"
#include "AliStack.h"
#include "AliHeader.h"
#include "AliGenEventHeader.h"
#include "AliPID.h"
-#else
+#endif
const Int_t kXiMinus = 3312;
const Int_t kOmegaMinus = 3334;
-#endif
TH1F* CreateHisto(const char* name, const char* title,
nRecCascades++;
}
- // loop over the PHOS clusters
+ // loop over the 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)->E());
+ for (Int_t iCluster=0; iCluster<esd->GetNumberOfCaloClusters(); iCluster++) {
+ AliESDCaloCluster * clust = esd->GetCaloCluster(iCluster);
+ if (clust->IsPHOS()) hEPHOS->Fill(clust->E());
+ if (clust->IsEMCAL()) hEEMCAL->Fill(clust->E());
+ }
}
- // 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)->E());
- }
}
// perform checks