1 /****************************************************************************
2 * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch *
3 ****************************************************************************/
7 #include "AliTPCParam.h"
8 #include "AliTPCtracker.h"
11 #include "TStopwatch.h"
14 Int_t AliTPCFindTracks(Int_t N=-1) {
16 cerr<<"Looking for tracks...\n";
20 delete gAlice->GetRunLoader();
25 rl = AliRunLoader::Open("galice.root");
28 cerr<<"Can not open session"<<endl;
31 tpcl = (AliTPCLoader*)rl->GetLoader("TPCLoader");
34 cerr<<"Can not get TPC Loader"<<endl;
40 cerr<<"Error occured while l"<<endl;
43 AliKalmanTrack::SetConvConst(1000/0.299792458/rl->GetAliRun()->Field()->SolenoidField());
45 AliTPCParam *dig=(AliTPCParam *)gDirectory->Get("75x40_100x60_150x60");
48 dig=(AliTPCParam *)gDirectory->Get("75x40_100x60");
51 cerr<<"TPC parameters have not been found !\n";
56 cout<<"TPC 75x40_100x60 geometry found"<<endl;
61 cout<<"TPC 75x40_100x60_150x60 geometry found"<<endl;
66 tpcl->LoadRecPoints("read");
67 tpcl->LoadTracks("recreate");
72 eventn = rl->GetNumberOfEvents();
80 for (Int_t i=0;i<eventn;i++){
81 printf("Processing event %d\n",i);
82 AliTPCtracker *tracker = new AliTPCtracker(dig,i);
83 rc=tracker->Clusters2Tracks();
86 timer.Stop(); timer.Print();
88 delete dig; //Thanks to Mariana Bondila