ITS/AliV0FindVertices.C

   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.015,// min. allowed negative daughter's impact parameter
  19                     0.015,// min. allowed positive daughter's impact parameter
  20                     0.060,// max. allowed DCA between the daughter tracks
  21                     0.997,// 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 }