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

Annotation of /DarthVader/src/DarthVader.cpp

Parent Directory Parent Directory | Revision Log Revision Log

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