Bugfix in AliL3Transformer Sin/Cos tables. Rotation to slice 0 was hidden by 35.
[u/mrichter/AliRoot.git] / ITS / AliCascadeFindVertices.C
CommitLineData
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
8Int_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}