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

Diff of /YodaProfiler/src/YodaProfiler.cpp

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

revision 1.14 by mocchiut, Tue Nov 28 13:30:39 2006 UTC revision 1.28 by mocchiut, Mon Mar 10 10:58:02 2008 UTC
# Line 5  Line 5 
5  //  //
6  #include <TSystem.h>  #include <TSystem.h>
7  //  //
8    #include <GLTables.h>
9  #include <PamelaDBOperations.h>  #include <PamelaDBOperations.h>
10  #include <YodaProfilerVerl2.h>  #include <YodaProfilerVerl2.h>
11  //  //
# Line 23  void usage(){ Line 24  void usage(){
24    printf(" -v | --verbose  be verbose [default]\n");    printf(" -v | --verbose  be verbose [default]\n");
25    printf(" -s | --silent   print nothing on STDOUT\n");    printf(" -s | --silent   print nothing on STDOUT\n");
26    printf(" -g | --debug    be very verbose [default: no]\n");    printf(" -g | --debug    be very verbose [default: no]\n");
27      printf(" -p | --pedantic be strict in checks, exit with error if any problem is found [default: no]\n");
28    printf(" -boot number    CPU boot number [default = taken from VarDump]\n");    printf(" -boot number    CPU boot number [default = taken from VarDump]\n");
29    printf(" -autoboot       if no VarDump found try to determine the BOOT number\n");    printf(" -autoboot       if no VarDump found try to determine the BOOT number\n");
30    printf("                 looking at timesync [default]\n");    printf("                 looking at timesync [default]\n");
31    printf(" -no-autoboot    disable previous check\n");    printf(" -no-autoboot    disable previous check\n");
32    printf(" -tsync number   timesync (s) [default = taken from data]\n");    printf(" -tsync number   timesync (s) [default = taken from data]\n");
33    printf(" -obt0 number    obt at timesync (ms) [default = taken from data]\n");    printf(" -obt0 number    obt at timesync (ms) [default = taken from data]\n");
34      printf(" -dworbit number number is the downlink orbit number (to be used with non-standard names)\n");
35      printf(" -static         use static path inside the DB instead of symbolic names\n");
36      //  printf(" -gpamela        shortcut to \"-clean 0 -obt0 1 -tsync 1 -no-autoboot -boot 1 -dworbit 1 -static\"\n");
37      printf(" -gpamela        shortcut to \"-clean 0 -obt0 1 -no-autoboot -dworbit 1 -static\"\n");
38    printf(" -clean number   number in seconds after which the fragment table\n");    printf(" -clean number   number in seconds after which the fragment table\n");
39    printf("                 can be cleaned and runs validated [default = -1 do not clean],\n");    printf("                 can be cleaned and runs validated [default = -1 do not clean],\n");
40    printf("                 if 0 force cleaning immediatly, if negative do not clean\n");    printf("                 if 0 force cleaning immediatly, if negative do not clean\n");
# Line 43  void usage(){ Line 49  void usage(){
49    printf("                 'same' can be used if in conjuction with -yodaFile\n");    printf("                 'same' can be used if in conjuction with -yodaFile\n");
50    printf(" -nofrag         do not leave runs in the fragment table and look for fragments\n");    printf(" -nofrag         do not leave runs in the fragment table and look for fragments\n");
51    printf("                 in the GL_RUN table.\n");    printf("                 in the GL_RUN table.\n");
52    printf(" -force          to be used to reprocess a file or to process a file\n");    printf(" -f | -force     to be used to reprocess a file or to process a file\n");
53    printf("                 when already validated the surroundings, it is equivalent to:\n");    printf("                 when already validated the surroundings, it is equivalent to:\n");
54    printf("                 -remove same -validate same -cleanfrag same -nofrag\n");    printf("                 -remove same -validate same -cleanfrag same -nofrag\n");
55      printf(" -check          performs some test on the GL_RUN table to check its consistency\n");
56      printf(" -checkfrom idrun performs some test on the GL_RUN table to check its consistency starting from run number idrun\n");
57      printf(" -checkto idrun   performs some test on the GL_RUN table to check its consistency arriving till run number idrun\n");
58    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");
59    printf(" -user           username for the DB connection [default = $PAM_DBUSER or \"anonymous\"] \n");    printf(" -user           username for the DB connection [default = $PAM_DBUSER or \"anonymous\"] \n");
60    printf(" -psw            password for the DB connection [default = $PAM_DBPSW or \"\"]\n");    printf(" -psw            password for the DB connection [default = $PAM_DBPSW or \"\"]\n");
61    printf(" -tle <file>     ascii file containing TLE obtained from celestrak.org or space-track.org [default = no]\n");    printf(" -tle <file>     ascii file containing TLE obtained from celestrak.org or space-track.org [default = no]\n");
62    printf("\n The order of input files and options does not matter. \n");    printf("\n The order of input files and options does not matter. \n");
63    printf("\n Example: \n");    printf("\n Example: \n");
64    printf(" YodaProfiler -yodaFile /path/to/raw/files/000_000_00000_cln2.pam -rawFile /path/to/filesfromyoda/000_000_00000_cln2.root -v \n\n");    printf(" YodaProfiler -yodaFile /path/to/raw/files/000_000_00000_cln2.root -rawFile /path/to/filesfromyoda/000_000_00000_cln2.pam -v \n\n");
65  };  };
66  //  //
67  int main(int numinp, char *inps[]){  int main(int numinp, char *inps[]){
# Line 64  int main(int numinp, char *inps[]){ Line 73  int main(int numinp, char *inps[]){
73    UInt_t boot = 0;    UInt_t boot = 0;
74    UInt_t tsync = 0;    UInt_t tsync = 0;
75    UInt_t obt0 = 0;    UInt_t obt0 = 0;
76      UInt_t dwinput = 0;
77      UInt_t chfrom = 0;
78      UInt_t chto = 0;
79    Long64_t olderthan = -1LL;      Long64_t olderthan = -1LL;  
80    //    //
81    //    //
# Line 87  int main(int numinp, char *inps[]){ Line 99  int main(int numinp, char *inps[]){
99    if ( strcmp(pamdbpsw,"") ) password = pamdbpsw;    if ( strcmp(pamdbpsw,"") ) password = pamdbpsw;
100    //    //
101    //        //    
102      Bool_t staticp = false;
103    Bool_t beverbose = true;    Bool_t beverbose = true;
104    Bool_t debug = false;    Bool_t debug = false;
105    Bool_t autoboot = true;    Bool_t autoboot = true;
106      Bool_t gpamela = false;
107    //    //
108    Bool_t remove = false;    Bool_t remove = false;
109    TString remfile = "";    TString remfile = "";
# Line 104  int main(int numinp, char *inps[]){ Line 118  int main(int numinp, char *inps[]){
118    //    //
119    Bool_t force = false;    Bool_t force = false;
120    //    //
121      Bool_t check = false;
122      //
123      Bool_t pedantic = false;
124      //
125    Int_t i = 0;    Int_t i = 0;
126    //    //
127    if ( numinp > 1 ){    if ( numinp > 1 ){
# Line 155  int main(int numinp, char *inps[]){ Line 173  int main(int numinp, char *inps[]){
173          validfile = (TString)inps[i+1];          validfile = (TString)inps[i+1];
174          validate = true;          validate = true;
175        };        };
176          if ( !strcmp(inps[i],"-gpamela") ) {
177            gpamela = true;
178          };
179          if ( !strcmp(inps[i],"-static") ) {
180            staticp = true;
181          };
182        if ( !strcmp(inps[i],"-autoboot") ) {        if ( !strcmp(inps[i],"-autoboot") ) {
183          autoboot = true;          autoboot = true;
184        };        };
# Line 168  int main(int numinp, char *inps[]){ Line 192  int main(int numinp, char *inps[]){
192          };          };
193          boot = atoi(inps[i+1]);          boot = atoi(inps[i+1]);
194        };        };
195          if ( !strcmp(inps[i],"-checkfrom") ) {
196            if ( numinp-1 < i+1 ){
197              usage();
198              exit(1);
199            };
200            chfrom = atoi(inps[i+1]);      
201            check = true;
202          };
203          if ( !strcmp(inps[i],"-checkto") ) {
204            if ( numinp-1 < i+1 ){
205              usage();
206              exit(1);
207            };
208            chto = atoi(inps[i+1]);
209            check = true;
210          };
211        if ( !strcmp(inps[i],"-tsync") ) {        if ( !strcmp(inps[i],"-tsync") ) {
212          if ( numinp-1 < i+1 ){          if ( numinp-1 < i+1 ){
213            usage();            usage();
# Line 182  int main(int numinp, char *inps[]){ Line 222  int main(int numinp, char *inps[]){
222          };          };
223          obt0 = (UInt_t)atoll(inps[i+1]);                  obt0 = (UInt_t)atoll(inps[i+1]);        
224        };        };
225          if ( !strcmp(inps[i],"-dworbit") ) {
226            if ( numinp-1 < i+1 ){
227              usage();
228              exit(1);
229            };
230            dwinput = (UInt_t)atoll(inps[i+1]);    
231          };
232        if ( !strcmp(inps[i],"-clean") ) {        if ( !strcmp(inps[i],"-clean") ) {
233          if ( numinp-1 < i+1 ){          if ( numinp-1 < i+1 ){
234            usage();            usage();
# Line 226  int main(int numinp, char *inps[]){ Line 273  int main(int numinp, char *inps[]){
273        //        //
274        if ( !strcmp(inps[i],"-f") || !strcmp(inps[i],"-force") ) force = true;        if ( !strcmp(inps[i],"-f") || !strcmp(inps[i],"-force") ) force = true;
275        //        //
276          if ( !strcmp(inps[i],"-p") || !strcmp(inps[i],"--pedantic") ) pedantic = true;
277          //
278        if ( !strcmp(inps[i],"-nofrag") ) nofrag = true;        if ( !strcmp(inps[i],"-nofrag") ) nofrag = true;
279        //        //
280          if ( !strcmp(inps[i],"-check") ) check = true;
281          //
282        i++;        i++;
283      };      };
284    } else {    } else {
# Line 272  int main(int numinp, char *inps[]){ Line 323  int main(int numinp, char *inps[]){
323      validfile = fcleanfile;      validfile = fcleanfile;
324    };    };
325    //    //
326    // Start:    const char *paml0 = gSystem->Getenv("PAM_L0");
327      const char *paml2 = gSystem->Getenv("PAM_L2");
328      const char *pamraw = gSystem->Getenv("PAM_RAW");
329      if ( !paml0 || !paml2 || !pamraw ){
330        printf(" WARNING: to use symbolic paths in DB you must set up PAM_L0, PAM_L2 and PAM_RAW env var!\n Using static paths\n");
331        staticp = true;
332      };
333      //
334      //
335      GL_TABLES *glt = new GL_TABLES(host,user,password);
336    //    //
337      if ( gpamela ){
338        //    //    -clean 0 -obt0 1 -tsync 1 -boot 1 -dworbit 1
339        //    -clean 0 -obt0 1 -dworbit 1
340        if ( olderthan < 0LL ) olderthan = 0LL;
341        if ( !obt0 ) obt0 = 1;
342        //    if ( !tsync ) tsync = 1;    
343        //    if ( !boot ) boot = 1;
344        if ( !dwinput ) dwinput = 1;
345        staticp = true;
346        autoboot = false;
347      };
348      //
349      if ( filerawname == "" && filerootname != "" ){
350        filerawname = "/not_given/"+(TString)gSystem->BaseName(filerootname.Data());
351        filerawname = filerawname.ReplaceAll(".root",5,".pam",4);
352        //    printf("filerawname %s \n",filerawname.Data());
353      };
354      //
355      // Start:
356    TString message;    TString message;
357    char *version = YodaProfilerInfo(false);    char *version = YodaProfilerInfo(false);
358    PamelaDBOperations *pamDB = 0;    PamelaDBOperations *pamDB = 0;
# Line 288  int main(int numinp, char *inps[]){ Line 367  int main(int numinp, char *inps[]){
367      // Create pamDB object and open SQL connection      // Create pamDB object and open SQL connection
368      //-------------------------------------------------------------------------------------------      //-------------------------------------------------------------------------------------------
369      if ( beverbose ) printf(" 1 => Initialize and open SQL connection \n");      if ( beverbose ) printf(" 1 => Initialize and open SQL connection \n");
370      pamDB = new PamelaDBOperations(host,user,password,filerawname,filerootname,boot,tsync,obt0,debug,tlefilename);      pamDB = new PamelaDBOperations(host,user,password,filerawname,filerootname,boot,tsync,obt0,debug,tlefilename,dwinput,staticp,gpamela);
371      pamDB->CheckConnection();      pamDB->CheckConnection();
372      //-------------------------------------------------------------------------------------------      //-------------------------------------------------------------------------------------------
373      //      //
374      pamDB->SetNoFrag(nofrag);      pamDB->SetNoFrag(nofrag);
375      pamDB->SetAutoBoot(autoboot);      pamDB->SetAutoBoot(autoboot);
376        pamDB->SetPedantic(pedantic);
377      //      //
378      if ( remove ){      if ( remove ){
379        //-------------------------------------------------------------------------------------------        //-------------------------------------------------------------------------------------------
# Line 353  int main(int numinp, char *inps[]){ Line 433  int main(int numinp, char *inps[]){
433        if ( beverbose ) printf(" 7 => Insert calorimeter calibrations in the GL_CALO_CALIB table\n");        if ( beverbose ) printf(" 7 => Insert calorimeter calibrations in the GL_CALO_CALIB table\n");
434        WAR[5] = pamDB->insertCALO_CALIB();        WAR[5] = pamDB->insertCALO_CALIB();
435        //-------------------------------------------------------------------------------------------            //-------------------------------------------------------------------------------------------    
436            
437          if ( !gpamela ){
438            //---------------------------------------------------------------------------------------------------------------------
439            //Insert in GL_CALOPULSE_CALIB calibration PULSE information records relative to a single unpack (only for real data!)
440            //---------------------------------------------------------------------------------------------------------------------
441            if ( beverbose ) printf(" 7b => Insert calorimeter pulse calibrations in the GL_CALOPULSE_CALIB table\n");
442            WAR[5] += pamDB->insertCALOPULSE_CALIB();
443            //-------------------------------------------------------------------------------------------    
444          };
445    
446        //-------------------------------------------------------------------------------------------        //-------------------------------------------------------------------------------------------
447        //Insert in GL_TRK_CALIB calibration information records relative to a single unpack        //Insert in GL_TRK_CALIB calibration information records relative to a single unpack
448        //-------------------------------------------------------------------------------------------        //-------------------------------------------------------------------------------------------
# Line 411  int main(int numinp, char *inps[]){ Line 500  int main(int numinp, char *inps[]){
500        if ( beverbose ) printf(" 12 => Insert TLE elements in the GL_TLE table from file %s \n",tlefilename.Data());        if ( beverbose ) printf(" 12 => Insert TLE elements in the GL_TLE table from file %s \n",tlefilename.Data());
501        WAR[10] = pamDB->populateTLE();        WAR[10] = pamDB->populateTLE();
502      }      }
503    
504        if ( check ){
505          //----------------------------------------------
506          //Check GL_RUN table consistency
507          //----------------------------------------------
508          if ( beverbose ) printf(" 13 => Check GL_RUN table consistency\n");
509          UInt_t ch = pamDB->Check(chfrom,chto);
510          if ( !ch ){
511            printf(" GL_RUN table seems ok!\n");
512          } else {
513            printf(" Problems found in the GL_RUN table!\n");
514          };
515        };
516                    
517    } catch (Int_t exc) {    } catch (Int_t exc) {
518      signal = exc;      signal = exc;
# Line 444  int main(int numinp, char *inps[]){ Line 546  int main(int numinp, char *inps[]){
546      case -27: message += " Cannot determine downlink orbit, wrong filename format?"; break;      case -27: message += " Cannot determine downlink orbit, wrong filename format?"; break;
547      case -28: message += " Cannot assign a RUN ID unique number!"; break;      case -28: message += " Cannot assign a RUN ID unique number!"; break;
548      case -29: message += " No VarDump and impossible to find a consistent BOOT number, use the -boot option to override"; break;      case -29: message += " No VarDump and impossible to find a consistent BOOT number, use the -boot option to override"; break;
549        case -30: message += " Old database selected, no GL_CALOPULSE_CALIB table in the DB!"; break;
550        case -31: message += " No CalibCalPulse1 tree in Level0 file"; break;
551        case -32: message += " No CalibCalPulse2 tree in Level0 file"; break;
552        case -66: message += " PEDANTIC: REPETITION"; break;
553        case -67: message += " PEDANTIC: TIMESYNC ALREADY INSERTED"; break;
554        case -68: message += " PEDANTIC: NO TIMESYNC IN THE FILE"; break;
555        case -69: message += " PEDANTIC: RUN FRAGMENT ALREADY INSERTED"; break;
556        case -70: message += " PEDANTIC: RUN ALREADY INSERTED"; break;
557        case -71: message += " PEDANTIC: OVERLAPPING RUNS (FIRST)"; break;
558        case -72: message += " PEDANTIC: OVERLAPPING RUNS (LAST)"; break;
559        case -73: message += " PEDANTIC: CALO CALIB ALREADY INSERTED"; break;
560        case -74: message += " PEDANTIC: CALO CALIB REPEATED"; break;
561        case -75: message += " PEDANTIC: CALOPULSE1 CALIB ALREADY INSERTED"; break;
562        case -76: message += " PEDANTIC: CALOPULSE1 CALIB REPEATED"; break;
563        case -77: message += " PEDANTIC: CALOPULSE2 CALIB ALREADY INSERTED"; break;
564        case -78: message += " PEDANTIC: CALOPULSE2 CALIB REPEATED"; break;
565        case -79: message += " PEDANTIC: TRK1 CALIB REPEATED"; break;
566        case -80: message += " PEDANTIC: TRK CALIB ALREADY INSERTED"; break;
567        case -81: message += " PEDANTIC: S4 CALIB ALREADY INSERTED"; break;
568        case -82: message += " PEDANTIC: S4 CALIB REPEATED"; break;
569        case -83: message += " PEDANTIC: CLEANING ERROR, RUN ALREADY INSERTED IN THE GL_RUN TABLE"; break;
570      default: message += " Unidentified error"; break;      default: message += " Unidentified error"; break;
571      };      };
572      printf("\n");      printf("\n");
# Line 476  int main(int numinp, char *inps[]){ Line 599  int main(int numinp, char *inps[]){
599              else if ( bit == 2 ) message += "=> No runheaders in the file\n";              else if ( bit == 2 ) message += "=> No runheaders in the file\n";
600              else if ( bit == 3 ) message += "=> No runtrailers in the file\n";              else if ( bit == 3 ) message += "=> No runtrailers in the file\n";
601              else if ( bit == 4 ) message += "=> No mcmd inclination in the file\n";              else if ( bit == 4 ) message += "=> No mcmd inclination in the file\n";
602                else if ( bit == 5 ) message += "=> Inconsistent PKT/OBT sequence\n";
603                else if ( bit == 6 ) message += "=> No physics events in the file\n";
604                else if ( bit == 7 ) message += "=> Less than 2 physics events in the file\n";
605              else  message += "=> Unidentified insertTimeSync warning\n";              else  message += "=> Unidentified insertTimeSync warning\n";
606            };            };
607          };          };
# Line 509  int main(int numinp, char *inps[]){ Line 635  int main(int numinp, char *inps[]){
635              if      ( bit == 0 ) message += "=> Inconsistent PKT/OBT sequence\n";              if      ( bit == 0 ) message += "=> Inconsistent PKT/OBT sequence\n";
636              else if ( bit == 1 ) message += "=> No physics events in the file\n";              else if ( bit == 1 ) message += "=> No physics events in the file\n";
637              else if ( bit == 2 ) message += "=> Less than 2 physics events in the file\n";              else if ( bit == 2 ) message += "=> Less than 2 physics events in the file\n";
638                else if ( bit == 3 ) message += "=> Not supported yet: run with no events, no runtrailer, no runheader\n";
639                else if ( bit == 4 ) message += "=> File with no events, no runtrailers, no runheaders\n";
640              else  message += "=> Unidentified insertPamelaRun warning\n";              else  message += "=> Unidentified insertPamelaRun warning\n";
641            };            };
642          };          };
# Line 598  int main(int numinp, char *inps[]){ Line 726  int main(int numinp, char *inps[]){
726    if ( beverbose ) printf(" 13 => Free objects and close SQL connection \n");    if ( beverbose ) printf(" 13 => Free objects and close SQL connection \n");
727    pamDB->Close();    pamDB->Close();
728    //    //
729      if ( debug ) printf(" Total number of queries through GLTABLES: %u \n",glt->GetNqueries());
730      delete glt;
731      //
732    printf("\n");    printf("\n");
733    printf(" Finished, exiting...\n");    printf(" Finished, exiting...\n");
734    printf("\n");    printf("\n");
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.23