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

Annotation of /DarthVader/src/DarthVader.cpp

Parent Directory Parent Directory | Revision Log Revision Log

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