3 #include "AliV0vertexer.h"
5 #include "TStopwatch.h"
8 Int_t AliV0FindVertices(Int_t nev=1) {
9 cerr<<"Looking for V0 vertices...\n";
11 TFile *out=TFile::Open("AliV0vertices.root","new");
12 if (!out->IsOpen()) {cerr<<"Delete old AliV0vertices.root !\n"; return 1;}
14 TFile *in=TFile::Open("AliITStracksV2.root");
15 if (!in->IsOpen()) {cerr<<"Can't open AliITStracksV2.root !\n"; return 2;}
17 Double_t cuts[]={33, // max. allowed chi2
18 0.16,// min. allowed negative daughter's impact parameter
19 0.05,// min. allowed positive daughter's impact parameter
20 0.080,// max. allowed DCA between the daughter tracks
21 0.998,// max. allowed cosine of V0's pointing angle
22 0.9, // min. radius of the fiducial volume
23 2.9 // max. radius of the fiducial volume
26 AliV0vertexer vtxer(cuts);
28 for (Int_t i=0; i<nev; i++) {
29 //Double_t vtx[3]={0.,0.,0.}; vtxer.SetVertex(vtx); // primary vertex (cm)
31 rc=vtxer.Tracks2V0vertices(in,out);
33 timer.Stop(); timer.Print();