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

Annotation of /DarthVader/src/DarthVader.cpp

Parent Directory Parent Directory | Revision Log Revision Log

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