X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=STEER%2FCheckESD.C;h=0784298c354b61bf4adfa07c7a4e3df12835c0ea;hb=40389866548fc06fcc5655701a58bb65f61e09a6;hp=9d823e3eeccfa79f9c613e338e807a88b2365093;hpb=cbc28ac7900dd4244096793957e903627cc3d725;p=u%2Fmrichter%2FAliRoot.git diff --git a/STEER/CheckESD.C b/STEER/CheckESD.C index 9d823e3eecc..0784298c354 100644 --- a/STEER/CheckESD.C +++ b/STEER/CheckESD.C @@ -12,16 +12,19 @@ #include "AliRunLoader.h" #include "AliLoader.h" -#include "AliESD.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, @@ -158,7 +161,7 @@ Bool_t CheckESD(const char* gAliceFileName = "galice.root", Error("CheckESD", "opening ESD file %s failed", esdFileName); return kFALSE; } - AliESD* esd = new AliESD; + AliESDEvent * esd = new AliESDEvent; TTree* tree = (TTree*) esdFile->Get("esdTree"); if (!tree) { Error("CheckESD", "no ESD tree found"); @@ -166,7 +169,7 @@ Bool_t CheckESD(const char* gAliceFileName = "galice.root", } esd->ReadFromTree(tree); - // efficienc and resolution histograms + // efficiency and resolution histograms Int_t nBinsPt = 15; Float_t minPt = 0.1; Float_t maxPt = 3.1; @@ -438,21 +441,15 @@ Bool_t CheckESD(const char* gAliceFileName = "galice.root", 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; iClusterFill(esd->GetCaloCluster(iCluster)->E()); + for (Int_t iCluster=0; iClusterGetNumberOfCaloClusters(); 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; iClusterFill(esd->GetCaloCluster(iCluster)->E()); - } } // perform checks