/[PAMELA software]/DarthVader/TrackerLevel2/src/TrkProcess.cpp
ViewVC logotype

Annotation of /DarthVader/TrackerLevel2/src/TrkProcess.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.14 - (hide annotations) (download)
Wed Mar 5 17:00:19 2008 UTC (16 years, 9 months ago) by pam-fi
Branch: MAIN
CVS Tags: v5r00, v9r00, v9r01, v6r01, v6r00
Changes since 1.13: +13 -1 lines
modified TrkSinglet, optimized DoTrack2, fixed bug in evaluation of effective angle

1 mocchiut 1.1 /**
2     * \file TrkProcess.cpp
3     * \author Elena Vannuccini
4     */
5     #include <TrkProcess.h>
6     #include <iostream>
7     using namespace std;
8 pam-fi 1.2 /**
9     * Create TrkProcess object, initializing the tracker processing variables at default values
10     */
11 mocchiut 1.1 TrkProcess::TrkProcess(){
12    
13     idrun = 0;
14 pam-fi 1.11 // dbg_mode.SetNone();
15 pam-fi 1.7 get1 = 0;
16     full1 = 0;
17     geth = 0;
18     get2 = 1;
19 mocchiut 1.1 standalone = false;
20     frame2 = "root";
21 pam-fi 1.6 frame1 = "root";
22     frameh = "root";
23 mocchiut 1.1 outdir = gSystem->WorkingDirectory();
24     pfolder = "/TrackerFolder";
25     if (!frame2.CompareTo("root", TString::kIgnoreCase)) ifroot2 = true;
26 pam-fi 1.2 else if (!frame2.CompareTo("hbook", TString::kIgnoreCase)) ifroot2 = false;
27 mocchiut 1.1
28     if (!frame1.CompareTo("root", TString::kIgnoreCase)) ifroot1 = true;
29 pam-fi 1.2 else if (!frame1.CompareTo("hbook", TString::kIgnoreCase)) ifroot1 = false;
30 mocchiut 1.1
31 pam-fi 1.6 if (!frameh.CompareTo("root", TString::kIgnoreCase)) ifrooth = true;
32     else if (!frameh.CompareTo("hbook", TString::kIgnoreCase)) ifrooth = false;
33    
34 mocchiut 1.1 file1 = "";
35     file2 = "";
36    
37 pam-fi 1.2 ostatus = 0;
38 pam-fi 1.10
39 pam-fi 1.14 npar = 0;
40    
41     nskip = 0;
42 pam-fi 1.11
43 pam-fi 1.13 // TrkParams::SetTrackingMode();
44     // TrkParams::SetPrecisionFactor();
45     // TrkParams::SetStepMin();
46     TrkParams::SetMiniDefault();
47 pam-fi 1.10 TrkParams::SetPFA();
48 mocchiut 1.1 };
49    
50 pam-fi 1.2 /**
51     * Create TrkProcess object, initializing the tracker processing variables according to input parameters.
52     * @param run Run ID
53     * @param f2 Pointer to output level2 file
54     */
55 mocchiut 1.1 TrkProcess::TrkProcess(ULong64_t run, TFile *f2){
56    
57     idrun = run;
58 pam-fi 1.2 // DEBUG = false;
59     // VERBOSE = false;
60 pam-fi 1.11 // dbg_mode.SetWarning();
61 pam-fi 1.7 get1 = 0;
62     full1 = 0;
63 mocchiut 1.1 get2 = 1;
64 pam-fi 1.6 geth = 0;
65 mocchiut 1.1 standalone = false;
66     frame2 = "root";
67 pam-fi 1.4 frame1 = "root";
68 pam-fi 1.6 frameh = "root";
69     pfolder = "/TrackerFolder";
70     file1 = "";
71 mocchiut 1.1 if(f2->IsOpen()){
72 pam-fi 1.6 file2 = f2->GetPath();
73     outdir = gSystem->DirName(gSystem->DirName(file2));
74 pam-fi 1.4 // check if the indicated output directory exists
75 pam-fi 1.6 FileStat_t t;
76     if( gSystem->GetPathInfo(outdir.Data(),t) )throw -12;
77 mocchiut 1.1 }else{
78 pam-fi 1.6 file2 = "";
79     outdir = gSystem->WorkingDirectory();
80 mocchiut 1.1 };
81    
82     if (!frame2.CompareTo("root", TString::kIgnoreCase)) ifroot2 = true;
83     else if (!frame2.CompareTo("hbook", TString::kIgnoreCase)) ifroot2 =false;
84    
85     if (!frame1.CompareTo("root", TString::kIgnoreCase)) ifroot1 = true;
86     else if (!frame1.CompareTo("hbook", TString::kIgnoreCase)) ifroot1 =false;
87    
88 pam-fi 1.6 if (!frameh.CompareTo("root", TString::kIgnoreCase)) ifrooth = true;
89     else if (!frameh.CompareTo("hbook", TString::kIgnoreCase)) ifrooth =false;
90    
91 pam-fi 1.2 ostatus = 0;
92 pam-fi 1.11
93     npar=0;
94 mocchiut 1.1
95 pam-fi 1.14 nskip = 0;
96    
97 pam-fi 1.13 // TrkParams::SetTrackingMode();
98     // TrkParams::SetPrecisionFactor();
99     // TrkParams::SetStepMin();
100     TrkParams::SetMiniDefault();
101 pam-fi 1.10 TrkParams::SetPFA();
102    
103 mocchiut 1.1 };
104 pam-fi 1.2 /**
105     * Set processing variables according to input custom parameters
106     */
107 pam-fi 1.12 int TrkProcess::HandleCustomPar(int ncustom, char *vcustom[]){
108 mocchiut 1.1
109     for (int i=0; i< ncustom; i++){
110 pam-fi 1.4
111 mocchiut 1.1 // -----------------------------------------------------//
112     if (!strcmp(vcustom[i], "-processFolder")){
113     if (++i >= ncustom) throw -3;
114     pfolder = vcustom[i];
115     continue;
116     }
117     // -----------------------------------------------------//
118 pam-fi 1.6 if (!strcmp(vcustom[i], "--get1") || !strcmp(vcustom[i], "--get-level1") ){
119 pam-fi 1.7 get1 = 1;
120     full1 = 0;
121     continue;
122     }
123     // -----------------------------------------------------//
124     if (!strcmp(vcustom[i], "--get1-full") || !strcmp(vcustom[i], "--get-level1-full") ){
125     get1 = 1;
126     full1 = 1;
127 mocchiut 1.1 continue;
128     }
129     // -----------------------------------------------------//
130 pam-fi 1.6 if (!strcmp(vcustom[i], "--geth") || !strcmp(vcustom[i], "--get-hough") ){
131     geth = 1;
132     continue;
133     }
134     // -----------------------------------------------------//
135     if (!strcmp(vcustom[i], "--dontget2") ||!strcmp(vcustom[i], "--dontget-level2") ){
136 mocchiut 1.1 get2 = 0;
137     continue;
138     }
139     // -----------------------------------------------------//
140     if (!strcmp(vcustom[i], "--standalone") || !strcmp(vcustom[i], "-s")){
141     standalone = true;
142     continue;
143     }
144     // -----------------------------------------------------//
145     if (!strcmp(vcustom[i], "-processFile1")){
146     if (++i >= ncustom) throw -3;
147     get1 = 1;
148     file1 = vcustom[i];
149     continue;
150     }
151     // -----------------------------------------------------//
152     if (!strcmp(vcustom[i], "-frame1")){
153     if (++i >= ncustom)throw -3;
154     get1 = 1;
155     frame1 = vcustom[i];
156     continue;
157     }
158     // -----------------------------------------------------//
159     if (!strcmp(vcustom[i], "-frame2")){
160     if (++i >= ncustom)throw -3;
161     get2 = 1;
162     frame2 = vcustom[i];
163     continue;
164     }
165     // -----------------------------------------------------//
166 pam-fi 1.12 if (!strcmp(vcustom[i], "-pfa")){
167     if (++i >= ncustom)throw -3;
168     TrkParams::SetPFA(atoi(vcustom[i]));
169     continue;
170     }
171     // -----------------------------------------------------//
172 pam-fi 1.14 if (!strcmp(vcustom[i], "-skip")){
173     if (++i >= ncustom)throw -3;
174     nskip = atoi(vcustom[i]);
175     continue;
176     }
177     // -----------------------------------------------------//
178 pam-fi 1.11 if ( atoi(vcustom[i]) < 0){
179     if (++i >= ncustom)throw -3;
180     partype[npar]=-1*atoi(vcustom[i-1]);
181     parpath[npar]=vcustom[i];
182     npar++;
183     continue;
184     }
185     // -----------------------------------------------------//
186 pam-fi 1.2 if (!strcmp(vcustom[i], "--verbose") || !strcmp(vcustom[i], "-v")){
187 pam-fi 1.11 TrkParams::SetVerboseMode();
188 pam-fi 1.2 continue;
189     }
190     // -----------------------------------------------------//
191 pam-fi 1.12 if (!strcmp(vcustom[i], "--warning") || !strcmp(vcustom[i], "-w")){
192     TrkParams::SetWarningMode();
193     continue;
194     }
195     // -----------------------------------------------------//
196     if (!strcmp(vcustom[i], "--help") || !strcmp(vcustom[i], "-h")){
197     PrintHelp();
198     return 1;
199     }
200     // -----------------------------------------------------//
201 pam-fi 1.2 else if (!strcmp(vcustom[i], "--debug") || !strcmp(vcustom[i], "-d")){
202 pam-fi 1.11 TrkParams::SetDebugMode();
203 mocchiut 1.1 continue;
204 pam-fi 1.2 };
205 mocchiut 1.1 }
206 pam-fi 1.4
207 mocchiut 1.1 if (!frame2.CompareTo("root", TString::kIgnoreCase)) ifroot2 = true;
208     else if (!frame2.CompareTo("hbook", TString::kIgnoreCase)) ifroot2 =false;
209     else throw -201;
210    
211     if (!frame1.CompareTo("root", TString::kIgnoreCase)) ifroot1 = true;
212     else if (!frame1.CompareTo("hbook", TString::kIgnoreCase)) ifroot1 =false;
213     else throw -201;
214 pam-fi 1.4
215 pam-fi 1.6 if (!frameh.CompareTo("root", TString::kIgnoreCase)) ifrooth = true;
216     else if (!frameh.CompareTo("hbook", TString::kIgnoreCase)) ifrooth =false;
217     else throw -201;
218 mocchiut 1.1
219 pam-fi 1.7 if(get1 && !get2)full1=true;
220    
221 pam-fi 1.2 ostatus = 0;
222 mocchiut 1.1
223 pam-fi 1.12 return 0;
224    
225 pam-fi 1.2 };
226     /**
227     * \brief Process Level0 event
228     * @param l0_event Pointer to an object of the tracker level0 class
229     * Starting from a Level0 object, this routine fills the level0 common, then calls
230     * the fortran routines, which process the event and fill level1 and level2 commons.
231     */
232     void TrkProcess::ProcessEvent(TrkLevel0 *l0_event){
233    
234     // fill Level0 common from Level0 object
235 pam-fi 1.8 l0_event->GetLevel0Struct();
236 pam-fi 1.2 // process the event Level0->Level1->Level2
237     int F77err = 0;
238     //
239     reductionflight_(&F77err);
240     if(F77err < 0)throw F77err;
241     //
242 pam-fi 1.10 if(get2) analysisflight_();
243 pam-fi 1.2 //
244     ostatus = F77err;
245    
246 mocchiut 1.1 };
247    
248     void TrkProcess::Dump(){
249 pam-fi 1.11 cout << endl;
250     cout << " |-------------------------------------| "<<endl;
251     cout << " tracker-processing input-parameters "<< endl;
252     cout << " |-------------------------------------| "<<endl;
253     cout << " idRun " << idrun << endl;
254 pam-fi 1.14 cout << " (skip " << nskip <<" events)"<< endl;
255 pam-fi 1.11 cout << " get level1 " << get1 << endl;
256     cout << " full level1 " << full1 << endl;
257     cout << " get level2 " << get2 << endl;
258     cout << " get hough " << geth << endl;
259     // cout << " frame1 " << frame1 << endl;
260     // cout << " frame2 " << frame2 << endl;
261     // cout << " frameh " << frame2 << endl;
262     // cout << " file1 " << file1 << endl;
263     cout << " output file " << file2 << endl;
264     // cout << " outdir " << outdir << endl;
265     // cout << " process folder " << pfolder << endl;
266     // cout << " standalone mode " << standalone << endl;
267     for(int i=0; i<npar; i++)
268     cout << " force loading parameters of type "<<partype[i]<<" from <--- "<<parpath[i]<<endl;
269     cout<<endl;
270    
271     // cout << " debug mode warning " << dbg_mode.warning << endl;
272     // cout << " debug mode verbose " << dbg_mode.verbose << endl;
273     // cout << " debug mode debug " << dbg_mode.debug << endl << endl;
274 mocchiut 1.1
275     }
276 pam-fi 1.12
277     void TrkProcess::PrintHelp(){
278    
279     printf( "\n\n +TRK [ options ] \n \n");
280     // printf( "\n --version : Print tracker software version and exit ");
281     printf( "Options:");
282     printf( "\n --help, -h : Print this help and exit ");
283     // printf( "\n -idRun RUN : ID number of the run to be processed (for reprocessing RUN=0) \n");
284     // printf( "\n -outDir OUTDIR : Path where to put the LEVEL2 output [default ./ ] ");
285     // printf( "\n -processFolder DIR : Directory (relative to OUTDIR) for other output files [default TrackerFolder/ ] ");
286     // printf( "\n -processFile FILE : Name of the LEVEL2 output file [default RUN.Level2.root]");
287     // printf( "\n -processFile1 FILE : Name of the LEVEL1 output file [default RUN.Level1.rz]");
288     // printf( "\n -frame1 FRAME : type of output for LEVEL1, root/hbook [default hbook ] ");
289     // printf( "\n -frame2 FRAME : type of output for LEVEL2, root/hbook [default root ] ");
290     printf( "\n --get1, --get-level1 : get LEVEL1 output (store only cluster associated to tracks) [default (no LEVEL1 output) ] ");
291     printf( "\n --get1-full, --get-level1-full : get LEVEL1 output (store all clusters) [default (no LEVEL1 output) ] ");
292     printf( "\n --dontget2, --dontget-level2 : do not get LEVEL2 output [default (get LEVEL2 output)]");
293     printf( "\n --geth,, --get-hough : get Hough-transform output [default (no Hough-t. output) ] ");
294     printf( "\n -pfa PFAID : type of output for LEVEL2, root/hbook [default 14 (COG4) ] ");
295     printf( "\n -N PATH : load parameter of type N=1,2... from PATH [default (from DB) ] ");
296     printf( "\n --verbose, -v : verbose mode ");
297     printf( "\n --debug, -d : debug mode ");
298     printf( "\n --warning, -w : warning mode ");
299 pam-fi 1.14 printf( "\n -skip NEV : skip first NEV events from each run ");
300 pam-fi 1.12 // printf( "\n --standalone, -s : standalone mode (without RunInfo) ... but it might not work \n");
301     // printf( "\n -host HOST : Name for the host [default mysql://localhost/pamelaprod ]");
302     // printf( "\n -user USER : Username for the DB [default anonymous] ");
303     // printf( "\n -psw PSW : Password for the DB [default (none)]\n \n \n");
304    
305    
306     }

  ViewVC Help
Powered by ViewVC 1.1.23