3 #include "AliCascadeVertexer.h"
5 #include "TStopwatch.h"
8 Int_t AliCascadeFindVertices() {
9 cerr<<"Looking for cascade vertices...\n";
11 TFile *out=TFile::Open("AliCascadeVertices.root","new");
13 cerr<<"Delete old AliCascadeVertices.root !\n"; return 1;
15 TFile *in=TFile::Open("AliITStracksV2.root");
16 if (!in->IsOpen()) {cerr<<"Can't open AliITStracksV2.root !\n"; return 2;}
18 TFile *file=TFile::Open("AliV0vertices.root");
19 if (!file->IsOpen()) {
20 cerr<<"Can't open AliV0vertices.root !\n";return 3;
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
32 AliCascadeVertexer *vertexer=new AliCascadeVertexer(cuts);
33 Int_t rc=vertexer->V0sTracks2CascadeVertices(in,out);
35 timer.Stop(); timer.Print();