#ifndef __CINT__
  #include <Riostream.h>
  #include "AliCascadeVertexer.h"
  #include "TFile.h"
  #include "TStopwatch.h"
#endif

Int_t AliCascadeFindVertices() {
   cerr<<"Looking for cascade vertices...\n";

   TFile *out=TFile::Open("AliCascadeVertices.root","new");
   if (!out->IsOpen()) {
      cerr<<"Delete old AliCascadeVertices.root !\n"; return 1;
   }
   TFile *in=TFile::Open("AliITStracksV2.root");
   if (!in->IsOpen()) {cerr<<"Can't open AliITStracksV2.root !\n"; return 2;}

   TFile *file=TFile::Open("AliV0vertices.root");
   if (!file->IsOpen()) {
      cerr<<"Can't open AliV0vertices.root !\n";return 3;
   }
   Double_t cuts[]={33.,    // max. allowed chi2
                    0.05,   // min. allowed V0 impact parameter 
                    0.008,  // window around the Lambda mass 
                    0.035,  // min. allowed bachelor's impact parameter 
                    0.10,   // max. allowed DCA between a V0 and a track
                    0.9985, // max. allowed cosine of the cascade pointing angle
                    0.9,    // min. radius of the fiducial volume
                    2.9     // max. radius of the fiducial volume
                   };
   TStopwatch timer;
   AliCascadeVertexer *vertexer=new AliCascadeVertexer(cuts);
   Int_t rc=vertexer->V0sTracks2CascadeVertices(in,out);
   delete vertexer;
   timer.Stop(); timer.Print();
    
   file->Close();
   in->Close();
   out->Close();

   return rc;
}