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

Annotation of /DarthVader/src/DarthVader.cpp

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.55 - (hide annotations) (download)
Tue Oct 14 14:07:51 2014 UTC (10 years, 1 month ago) by mocchiut
Branch: MAIN
CVS Tags: v10RED, v10REDr01
Changes since 1.54: +4 -3 lines 
10RED: lost sync bug fixed
1 mocchiut 1.1 //
2     // C/C++ headers
3     //
4     #include <iostream>
5     #include <sstream>
6     //
7     // ROOT headers
8     //
9 mocchiut 1.52 #include <TTimeStamp.h>
10 mocchiut 1.1 #include <TString.h>
11     #include <TSQLServer.h>
12     #include <TFile.h>
13     #include <TSystem.h>
14 mocchiut 1.19 #include <TStopwatch.h>
15 mocchiut 1.27 #include <TObjectTable.h>
16     //
17     #include <GLTables.h>
18 mocchiut 1.1 //
19     // Detector's package headers
20     //
21     #include <RunInfoCore.h>
22     #include <TrkCore.h>
23     #include <TrigCore.h>
24     #include <ToFCore.h>
25     #include <CaloCore.h>
26     #include <AcCore.h>
27     #include <S4Core.h>
28     #include <NDCore.h>
29     #include <OrbitalInfoCore.h>
30 mocchiut 1.4 #include <RunInfo.h>
31 mocchiut 1.1 //
32     using namespace std;
33     //
34 mocchiut 1.33 GL_TABLES *glt = NULL;
35 mocchiut 1.1 //
36     //
37     #include <DarthVaderVerl2.h>
38     //
39     // Usage subroutine
40     //
41     void usage(){
42     printf("\nUsage:\n");
43 mocchiut 1.19 printf("\n DarthVader [ options ] -idRun ID_RUN [+-all] [+-detector [ detector options ] ] \n");
44 mocchiut 1.22 printf("\n Options are: \n\n");
45     printf(" --version print informations about compilation and exit\n");
46 mocchiut 1.28 printf(" -h || --help print this help and exit \n");
47     printf(" -v || --verbose be verbose [default]\n");
48     printf(" -s || --silent print nothing on STDOUT\n");
49     printf(" -c || --clean remove file if exiting with errors\n");
50     printf(" -b || --benchmark perform and print a benchmark test\n");
51 mocchiut 1.47 printf(" -r || --reprocess force running on existing file [default: exit with error, NB: DEFAULT CHANGED since 10RED!!]\n");
52     printf(" -n || --new-fit use new fitting algorithm [default]\n");
53 mocchiut 1.53 printf(" --no-new-fit use only standard (up to 9RED) fitting algorithm\n");
54     printf(" --force-TRK force TRK and CAL/TOF simultaneous reprocessing (use only if you know the meaning! could crash!)\n");
55 mocchiut 1.34 printf(" -auto || -AUTO exclude from processing detector which are NOT in the acquisition\n");
56 mocchiut 1.22 printf(" -zerofill if a detector is not in the acquisition the routine is called anyway \n");
57 mocchiut 1.34 printf(" but all detector's data will be marked as bad [default]\n");
58 mocchiut 1.22 printf(" -tedious exit with error if a detector is not in the acquisition and \n");
59     printf(" it has not been excluded from processing\n");
60     printf(" -host name of the DB host [default = $PAM_DBHOST or mysql://localhost/pamelaprod]\n");
61     printf(" -user username for the DB connection [default = $PAM_DBUSER or \"anonymous\"] \n");
62     printf(" -psw password for the DB connection [default = $PAM_DBPSW or \"\"]\n");
63     printf(" -idRun ID_RUN ID number (from the DB) of the run to be processed \n");
64     printf(" -processFile file output filename [default ID_RUN.Level2.root]\n");
65 mocchiut 1.28 printf(" +all || +ALL call all detectors software [default]\n");
66     printf(" -all || -ALL call nothing\n");
67 mocchiut 1.22 printf(" +detector process detector; detector can be: TOF,TRK,CAL,TRG,ORB,S4,ND,AC,RUN\n");
68     printf(" -detector do not process detector (as above)\n");
69     printf(" detector options must be included in square parenthesis with spaces, for example:\n");
70     printf(" +CAL [ --verbose -g ] +TRK [ -v --level1 ] \n");
71     printf("\n Examples: \n");
72     printf(" Standard call:\n DarthVader -idRun 1085 \n");
73     printf(" Process only RunInfo and Tracker (be verbose for tracker):\n DarthVader -idRun 1085 -all +RUN +TRK [ --verbose ] \n");
74     printf(" Process all and be verbose for calorimeter:\n DarthVader -idRun 1085 +CAL [ --verbose ] \n\n");
75 mocchiut 1.46 printf(" Set Tracker to work with simulated data:\n DarthVader -idRun 1085 +TRK [ --simu ] \n\n");
76 mocchiut 1.1 };
77 mocchiut 1.26
78 mocchiut 1.1 //
79     // Here the main
80     //
81     int main(int numinp, char *inps[]){
82     //
83     // Variables booking
84     //
85 mocchiut 1.19 //
86     // benchmarck variables
87     //
88     TStopwatch timer;
89     TStopwatch dvtimer;
90     dvtimer.Start(kTRUE);
91     UInt_t nevents = 0;
92     UInt_t nruns = 0;
93     Double_t runtime = 0.;
94     Double_t cruntime = 0.;
95     Double_t trktime = 0.;
96     Double_t ctrktime = 0.;
97     Double_t caltime = 0.;
98     Double_t ccaltime = 0.;
99 mocchiut 1.47 Double_t caltime1 = 0.;
100     Double_t ccaltime1 = 0.;
101 mocchiut 1.19 Double_t toftime = 0.;
102     Double_t ctoftime = 0.;
103 mocchiut 1.47 Double_t toftime1 = 0.;
104     Double_t ctoftime1 = 0.;
105 mocchiut 1.19 Double_t trgtime = 0.;
106     Double_t ctrgtime = 0.;
107     Double_t actime = 0.;
108     Double_t cactime = 0.;
109     Double_t s4time = 0.;
110     Double_t cs4time = 0.;
111     Double_t ndtime = 0.;
112     Double_t cndtime = 0.;
113     Double_t orbtime = 0.;
114     Double_t corbtime = 0.;
115     Double_t dvtime = 0.;
116     Double_t cdvtime = 0.;
117     //
118 mocchiut 1.1 TString message;
119     int nul = 0;
120     Int_t error = 0;
121     //
122     Int_t CALSGN = 0;
123 mocchiut 1.47 Int_t CALSGN1 = 0;
124 mocchiut 1.1 Int_t TRKSGN = 0;
125     Int_t TRGSGN = 0;
126     Int_t TOFSGN = 0;
127 mocchiut 1.47 Int_t TOFSGN1 = 0;
128 mocchiut 1.1 Int_t RUNSGN = 0;
129     Int_t ORBSGN = 0;
130     Int_t ACSGN = 0;
131     Int_t S4SGN = 0;
132     Int_t NDSGN = 0;
133 mocchiut 1.12 Int_t DVSGN = 0;
134 mocchiut 1.1 //
135 mocchiut 1.27 UInt_t NQRUN = 0;
136     UInt_t NQTRK = 0;
137     UInt_t NQCAL = 0;
138 mocchiut 1.47 UInt_t NQCAL1 = 0;
139 mocchiut 1.27 UInt_t NQTOF = 0;
140 mocchiut 1.47 UInt_t NQTOF1 = 0;
141 mocchiut 1.27 UInt_t NQORB = 0;
142     UInt_t NQTRG = 0;
143     UInt_t NQAC = 0;
144     UInt_t NQND = 0;
145     UInt_t NQS4 = 0;
146     UInt_t NQTOT = 0;
147     //
148 mocchiut 1.36 Bool_t autom = false;
149 mocchiut 1.34 Bool_t zerofill = true;
150 mocchiut 1.17 Bool_t tedious = false;
151 mocchiut 1.15 Bool_t remfile = false;
152 mocchiut 1.1 Bool_t debug = false;
153 mocchiut 1.10 Bool_t beverbose = true;
154 mocchiut 1.1 Bool_t givenid = false;
155 mocchiut 1.19 Bool_t bench = false;
156 mocchiut 1.47 Bool_t reprocess = false;
157     Bool_t newfit = true;
158 mocchiut 1.53 Bool_t forcetrk = false;
159 mocchiut 1.1 Bool_t CAL = true;
160     Bool_t TRK = true;
161     Bool_t TRG = true;
162     Bool_t TOF = true;
163     Bool_t S4 = true;
164     Bool_t ND = true;
165     Bool_t AC = true;
166     Bool_t ORB = true;
167 mocchiut 1.19 Bool_t RUN = true;
168 mocchiut 1.1 //
169     Int_t calargc = 0;
170     char *calargv[50];
171     Int_t trkargc = 0;
172     char *trkargv[50];
173     Int_t tofargc = 0;
174     char *tofargv[50];
175     Int_t trgargc = 0;
176     char *trgargv[50];
177     Int_t orbargc = 0;
178     char *orbargv[50];
179     Int_t runargc = 0;
180     char *runargv[50];
181     Int_t acargc = 0;
182     char *acargv[50];
183     Int_t s4argc = 0;
184     char *s4argv[50];
185     Int_t ndargc = 0;
186     char *ndargv[50];
187     //
188     //
189     //
190 mocchiut 1.8 UInt_t run = 0;
191 mocchiut 1.1 //
192     TString filename;
193     TString outDir = gSystem->WorkingDirectory();
194     //
195     TSQLServer *dbc = 0;
196     TString host = "mysql://localhost/pamelaprod";
197     TString user = "anonymous";
198     TString psw = "";
199 mocchiut 1.27 //
200 mocchiut 1.10 //
201     const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");
202     const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");
203     const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");
204 mocchiut 1.12 if ( !pamdbhost ) pamdbhost = "";
205     if ( !pamdbuser ) pamdbuser = "";
206     if ( !pamdbpsw ) pamdbpsw = "";
207 mocchiut 1.10 if ( strcmp(pamdbhost,"") ) host = pamdbhost;
208     if ( strcmp(pamdbuser,"") ) user = pamdbuser;
209     if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;
210     //
211     //
212 mocchiut 1.1 TFile *processFile = 0;
213     //
214     // Checking input parameters
215     //
216     Int_t i = 0;
217     try {
218     if ( numinp > 1 ){
219     while ( i < numinp ){
220 mocchiut 1.8 Bool_t found = false;
221 mocchiut 1.1 if ( !strcmp(inps[i],"--version") ){
222     DarthVaderInfo(true);
223     exit(0);
224     };
225     if ( !strcmp(inps[i],"-h") || !strcmp(inps[i],"--help") ){
226     usage();
227     exit(0);
228     };
229     if ( !strcmp(inps[i],"-idRun") ) {
230     if ( numinp-1 < i+1 ) {
231     usage();
232     throw -3;
233     };
234 mocchiut 1.8 i++;
235     found = true;
236 mocchiut 1.1 givenid = true;
237 mocchiut 1.8 run = atoll(inps[i]);
238 mocchiut 1.1 };
239     if ( !strcmp(inps[i],"-processFile") ) {
240     if ( numinp-1 < i+1 ){
241     usage();
242     throw -3;
243     };
244 mocchiut 1.8 i++;
245     found = true;
246     filename = (TString)inps[i];
247 mocchiut 1.1 };
248     if ( !strcmp(inps[i],"-outDir") ) {
249     if ( numinp-1 < i+1 ){
250     usage();
251     throw -3;
252     };
253 mocchiut 1.8 i++;
254     found = true;
255     outDir = (TString)inps[i];
256 mocchiut 1.1 };
257     if ( !strcmp(inps[i],"-host") ) {
258     if ( numinp-1 < i+1 ){
259     usage();
260     throw -3;
261     };
262 mocchiut 1.8 i++;
263     found = true;
264     host = (TString)inps[i];
265 mocchiut 1.1 };
266     if ( !strcmp(inps[i],"-user") ) {
267     if ( numinp-1 < i+1 ){
268     usage();
269     throw -3;
270     };
271 mocchiut 1.8 i++;
272     found = true;
273     user = (TString)inps[i];
274 mocchiut 1.1 };
275     if ( !strcmp(inps[i],"-psw") ) {
276     if ( numinp-1 < i+1 ){
277     usage();
278     throw -3;
279     };
280 mocchiut 1.8 i++;
281     found = true;
282     psw = (TString)inps[i];
283     };
284     if ( !strcmp(inps[i],"-v") || !strcmp(inps[i],"--verbose") ){
285     found = true;
286     beverbose = true;
287 mocchiut 1.1 };
288 mocchiut 1.15 if ( !strcmp(inps[i],"-c") || !strcmp(inps[i],"--clean")){
289     found = true;
290     remfile = true;
291     };
292 mocchiut 1.10 if ( !strcmp(inps[i],"-s") || !strcmp(inps[i],"--silent") ){
293     found = true;
294     beverbose = false;
295     };
296 mocchiut 1.8 if ( !strcmp(inps[i],"-g") || !strcmp(inps[i],"--debug") ){
297     debug = true;
298     found = true;
299     };
300 mocchiut 1.19 if ( !strcmp(inps[i],"-b") || !strcmp(inps[i],"--benchmark") ){
301     bench = true;
302     found = true;
303     };
304 mocchiut 1.47 if ( !strcmp(inps[i],"-r") || !strcmp(inps[i],"--reprocess") ){
305     reprocess = true;
306     found = true;
307     };
308     if ( !strcmp(inps[i],"-n") || !strcmp(inps[i],"--new-fit") ){
309     newfit = true;
310     found = true;
311     };
312     if ( !strcmp(inps[i],"--no-new-fit") ){
313     newfit = false;
314     found = true;
315     };
316 mocchiut 1.53 if ( !strcmp(inps[i],"--force-TRK") ){
317     forcetrk = true;
318     found = true;
319     };
320 mocchiut 1.17 if ( !strcmp(inps[i],"-auto") || !strcmp(inps[i],"-AUTO") ){
321     autom = true;
322     tedious = false;
323     zerofill = false;
324     found = true;
325     };
326     if ( !strcmp(inps[i],"-zerofill") ){
327     tedious = false;
328     autom = false;
329     zerofill = true;
330     found = true;
331     };
332     if ( !strcmp(inps[i],"-tedious") ){
333     tedious = true;
334     zerofill = false;
335     autom = false;
336     found = true;
337     };
338 mocchiut 1.1 //
339     if ( !strcmp(inps[i],"-CAL") ) {
340 mocchiut 1.8 found = true;
341 mocchiut 1.1 CAL = false;
342     };
343     if ( !strcmp(inps[i],"-TRK") ) {
344 mocchiut 1.8 found = true;
345 mocchiut 1.1 TRK = false;
346     };
347     if ( !strcmp(inps[i],"-TOF") ) {
348 mocchiut 1.8 found = true;
349 mocchiut 1.1 TOF = false;
350     };
351     if ( !strcmp(inps[i],"-TRG") ) {
352 mocchiut 1.8 found = true;
353 mocchiut 1.1 TRG = false;
354     };
355     if ( !strcmp(inps[i],"-S4") ) {
356 mocchiut 1.8 found = true;
357 mocchiut 1.1 S4 = false;
358     };
359     if ( !strcmp(inps[i],"-ND") ) {
360 mocchiut 1.8 found = true;
361 mocchiut 1.1 ND = false;
362     };
363     if ( !strcmp(inps[i],"-AC") ) {
364 mocchiut 1.8 found = true;
365 mocchiut 1.1 AC = false;
366     };
367     if ( !strcmp(inps[i],"-RUN") ) {
368 mocchiut 1.8 found = true;
369 mocchiut 1.1 RUN = false;
370     };
371     if ( !strcmp(inps[i],"-ORB") ) {
372 mocchiut 1.8 found = true;
373 mocchiut 1.1 ORB = false;
374     };
375     //
376 mocchiut 1.8 if ( !strcmp(inps[i],"-all") || !strcmp(inps[i],"-ALL") ) {
377 mocchiut 1.1 CAL = false;
378     ORB = false;
379     TRK = false;
380     TRG = false;
381     TOF = false;
382     S4 = false;
383     ND = false;
384     AC = false;
385     RUN = false;
386 mocchiut 1.8 found = true;
387 mocchiut 1.1 };
388     //
389 mocchiut 1.8 if ( !strcmp(inps[i],"+all") || !strcmp(inps[i],"+ALL") ) {
390 mocchiut 1.1 CAL = true;
391     ORB = true;
392     TRK = true;
393     TRG = true;
394     TOF = true;
395     S4 = true;
396     ND = true;
397     AC = true;
398     RUN = true;
399 mocchiut 1.8 found = true;
400 mocchiut 1.1 };
401     //
402     if ( !strcmp(inps[i],"+CAL") ) {
403 mocchiut 1.8 found = true;
404 mocchiut 1.1 CAL = true;
405     if ( numinp >= i+2 && !strcmp(inps[i+1],"[") ){
406     if ( numinp < i+2 ){
407     usage();
408     throw -3;
409     };
410     i += 2;
411     calargc = 0;
412     while ( strcmp(inps[i],"]") ){
413     calargv[calargc] = inps[i];
414     calargc++;
415     i++;
416     if ( i > numinp-1 ){
417     usage();
418     throw -3;
419     };
420     };
421     };
422     };
423     if ( !strcmp(inps[i],"+TRK") ) {
424 mocchiut 1.8 found = true;
425 mocchiut 1.1 TRK = true;
426     if ( numinp >= i+2 && !strcmp(inps[i+1],"[") ){
427     if ( numinp-1 < i+2 ){
428     usage();
429     throw -3;
430     };
431     i += 2;
432     trkargc = 0;
433     while ( strcmp(inps[i],"]") ){
434     trkargv[trkargc] = inps[i];
435     trkargc++;
436     i++;
437     if ( i > numinp-1 ){
438     usage();
439     throw -3;
440     };
441     };
442     };
443     };
444     if ( !strcmp(inps[i],"+TOF") ) {
445 mocchiut 1.8 found = true;
446 mocchiut 1.1 TOF = true;
447     if ( numinp >= i+2 && !strcmp(inps[i+1],"[") ){
448     i += 2;
449     tofargc = 0;
450     while ( strcmp(inps[i],"]") ){
451     tofargv[tofargc] = inps[i];
452     tofargc++;
453     i++;
454     if ( i > numinp-1 ){
455     usage();
456     throw -3;
457     };
458     };
459     };
460     };
461     if ( !strcmp(inps[i],"+TRG") ) {
462 mocchiut 1.8 found = true;
463 mocchiut 1.1 TRG = true;
464     if ( numinp >= i+2 && !strcmp(inps[i+1],"[") ){
465     i += 2;
466     trgargc = 0;
467     while ( strcmp(inps[i],"]") ){
468     trgargv[trgargc] = inps[i];
469     trgargc++;
470     i++;
471     if ( i > numinp-1 ){
472     usage();
473     throw -3;
474     };
475     };
476     };
477     };
478     if ( !strcmp(inps[i],"+ORB") ) {
479 mocchiut 1.8 found = true;
480 mocchiut 1.1 ORB = true;
481     if ( numinp >= i+2 && !strcmp(inps[i+1],"[") ){
482     i += 2;
483     orbargc = 0;
484     while ( strcmp(inps[i],"]") ){
485     orbargv[orbargc] = inps[i];
486     orbargc++;
487     i++;
488     if ( i > numinp-1 ){
489     usage();
490     throw -3;
491     };
492     };
493     };
494     };
495     if ( !strcmp(inps[i],"+RUN") ) {
496 mocchiut 1.8 found = true;
497 mocchiut 1.1 RUN = true;
498     if ( numinp >= i+2 && !strcmp(inps[i+1],"[") ){
499     i += 2;
500     runargc = 0;
501     while ( strcmp(inps[i],"]") ){
502     runargv[runargc] = inps[i];
503     runargc++;
504     i++;
505     if ( i > numinp-1 ){
506     usage();
507     throw -3;
508     };
509     };
510     };
511     };
512     if ( !strcmp(inps[i],"+AC") ) {
513 mocchiut 1.8 found = true;
514 mocchiut 1.1 AC = true;
515     if ( numinp >= i+2 && !strcmp(inps[i+1],"[") ){
516     i += 2;
517     acargc = 0;
518     while ( strcmp(inps[i],"]") ){
519     acargv[acargc] = inps[i];
520     acargc++;
521     i++;
522     if ( i > numinp-1 ){
523     usage();
524     throw -3;
525     };
526     };
527     };
528     };
529     if ( !strcmp(inps[i],"+S4") ) {
530 mocchiut 1.8 found = true;
531 mocchiut 1.1 S4 = true;
532     if ( numinp >= i+2 && !strcmp(inps[i+1],"[") ){
533     i += 2;
534     s4argc = 0;
535     while ( strcmp(inps[i],"]") ){
536     s4argv[s4argc] = inps[i];
537     s4argc++;
538     i++;
539     if ( i > numinp-1 ){
540     usage();
541     throw -3;
542     };
543     };
544     };
545     };
546     if ( !strcmp(inps[i],"+ND") ) {
547 mocchiut 1.8 found = true;
548 mocchiut 1.1 ND = true;
549     if ( numinp >= i+2 && !strcmp(inps[i+1],"[") ){
550     i += 2;
551     ndargc = 0;
552     while ( strcmp(inps[i],"]") ){
553     ndargv[ndargc] = inps[i];
554     ndargc++;
555     i++;
556     if ( i > numinp-1 ){
557     usage();
558     throw -3;
559     };
560     };
561     };
562     };
563     //
564 mocchiut 1.8 if ( !found && i > 0 ){
565     usage();
566     printf(" Unknown input number %i, that is \"%s\" \n",i,inps[i]);
567     throw -22;
568     };
569     //
570 mocchiut 1.1 i++;
571     };
572     //
573     } else {
574     //
575     // no input parameters exit with error, we need at least the run id.
576     //
577     throw -1;
578     };
579     //
580     // If not in verbose mode redirect to /dev/null the stdout and stderr
581     //
582     if ( !beverbose ){
583 pam-fi 1.42 nul = open("/dev/null", O_CREAT | O_RDWR,S_IRUSR | S_IWUSR);
584 mocchiut 1.1 dup2(nul,1);
585     dup2(nul,2);
586     };
587     //
588     // Check that an input run number has been given
589     //
590     if ( !givenid ) throw -1;
591     //
592 mocchiut 1.15 TString version = DarthVaderInfo(false);
593 mocchiut 1.1 //
594     // Start:
595     //
596 mocchiut 1.50 printf("\n Welcome to the PAMELA LEVEL2 flight software, version %s (The Ultimate Warrior)\n\n",version.Data());
597 mocchiut 1.51 if ( run ){
598     printf("\n Processing run number %u \n",run);
599     } else {
600     printf("\n Re-processing all runs?\n");
601     }
602 mocchiut 1.1 //
603 mocchiut 1.53 // check tracker reprocessing
604     //
605     if ( TRK && reprocess ){
606     printf(" WARNING: Tracker reprocessing can be incomplete and/or not working properly! \n");
607     if ( ( CAL && TRK && reprocess ) || ( TOF && TRK && reprocess ) ){
608     printf(" ERROR: cannot reprocess Tracker and Calorimeter or Tracker and ToF simultaneously! \n");
609     printf(" If CaloLevel1 and ToF level0 remain unchanged: reprocess first tracker only and the CAL and/or TOF \n");
610     printf(" If CaloLevel1 and ToF level0 could change: reprocess first CAL and/or TOF, then reprocess TRK only and finally reprocess CAL and/or TOF \n");
611     printf(" IF you are crazy you can override this error using flag: --force-TRK (i.e. DarthVader -v -b -idRun 0 -processFile myfile.root --force-TRK -ALL +TRK ) \n");
612     }
613     if ( !forcetrk ) throw -38;
614     }
615    
616     //
617 mocchiut 1.1 // Connect to the DB
618     //
619     if ( debug ) printf("\nConnecting to database... \n");
620     //
621     dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());
622     if( !dbc ) throw -2;
623     //
624 mocchiut 1.18 Bool_t connect = dbc->IsConnected();
625 mocchiut 1.1 //
626     if( !connect ) throw -2;
627     //
628     if ( debug ) printf("...connected! \n\n");
629     //
630 mocchiut 1.33 glt = new GL_TABLES(host,user,psw);
631     //GL_TABLES *glt = new GL_TABLES(host,user,psw);
632 mocchiut 1.27 //
633 mocchiut 1.51 printf("\n DB settings:\n SQL: %s Version: %s Host %s Port %i \n",dbc->GetDBMS(),dbc->ServerInfo(),dbc->GetHost(),dbc->GetPort());
634     printf(" DB %s --- User %s \n\n",host.Data(),user.Data());
635 mocchiut 1.19 if ( debug ) printf("\n DB INFORMATIONS:\n SQL: %s Version: %s Host %s Port %i \n\n",dbc->GetDBMS(),dbc->ServerInfo(),dbc->GetHost(),dbc->GetPort());
636     //
637 mocchiut 1.15 // Use UTC in the DB
638     //
639 mocchiut 1.14 stringstream myquery;
640     myquery.str("");
641 mocchiut 1.55 myquery << "SET time_zone='+0:00';";
642     delete dbc->Query(myquery.str().c_str());
643 mocchiut 1.27 myquery.str("");
644     myquery << "SET wait_timeout=173000;";
645 mocchiut 1.55 delete dbc->Query(myquery.str().c_str());
646     delete dbc->Query("SET sql_mode = 'NO_UNSIGNED_SUBTRACTION';");
647 mocchiut 1.14 //
648     //
649 mocchiut 1.1 // Create LEVEL2 filename and open it in update mode
650     //
651     if ( filename.IsNull() ){
652     stringstream strun;
653     strun.str("");
654     strun << run;
655 mocchiut 1.3 filename += outDir;
656     filename += "/";
657 mocchiut 1.1 filename += strun.str();
658     filename += ".Level2.root";
659     };
660 mocchiut 1.47 //
661     // check if file exists and check reprocess flag
662     //
663     if ( !reprocess ){
664     if ( !gSystem->GetPathInfo(filename.Data(),NULL,(Long_t*)NULL,NULL,NULL) ) {
665     if ( beverbose ) printf(" File %s already exists!\n",filename.Data());
666     throw -37;
667     }
668     }
669    
670 mocchiut 1.1 processFile = new TFile(filename.Data(),"UPDATE");
671     if ( !processFile->IsOpen() ) throw -15;
672     //
673 mocchiut 1.23 Long64_t maxsize = 10000000000LL;
674     TTree::SetMaxTreeSize(maxsize);
675 mocchiut 1.24 processFile->SetCompressionLevel(1);
676 mocchiut 1.23 //
677 mocchiut 1.52 // ok, we are going to do some kinf of processing, the file is opened and ready. Save processing infos first
678     //
679     ProcInfo *procinfo = new ProcInfo();
680     procinfo->runID = run;
681     TTimeStamp *dt = new TTimeStamp();
682     procinfo->date = dt->AsString();
683     delete dt;
684     for ( Int_t icl = 0; icl<numinp; icl++){
685     procinfo->commandLine += Form("%s ",inps[icl]);
686     }
687     procinfo->outputFilename = filename;
688     procinfo->localDir = gSystem->WorkingDirectory();
689     procinfo->uname = gSystem->GetFromPipe("uname -a");
690     procinfo->DB = host;
691    
692     TTree *pinfo = 0;
693     pinfo = (TTree*)processFile->Get("ProcessingInfo");
694     // ProcessingInfo tree does not exist, crating proc info
695     if ( !pinfo ){
696     if ( debug ) printf("ProcessingInfo tree does not exist, crating proc info\n");
697     pinfo = new TTree("ProcessingInfo","Log of data processing");
698     pinfo->Branch("ProcInfo","ProcInfo",&procinfo);
699     } else {
700     pinfo->SetBranchAddress("ProcInfo",&procinfo);
701     }
702     pinfo->Fill();
703     processFile->cd();
704     pinfo->Write("ProcessingInfo",TObject::kOverwrite);
705     if ( procinfo ) delete procinfo;
706     if ( pinfo ) pinfo->Delete();
707 mocchiut 1.23 //
708 mocchiut 1.1 // Run the core program, put any output error in the "error" variable
709     //
710 mocchiut 1.19 if ( debug ) printf("\n\n Pre-processing:\n\n");
711 mocchiut 1.15 //
712 mocchiut 1.19 timer.Start(kTRUE);
713 mocchiut 1.1 if ( RUN ) {
714 mocchiut 1.27 glt->ResetCounters();
715 mocchiut 1.19 if ( debug ) printf(" Retrieve, if the case, the RUN informations from the DB...\n");
716 mocchiut 1.15 printf(" RunInfo called\n");
717     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start RunInfoCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
718     RUNSGN = RunInfoCore(run,processFile,dbc,version,runargc,runargv);
719     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end RunInfoCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
720 mocchiut 1.19 if ( debug ) printf(" ...done\n");
721 mocchiut 1.27 NQRUN = glt->GetNqueries();
722 mocchiut 1.1 };
723 mocchiut 1.19 timer.Stop();
724     runtime = timer.RealTime();
725     cruntime = timer.CpuTime();
726 mocchiut 1.4 //
727     // From the Run Infos extract acq_var_info to determine if detectors are in the acquisition or not
728     //
729 mocchiut 1.19 if ( debug ) printf(" Checking if requested detectors are in the acquisition\n");
730 mocchiut 1.4 ItoRunInfo *runinfo = new ItoRunInfo(processFile);
731 mocchiut 1.15 Int_t sgnl = runinfo->Read(run);
732     if ( sgnl != 0 ) throw -17;
733 mocchiut 1.19 nevents = runinfo->GetEntries();
734 mocchiut 1.53 if ( run == 0 ) nevents = runinfo->GetFileEntries();
735 mocchiut 1.19 nruns = runinfo->GetNoRun();
736     if ( debug ) printf(" => ACQ_VAR_INFO = %i \n",runinfo->ACQ_VAR_INFO);
737 mocchiut 1.17 if ( TRK && !(runinfo->ACQ_VAR_INFO & (1 << 4)) ){
738     DVSGN += 1;
739     if ( tedious ) throw -30;
740     if ( autom ){
741     TRK = 0;
742 mocchiut 1.19 if ( debug ) printf(" TRK excluded from processing \n");
743 mocchiut 1.17 };
744 mocchiut 1.19 if ( zerofill && debug ) printf(" TRK is not in the acquisition! \n");
745 mocchiut 1.17 };
746     if ( TOF && !(runinfo->ACQ_VAR_INFO & (1 << 0)) ){
747     DVSGN += 2;
748     if ( tedious ) throw -31;
749     if ( autom ){
750     TOF = 0;
751 mocchiut 1.19 if ( debug ) printf(" TOF excluded from processing \n");
752 mocchiut 1.17 };
753 mocchiut 1.19 if ( zerofill && debug ) printf(" TOF is not in the acquisition! \n");
754 mocchiut 1.17 };
755     if ( CAL && !(runinfo->ACQ_VAR_INFO & (1 << 3)) ){
756     DVSGN += 4;
757     if ( tedious ) throw -32;
758     if ( autom ){
759     CAL = 0;
760 mocchiut 1.19 if ( debug ) printf(" CAL excluded from processing \n");
761 mocchiut 1.17 };
762 mocchiut 1.19 if ( zerofill && debug ) printf(" CAL is not in the acquisition! \n");
763 mocchiut 1.17 };
764     if ( AC && (!(runinfo->ACQ_VAR_INFO & (1 << 1)) || !(runinfo->ACQ_VAR_INFO & (1 << 2))) ){
765     DVSGN += 8;
766     if ( tedious ) throw -33;
767     if ( autom ){
768     AC = 0;
769 mocchiut 1.19 if ( debug ) printf(" AC excluded from processing \n");
770 mocchiut 1.17 };
771 mocchiut 1.19 if ( zerofill && debug ) printf(" AC is not in the acquisition! \n");
772 mocchiut 1.17 };
773     if ( S4 && !(runinfo->ACQ_VAR_INFO & (1 << 5)) ){
774     DVSGN += 16;
775     if ( tedious ) throw -34;
776     if ( autom ){
777     S4 = 0;
778 mocchiut 1.19 if ( debug ) printf(" S4 excluded from processing \n");
779 mocchiut 1.17 };
780 mocchiut 1.19 if ( zerofill && debug ) printf(" S4 is not in the acquisition! \n");
781 mocchiut 1.17 };
782     if ( ND && !(runinfo->ACQ_VAR_INFO & (1 << 9)) ){
783     DVSGN += 32;
784     if ( tedious ) throw -35;
785     if ( autom ){
786     ND = 0;
787 mocchiut 1.19 if ( debug ) printf(" ND excluded from processing \n");
788 mocchiut 1.17 };
789 mocchiut 1.19 if ( zerofill && debug ) printf(" ND is not in the acquisition! \n");
790 mocchiut 1.17 };
791 mocchiut 1.15 //
792     if ( !DVSGN ){
793 mocchiut 1.19 if ( debug ) printf(" OK! Start processing detector's data. \n");
794 mocchiut 1.15 } else {
795 mocchiut 1.19 if ( debug ) printf(" WARNING! missing detector(s)! Start anyway processing detector's data. \n");
796 mocchiut 1.15 };
797 mocchiut 1.19 if ( debug ) printf("\n End pre-processing \n\n");
798 mocchiut 1.4 //
799 mocchiut 1.33 dbc->Close();
800     //
801 mocchiut 1.43 //
802     timer.Start(kTRUE);
803     if ( TRG ) {
804     glt->ResetCounters();
805     printf(" TriggerLevel2 called\n");
806     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start TrigCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
807     TRGSGN = TrigCore(run,processFile,glt,trgargc,trgargv);
808     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end TrigCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
809     NQTRG = glt->GetNqueries();
810     };
811     timer.Stop();
812     trgtime = timer.RealTime();
813     ctrgtime = timer.CpuTime();
814     //
815 mocchiut 1.47 //
816     // new fitting algorithm require calorimeter level1 and tof XXX to be processed first
817     //
818     if ( newfit ) {
819     timer.Start(kTRUE);
820 mocchiut 1.53 if ( TOF && TRK ) {
821 mocchiut 1.47 glt->ResetCounters();
822     printf(" ToFLevel2 called, 1st call\n");
823     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
824 mocchiut 1.48 string tstring = "--level1-only";
825     tofargv[tofargc] = (char *)tstring.c_str();
826     tofargc++;
827 mocchiut 1.47 TOFSGN1 = ToFCore(run,processFile,glt,tofargc,tofargv);
828     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
829     NQTOF1 = glt->GetNqueries();
830     };
831     timer.Stop();
832     toftime1 = timer.RealTime();
833     ctoftime1 = timer.CpuTime();
834     //
835     timer.Start(kTRUE);
836 mocchiut 1.53 if ( CAL && TRK ) {
837 mocchiut 1.47 glt->ResetCounters();
838     printf(" CalorimeterLevel2 called, 1st call\n");
839     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
840     string cstring = "--level1-only";
841     calargv[calargc] = (char *)cstring.c_str();
842     calargc++;
843     CALSGN1 = CaloCore(run,processFile,glt,calargc,calargv);
844     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
845     NQCAL1 = glt->GetNqueries();
846     };
847     timer.Stop();
848     caltime1 = timer.RealTime();
849     ccaltime1 = timer.CpuTime();
850     }
851     //
852     //
853 mocchiut 1.19 timer.Start(kTRUE);
854 mocchiut 1.1 if ( TRK ) {
855 mocchiut 1.27 glt->ResetCounters();
856 mocchiut 1.15 printf(" TrackerLevel2 called\n");
857     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start TrkCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
858 mocchiut 1.33 TRKSGN = TrkCore(run,processFile,glt,trkargc,trkargv);
859 mocchiut 1.11 gSystem->Unlink("TrackerFolder"); //patch
860 mocchiut 1.15 if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end TrkCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
861 mocchiut 1.27 NQTRK = glt->GetNqueries();
862 mocchiut 1.1 };
863 mocchiut 1.19 timer.Stop();
864     trktime = timer.RealTime();
865     ctrktime = timer.CpuTime();
866     //
867     timer.Start(kTRUE);
868 mocchiut 1.1 if ( TOF ) {
869 mocchiut 1.27 glt->ResetCounters();
870 mocchiut 1.15 printf(" ToFLevel2 called\n");
871     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
872 mocchiut 1.48 if ( !reprocess ){
873     string tstring = "--delete-tree";
874     tofargv[tofargc] = (char *)tstring.c_str();
875     } else {
876 mocchiut 1.53 if ( TRK ) tofargc--;
877 mocchiut 1.48 }
878 mocchiut 1.33 TOFSGN = ToFCore(run,processFile,glt,tofargc,tofargv);
879 mocchiut 1.15 if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
880 mocchiut 1.27 NQTOF = glt->GetNqueries();
881 mocchiut 1.1 };
882 mocchiut 1.19 timer.Stop();
883     toftime = timer.RealTime();
884     ctoftime = timer.CpuTime();
885     //
886     timer.Start(kTRUE);
887 mocchiut 1.1 if ( CAL ) {
888 mocchiut 1.27 glt->ResetCounters();
889 mocchiut 1.15 printf(" CalorimeterLevel2 called\n");
890     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
891 mocchiut 1.47 if ( !reprocess ){
892     string cstring = "--delete-calo-tree";
893     calargv[calargc] = (char *)cstring.c_str();
894     } else {
895 mocchiut 1.53 if ( TRK ) calargc--;
896 mocchiut 1.47 }
897 mocchiut 1.33 CALSGN = CaloCore(run,processFile,glt,calargc,calargv);
898 mocchiut 1.15 if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
899 mocchiut 1.27 NQCAL = glt->GetNqueries();
900 mocchiut 1.1 };
901 mocchiut 1.19 timer.Stop();
902     caltime = timer.RealTime();
903     ccaltime = timer.CpuTime();
904     //
905     timer.Start(kTRUE);
906 mocchiut 1.1 if ( AC ) {
907 mocchiut 1.27 glt->ResetCounters();
908 mocchiut 1.15 printf(" AnticounterLevel2 called\n");
909     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start AcCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
910 mocchiut 1.33 ACSGN = AcCore(run,processFile,glt,acargc,acargv);
911 mocchiut 1.15 if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end AcCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
912 mocchiut 1.27 NQAC = glt->GetNqueries();
913 mocchiut 1.1 };
914 mocchiut 1.19 timer.Stop();
915     actime = timer.RealTime();
916     cactime = timer.CpuTime();
917     //
918     timer.Start(kTRUE);
919 mocchiut 1.1 if ( S4 ) {
920 mocchiut 1.27 glt->ResetCounters();
921 mocchiut 1.15 printf(" S4Level2 called\n");
922     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start S4Core STDOUT |<<<<<<<<<<<<<<<<<<<\n");
923 mocchiut 1.33 S4SGN = S4Core(run,processFile,glt,s4argc,s4argv);
924 mocchiut 1.15 if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end S4Core STDOUT |<<<<<<<<<<<<<<<<<<<\n");
925 mocchiut 1.27 NQS4 = glt->GetNqueries();
926 mocchiut 1.1 };
927 mocchiut 1.19 timer.Stop();
928     s4time = timer.RealTime();
929     cs4time = timer.CpuTime();
930     //
931     timer.Start(kTRUE);
932 mocchiut 1.1 if ( ND ) {
933 mocchiut 1.27 glt->ResetCounters();
934 mocchiut 1.15 printf(" NDLevel2 called\n");
935     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start NDCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
936 mocchiut 1.33 NDSGN = NDCore(run,processFile,glt,ndargc,ndargv);
937 mocchiut 1.15 if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end NDCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
938 mocchiut 1.27 NQND = glt->GetNqueries();
939 mocchiut 1.1 };
940 mocchiut 1.19 timer.Stop();
941     ndtime = timer.RealTime();
942     cndtime = timer.CpuTime();
943     //
944     timer.Start(kTRUE);
945 mocchiut 1.1 if ( ORB ) {
946 mocchiut 1.27 glt->ResetCounters();
947 mocchiut 1.15 printf(" OrbitalInfo called\n");
948     if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start OrbitalInfoCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
949 mocchiut 1.33 ORBSGN = OrbitalInfoCore(run,processFile,glt,orbargc,orbargv);
950 mocchiut 1.15 if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end OrbitalInfoCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
951 mocchiut 1.27 NQORB = glt->GetNqueries();
952 mocchiut 1.1 };
953 mocchiut 1.19 timer.Stop();
954     orbtime = timer.RealTime();
955     corbtime = timer.CpuTime();
956 mocchiut 1.1 //
957 mocchiut 1.47 NQTOT = 2 + NQRUN + NQTRK + NQCAL + NQND + NQS4 + NQAC + NQTRG + NQTOF + NQORB + NQCAL1 + NQTOF1;
958 mocchiut 1.27 //
959     delete glt;
960     delete runinfo;
961     //
962 mocchiut 1.1 } catch(Int_t signal) {
963     error = signal;
964     switch(signal){
965     case 51: message += " GLTABLES - TO_TIME < time when querying tables"; break;
966     case -1: message += " Missing/wrong run ID input parameter"; break;
967     case -2: message += " DB connection failure"; break;
968     case -3: message += " Error in input parameters (check format)"; break;
969     case -4: message += " Request reprocessing of all runs (idRun = 0) but processFile is missing"; break;
970 mocchiut 1.38 case -5: message += " RUNINFO - ERROR: no run with this ID_RUN (ID_RUN mismatch?) "; break;
971 mocchiut 1.1 case -6: message += " No LEVEL0 file "; break;
972     case -7: message += " No Physics tree in LEVEL0 file"; break;
973     case -8: message += " No Header branch in LEVEL0 Physics tree"; break;
974     case -9: message += " No Registry branch in LEVEL0 Physics tree"; break;
975     case -11: message += " LEVEL0 Physics tree is empty"; break;
976 mocchiut 1.9 case -12: message += " Too few entries in the tree"; break;
977 mocchiut 1.1 case -13: message += " Cannot create processFolder directory"; break;
978     case -14: message += " Error querying the DB"; break;
979     case -15: message += " Cannot open file for writing"; break;
980 mocchiut 1.15 case -17: message += " Error during pre-processing"; break;
981 mocchiut 1.8 case -22: message += " Unknown input or wrong syntax in input paramters!"; break;
982 mocchiut 1.17 case -30: message += " No TRK in the acquisition"; break;
983     case -31: message += " No TOF in the acquisition"; break;
984     case -32: message += " No CAL in the acquisition"; break;
985     case -33: message += " No AC in the acquisition"; break;
986     case -34: message += " No S4 in the acquisition"; break;
987     case -35: message += " No ND in the acquisition"; break;
988 mocchiut 1.41 case -36: message += " I/O error or missing entry"; break;
989 mocchiut 1.47 case -37: message += " Cannot override output file, see help and use -r option"; break;
990 mocchiut 1.53 case -38: message += " Cannot process TRK and CAL or TOF simultaneously"; break;
991 mocchiut 1.41
992 mocchiut 1.1 //
993     case -50: message += " GLTABLES - No entries matching GL_RUN query"; break;
994     case -51: message += " GLTABLES - No entries matching GL_ROOT query"; break;
995     case -52: message += " GLTABLES - No entries matching GL_PARAM query"; break;
996     case -53: message += " GLTABLES - No entries matching GL_TRK_CALIB query"; break;
997     case -54: message += " GLTABLES - No entries matching GL_CALO_CALIB query"; break;
998 mocchiut 1.30 case -55: message += " GLTABLES - No entries matching GL_S4_CALIB query"; break;
999 mocchiut 1.25 case -56: message += " GLTABLES - No entries matching GL_TLE query"; break;
1000 mocchiut 1.27 case -57: message += " GLTABLES - DB connection gone and not able to reconnect"; break;
1001 mocchiut 1.1 //
1002     case -100: message += " CALORIMETERLEVEL2 - No Level2 input file"; break;
1003     case -101: message += " CALORIMETERLEVEL2 - Cannot open Level2 file"; break;
1004     case -102: message += " CALORIMETERLEVEL2 - No Tracker TTree in Level2 file"; break;
1005     case -103: message += " CALORIMETERLEVEL2 - No Calorimeter TBranch in Level0 file"; break;
1006     case -104: message += " CALORIMETERLEVEL2 - No Trigger TTree in Level0 file"; break;
1007     case -105: message += " CALORIMETERLEVEL2 - No Calorimeter ADC2MIP conversion file"; break;
1008     case -106: message += " CALORIMETERLEVEL2 - No Calorimeter alignement file"; break;
1009     case -107: message += " CALORIMETERLEVEL2 - Cannot find Level0 file needed for the calibration"; break;
1010     case -108: message += " CALORIMETERLEVEL2 - Cannot open Level0 file needed for the calibration"; break;
1011     case -109: message += " CALORIMETERLEVEL2 - No CalibCalPed TTree in Level0 file needed for the calibration"; break;
1012     case -110: message += " CALORIMETERLEVEL2 - No calibrations in Level0 file needed for the calibration"; break;
1013     case -111: message += " CALORIMETERLEVEL2 - Corrupted calibration"; break;
1014     case -112: message += " CALORIMETERLEVEL2 - No physics event related to registry entry in Level0 file"; break;
1015     case -113: message += " CALORIMETERLEVEL2 - No tracker event related to registry entry in Level2 file"; break;
1016 mocchiut 1.17 case -114: message += " CALORIMETERLEVEL2 - Help called"; break;
1017     case -115: message += " CALORIMETERLEVEL2 - No Calorimeter bad strip offline mask file"; break;
1018 mocchiut 1.18 case -116: message += " CALORIMETERLEVEL2 - DB connection problems"; break;
1019 mocchiut 1.20 case -117: message += " CALORIMETERLEVEL2 - Cannot reprocess file with Level1 data without calling Level1 routine"; break;
1020     case -118: message += " CALORIMETERLEVEL2 - Cannot reprocess file without Level1 data calling Level1 routine"; break;
1021 mocchiut 1.35 case -119: message += " CALORIMETERLEVEL2 - No CalibCalPulse2 TTree in Level0 file needed for the calibration"; break;
1022     case -120: message += " CALORIMETERLEVEL2 - No CalibCalPulse1 TTree in Level0 file needed for the calibration"; break;
1023 mocchiut 1.37 case -121: message += " CALORIMETERLEVEL2 - Cannot open calorimeter neighbour crosstalk correction table file"; break;
1024     case -122: message += " CALORIMETERLEVEL2 - Cannot open calorimeter second neighbour crosstalk correction table file"; break;
1025     case -123: message += " CALORIMETERLEVEL2 - Cannot open calorimeter special calibration file"; break;
1026     case -124: message += " CALORIMETERLEVEL2 - Cannot open calorimeter max rms file"; break;
1027     case -125: message += " CALORIMETERLEVEL2 - Cannot open calorimeter silicon crosstalk correction table file"; break;
1028 mocchiut 1.1 //
1029     case -200: message += " TRACKERLEVEL2 - LEVEL1 framework unknown (HBOOK/ROOT)"; break;
1030     case -201: message += " TRACKERLEVEL2 - LEVEL2 framework unknown (HBOOK/ROOT)"; break;
1031     case -202: message += " TRACKERLEVEL2 - Neither LEVEL1 nor LEVEL2 output requested"; break;
1032     case -203: message += " TRACKERLEVEL2 - No Tracker branch in LEVEL0 Physics tree"; break;
1033     case -204: message += " TRACKERLEVEL2 - No reprocessing implemented for LEVEL1 output"; break;
1034     case -205: message += " TRACKERLEVEL2 - Error accessing RunInfo "; break;
1035 pam-fi 1.49 case -206: message += " TRACKERLEVEL2 - Extended-tracking mode - no Calorimeter and/or ToF trees found"; break;
1036 mocchiut 1.1 case -210: message += " TRACKERLEVEL2 - Error opening/reading trk mask GL_PARAM parameters "; break;
1037     case -211: message += " TRACKERLEVEL2 - Error opening/reading trk alignment GL_PARAM parameters"; break;
1038     case -212: message += " TRACKERLEVEL2 - Error opening/reading trk mip GL_PARAM parameters"; break;
1039     case -213: message += " TRACKERLEVEL2 - Error opening/reading trk charge GL_PARAM parameters"; break;
1040     case -214: message += " TRACKERLEVEL2 - Error opening/reading trk pfa GL_PARAM parameters"; break;
1041     case -215: message += " TRACKERLEVEL2 - Error opening/reading field GL_PARAM parameters"; break;
1042     case -216: message += " TRACKERLEVEL2 - Error opening/reading default calibration GL_PARAM parameters"; break;
1043     case -298: message += " TRACKERLEVEL2 - Reprocessing not implemented in standalone mode"; break;
1044     case -299: message += " TRACKERLEVEL2 - Not yet implemented"; break;
1045     //
1046     case -300: message += " TOFLEVEL2 - No Trigger branch in Level0 tree"; break;
1047     case -301: message += " TOFLEVEL2 - Cannot open file for writing"; break;
1048     case -302: message += " TOFLEVEL2 - No tracker tree in Level2 file"; break;
1049     case -303: message += " TOFLEVEL2 - No Tof branch in Level0 file"; break;
1050 mocchiut 1.43 case -304: message += " TOFLEVEL2 - No trigger tree in Level2 file"; break;
1051 mocchiut 1.1 case -313: message += " TOFLEVEL2 - No more tracker events in Level2 file"; break;
1052 mocchiut 1.18 case -314: message += " TOFLEVEL2 - DB connection problems"; break;
1053 mocchiut 1.40 case -315: message += " TOFLEVEL2 - Problems with dE/dx II order correction file parameter"; break;
1054     case -316: message += " TOFLEVEL2 - Problems with dE/dx II order correction, dividing by zero!"; break;
1055     case -317: message += " TOFLEVEL2 - Problems with dE/dx II order correction, outside time limits!"; break;
1056     case -318: message += " TOFLEVEL2 - Problems with dE/dx II order correction, too many time intervals!"; break;
1057 mocchiut 1.43 case -319: message += " TOFLEVEL2 - No more trigger events in Level2 file"; break;
1058 mocchiut 1.1 //
1059     case -401: message += " TRIGGERLEVEL2 - Cannot open file for writing"; break;
1060     case -402: message += " TRIGGERLEVEL2 - No Trigger branch in Level0 tree"; break;
1061 mocchiut 1.18 case -403: message += " TRIGGERLEVEL2 - DB connection problems"; break;
1062 mocchiut 1.1 //
1063     case -500: message += " S4LEVEL2 - No level2 file"; break;
1064     case -501: message += " S4LEVEL2 - Cannot open file for writing"; break;
1065     case -502: message += " S4LEVEL2 - No result from GL_S4_CALIB"; break;
1066     case -503: message += " S4LEVEL2 - No S4 branch in Level0 tree"; break;
1067 mocchiut 1.18 case -504: message += " S4LEVEL2 - DB connection problems"; break;
1068 mocchiut 1.1 //
1069     case -600: message += " NDLEVEL2 - No level2 file"; break;
1070     case -601: message += " NDLEVEL2 - Cannot open file for writing"; break;
1071     case -603: message += " NDLEVEL2 - No S4Level2 branch in Level0 tree"; break;
1072 mocchiut 1.18 case -604: message += " NDLEVEL2 - DB connection problems"; break;
1073 mocchiut 1.1 //
1074 mocchiut 1.17 case -701: message += " ACLEVEL2 - Cannot open file for writing"; break;
1075     case -704: message += " ACLEVEL2 - No Anticounter branch in Level0 tree"; break;
1076 mocchiut 1.18 case -705: message += " ACLEVEL2 - DB connection problems"; break;
1077 mocchiut 1.1 //
1078     case -800: message += " RUNINFO - No such run in the RunInfo list"; break;
1079     case -801: message += " RUNINFO - No RunInfo tree in Level2 file"; break;
1080     case -802: message += " RUNINFO - Cannot open file for writing"; break;
1081     case -803: message += " RUNINFO - Updating but no RunInfo tree in Level2 file"; break;
1082     case -804: message += " RUNINFO - Unknown detector"; break;
1083 mocchiut 1.2 case -805: message += " RUNINFO - Reprocessing data but no RunInfo tree in Level2 file"; break;
1084 mocchiut 1.7 case -806: message += " RUNINFO - Can not handle more than 500 runs"; break;
1085 mocchiut 1.18 case -807: message += " RUNINFO - DB connection problems"; break;
1086 mocchiut 1.39 //
1087     case -900: message += " OrbitalInfo - no ToF tree"; break;
1088 mocchiut 1.50 case -901: message += " OrbitalInfo - Mismatch between tracker and tof tree branches concerning extended tracking algorithm(s)"; break;
1089 mocchiut 1.18 //
1090 mocchiut 1.1 default: message += "Unidentified error or warning"; break;
1091     };
1092     printf("\n");
1093     if ( signal < 0 ) cout << " ERROR ("<< signal << ") "<< message <<endl;
1094     }
1095 mocchiut 1.12 //
1096 mocchiut 1.15 // Warnings from XCore routines and from DV:
1097     //
1098     if ( debug && DVSGN ) printf(" DVSGN = %i \n",DVSGN);
1099     if ( DVSGN ) printf("\n WARNING DarthVader: \n");
1100     if ( DVSGN & (1 << 0) ) printf(" - No tracker in the acquisition\n");
1101     if ( DVSGN & (1 << 1) ) printf(" - No ToF in the acquisition\n");
1102     if ( DVSGN & (1 << 2) ) printf(" - No calorimeter in the acquisition\n");
1103     if ( DVSGN & (1 << 3) ) printf(" - No anticounters in the acquisition\n");
1104     if ( DVSGN & (1 << 4) ) printf(" - No S4 in the acquisition\n");
1105     if ( DVSGN & (1 << 5) ) printf(" - No neutron detector in the acquisition\n");
1106     //
1107 mocchiut 1.47 // CaloCore, first call:
1108     //
1109     switch(CALSGN1){
1110     case 100: printf("\n WARNING CALORIMETER - Data with no associated calibration\n");
1111     case 101: printf("\n WARNING CALORIMETER - No tracks or good events in this run\n");
1112     }
1113     //
1114 mocchiut 1.15 // CaloCore:
1115 mocchiut 1.12 //
1116 mocchiut 1.1 switch(CALSGN){
1117 mocchiut 1.13 case 100: printf("\n WARNING CALORIMETER - Data with no associated calibration\n");
1118     case 101: printf("\n WARNING CALORIMETER - No tracks or good events in this run\n");
1119 mocchiut 1.1 };
1120     //
1121 mocchiut 1.29 // OrbitalInfoCore:
1122     //
1123     switch(ORBSGN){
1124     case 900: printf("\n WARNING ORBITALINFO - No inclination MCMDs\n");
1125     };
1126     //
1127 mocchiut 1.15 // no other at the moment
1128     //
1129     //
1130 mocchiut 1.1 // Close the DB connection
1131     //
1132     if ( dbc ){
1133     if ( debug ) printf("\nClose the connection to the database... \n");
1134     dbc->Close();
1135     if ( debug ) printf("...connection terminated!\n\n");
1136     };
1137     if ( processFile ){
1138     processFile->cd();
1139 mocchiut 1.24 processFile->WriteStreamerInfo();
1140 mocchiut 1.45 if ( debug ) processFile->ls();
1141 mocchiut 1.1 processFile->Close();
1142     };
1143     //
1144 mocchiut 1.19 // print benchmark results (if the case)
1145     //
1146     dvtimer.Stop();
1147     dvtime = dvtimer.RealTime();
1148     cdvtime = dvtimer.CpuTime();
1149     if ( bench ){
1150     Float_t runrt = 0.;
1151     Float_t trkrt = 0.;
1152     Float_t calrt = 0.;
1153     Float_t tofrt = 0.;
1154 mocchiut 1.47 Float_t calrt1 = 0.;
1155     Float_t tofrt1 = 0.;
1156 mocchiut 1.19 Float_t trgrt = 0.;
1157     Float_t acrt = 0.;
1158     Float_t s4rt = 0.;
1159     Float_t ndrt = 0.;
1160     Float_t orbrt = 0.;
1161 mocchiut 1.31 Float_t dvrt = (nevents+1)/dvtime;
1162 mocchiut 1.19 if ( RUN ) runrt = nruns/runtime;
1163 mocchiut 1.31 if ( TRK ) trkrt = (nevents+1)/trktime;
1164     if ( CAL ) calrt = (nevents+1)/caltime;
1165 mocchiut 1.53 if ( TOF && TRK && newfit ) tofrt1 = (nevents+1)/toftime1;
1166     if ( CAL && TRK && newfit ) calrt1 = (nevents+1)/caltime1;
1167 mocchiut 1.31 if ( TOF ) tofrt = (nevents+1)/toftime;
1168     if ( TRG ) trgrt = (nevents+1)/trgtime;
1169     if ( AC ) acrt = (nevents+1)/actime;
1170     if ( S4 ) s4rt = (nevents+1)/s4time;
1171     if ( ND ) ndrt = (nevents+1)/ndtime;
1172     if ( ORB ) orbrt = (nevents+1)/orbtime;
1173 mocchiut 1.19 //
1174 mocchiut 1.27 printf("\n\n###########################################################################################\n");
1175 mocchiut 1.31 printf("# Benchmark results: nevents = %10u runs = %3u #\n",(nevents+1),nruns);
1176 mocchiut 1.27 printf("###########################################################################################\n");
1177     printf("# Detector # Core routine called # Queries # Real Time # CPU time # Events/s #\n");
1178     printf("###########################################################################################\n");
1179     printf("# RUN # %i # %3u # %8.2f # %8.2f # %8.2f #\n",RUN,NQRUN,fabs(runtime),fabs(cruntime),runrt);
1180 mocchiut 1.47 printf("# TRG # %i # %3u # %8.2f # %8.2f # %8.2f #\n",TRG,NQTRG,fabs(trgtime),fabs(ctrgtime),trgrt);
1181     if ( newfit ){
1182 mocchiut 1.53 printf("# TOF L1 # %i # %3u # %8.2f # %8.2f # %8.2f #\n",(TOF&&TRK),NQTOF1,fabs(toftime1),fabs(ctoftime1),tofrt1);
1183 mocchiut 1.47 printf("# CAL L1 # %i # %3u # %8.2f # %8.2f # %8.2f #\n",CAL,NQCAL1,fabs(caltime1),fabs(ccaltime1),calrt1);
1184     }
1185 mocchiut 1.27 printf("# TRK # %i # %3u # %8.2f # %8.2f # %8.2f #\n",TRK,NQTRK,fabs(trktime),fabs(ctrktime),trkrt);
1186 mocchiut 1.47 printf("# TOF # %i # %3u # %8.2f # %8.2f # %8.2f #\n",TOF,NQTOF,fabs(toftime),fabs(ctoftime),tofrt);
1187 mocchiut 1.27 printf("# CAL # %i # %3u # %8.2f # %8.2f # %8.2f #\n",CAL,NQCAL,fabs(caltime),fabs(ccaltime),calrt);
1188     printf("# AC # %i # %3u # %8.2f # %8.2f # %8.2f #\n",AC,NQAC,fabs(actime),fabs(cactime),acrt);
1189     printf("# S4 # %i # %3u # %8.2f # %8.2f # %8.2f #\n",S4,NQS4,fabs(s4time),fabs(cs4time),s4rt);
1190     printf("# ND # %i # %3u # %8.2f # %8.2f # %8.2f #\n",ND,NQND,fabs(ndtime),fabs(cndtime),ndrt);
1191     printf("# ORB # %i # %3u # %8.2f # %8.2f # %8.2f #\n",ORB,NQORB,fabs(orbtime),fabs(corbtime),orbrt);
1192     printf("###########################################################################################\n");
1193     printf("# Total # # %4u # %8.2f # %8.2f # %8.2f #\n",NQTOT,fabs(dvtime),fabs(cdvtime),dvrt);
1194     printf("###########################################################################################\n");
1195 mocchiut 1.19 };
1196     //
1197     //
1198 mocchiut 1.1 if ( error != 0 ) printf("\n\n WARNING: exiting with signal %i \n\n",error);
1199 mocchiut 1.15 //
1200     //
1201     //
1202     if ( remfile && error < 0 ){
1203     printf(" -c or --clean option used, deleting file %s \n\n",filename.Data());
1204     gSystem->Unlink(filename.Data());
1205     };
1206     //
1207 mocchiut 1.47 if ( !error && (CALSGN1 || TOFSGN1 || CALSGN || TRKSGN || TRGSGN || TOFSGN || RUNSGN || ORBSGN || ACSGN || S4SGN || NDSGN || DVSGN) ) error = 1;
1208 mocchiut 1.34 //
1209 mocchiut 1.15 printf("\n Finished, exiting...\n\n");
1210 mocchiut 1.1 //
1211     // Close redirection if the case.
1212     //
1213     if ( !beverbose ) close(nul);
1214     //
1215     //
1216 mocchiut 1.34 if ( !error ) exit(0); // no errors
1217     if ( error == 1 ) exit(255); // warnings
1218     exit(1); // errors
1219 mocchiut 1.1 }
  ViewVC Help
Powered by ViewVC 1.1.23