#if !defined(__CINT__) || defined(__MAKECINT__)
+#include <TROOT.h>
#include <TClassTable.h>
+#include <TStopwatch.h>
#include <TFile.h>
#include <TTree.h>
#include <Riostream.h>
#include <AliRunLoader.h>
#include <AliITSVertexerIons.h>
#include <AliITSLoader.h>
+#include <AliGenHijingEventHeader.h>
#include <unistd.h>
#endif
// Open output file for vertices (default name: ITS.Vertex.root
// and Create vertexer
- AliITSVertexerIons *vertexer = new AliITSVertexerIons("default");
+ AliITSVertexerIons *vertexer = new AliITSVertexerIons();
+ vertexer->Init("default");
//vertexer->SetDebug(1);
AliESDVertex *V;
// Loop over events
+ AliITSLoader* itsloader = (AliITSLoader*) rl->GetLoader("ITSLoader");
+ itsloader->LoadRecPoints("read");
+
for (int nev=evNumber1; nev< evNumber2; nev++) {
cout<<"=============================================================\n";
cout<<" Processing event "<<nev<<endl;
TStopwatch timer;
timer.Start();
- V=vertexer->FindVertexForCurrentEvent(nev);
+ TTree* cltree = itsloader->TreeR();
+ V=vertexer->FindVertexForCurrentEvent(cltree);
TVector3 vtrue(primaryVertex[0],primaryVertex[1],primaryVertex[2]);
TVector3 vfound(V->GetXv(),V->GetYv(),V->GetZv());
cout << "True vertex coordinates (cm) = " << vtrue.X() << " " << vtrue.Y() << " " << vtrue.Z() << endl;
cout << "Found vertex coordinates (cm) = " << vfound.X() << " " << vfound.Y() << " " << vfound.Z() << endl; cout << "Difference true - found (cm) = " << dif.Mag() << " " << dif.X() << " " << dif.Y() << " " << dif.Z() << endl;
- if(V){
- Double_t pos[3];
- for(Int_t kk=0;kk<3;kk++)pos[kk]=(Double_t)primaryVertex[kk];
- V->SetTruePos(pos);
- }
timer.Stop();
timer.Print();