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();
char name[30];
sprintf(name,"Vertex_%d",evnumb);
fCurrentVertex = new AliESDVertex(vrtx,fSmear,name);
- fCurrentVertex->SetTruePos(vrttrue);
return fCurrentVertex;
}
}
}
if(vert){
- Double_t pos[3];
- for(Int_t kk=0;kk<3;kk++)pos[kk]=(Double_t)primaryVertex[kk];
- vert->SetTruePos(pos);
Float_t found = vert->GetZv();
diff2->Fill(primaryVertex[2],found);
found = 10000.*(found-primaryVertex[2]);