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"); |
70 |
printf(" Standard call:\n DarthVader -idRun 1085 \n"); |
printf(" Standard call:\n DarthVader -idRun 1085 \n"); |
71 |
printf(" Process only RunInfo and Tracker (be verbose for tracker):\n DarthVader -idRun 1085 -all +RUN +TRK [ --verbose ] \n"); |
printf(" Process only RunInfo and Tracker (be verbose for tracker):\n DarthVader -idRun 1085 -all +RUN +TRK [ --verbose ] \n"); |
72 |
printf(" Process all and be verbose for calorimeter:\n DarthVader -idRun 1085 +CAL [ --verbose ] \n\n"); |
printf(" Process all and be verbose for calorimeter:\n DarthVader -idRun 1085 +CAL [ --verbose ] \n\n"); |
73 |
|
printf(" Set Tracker to work with simulated data:\n DarthVader -idRun 1085 +TRK [ --simu ] \n\n"); |
74 |
}; |
}; |
75 |
|
|
76 |
// |
// |
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.; |
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; |
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; |
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; |
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; |
586 |
// |
// |
587 |
// Start: |
// Start: |
588 |
// |
// |
589 |
printf("\n Welcome to the PAMELA LEVEL2 flight software, version %s \n\n",version.Data()); |
printf("\n Welcome to the PAMELA LEVEL2 flight software, version %s (The Ultimate Warrior)\n\n",version.Data()); |
590 |
if ( run ) printf("\n Processing run number %u \n\n",run); |
if ( run ) printf("\n Processing run number %u \n\n",run); |
591 |
// |
// |
592 |
// Connect to the DB |
// Connect to the DB |
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 |
// |
// |
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 ) { |
762 |
|
glt->ResetCounters(); |
763 |
|
printf(" ToFLevel2 called, 1st call\n"); |
764 |
|
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
765 |
|
string tstring = "--level1-only"; |
766 |
|
tofargv[tofargc] = (char *)tstring.c_str(); |
767 |
|
tofargc++; |
768 |
|
TOFSGN1 = ToFCore(run,processFile,glt,tofargc,tofargv); |
769 |
|
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
770 |
|
NQTOF1 = glt->GetNqueries(); |
771 |
|
}; |
772 |
|
timer.Stop(); |
773 |
|
toftime1 = timer.RealTime(); |
774 |
|
ctoftime1 = timer.CpuTime(); |
775 |
|
// |
776 |
|
timer.Start(kTRUE); |
777 |
|
if ( CAL ) { |
778 |
|
glt->ResetCounters(); |
779 |
|
printf(" CalorimeterLevel2 called, 1st call\n"); |
780 |
|
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
781 |
|
string cstring = "--level1-only"; |
782 |
|
calargv[calargc] = (char *)cstring.c_str(); |
783 |
|
calargc++; |
784 |
|
CALSGN1 = CaloCore(run,processFile,glt,calargc,calargv); |
785 |
|
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
786 |
|
NQCAL1 = glt->GetNqueries(); |
787 |
|
}; |
788 |
|
timer.Stop(); |
789 |
|
caltime1 = timer.RealTime(); |
790 |
|
ccaltime1 = timer.CpuTime(); |
791 |
|
} |
792 |
|
// |
793 |
|
// |
794 |
timer.Start(kTRUE); |
timer.Start(kTRUE); |
795 |
if ( TRK ) { |
if ( TRK ) { |
796 |
glt->ResetCounters(); |
glt->ResetCounters(); |
810 |
glt->ResetCounters(); |
glt->ResetCounters(); |
811 |
printf(" ToFLevel2 called\n"); |
printf(" ToFLevel2 called\n"); |
812 |
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
813 |
|
if ( !reprocess ){ |
814 |
|
string tstring = "--delete-tree"; |
815 |
|
tofargv[tofargc] = (char *)tstring.c_str(); |
816 |
|
} else { |
817 |
|
tofargc--; |
818 |
|
} |
819 |
TOFSGN = ToFCore(run,processFile,glt,tofargc,tofargv); |
TOFSGN = ToFCore(run,processFile,glt,tofargc,tofargv); |
820 |
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end ToFCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
821 |
NQTOF = glt->GetNqueries(); |
NQTOF = glt->GetNqueries(); |
829 |
glt->ResetCounters(); |
glt->ResetCounters(); |
830 |
printf(" CalorimeterLevel2 called\n"); |
printf(" CalorimeterLevel2 called\n"); |
831 |
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| start CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
832 |
|
if ( !reprocess ){ |
833 |
|
string cstring = "--delete-calo-tree"; |
834 |
|
calargv[calargc] = (char *)cstring.c_str(); |
835 |
|
} else { |
836 |
|
calargc--; |
837 |
|
} |
838 |
CALSGN = CaloCore(run,processFile,glt,calargc,calargv); |
CALSGN = CaloCore(run,processFile,glt,calargc,calargv); |
839 |
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
if ( debug ) printf(" >>>>>>>>>>>>>>>>>>>| end CaloCore STDOUT |<<<<<<<<<<<<<<<<<<<\n"); |
840 |
NQCAL = glt->GetNqueries(); |
NQCAL = glt->GetNqueries(); |
895 |
orbtime = timer.RealTime(); |
orbtime = timer.RealTime(); |
896 |
corbtime = timer.CpuTime(); |
corbtime = timer.CpuTime(); |
897 |
// |
// |
898 |
NQTOT = 2 + NQRUN + NQTRK + NQCAL + NQND + NQS4 + NQAC + NQTRG + NQTOF + NQORB; |
NQTOT = 2 + NQRUN + NQTRK + NQCAL + NQND + NQS4 + NQAC + NQTRG + NQTOF + NQORB + NQCAL1 + NQTOF1; |
899 |
// |
// |
900 |
delete glt; |
delete glt; |
901 |
delete runinfo; |
delete runinfo; |
927 |
case -34: message += " No S4 in the acquisition"; break; |
case -34: message += " No S4 in the acquisition"; break; |
928 |
case -35: message += " No ND in the acquisition"; break; |
case -35: message += " No ND in the acquisition"; break; |
929 |
case -36: message += " I/O error or missing entry"; break; |
case -36: message += " I/O error or missing entry"; break; |
930 |
|
case -37: message += " Cannot override output file, see help and use -r option"; break; |
931 |
|
|
932 |
// |
// |
933 |
case -50: message += " GLTABLES - No entries matching GL_RUN query"; break; |
case -50: message += " GLTABLES - No entries matching GL_RUN query"; break; |
972 |
case -203: message += " TRACKERLEVEL2 - No Tracker branch in LEVEL0 Physics tree"; break; |
case -203: message += " TRACKERLEVEL2 - No Tracker branch in LEVEL0 Physics tree"; break; |
973 |
case -204: message += " TRACKERLEVEL2 - No reprocessing implemented for LEVEL1 output"; break; |
case -204: message += " TRACKERLEVEL2 - No reprocessing implemented for LEVEL1 output"; break; |
974 |
case -205: message += " TRACKERLEVEL2 - Error accessing RunInfo "; break; |
case -205: message += " TRACKERLEVEL2 - Error accessing RunInfo "; break; |
975 |
|
case -206: message += " TRACKERLEVEL2 - Extended-tracking mode - no Calorimeter and/or ToF trees found"; break; |
976 |
case -210: message += " TRACKERLEVEL2 - Error opening/reading trk mask GL_PARAM parameters "; break; |
case -210: message += " TRACKERLEVEL2 - Error opening/reading trk mask GL_PARAM parameters "; break; |
977 |
case -211: message += " TRACKERLEVEL2 - Error opening/reading trk alignment GL_PARAM parameters"; break; |
case -211: message += " TRACKERLEVEL2 - Error opening/reading trk alignment GL_PARAM parameters"; break; |
978 |
case -212: message += " TRACKERLEVEL2 - Error opening/reading trk mip GL_PARAM parameters"; break; |
case -212: message += " TRACKERLEVEL2 - Error opening/reading trk mip GL_PARAM parameters"; break; |
1025 |
case -807: message += " RUNINFO - DB connection problems"; break; |
case -807: message += " RUNINFO - DB connection problems"; break; |
1026 |
// |
// |
1027 |
case -900: message += " OrbitalInfo - no ToF tree"; break; |
case -900: message += " OrbitalInfo - no ToF tree"; break; |
1028 |
|
case -901: message += " OrbitalInfo - Mismatch between tracker and tof tree branches concerning extended tracking algorithm(s)"; break; |
1029 |
// |
// |
1030 |
default: message += "Unidentified error or warning"; break; |
default: message += "Unidentified error or warning"; break; |
1031 |
}; |
}; |
1044 |
if ( DVSGN & (1 << 4) ) printf(" - No S4 in the acquisition\n"); |
if ( DVSGN & (1 << 4) ) printf(" - No S4 in the acquisition\n"); |
1045 |
if ( DVSGN & (1 << 5) ) printf(" - No neutron detector in the acquisition\n"); |
if ( DVSGN & (1 << 5) ) printf(" - No neutron detector in the acquisition\n"); |
1046 |
// |
// |
1047 |
|
// CaloCore, first call: |
1048 |
|
// |
1049 |
|
switch(CALSGN1){ |
1050 |
|
case 100: printf("\n WARNING CALORIMETER - Data with no associated calibration\n"); |
1051 |
|
case 101: printf("\n WARNING CALORIMETER - No tracks or good events in this run\n"); |
1052 |
|
} |
1053 |
|
// |
1054 |
// CaloCore: |
// CaloCore: |
1055 |
// |
// |
1056 |
switch(CALSGN){ |
switch(CALSGN){ |
1077 |
if ( processFile ){ |
if ( processFile ){ |
1078 |
processFile->cd(); |
processFile->cd(); |
1079 |
processFile->WriteStreamerInfo(); |
processFile->WriteStreamerInfo(); |
1080 |
processFile->ls(); |
if ( debug ) processFile->ls(); |
1081 |
processFile->Close(); |
processFile->Close(); |
1082 |
}; |
}; |
1083 |
// |
// |
1091 |
Float_t trkrt = 0.; |
Float_t trkrt = 0.; |
1092 |
Float_t calrt = 0.; |
Float_t calrt = 0.; |
1093 |
Float_t tofrt = 0.; |
Float_t tofrt = 0.; |
1094 |
|
Float_t calrt1 = 0.; |
1095 |
|
Float_t tofrt1 = 0.; |
1096 |
Float_t trgrt = 0.; |
Float_t trgrt = 0.; |
1097 |
Float_t acrt = 0.; |
Float_t acrt = 0.; |
1098 |
Float_t s4rt = 0.; |
Float_t s4rt = 0.; |
1102 |
if ( RUN ) runrt = nruns/runtime; |
if ( RUN ) runrt = nruns/runtime; |
1103 |
if ( TRK ) trkrt = (nevents+1)/trktime; |
if ( TRK ) trkrt = (nevents+1)/trktime; |
1104 |
if ( CAL ) calrt = (nevents+1)/caltime; |
if ( CAL ) calrt = (nevents+1)/caltime; |
1105 |
|
if ( TOF && newfit ) tofrt1 = (nevents+1)/toftime1; |
1106 |
|
if ( CAL && newfit ) calrt1 = (nevents+1)/caltime1; |
1107 |
if ( TOF ) tofrt = (nevents+1)/toftime; |
if ( TOF ) tofrt = (nevents+1)/toftime; |
1108 |
if ( TRG ) trgrt = (nevents+1)/trgtime; |
if ( TRG ) trgrt = (nevents+1)/trgtime; |
1109 |
if ( AC ) acrt = (nevents+1)/actime; |
if ( AC ) acrt = (nevents+1)/actime; |
1117 |
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"); |
1118 |
printf("###########################################################################################\n"); |
printf("###########################################################################################\n"); |
1119 |
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); |
1120 |
|
printf("# TRG # %i # %3u # %8.2f # %8.2f # %8.2f #\n",TRG,NQTRG,fabs(trgtime),fabs(ctrgtime),trgrt); |
1121 |
|
if ( newfit ){ |
1122 |
|
printf("# TOF L1 # %i # %3u # %8.2f # %8.2f # %8.2f #\n",TOF,NQTOF1,fabs(toftime1),fabs(ctoftime1),tofrt1); |
1123 |
|
printf("# CAL L1 # %i # %3u # %8.2f # %8.2f # %8.2f #\n",CAL,NQCAL1,fabs(caltime1),fabs(ccaltime1),calrt1); |
1124 |
|
} |
1125 |
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); |
|
1126 |
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); |
1127 |
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); |
1128 |
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); |
1129 |
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); |
1130 |
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); |
1144 |
gSystem->Unlink(filename.Data()); |
gSystem->Unlink(filename.Data()); |
1145 |
}; |
}; |
1146 |
// |
// |
1147 |
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; |
1148 |
// |
// |
1149 |
printf("\n Finished, exiting...\n\n"); |
printf("\n Finished, exiting...\n\n"); |
1150 |
// |
// |