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

Diff of /DarthVader/src/DarthVader.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.27 by mocchiut, Wed Jan 31 16:15:57 2007 UTC revision 1.35 by mocchiut, Thu Oct 25 13:24:22 2007 UTC
# Line 30  Line 30 
30  //  //
31  using namespace std;  using namespace std;
32  //  //
33    GL_TABLES *glt = NULL;
34  //  //
35  //  //
36  #include <DarthVaderVerl2.h>  #include <DarthVaderVerl2.h>
# Line 41  void usage(){ Line 42  void usage(){
42    printf("\n DarthVader [ options ] -idRun ID_RUN [+-all] [+-detector [ detector options ] ] \n");    printf("\n DarthVader [ options ] -idRun ID_RUN [+-all] [+-detector [ detector options ] ] \n");
43    printf("\n Options are: \n\n");    printf("\n Options are: \n\n");
44    printf(" --version         print informations about compilation and exit\n");    printf(" --version         print informations about compilation and exit\n");
45    printf(" -h | --help       print this help and exit \n");    printf(" -h || --help      print this help and exit \n");
46    printf(" -v | --verbose    be verbose [default]\n");    printf(" -v || --verbose   be verbose [default]\n");
47    printf(" -s | --silent     print nothing on STDOUT\n");    printf(" -s || --silent    print nothing on STDOUT\n");
48    printf(" -c | --clean      remove file if exiting with errors\n");    printf(" -c || --clean     remove file if exiting with errors\n");
49    printf(" -b | --benchmark  perform and print a benchmark test\n");    printf(" -b || --benchmark perform and print a benchmark test\n");
50    printf(" -auto | -AUTO     exclude from processing detector which are NOT in the acquisition [default]\n");    printf(" -auto || -AUTO    exclude from processing detector which are NOT in the acquisition\n");
51    printf(" -zerofill         if a detector is not in the acquisition the routine is called anyway \n");    printf(" -zerofill         if a detector is not in the acquisition the routine is called anyway \n");
52    printf("                   but all detector's data will be marked as bad\n");    printf("                   but all detector's data will be marked as bad [default]\n");
53    printf(" -tedious          exit with error if a detector is not in the acquisition and \n");    printf(" -tedious          exit with error if a detector is not in the acquisition and \n");
54    printf("                   it has not been excluded from processing\n");    printf("                   it has not been excluded from processing\n");
55    printf(" -host             name of the DB host [default = $PAM_DBHOST or mysql://localhost/pamelaprod]\n");    printf(" -host             name of the DB host [default = $PAM_DBHOST or mysql://localhost/pamelaprod]\n");
# Line 56  void usage(){ Line 57  void usage(){
57    printf(" -psw              password for the DB connection [default = $PAM_DBPSW or \"\"]\n");    printf(" -psw              password for the DB connection [default = $PAM_DBPSW or \"\"]\n");
58    printf(" -idRun ID_RUN     ID number (from the DB) of the run to be processed \n");    printf(" -idRun ID_RUN     ID number (from the DB) of the run to be processed \n");
59    printf(" -processFile file output filename [default ID_RUN.Level2.root]\n");    printf(" -processFile file output filename [default ID_RUN.Level2.root]\n");
60    printf(" +all | +ALL       call all detectors software [default]\n");    printf(" +all || +ALL      call all detectors software [default]\n");
61    printf(" -all | -ALL       call nothing\n");    printf(" -all || -ALL      call nothing\n");
62    printf(" +detector         process detector; detector can be: TOF,TRK,CAL,TRG,ORB,S4,ND,AC,RUN\n");    printf(" +detector         process detector; detector can be: TOF,TRK,CAL,TRG,ORB,S4,ND,AC,RUN\n");
63    printf(" -detector         do not process detector (as above)\n");    printf(" -detector         do not process detector (as above)\n");
64    printf("                   detector options must be included in square parenthesis with spaces, for example:\n");    printf("                   detector options must be included in square parenthesis with spaces, for example:\n");
# Line 131  int main(int numinp, char *inps[]){ Line 132  int main(int numinp, char *inps[]){
132    UInt_t NQTOT = 0;    UInt_t NQTOT = 0;
133    //    //
134    Bool_t autom = true;    Bool_t autom = true;
135    Bool_t zerofill = false;    Bool_t zerofill = true;
136    Bool_t tedious = false;    Bool_t tedious = false;
137    Bool_t remfile = false;    Bool_t remfile = false;
138    Bool_t debug = false;    Bool_t debug = false;
# Line 575  int main(int numinp, char *inps[]){ Line 576  int main(int numinp, char *inps[]){
576      //      //
577      if ( debug ) printf("...connected! \n\n");      if ( debug ) printf("...connected! \n\n");
578      //      //
579      GL_TABLES *glt = new GL_TABLES(host,user,psw);      glt = new GL_TABLES(host,user,psw);
580        //GL_TABLES *glt = new GL_TABLES(host,user,psw);
581      //      //
582      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());      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());
583      //      //
# Line 699  int main(int numinp, char *inps[]){ Line 701  int main(int numinp, char *inps[]){
701      };      };
702      if ( debug ) printf("\n End pre-processing \n\n");      if ( debug ) printf("\n End pre-processing \n\n");
703      //      //
704        dbc->Close();
705        //
706      timer.Start(kTRUE);      timer.Start(kTRUE);
707      if ( TRK ) {      if ( TRK ) {
708        glt->ResetCounters();        glt->ResetCounters();
709        printf(" TrackerLevel2 called\n");                        printf(" TrackerLevel2 called\n");                
710        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start TrkCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start TrkCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
711        TRKSGN = TrkCore(run,processFile,dbc,trkargc,trkargv);        TRKSGN = TrkCore(run,processFile,glt,trkargc,trkargv);
712        gSystem->Unlink("TrackerFolder"); //patch        gSystem->Unlink("TrackerFolder"); //patch
713        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end TrkCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end TrkCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
714        NQTRK = glt->GetNqueries();        NQTRK = glt->GetNqueries();
# Line 718  int main(int numinp, char *inps[]){ Line 722  int main(int numinp, char *inps[]){
722        glt->ResetCounters();        glt->ResetCounters();
723        printf(" ToFLevel2 called\n");        printf(" ToFLevel2 called\n");
724        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
725        TOFSGN = ToFCore(run,processFile,dbc,tofargc,tofargv);        TOFSGN = ToFCore(run,processFile,glt,tofargc,tofargv);
726        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end ToFCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end ToFCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
727        NQTOF = glt->GetNqueries();        NQTOF = glt->GetNqueries();
728      };      };
# Line 731  int main(int numinp, char *inps[]){ Line 735  int main(int numinp, char *inps[]){
735        glt->ResetCounters();        glt->ResetCounters();
736        printf(" CalorimeterLevel2 called\n");        printf(" CalorimeterLevel2 called\n");
737        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
738        CALSGN = CaloCore(run,processFile,dbc,calargc,calargv);        CALSGN = CaloCore(run,processFile,glt,calargc,calargv);
739        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end CaloCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end CaloCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
740        NQCAL = glt->GetNqueries();        NQCAL = glt->GetNqueries();
741      };      };
# Line 744  int main(int numinp, char *inps[]){ Line 748  int main(int numinp, char *inps[]){
748        glt->ResetCounters();        glt->ResetCounters();
749        printf(" TriggerLevel2 called\n");        printf(" TriggerLevel2 called\n");
750        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start TrigCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start TrigCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
751        TRGSGN = TrigCore(run,processFile,dbc,trgargc,trgargv);        TRGSGN = TrigCore(run,processFile,glt,trgargc,trgargv);
752        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end TrigCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end TrigCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
753        NQTRG = glt->GetNqueries();        NQTRG = glt->GetNqueries();
754      };      };
# Line 757  int main(int numinp, char *inps[]){ Line 761  int main(int numinp, char *inps[]){
761        glt->ResetCounters();        glt->ResetCounters();
762        printf(" AnticounterLevel2 called\n");        printf(" AnticounterLevel2 called\n");
763        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start AcCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start AcCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
764        ACSGN = AcCore(run,processFile,dbc,acargc,acargv);        ACSGN = AcCore(run,processFile,glt,acargc,acargv);
765        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end AcCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end AcCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
766        NQAC = glt->GetNqueries();        NQAC = glt->GetNqueries();
767      };      };
# Line 770  int main(int numinp, char *inps[]){ Line 774  int main(int numinp, char *inps[]){
774        glt->ResetCounters();        glt->ResetCounters();
775        printf(" S4Level2 called\n");        printf(" S4Level2 called\n");
776        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start S4Core STDOUT  |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start S4Core STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
777        S4SGN = S4Core(run,processFile,dbc,s4argc,s4argv);        S4SGN = S4Core(run,processFile,glt,s4argc,s4argv);
778        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end S4Core STDOUT   |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end S4Core STDOUT   |<<<<<<<<<<<<<<<<<<<\n");
779        NQS4 = glt->GetNqueries();        NQS4 = glt->GetNqueries();
780      };      };
# Line 783  int main(int numinp, char *inps[]){ Line 787  int main(int numinp, char *inps[]){
787        glt->ResetCounters();        glt->ResetCounters();
788        printf(" NDLevel2 called\n");        printf(" NDLevel2 called\n");
789        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start NDCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start NDCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
790        NDSGN = NDCore(run,processFile,dbc,ndargc,ndargv);        NDSGN = NDCore(run,processFile,glt,ndargc,ndargv);
791        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end NDCore STDOUT   |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end NDCore STDOUT   |<<<<<<<<<<<<<<<<<<<\n");
792        NQND = glt->GetNqueries();        NQND = glt->GetNqueries();
793      };      };
# Line 796  int main(int numinp, char *inps[]){ Line 800  int main(int numinp, char *inps[]){
800        glt->ResetCounters();        glt->ResetCounters();
801        printf(" OrbitalInfo called\n");        printf(" OrbitalInfo called\n");
802        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start OrbitalInfoCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start OrbitalInfoCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
803        ORBSGN = OrbitalInfoCore(run,processFile,dbc,orbargc,orbargv);        ORBSGN = OrbitalInfoCore(run,processFile,glt,orbargc,orbargv);
804        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end OrbitalInfoCore STDOUT   |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end OrbitalInfoCore STDOUT   |<<<<<<<<<<<<<<<<<<<\n");
805        NQORB = glt->GetNqueries();        NQORB = glt->GetNqueries();
806      };          };    
# Line 840  int main(int numinp, char *inps[]){ Line 844  int main(int numinp, char *inps[]){
844      case -52:  message += " GLTABLES - No entries matching GL_PARAM query"; break;      case -52:  message += " GLTABLES - No entries matching GL_PARAM query"; break;
845      case -53:  message += " GLTABLES - No entries matching GL_TRK_CALIB query"; break;      case -53:  message += " GLTABLES - No entries matching GL_TRK_CALIB query"; break;
846      case -54:  message += " GLTABLES - No entries matching GL_CALO_CALIB query"; break;      case -54:  message += " GLTABLES - No entries matching GL_CALO_CALIB query"; break;
847      case -55:  message += " GLTABLES - No entries matching GL_CALO_CALIB query"; break;      case -55:  message += " GLTABLES - No entries matching GL_S4_CALIB query"; break;
848      case -56:  message += " GLTABLES - No entries matching GL_TLE query"; break;      case -56:  message += " GLTABLES - No entries matching GL_TLE query"; break;
849      case -57:  message += " GLTABLES - DB connection gone and not able to reconnect"; break;      case -57:  message += " GLTABLES - DB connection gone and not able to reconnect"; break;
850        //              //      
# Line 863  int main(int numinp, char *inps[]){ Line 867  int main(int numinp, char *inps[]){
867      case -116: message += " CALORIMETERLEVEL2 -  DB connection problems"; break;      case -116: message += " CALORIMETERLEVEL2 -  DB connection problems"; break;
868      case -117: message += " CALORIMETERLEVEL2 -  Cannot reprocess file with Level1 data without calling Level1 routine"; break;      case -117: message += " CALORIMETERLEVEL2 -  Cannot reprocess file with Level1 data without calling Level1 routine"; break;
869      case -118: message += " CALORIMETERLEVEL2 -  Cannot reprocess file without Level1 data calling Level1 routine"; break;      case -118: message += " CALORIMETERLEVEL2 -  Cannot reprocess file without Level1 data calling Level1 routine"; break;
870        case -119: message += " CALORIMETERLEVEL2 -  No CalibCalPulse2 TTree in Level0 file needed for the calibration"; break;
871        case -120: message += " CALORIMETERLEVEL2 -  No CalibCalPulse1 TTree in Level0 file needed for the calibration"; break;
872        //        //
873      case -200: message += " TRACKERLEVEL2 - LEVEL1 framework unknown (HBOOK/ROOT)"; break;      case -200: message += " TRACKERLEVEL2 - LEVEL1 framework unknown (HBOOK/ROOT)"; break;
874      case -201: message += " TRACKERLEVEL2 - LEVEL2 framework unknown (HBOOK/ROOT)"; break;      case -201: message += " TRACKERLEVEL2 - LEVEL2 framework unknown (HBOOK/ROOT)"; break;
# Line 939  int main(int numinp, char *inps[]){ Line 945  int main(int numinp, char *inps[]){
945    case 101: printf("\n WARNING CALORIMETER - No tracks or good events in this run\n");    case 101: printf("\n WARNING CALORIMETER - No tracks or good events in this run\n");
946    };    };
947    //    //
948      // OrbitalInfoCore:
949      //
950      switch(ORBSGN){
951      case 900: printf("\n WARNING ORBITALINFO - No inclination MCMDs\n");
952      };
953      //
954    // no other at the moment    // no other at the moment
955    //    //
956    //    //
# Line 970  int main(int numinp, char *inps[]){ Line 982  int main(int numinp, char *inps[]){
982      Float_t s4rt = 0.;      Float_t s4rt = 0.;
983      Float_t ndrt = 0.;      Float_t ndrt = 0.;
984      Float_t orbrt = 0.;      Float_t orbrt = 0.;
985      Float_t dvrt = nevents/dvtime;      Float_t dvrt = (nevents+1)/dvtime;
986      if ( RUN ) runrt = nruns/runtime;      if ( RUN ) runrt = nruns/runtime;
987      if ( TRK ) trkrt = nevents/trktime;      if ( TRK ) trkrt = (nevents+1)/trktime;
988      if ( CAL ) calrt = nevents/caltime;      if ( CAL ) calrt = (nevents+1)/caltime;
989      if ( TOF ) tofrt = nevents/toftime;      if ( TOF ) tofrt = (nevents+1)/toftime;
990      if ( TRG ) trgrt = nevents/trgtime;      if ( TRG ) trgrt = (nevents+1)/trgtime;
991      if ( AC ) acrt = nevents/actime;      if ( AC ) acrt = (nevents+1)/actime;
992      if ( S4 ) s4rt = nevents/s4time;      if ( S4 ) s4rt = (nevents+1)/s4time;
993      if ( ND ) ndrt = nevents/ndtime;      if ( ND ) ndrt = (nevents+1)/ndtime;
994      if ( ORB ) orbrt = nevents/orbtime;      if ( ORB ) orbrt = (nevents+1)/orbtime;
995      //      //
996      printf("\n\n###########################################################################################\n");      printf("\n\n###########################################################################################\n");
997      printf("#              Benchmark results:      nevents = %10u     runs = %3u                #\n",nevents,nruns);      printf("#              Benchmark results:      nevents = %10u     runs = %3u                #\n",(nevents+1),nruns);
998      printf("###########################################################################################\n");      printf("###########################################################################################\n");
999      printf("# Detector # Core routine called # Queries  #     Real Time   #   CPU time   #  Events/s  #\n");      printf("# Detector # Core routine called # Queries  #     Real Time   #   CPU time   #  Events/s  #\n");
1000      printf("###########################################################################################\n");      printf("###########################################################################################\n");
# Line 1010  int main(int numinp, char *inps[]){ Line 1022  int main(int numinp, char *inps[]){
1022      gSystem->Unlink(filename.Data());      gSystem->Unlink(filename.Data());
1023    };    };
1024    //    //
1025      if ( !error && (CALSGN || TRKSGN || TRGSGN || TOFSGN || RUNSGN || ORBSGN || ACSGN || S4SGN || NDSGN || DVSGN) ) error = 1;
1026      //
1027    printf("\n Finished, exiting...\n\n");    printf("\n Finished, exiting...\n\n");
1028    //    //
1029    // Close redirection if the case.    // Close redirection if the case.
# Line 1017  int main(int numinp, char *inps[]){ Line 1031  int main(int numinp, char *inps[]){
1031    if ( !beverbose ) close(nul);    if ( !beverbose ) close(nul);
1032    //    //
1033    //    //
1034    exit(error);    if ( !error ) exit(0);       // no errors
1035      if ( error == 1 ) exit(255);  // warnings
1036      exit(1);                     // errors
1037  }  }
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.23