Adding new file also for the new Makefile
[u/mrichter/AliRoot.git] / ITS / AliV0FindVertices.C
CommitLineData
ca28c5f5 1#ifndef __CINT__
2 #include <iostream.h>
3 #include "AliV0vertexer.h"
4 #include "TFile.h"
5 #include "TStopwatch.h"
6#endif
7
8Int_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.12,// min. allowed negative daughter's impact parameter
19 0.06,// min. allowed positive daughter's impact parameter
20 0.080,// max. allowed DCA between the daughter tracks
21 0.998,// 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}