ITS/AliITSrefitV2.C

   1 #ifndef __CINT__
   2   #include <Riostream.h>
   3   #include "AliITSgeom.h"
   4   #include "AliITStrackerV2.h"
   5
   6   #include "TFile.h"
   7   #include "TStopwatch.h"
   8 #endif
   9
  10 Int_t AliITSrefitV2(Int_t nev=1) {
  11    cerr<<"Propagating tracks inward through the ITS...\n";
  12
  13    TFile *in=TFile::Open("AliTPCrefited.root");
  14    if (!in->IsOpen()) {cerr<<"Can't open AliTPCrefited.root !\n"; return 1;}
  15
  16    TFile *out=TFile::Open("AliITStracksV2.root","update");
  17    if (!out->IsOpen()) {
  18       cerr<<"Can't open AliITStracksV2.root !\n"; return 2;
  19    }
  20    TFile *file=TFile::Open("AliITSclustersV2.root");
  21    if (!file->IsOpen()) {
  22       cerr<<"Can't open AliITSclustersV2.root !\n";return 3;
  23    }
  24    AliITSgeom *geom=(AliITSgeom*)file->Get("AliITSgeom");
  25
  26    Int_t rc=0;
  27    TStopwatch timer;
  28    AliITStrackerV2 tracker(geom);
  29    for (Int_t i=0; i<nev; i++) {
  30      cerr<<"Processing event number : "<<i<<endl;
  31      tracker.SetEventNumber(i);
  32      rc=tracker.RefitInward(in,out);
  33    }
  34    timer.Stop(); timer.Print();
  35
  36    delete geom;
  37
  38    file->Close();
  39    in->Close();
  40    out->Close();
  41
  42    return rc;
  43 }