]>
Commit | Line | Data |
---|---|---|
1c53abe2 | 1 | /**************************************************************************** |
2 | * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch * | |
3 | ****************************************************************************/ | |
4 | ||
5 | #ifndef __CINT__ | |
434772a9 | 6 | #include <iostream.h> |
7 | #include "AliTPCParam.h" | |
8 | #include "AliTPCtracker.h" | |
9 | ||
10 | #include "TFile.h" | |
11 | #include "TStopwatch.h" | |
1c53abe2 | 12 | #endif |
13 | ||
434772a9 | 14 | Int_t AliTPCFindTracksMI(Int_t N=-1) { |
15 | ||
1c53abe2 | 16 | cerr<<"Looking for tracks...\n"; |
1c53abe2 | 17 | |
434772a9 | 18 | if (gAlice) |
19 | { | |
20 | delete gAlice->GetRunLoader(); | |
21 | delete gAlice; | |
22 | gAlice = 0x0; | |
23 | } | |
24 | ||
25 | rl = AliRunLoader::Open("galice.root"); | |
26 | if (rl == 0x0) | |
27 | { | |
28 | cerr<<"Can not open session"<<endl; | |
29 | return 1; | |
30 | } | |
31 | tpcl = (AliTPCLoader*)rl->GetLoader("TPCLoader"); | |
32 | if (tpcl == 0x0) | |
33 | { | |
34 | cerr<<"Can not get TPC Loader"<<endl; | |
35 | return 1; | |
36 | } | |
37 | ||
38 | if (rl->LoadgAlice()) | |
39 | { | |
40 | cerr<<"Error occured while l"<<endl; | |
41 | return 1; | |
42 | } | |
43 | AliKalmanTrack::SetConvConst(1000/0.299792458/rl->GetAliRun()->Field()->SolenoidField()); | |
44 | rl->UnloadgAlice(); | |
45 | ||
46 | rl->CdGAFile(); | |
47 | ||
48 | AliTPCParam *dig=(AliTPCParam *)gDirectory->Get("75x40_100x60_150x60"); | |
49 | if (!dig) | |
50 | { | |
51 | dig=(AliTPCParam *)gDirectory->Get("75x40_100x60"); | |
52 | if (!param) | |
53 | { | |
54 | cerr<<"TPC parameters have not been found !\n"; | |
55 | return 1; | |
56 | } | |
57 | else | |
58 | { | |
59 | cout<<"TPC 75x40_100x60 geometry found"<<endl; | |
60 | } | |
61 | } | |
62 | else | |
63 | { | |
64 | cout<<"TPC 75x40_100x60_150x60 geometry found"<<endl; | |
65 | } | |
1c53abe2 | 66 | |
434772a9 | 67 | |
68 | tpcl->LoadTracks("recreate"); | |
69 | ||
70 | Int_t eventn; | |
71 | if (N<=0) | |
72 | { | |
73 | eventn = rl->GetNumberOfEvents(); | |
74 | rl->UnloadHeader(); | |
75 | } | |
76 | else | |
77 | eventn = N; | |
78 | ||
1c53abe2 | 79 | TStopwatch timer; |
1c53abe2 | 80 | Int_t rc=0; |
434772a9 | 81 | for (Int_t i=0;i<eventn;i++) |
82 | { | |
83 | TTree * input = tpcl->TreeR(); | |
84 | if (input == 0x0) | |
85 | { | |
86 | tpcl->LoadRecPoints("read"); | |
87 | input = tpcl->TreeR(); | |
88 | if (input == 0x0) | |
89 | { | |
90 | cerr << "Problems with input tree (TreeR) for event " << i <<endl; | |
91 | continue; | |
92 | } | |
93 | } | |
94 | TTree * output = tpcl->TreeT(); | |
95 | if (output == 0x0) | |
96 | { | |
97 | tpcl->MakeTree("T"); | |
98 | output = tpcl->TreeT(); | |
99 | if (output == 0x0) | |
100 | { | |
101 | cerr << "Problems with output tree (TreeT) for event " << i <<endl; | |
102 | continue; | |
103 | } | |
104 | } | |
105 | ||
106 | printf("Processing event %d\n",i); | |
107 | AliTPCtrackerMI *tracker = new AliTPCtrackerMI(dig); | |
108 | rc=tracker->Clusters2Tracks(); | |
109 | delete tracker; | |
110 | } | |
1c53abe2 | 111 | timer.Stop(); timer.Print(); |
112 | ||
434772a9 | 113 | delete dig; //Thanks to Mariana Bondila |
114 | delete rl; | |
1c53abe2 | 115 | return rc; |
116 | } |