TPC/AliTPCFindTracks.C

   1 Int_t AliTPCFindTracks() {
   2    cerr<<"Looking for tracks...\n";
   3
   4    TFile *out=TFile::Open("AliTPCtracks.root","new");
   5    if (!out->IsOpen()) {cerr<<"Delete old AliTPCtracks.root !\n"; return 1;}
   6
   7    TFile *in=TFile::Open("AliTPCclusters.root");
   8    if (!in->IsOpen()) {cerr<<"Can't open AliTPCclusters.root !\n"; return 2;}
   9
  10    AliTPCv2 TPC;
  11    AliTPCParam *digp= (AliTPCParam*)in->Get("75x40_100x60");
  12    if (!digp) {cerr<<"TPC parameters have not been found !\n"; return 3;}
  13    TPC.SetParam(digp);
  14
  15    TStopwatch timer;
  16    TPC.Clusters2Tracks(out);
  17    timer.Stop(); timer.Print();
  18
  19    in->Close();
  20    out->Close();
  21
  22    return 0;
  23 }