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

Annotation of /DarthVader/src/DarthVader.cpp

Parent Directory Parent Directory | Revision Log Revision Log

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