ITS/AliV0FindVertices.C

   1 #ifndef __CINT__
   2   #include "Riostream.h"
   3   #include "AliV0vertexer.h"
   4   #include "TFile.h"
   5   #include "TStopwatch.h"
   6 #endif
   7
   8 Int_t AliV0FindVertices(Int_t nev=1) {
   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 vtxer(cuts);
  27    Int_t rc=0;
  28    for (Int_t i=0; i<nev; i++) {
  29      //Double_t vtx[3]={0.,0.,0.}; vtxer.SetVertex(vtx); // primary vertex (cm)
  30      vtxer.SetEvent(i);
  31      rc=vtxer.Tracks2V0vertices(in,out);
  32    }
  33    timer.Stop(); timer.Print();
  34
  35    in->Close();
  36    out->Close();
  37
  38    return rc;
  39 }