Updated code for tracking V2 (from Y. Belikov)
[u/mrichter/AliRoot.git] / ITS / AliCascadeFindVertices.C
1 #ifndef __CINT__
2   #include <iostream.h>
3   #include "AliCascadeVertexer.h"
4   #include "TFile.h"
5   #include "TStopwatch.h"
6 #endif
7
8 Int_t AliCascadeFindVertices() {
9    cerr<<"Looking for cascade vertices...\n";
10
11    TFile *out=TFile::Open("AliCascadeVertices.root","new");
12    if (!out->IsOpen()) {
13       cerr<<"Delete old AliCascadeVertices.root !\n"; return 1;
14    }
15    TFile *in=TFile::Open("AliITStracksV2.root");
16    if (!in->IsOpen()) {cerr<<"Can't open AliITStracksV2.root !\n"; return 2;}
17
18    TFile *file=TFile::Open("AliV0vertices.root");
19    if (!file->IsOpen()) {
20       cerr<<"Can't open AliV0vertices.root !\n";return 3;
21    }
22    Double_t cuts[]={33.,  // max. allowed chi2
23                     0.015,// min. allowed V0 impact parameter 
24                     0.05, // window around the Lambda mass 
25                     0.015,// min. allowed track impact parameter 
26                     0.060,// max. allowed DCA between a V0 and a track
27                     0.997,// max. allowed cosine of the cascade pointing angle
28                     0.9,  // min. radius of the fiducial volume
29                     2.9   // max. radius of the fiducial volume
30                    };
31    TStopwatch timer;
32    AliCascadeVertexer *vertexer=new AliCascadeVertexer(cuts);
33    Int_t rc=vertexer->V0sTracks2CascadeVertices(in,out);
34    delete vertexer;
35    timer.Stop(); timer.Print();
36     
37    file->Close();
38    in->Close();
39    out->Close();
40
41    return rc;
42 }