TPC/AliTPCrefit.C

   1 /****************************************************************************
   2  *           Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch                 *
   3  ****************************************************************************/
   4
   5 #ifndef __CINT__
   6   #include <Riostream.h>
   7   #include "AliTPCtracker.h"
   8
   9   #include "TFile.h"
  10   #include "TStopwatch.h"
  11 #endif
  12
  13 Int_t AliTPCrefit() {
  14
  15   cout << "Propagating tracks inward through the TPC..." << endl;
  16
  17   AliKalmanTrack::SetConvConst(1000/0.299792458/4.);
  18
  19   TFile *in = new TFile("AliTPCBackTracks.root", "READ");
  20   TFile *out= new TFile("AliTPCrefited.root", "RECREATE");
  21   TFile *file = new TFile("AliTPCclusters.root");
  22
  23   AliTPCParam *param=(AliTPCParam*)file->Get("75x40_100x60_150x60");
  24   if (!param) {cerr<<"Can't get TPC parameters !\n"; return 4;}
  25
  26   TStopwatch timer;
  27   AliTPCtracker *tracker = new AliTPCtracker(param);
  28   Int_t rc=tracker->RefitInward(in,out);
  29   delete tracker;
  30   timer.Stop(); timer.Print();
  31
  32   file->Close();
  33   in->Close();
  34   out->Close();
  35
  36   return rc;
  37 }