* Origin: Christian Kuhn, IReS, Strasbourg, christian.kuhn@ires.in2p3.fr *
****************************************************************************/
+#include <TPDGCode.h>
#if !defined(__CINT__) || defined(__MAKECINT__)
#include <TMath.h>
#include <TError.h>
#include "AliTrackReference.h"
#include "AliRunLoader.h"
#include "AliRun.h"
- #include "AliESD.h"
-#else
-const Int_t kXiMinus = 3312;
-const Int_t kXiPlusBar = -3312;
-const Int_t kOmegaMinus = 3334;
-const Int_t kOmegaPlusBar = -3334;
+ #include "AliESDEvent.h"
+ #include "AliESDcascade.h"
#endif
Int_t GoodCascades(const Char_t *dir=".");
return 5;
}
}
- AliESD* event = new AliESD;
+ AliESDEvent* event = new AliESDEvent();
TTree* esdTree = (TTree*) ef->Get("esdTree");
if (!esdTree) {
::Error("AliCascadeComparison.C", "no ESD tree found");
return 6;
}
- esdTree->SetBranchAddress("ESD", &event);
+ event->ReadFromTree(esdTree);
//******* Loop over events *********
} //**** End of the loop over events
delete event;
+ delete esdTree;
ef->Close();
delete csTree;
refFile->Close();
- Stat_t ng=hgood->GetEntries(), nf=hfound->GetEntries();
- if (ng!=0) cout<<"Integral efficiency is about "<<nf/ng*100.<<" %\n";
+ Stat_t ngg=hgood->GetEntries(), nf=hfound->GetEntries();
+ if (ngg!=0) cout<<"Integral efficiency is about "<<nf/ngg*100.<<" %\n";
cout<<
"Total number of found cascades: "<<allfound<<" ("<<nf<<" in the peak)\n";
cout<<"Total number of \"good\" cascades: "<<allgood<<endl;
Int_t GoodCascades(const Char_t *dir) {
if (gAlice) {
- delete gAlice->GetRunLoader();
+ delete AliRunLoader::Instance();
delete gAlice;
gAlice=0;
}