Updated code for tracking V2 (from Y. Belikov)
[u/mrichter/AliRoot.git] / ITS / AliV0FindVertices.C
1 #ifndef __CINT__
2   #include <iostream.h>
3   #include "AliV0vertexer.h"
4   #include "TFile.h"
5   #include "TStopwatch.h"
6 #endif
7
8 Int_t AliV0FindVertices() {
9    cerr<<"Looking for V0 vertices...\n";
10
11    TFile *out=TFile::Open("AliV0vertices.root","new");
12    if (!out->IsOpen()) {cerr<<"Delete old AliV0vertices.root !\n"; return 1;}
13
14    TFile *in=TFile::Open("AliITStracksV2.root");
15    if (!in->IsOpen()) {cerr<<"Can't open AliITStracksV2.root !\n"; return 2;}
16
17    Double_t cuts[]={33,  // max. allowed chi2
18                     0.02,// min. allowed negative daughter's impact parameter 
19                     0.02,// min. allowed positive daughter's impact parameter 
20                     0.080,// max. allowed DCA between the daughter tracks
21                     0.9985,// 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
24                    };
25    TStopwatch timer;
26    AliV0vertexer *vertexer=new AliV0vertexer(cuts);
27    Int_t rc=vertexer->Tracks2V0vertices(in,out);
28    delete vertexer;
29    timer.Stop(); timer.Print();
30     
31    in->Close();
32    out->Close();
33
34    return rc;
35 }