AliGeomManager::LoadGeometry("geometry.root");
- AliITSVertexerZ *dovert = new AliITSVertexerZ("default",initx,inity);
+ 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";
cout << "Event number: "<<i<<") Z Vertex:"<<endl;
if(vert){
- cout<<"FOUND: "<<vert->GetZv()<<"; ";
+ cout<<"FOUND: "<<vert->GetZ()<<"; ";
cout<<vert->GetZRes()<<"; "<<vert->GetNContributors()<<endl;
cout <<" True Z position "<<primaryVertex[2]<<", diff= ";
- cout<<(primaryVertex[2]-vert->GetZv())*10000.<<endl;
+ cout<<(primaryVertex[2]-vert->GetZ())*10000.<<endl;
} else {
cout<<"NOT FOUND"<<endl;
}
}
if(vert){
- Float_t found = vert->GetZv();
+ Float_t found = vert->GetZ();
diff2->Fill(primaryVertex[2],found);
found = 10000.*(found-primaryVertex[2]);
if(vert->GetZRes()!=0){