ITS/AliCascadeFindVertices.C

   1 #ifndef __CINT__
   2   #include <iostream.h>
   3   #include "AliCascadeVertexer.h"
   4   #include "TFile.h"
   5   #include "TStopwatch.h"
   6 #endif
   7
   8 Int_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    }
  22    Double_t cuts[]={33.,  // max. allowed chi2
  23                     0.015,// min. allowed V0 impact parameter
  24                     0.05, // window around the Lambda mass
  25                     0.015,// min. allowed track impact parameter
  26                     0.060,// max. allowed DCA between a V0 and a track
  27                     0.997,// 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
  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 }