/[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.46 by mocchiut, Wed Mar 6 14:20:27 2013 UTC revision 1.47 by mocchiut, Tue Nov 5 14:02:34 2013 UTC
# Line 47  void usage(){ Line 47  void usage(){
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(" -r || --reprocess force running on existing file [default: exit with error, NB: DEFAULT CHANGED since 10RED!!]\n");
51      printf(" -n || --new-fit   use new fitting algorithm [default]\n");
52      printf(" --no-new-fit      use standard (up to 9RED) fitting algorithm\n");
53    printf(" -auto || -AUTO    exclude from processing detector which are NOT in the acquisition\n");    printf(" -auto || -AUTO    exclude from processing detector which are NOT in the acquisition\n");
54    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");
55    printf("                   but all detector's data will be marked as bad [default]\n");    printf("                   but all detector's data will be marked as bad [default]\n");
# Line 91  int main(int numinp, char *inps[]){ Line 94  int main(int numinp, char *inps[]){
94    Double_t ctrktime = 0.;    Double_t ctrktime = 0.;
95    Double_t caltime = 0.;    Double_t caltime = 0.;
96    Double_t ccaltime = 0.;    Double_t ccaltime = 0.;
97      Double_t caltime1 = 0.;
98      Double_t ccaltime1 = 0.;
99    Double_t toftime = 0.;    Double_t toftime = 0.;
100    Double_t ctoftime = 0.;    Double_t ctoftime = 0.;
101      Double_t toftime1 = 0.;
102      Double_t ctoftime1 = 0.;
103    Double_t trgtime = 0.;    Double_t trgtime = 0.;
104    Double_t ctrgtime = 0.;    Double_t ctrgtime = 0.;
105    Double_t actime = 0.;    Double_t actime = 0.;
# Line 111  int main(int numinp, char *inps[]){ Line 118  int main(int numinp, char *inps[]){
118    Int_t error = 0;    Int_t error = 0;
119    //    //
120    Int_t CALSGN = 0;    Int_t CALSGN = 0;
121      Int_t CALSGN1 = 0;
122    Int_t TRKSGN = 0;    Int_t TRKSGN = 0;
123    Int_t TRGSGN = 0;    Int_t TRGSGN = 0;
124    Int_t TOFSGN = 0;    Int_t TOFSGN = 0;
125      Int_t TOFSGN1 = 0;
126    Int_t RUNSGN = 0;    Int_t RUNSGN = 0;
127    Int_t ORBSGN = 0;    Int_t ORBSGN = 0;
128    Int_t ACSGN = 0;    Int_t ACSGN = 0;
# Line 124  int main(int numinp, char *inps[]){ Line 133  int main(int numinp, char *inps[]){
133    UInt_t NQRUN = 0;    UInt_t NQRUN = 0;
134    UInt_t NQTRK = 0;    UInt_t NQTRK = 0;
135    UInt_t NQCAL = 0;    UInt_t NQCAL = 0;
136      UInt_t NQCAL1 = 0;
137    UInt_t NQTOF = 0;    UInt_t NQTOF = 0;
138      UInt_t NQTOF1 = 0;
139    UInt_t NQORB = 0;    UInt_t NQORB = 0;
140    UInt_t NQTRG = 0;    UInt_t NQTRG = 0;
141    UInt_t NQAC = 0;    UInt_t NQAC = 0;
# Line 140  int main(int numinp, char *inps[]){ Line 151  int main(int numinp, char *inps[]){
151    Bool_t beverbose = true;    Bool_t beverbose = true;
152    Bool_t givenid = false;    Bool_t givenid = false;
153    Bool_t bench = false;    Bool_t bench = false;
154      Bool_t reprocess = false;
155      Bool_t newfit = true;
156    Bool_t CAL = true;    Bool_t CAL = true;
157    Bool_t TRK = true;    Bool_t TRK = true;
158    Bool_t TRG = true;    Bool_t TRG = true;
# Line 285  int main(int numinp, char *inps[]){ Line 298  int main(int numinp, char *inps[]){
298            bench = true;            bench = true;
299            found = true;            found = true;
300          };          };
301            if ( !strcmp(inps[i],"-r") || !strcmp(inps[i],"--reprocess") ){
302              reprocess = true;
303              found = true;
304            };
305            if ( !strcmp(inps[i],"-n") || !strcmp(inps[i],"--new-fit") ){
306              newfit = true;
307              found = true;
308            };
309            if ( !strcmp(inps[i],"--no-new-fit") ){
310              newfit = false;
311              found = true;
312            };
313          if ( !strcmp(inps[i],"-auto") || !strcmp(inps[i],"-AUTO") ){          if ( !strcmp(inps[i],"-auto") || !strcmp(inps[i],"-AUTO") ){
314            autom = true;            autom = true;
315            tedious = false;            tedious = false;
# Line 604  int main(int numinp, char *inps[]){ Line 629  int main(int numinp, char *inps[]){
629        filename += strun.str();        filename += strun.str();
630        filename += ".Level2.root";        filename += ".Level2.root";
631      };      };
632        //
633        // check if file exists and check reprocess flag
634        //
635        if ( !reprocess ){
636          if ( !gSystem->GetPathInfo(filename.Data(),NULL,(Long_t*)NULL,NULL,NULL) ) {
637            if ( beverbose ) printf(" File %s already exists!\n",filename.Data());
638            throw -37;
639          }
640        }
641        
642      processFile = new TFile(filename.Data(),"UPDATE");      processFile = new TFile(filename.Data(),"UPDATE");
643      if ( !processFile->IsOpen() ) throw -15;          if ( !processFile->IsOpen() ) throw -15;    
644      //      //
# Line 718  int main(int numinp, char *inps[]){ Line 753  int main(int numinp, char *inps[]){
753      trgtime = timer.RealTime();      trgtime = timer.RealTime();
754      ctrgtime = timer.CpuTime();      ctrgtime = timer.CpuTime();
755      //      //
756        //
757        // new fitting algorithm require calorimeter level1 and tof XXX to be processed first
758        //
759        if ( newfit ) {
760          timer.Start(kTRUE);
761          if ( TOF && false) {
762            glt->ResetCounters();
763            printf(" ToFLevel2 called, 1st call\n");
764            if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
765            TOFSGN1 = ToFCore(run,processFile,glt,tofargc,tofargv);
766            if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end ToFCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
767            NQTOF1 = glt->GetNqueries();
768          };
769          timer.Stop();
770          toftime1 = timer.RealTime();
771          ctoftime1 = timer.CpuTime();
772          //
773          timer.Start(kTRUE);
774          if ( CAL ) {
775            glt->ResetCounters();
776            printf(" CalorimeterLevel2 called, 1st call\n");
777            if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
778            string cstring = "--level1-only";
779            calargv[calargc] = (char *)cstring.c_str();
780            calargc++;
781            CALSGN1 = CaloCore(run,processFile,glt,calargc,calargv);
782            if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end CaloCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
783            NQCAL1 = glt->GetNqueries();
784          };
785          timer.Stop();
786          caltime1 = timer.RealTime();
787          ccaltime1 = timer.CpuTime();    
788        }
789        //
790        //
791      timer.Start(kTRUE);      timer.Start(kTRUE);
792      if ( TRK ) {      if ( TRK ) {
793        glt->ResetCounters();        glt->ResetCounters();
# Line 750  int main(int numinp, char *inps[]){ Line 820  int main(int numinp, char *inps[]){
820        glt->ResetCounters();        glt->ResetCounters();
821        printf(" CalorimeterLevel2 called\n");        printf(" CalorimeterLevel2 called\n");
822        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n");
823          if ( !reprocess ){
824            string cstring = "--delete-calo-tree";
825            calargv[calargc] = (char *)cstring.c_str();
826          } else {
827            calargc--;
828          }
829        CALSGN = CaloCore(run,processFile,glt,calargc,calargv);        CALSGN = CaloCore(run,processFile,glt,calargc,calargv);
830        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end CaloCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");        if ( debug ) printf("                      >>>>>>>>>>>>>>>>>>>|  end CaloCore STDOUT  |<<<<<<<<<<<<<<<<<<<\n");
831        NQCAL = glt->GetNqueries();        NQCAL = glt->GetNqueries();
# Line 810  int main(int numinp, char *inps[]){ Line 886  int main(int numinp, char *inps[]){
886      orbtime = timer.RealTime();      orbtime = timer.RealTime();
887      corbtime = timer.CpuTime();      corbtime = timer.CpuTime();
888      //      //
889      NQTOT = 2 + NQRUN + NQTRK + NQCAL + NQND + NQS4 + NQAC + NQTRG + NQTOF + NQORB;      NQTOT = 2 + NQRUN + NQTRK + NQCAL + NQND + NQS4 + NQAC + NQTRG + NQTOF + NQORB + NQCAL1 + NQTOF1;
890      //      //
891      delete glt;      delete glt;
892      delete runinfo;      delete runinfo;
# Line 842  int main(int numinp, char *inps[]){ Line 918  int main(int numinp, char *inps[]){
918      case -34:  message += " No S4 in the acquisition"; break;      case -34:  message += " No S4 in the acquisition"; break;
919      case -35:  message += " No ND in the acquisition"; break;      case -35:  message += " No ND in the acquisition"; break;
920      case -36:  message += " I/O error or missing entry"; break;      case -36:  message += " I/O error or missing entry"; break;
921        case -37:  message += " Cannot override output file, see help and use -r option"; break;
922        
923        //        //
924      case -50:  message += " GLTABLES - No entries matching GL_RUN query"; break;      case -50:  message += " GLTABLES - No entries matching GL_RUN query"; break;
# Line 956  int main(int numinp, char *inps[]){ Line 1033  int main(int numinp, char *inps[]){
1033    if ( DVSGN & (1 << 4) ) printf(" - No S4 in the acquisition\n");    if ( DVSGN & (1 << 4) ) printf(" - No S4 in the acquisition\n");
1034    if ( DVSGN & (1 << 5) ) printf(" - No neutron detector in the acquisition\n");    if ( DVSGN & (1 << 5) ) printf(" - No neutron detector in the acquisition\n");
1035    //    //
1036      // CaloCore, first call:
1037      //
1038      switch(CALSGN1){
1039      case 100: printf("\n WARNING CALORIMETER - Data with no associated calibration\n");
1040      case 101: printf("\n WARNING CALORIMETER - No tracks or good events in this run\n");
1041      }
1042      //
1043    // CaloCore:    // CaloCore:
1044    //    //
1045    switch(CALSGN){    switch(CALSGN){
# Line 996  int main(int numinp, char *inps[]){ Line 1080  int main(int numinp, char *inps[]){
1080      Float_t trkrt = 0.;      Float_t trkrt = 0.;
1081      Float_t calrt = 0.;      Float_t calrt = 0.;
1082      Float_t tofrt = 0.;      Float_t tofrt = 0.;
1083        Float_t calrt1 = 0.;
1084        Float_t tofrt1 = 0.;
1085      Float_t trgrt = 0.;      Float_t trgrt = 0.;
1086      Float_t acrt = 0.;      Float_t acrt = 0.;
1087      Float_t s4rt = 0.;      Float_t s4rt = 0.;
# Line 1005  int main(int numinp, char *inps[]){ Line 1091  int main(int numinp, char *inps[]){
1091      if ( RUN ) runrt = nruns/runtime;      if ( RUN ) runrt = nruns/runtime;
1092      if ( TRK ) trkrt = (nevents+1)/trktime;      if ( TRK ) trkrt = (nevents+1)/trktime;
1093      if ( CAL ) calrt = (nevents+1)/caltime;      if ( CAL ) calrt = (nevents+1)/caltime;
1094        if ( TOF && newfit ) tofrt1 = (nevents+1)/toftime1;
1095        if ( CAL && newfit ) calrt1 = (nevents+1)/caltime1;
1096      if ( TOF ) tofrt = (nevents+1)/toftime;      if ( TOF ) tofrt = (nevents+1)/toftime;
1097      if ( TRG ) trgrt = (nevents+1)/trgtime;      if ( TRG ) trgrt = (nevents+1)/trgtime;
1098      if ( AC ) acrt = (nevents+1)/actime;      if ( AC ) acrt = (nevents+1)/actime;
# Line 1018  int main(int numinp, char *inps[]){ Line 1106  int main(int numinp, char *inps[]){
1106      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");
1107      printf("###########################################################################################\n");      printf("###########################################################################################\n");
1108      printf("# RUN      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",RUN,NQRUN,fabs(runtime),fabs(cruntime),runrt);      printf("# RUN      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",RUN,NQRUN,fabs(runtime),fabs(cruntime),runrt);
1109        printf("# TRG      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",TRG,NQTRG,fabs(trgtime),fabs(ctrgtime),trgrt);
1110        if ( newfit ){
1111        printf("# TOF L1   #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",TOF,NQTOF1,fabs(toftime1),fabs(ctoftime1),tofrt1);    
1112        printf("# CAL L1   #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",CAL,NQCAL1,fabs(caltime1),fabs(ccaltime1),calrt1);
1113        }
1114      printf("# TRK      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",TRK,NQTRK,fabs(trktime),fabs(ctrktime),trkrt);      printf("# TRK      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",TRK,NQTRK,fabs(trktime),fabs(ctrktime),trkrt);
     printf("# CAL      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",CAL,NQCAL,fabs(caltime),fabs(ccaltime),calrt);  
1115      printf("# TOF      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",TOF,NQTOF,fabs(toftime),fabs(ctoftime),tofrt);      printf("# TOF      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",TOF,NQTOF,fabs(toftime),fabs(ctoftime),tofrt);
1116      printf("# TRG      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",TRG,NQTRG,fabs(trgtime),fabs(ctrgtime),trgrt);      printf("# CAL      #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",CAL,NQCAL,fabs(caltime),fabs(ccaltime),calrt);
1117      printf("# AC       #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",AC,NQAC,fabs(actime),fabs(cactime),acrt);      printf("# AC       #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",AC,NQAC,fabs(actime),fabs(cactime),acrt);
1118      printf("# S4       #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",S4,NQS4,fabs(s4time),fabs(cs4time),s4rt);      printf("# S4       #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",S4,NQS4,fabs(s4time),fabs(cs4time),s4rt);
1119      printf("# ND       #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",ND,NQND,fabs(ndtime),fabs(cndtime),ndrt);      printf("# ND       #         %i           #   %3u    #      %8.2f   #    %8.2f  #  %8.2f  #\n",ND,NQND,fabs(ndtime),fabs(cndtime),ndrt);
# Line 1041  int main(int numinp, char *inps[]){ Line 1133  int main(int numinp, char *inps[]){
1133      gSystem->Unlink(filename.Data());      gSystem->Unlink(filename.Data());
1134    };    };
1135    //    //
1136    if ( !error && (CALSGN || TRKSGN || TRGSGN || TOFSGN || RUNSGN || ORBSGN || ACSGN || S4SGN || NDSGN || DVSGN) ) error = 1;    if ( !error && (CALSGN1 || TOFSGN1 || CALSGN || TRKSGN || TRGSGN || TOFSGN || RUNSGN || ORBSGN || ACSGN || S4SGN || NDSGN || DVSGN) ) error = 1;
1137    //    //
1138    printf("\n Finished, exiting...\n\n");    printf("\n Finished, exiting...\n\n");
1139    //    //
Legend:
Removed from v.1.46  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.23