| 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.16,// min. allowed negative daughter's impact parameter |
| 19 | 0.05,// 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 | } |