]>
Commit | Line | Data |
---|---|---|
ca28c5f5 | 1 | #ifndef __CINT__ |
ddae8318 | 2 | #include <Riostream.h> |
ca28c5f5 | 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 | } | |
04b2a5f1 | 22 | Double_t cuts[]={33., // max. allowed chi2 |
23 | 0.05, // min. allowed V0 impact parameter | |
24 | 0.008, // window around the Lambda mass | |
25 | 0.035, // min. allowed bachelor's impact parameter | |
26 | 0.10, // max. allowed DCA between a V0 and a track | |
27 | 0.9985, // 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 | |
ca28c5f5 | 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 | } |