1 /****************************************************************************
2 * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch *
3 ****************************************************************************/
5 #if !defined(__CINT__) || defined(__MAKECINT__)
7 #include "AliTPCParam.h"
8 #include "AliTPCtracker.h"
9 #include "AliTPCpidESD.h"
12 #include "AliRunLoader.h"
13 #include "AliTPCLoader.h"
14 #include "AliESDpid.h"
18 #include "TStopwatch.h"
21 extern AliRun *gAlice;
23 Int_t AliTPCFindTracks(Int_t nev=5) {
25 cerr<<"Looking for tracks...\n";
28 delete gAlice->GetRunLoader();
33 AliRunLoader *rl = AliRunLoader::Open("galice.root");
35 cerr<<"Can not open session"<<endl;
39 AliTPCLoader *tpcl = (AliTPCLoader*)rl->GetLoader("TPCLoader");
41 cerr<<"Can not get TPC Loader"<<endl;
45 if (rl->LoadgAlice()) {
46 cerr<<"Error occured while loading gAlice"<<endl;
50 AliKalmanTrack::SetConvConst(
51 1000/0.299792458/rl->GetAliRun()->Field()->SolenoidField()
55 AliTPCParam *dig=(AliTPCParam *)gDirectory->Get("75x40_100x60_150x60");
57 cerr<<"TPC parameters have not been found !\n";
63 tpcl->LoadRecPoints("read");
65 if (nev>rl->GetNumberOfEvents()) nev=rl->GetNumberOfEvents();
67 TFile *ef=TFile::Open("AliESDtpc.root","RECREATE");
68 if ((!ef)||(!ef->IsOpen())) {
69 cerr<<"Can't AliESDtpc.root !\n"; return 1;
73 AliTPCtracker tracker(dig);
74 //An instance of the TPC PID maker
75 Double_t parTPC[]={45.0,0.08,10.}; // normalization constants,
76 // see AliTPCpidESD class
77 AliTPCpidESD tpcPID(parTPC);
78 for (Int_t i=0;i<nev;i++){
79 printf("Processing event %d\n",i);
80 AliESD *event=new AliESD();
83 TTree *in=tpcl->TreeR();
85 cerr<<"Can't get clusters tree !\n";
89 tracker.LoadClusters(in);
90 rc=tracker.Clusters2Tracks(event);
91 tpcPID.MakePID(event); // Preliminary PID
92 AliESDpid::MakePID(event); // for the ITS tracker
93 tracker.UnloadClusters();
97 sprintf(ename,"%d",i);
99 if (!event->Write(ename)) rc++;
102 cerr<<"Something bad happened...\n";
107 timer.Stop(); timer.Print();
111 delete dig; //Thanks to Mariana Bondila