4c3d32ed2bc614f7df31e77bd8106413e5185060
[u/mrichter/AliRoot.git] / TPC / AliTPCFindTracks.C
1 #ifndef __CINT__
2   #include <iostream.h>
3   #include "AliTPCParam.h"
4   #include "AliTPCtracker.h"
5
6   #include "TFile.h"
7   #include "TStopwatch.h"
8 #endif
9
10 Int_t AliTPCFindTracks(Int_t eventn=1) { 
11    cerr<<"Looking for tracks...\n";
12
13    TFile *out=TFile::Open("AliTPCtracks.root","new");
14    if (!out->IsOpen()) {cerr<<"Delete old AliTPCtracks.root !\n"; return 1;}
15
16    TFile *in=TFile::Open("AliTPCclusters.root");
17    if (!in->IsOpen()) {cerr<<"Can't open AliTPCclusters.root !\n"; return 2;}
18
19    AliTPCParam *par=(AliTPCParam*)in->Get("75x40_100x60_150x60");
20    if (!par) {cerr<<"Can't get TPC parameters !\n"; return 3;}
21  
22    TStopwatch timer;
23
24    Int_t rc=0;
25    for (Int_t i=0;i<eventn;i++){
26      printf("Processing event %d\n",i);
27      AliTPCtracker *tracker = new AliTPCtracker(par,i);
28      //Double_t xyz[]={0.,0.,0.}; tracker->SetVertex(xyz); //primary vertex
29      rc=tracker->Clusters2Tracks(0,out);
30      delete tracker;
31    }
32    timer.Stop(); timer.Print();
33  
34    delete par; //Thanks to Mariana Bondila
35
36    in->Close();
37    out->Close();
38
39    return rc;
40 }