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 AliTPCFindTracksMI(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());
48 AliTPCParam *dig=(AliTPCParam *)gDirectory->Get("75x40_100x60_150x60");
51 dig=(AliTPCParam *)gDirectory->Get("75x40_100x60");
54 cerr<<"TPC parameters have not been found !\n";
59 cout<<"TPC 75x40_100x60 geometry found"<<endl;
64 cout<<"TPC 75x40_100x60_150x60 geometry found"<<endl;
68 tpcl->LoadTracks("recreate");
73 eventn = rl->GetNumberOfEvents();
81 for (Int_t i=0;i<eventn;i++)
83 TTree * input = tpcl->TreeR();
86 tpcl->LoadRecPoints("read");
87 input = tpcl->TreeR();
90 cerr << "Problems with input tree (TreeR) for event " << i <<endl;
94 TTree * output = tpcl->TreeT();
98 output = tpcl->TreeT();
101 cerr << "Problems with output tree (TreeT) for event " << i <<endl;
106 printf("Processing event %d\n",i);
107 AliTPCtrackerMI *tracker = new AliTPCtrackerMI(dig);
108 rc=tracker->Clusters2Tracks();
111 timer.Stop(); timer.Print();
113 delete dig; //Thanks to Mariana Bondila