#include <AliRun.h>
#include <AliHeader.h>
#include <AliGenEventHeader.h>
+#include <AliGeomManager.h>
#include <AliITSVertexerZ.h>
#include <AliRunLoader.h>
#include <AliITSLoader.h>
cerr<<"AliITSVertexerZTest.C : LoadKinematics returned error"<<endl;
return;
}
- TGeoManager::Import("geometry.root");
- AliITSLoader* ITSloader = (AliITSLoader*) rl->GetLoader("ITSLoader");
+
+ AliGeomManager::LoadGeometry("geometry.root");
- if(!ITSloader){
- cerr<<"AliITSVertexerZTest.C : ITS loader not found"<<endl;
- return;
- }
- // ITSloader->LoadRecPoints("read");
- // TFile *fo = new TFile("vertici.root","recreate");
- // AliITSVertexerPPZ *dovert = new AliITSVertexerPPZ("default",initx,inity);
- AliITSVertexerZ *dovert = new AliITSVertexerZ("default",initx,inity);
- dovert->SetDebug(0);
+ AliITSVertexerZ *dovert = new AliITSVertexerZ(initx,inity);
+ dovert->Init("default");
+ //dovert->SetDebug(0);
// dovert->SetDiffPhiMax(delphi);
// dovert->SetWindow(window);
dovert->PrintStatus();
Int_t sigmazero=0;
AliESDVertex *vert = 0;
+ AliITSLoader* itsloader = (AliITSLoader*) rl->GetLoader("ITSLoader");
+ itsloader->LoadRecPoints("read");
for(Int_t i=0; i<rl->TreeE()->GetEntries(); i++){
rl->GetEvent(i);
// The true Z coord. is fetched for comparison
AliGenEventHeader* genEventHeader = header->GenEventHeader();
TArrayF primaryVertex(3);
genEventHeader->PrimaryVertex(primaryVertex);
- vert = dovert->FindVertexForCurrentEvent(i);
+ TTree* cltree = itsloader->TreeR();
+ vert = dovert->FindVertexForCurrentEvent(cltree);
if(kDebug>0){
// Prints the results
cout <<"========================================================\n";