]>
Commit | Line | Data |
---|---|---|
b9de75e1 | 1 | #ifndef __CINT__ |
2 | #include <iostream.h> | |
7f6ddf58 | 3 | #include "AliTPCParam.h" |
b9de75e1 | 4 | #include "AliTPCtracker.h" |
5 | ||
6 | #include "TFile.h" | |
7 | #include "TStopwatch.h" | |
8 | #endif | |
9 | ||
afc42102 | 10 | Int_t AliTPCFindTracks(Int_t eventn=1) { |
73042f01 | 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 | ||
b9de75e1 | 19 | AliTPCParam *par=(AliTPCParam*)in->Get("75x40_100x60"); |
20 | if (!par) {cerr<<"Can't get TPC parameters !\n"; return 3;} | |
21 | ||
73042f01 | 22 | TStopwatch timer; |
afc42102 | 23 | |
7f6ddf58 | 24 | Int_t rc=0; |
afc42102 | 25 | for (Int_t i=0;i<eventn;i++){ |
26 | printf("Processing event %d\n",i); | |
27 | AliTPCtracker *tracker = new AliTPCtracker(par,i); | |
7f6ddf58 | 28 | //Double_t xyz[]={0.,0.,0.}; tracker->SetVertex(xyz); //primary vertex |
29 | rc=tracker->Clusters2Tracks(0,out); | |
ec48216e | 30 | delete tracker; |
afc42102 | 31 | } |
73042f01 | 32 | timer.Stop(); timer.Print(); |
b9de75e1 | 33 | |
f38c8ae5 | 34 | delete par; //Thanks to Mariana Bondila |
35 | ||
73042f01 | 36 | in->Close(); |
37 | out->Close(); | |
38 | ||
7f6ddf58 | 39 | return rc; |
73042f01 | 40 | } |