]>
Commit | Line | Data |
---|---|---|
7f6ddf58 | 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.015,// min. allowed negative daughter's impact parameter | |
19 | 0.015,// min. allowed positive daughter's impact parameter | |
20 | 0.060,// max. allowed DCA between the daughter tracks | |
21 | 0.997,// 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 | } |