1 |
// |
2 |
// C/C++ headers |
3 |
// |
4 |
#include <iostream> |
5 |
#include <iomanip> |
6 |
#include <fstream> |
7 |
#include <sstream> |
8 |
// |
9 |
// ROOT headers |
10 |
// |
11 |
#include <TString.h> |
12 |
#include <TTimeStamp.h> |
13 |
#include <TSQLServer.h> |
14 |
#include <TFile.h> |
15 |
#include <TSystem.h> |
16 |
// |
17 |
// Detector's package headers |
18 |
// |
19 |
#include <GLTables.h> |
20 |
// |
21 |
using namespace std; |
22 |
// |
23 |
// |
24 |
// |
25 |
#include <YodaProfilerVerl2.h> |
26 |
// |
27 |
// Usage subroutine |
28 |
// |
29 |
void r2d2usage(){ |
30 |
printf("\nUsage:\n"); |
31 |
printf("\n R2-D2 [ options ]\n"); |
32 |
printf("\n Options are:\n\n"); |
33 |
printf(" --version print informations about compilation and exit\n"); |
34 |
printf(" -h | --help print this help and exit \n"); |
35 |
printf(" -v | --verbose be verbose [default]\n"); |
36 |
printf(" -n | --neat used with '-filename' returns only the run ID number\n"); |
37 |
printf(" -idRun run ID_RUN: ID number of the run \n"); |
38 |
printf(" -filename file output yoda filename \n"); |
39 |
printf(" -l2filename file output amidala filename \n"); |
40 |
printf(" -host name of the DB host [default = $PAM_DBHOST or mysql://localhost/pamelaprod]\n"); |
41 |
printf(" -user username for the DB connection [default = $PAM_DBUSER or \"anonymous\"] \n"); |
42 |
printf(" -psw password for the DB connection [default = $PAM_DBPSW or \"\"]\n"); |
43 |
printf(" -deps file shows dependencies of the given file\n"); |
44 |
printf(" -splitat file shows dependenices of the given file and the following files (given by -nsplit)\n"); |
45 |
printf(" -nsplit number number of files to be shown by -splitat [default 20]\n"); |
46 |
printf(" -tzone timezone the time zone: UTC,GMT,MSK,MSD,CET,CEST are accepted \n"); |
47 |
printf(" -convert dbtime convert the dbtime given in seconds (from the DB) to a string\n"); |
48 |
printf(" -runat date returns run number which contains the given date,\n"); |
49 |
printf(" for date use the SQL format \"yyyy-mm-dd hh:mm:ss\" \n"); |
50 |
printf(" -runatDB time returns run number which contains the given DB time\n"); |
51 |
printf(" -tsfile file yoda filename for the time sync (to be used with -obt)\n"); |
52 |
printf(" -obt OBT OBT in ms returns a date (to be used with -tsfile)\n"); |
53 |
printf(" -getRTime OBT OBT in ms returns Resurs time (to be used with -tsfile)\n"); |
54 |
printf(" -dumpTLEfor date save into file tle.txt the TLE for the given date,\n"); |
55 |
printf(" for date use the SQL format \"yyyy-mm-dd hh:mm:ss\" \n"); |
56 |
printf("\nExamples: \n"); |
57 |
printf(" R2-D2 -idRun 1085 \n"); |
58 |
printf(" R2-D2 -filename DW_050208_00900.root \n"); |
59 |
printf(" R2-D2 -idRun 1085 -filename DW_050208_00900.root \n\n"); |
60 |
}; |
61 |
// |
62 |
// Here the main |
63 |
// |
64 |
int main(int numinp, char *inps[]){ |
65 |
// |
66 |
// Variables booking |
67 |
// |
68 |
TString message; |
69 |
Int_t error = 0; |
70 |
// |
71 |
UInt_t run = 0ULL; |
72 |
// |
73 |
TString filename = ""; |
74 |
TString l2filename = ""; |
75 |
TString splitat = ""; |
76 |
UInt_t nsplit = 20; |
77 |
// |
78 |
TSQLServer *dbc = 0; |
79 |
TString host = "mysql://localhost/pamelaprod"; |
80 |
TString user = "anonymous"; |
81 |
TString psw = ""; |
82 |
// |
83 |
const char *pamdbhost=gSystem->Getenv("PAM_DBHOST"); |
84 |
const char *pamdbuser=gSystem->Getenv("PAM_DBUSER"); |
85 |
const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW"); |
86 |
if ( !pamdbhost ) pamdbhost = ""; |
87 |
if ( !pamdbuser ) pamdbuser = ""; |
88 |
if ( !pamdbpsw ) pamdbpsw = ""; |
89 |
if ( strcmp(pamdbhost,"") ) host = pamdbhost; |
90 |
if ( strcmp(pamdbuser,"") ) user = pamdbuser; |
91 |
if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw; |
92 |
// |
93 |
// printf(" host %s user %s psw %s \n",host.Data(),user.Data(),psw.Data()); |
94 |
// |
95 |
TString tzone = "UTC"; |
96 |
TString runtime = "1970-01-01 00:00:00"; |
97 |
TString tletime = "1970-01-01 00:00:00"; |
98 |
UInt_t dbti = 0; |
99 |
Bool_t convert = false; |
100 |
Bool_t convres = false; |
101 |
Bool_t ruti = false; |
102 |
Bool_t dtle = false; |
103 |
Bool_t ruti2 = false; |
104 |
Bool_t convobt = false; |
105 |
Bool_t convobtts = false; |
106 |
Bool_t neat = false; |
107 |
// |
108 |
UInt_t runtime2 = 0; |
109 |
UInt_t obt = 0; |
110 |
UInt_t res = 0; |
111 |
TString tsfile = ""; |
112 |
// |
113 |
TSQLResult *pResult; |
114 |
TSQLRow *Row; |
115 |
TSQLResult *pResult2 = 0; |
116 |
TSQLRow *Row2; |
117 |
int t; |
118 |
int r; |
119 |
stringstream myquery; |
120 |
// |
121 |
// Checking input parameters |
122 |
// |
123 |
Int_t i = 0; |
124 |
if ( numinp > 1 ){ |
125 |
while ( i < numinp ){ |
126 |
if ( !strcmp(inps[i],"--version") ){ |
127 |
YodaProfilerInfo(true); |
128 |
exit(0); |
129 |
}; |
130 |
if ( !strcmp(inps[i],"-h") || !strcmp(inps[i],"--help") ){ |
131 |
r2d2usage(); |
132 |
exit(0); |
133 |
}; |
134 |
if ( !strcmp(inps[i],"-n") || !strcmp(inps[i],"--neat") ){ |
135 |
neat = true; |
136 |
}; |
137 |
if ( !strcmp(inps[i],"-idRun") ) { |
138 |
if ( numinp-1 < i+1 ) { |
139 |
r2d2usage(); |
140 |
exit(-3); |
141 |
}; |
142 |
run = (UInt_t)atoll(inps[i+1]); |
143 |
}; |
144 |
if ( !strcmp(inps[i],"-filename") ) { |
145 |
if ( numinp-1 < i+1 ){ |
146 |
r2d2usage(); |
147 |
exit(-3); |
148 |
}; |
149 |
filename = (TString)inps[i+1]; |
150 |
}; |
151 |
if ( !strcmp(inps[i],"-splitat") ) { |
152 |
if ( numinp-1 < i+1 ){ |
153 |
r2d2usage(); |
154 |
exit(-3); |
155 |
}; |
156 |
splitat = (TString)inps[i+1]; |
157 |
}; |
158 |
// |
159 |
if ( !strcmp(inps[i],"-nsplit") ) { |
160 |
if ( numinp-1 < i+1 ){ |
161 |
r2d2usage(); |
162 |
exit(-3); |
163 |
}; |
164 |
nsplit = (UInt_t)atoll(inps[i+1]); |
165 |
}; |
166 |
if ( !strcmp(inps[i],"-deps") ) { |
167 |
if ( numinp-1 < i+1 ){ |
168 |
r2d2usage(); |
169 |
exit(-3); |
170 |
}; |
171 |
nsplit = 1; |
172 |
splitat = (TString)inps[i+1]; |
173 |
}; |
174 |
if ( !strcmp(inps[i],"-l2filename") ) { |
175 |
if ( numinp-1 < i+1 ){ |
176 |
r2d2usage(); |
177 |
exit(-3); |
178 |
}; |
179 |
l2filename = (TString)inps[i+1]; |
180 |
}; |
181 |
if ( !strcmp(inps[i],"-host") ) { |
182 |
if ( numinp-1 < i+1 ){ |
183 |
r2d2usage(); |
184 |
exit(-3); |
185 |
}; |
186 |
host = (TString)inps[i+1]; |
187 |
}; |
188 |
if ( !strcmp(inps[i],"-user") ) { |
189 |
if ( numinp-1 < i+1 ){ |
190 |
r2d2usage(); |
191 |
exit(-3); |
192 |
}; |
193 |
user = (TString)inps[i+1]; |
194 |
}; |
195 |
if ( !strcmp(inps[i],"-psw") ) { |
196 |
if ( numinp-1 < i+1 ){ |
197 |
r2d2usage(); |
198 |
exit(-3); |
199 |
}; |
200 |
psw = (TString)inps[i+1]; |
201 |
}; |
202 |
// |
203 |
if ( !strcmp(inps[i],"-tsfile") ) { |
204 |
convobtts = true; |
205 |
if ( numinp-1 < i+1 ){ |
206 |
r2d2usage(); |
207 |
exit(-3); |
208 |
}; |
209 |
tsfile = (TString)inps[i+1]; |
210 |
}; |
211 |
// |
212 |
if ( !strcmp(inps[i],"-obt") ) { |
213 |
convobt = true; |
214 |
if ( numinp-1 < i+1 ){ |
215 |
r2d2usage(); |
216 |
exit(-3); |
217 |
}; |
218 |
obt = (UInt_t)atoll(inps[i+1]); |
219 |
}; |
220 |
// |
221 |
if ( !strcmp(inps[i],"-getRTime") ) { |
222 |
convres = true; |
223 |
if ( numinp-1 < i+1 ){ |
224 |
r2d2usage(); |
225 |
exit(-3); |
226 |
}; |
227 |
res = (UInt_t)atoll(inps[i+1]); |
228 |
}; |
229 |
// |
230 |
// |
231 |
if ( !strcmp(inps[i],"-tzone") ) { |
232 |
if ( numinp-1 < i+1 ){ |
233 |
r2d2usage(); |
234 |
exit(-3); |
235 |
}; |
236 |
tzone = (TString)inps[i+1]; |
237 |
}; |
238 |
// |
239 |
if ( !strcmp(inps[i],"-convert") ) { |
240 |
convert = true; |
241 |
if ( numinp-1 < i+1 ){ |
242 |
r2d2usage(); |
243 |
exit(-3); |
244 |
}; |
245 |
dbti = (UInt_t)atoll(inps[i+1]); |
246 |
}; |
247 |
// |
248 |
if ( !strcmp(inps[i],"-runat") ) { |
249 |
ruti = true; |
250 |
if ( numinp-1 < i+1 ){ |
251 |
r2d2usage(); |
252 |
exit(-3); |
253 |
}; |
254 |
runtime = (TString)inps[i+1]; |
255 |
}; |
256 |
// |
257 |
if ( !strcmp(inps[i],"-dumpTLEfor") ) { |
258 |
dtle = true; |
259 |
if ( numinp-1 < i+1 ){ |
260 |
r2d2usage(); |
261 |
exit(-3); |
262 |
}; |
263 |
tletime = (TString)inps[i+1]; |
264 |
}; |
265 |
// |
266 |
// |
267 |
if ( !strcmp(inps[i],"-runatDB") ) { |
268 |
ruti2 = true; |
269 |
if ( numinp-1 < i+1 ){ |
270 |
r2d2usage(); |
271 |
exit(-3); |
272 |
}; |
273 |
runtime2 = (UInt_t)atoll(inps[i+1]); |
274 |
}; |
275 |
// |
276 |
i++; |
277 |
}; |
278 |
// |
279 |
} else { |
280 |
// |
281 |
// no input parameters exit with error, we need at least the run id. |
282 |
// |
283 |
r2d2usage(); |
284 |
exit(-2); |
285 |
}; |
286 |
// |
287 |
// Connect to the DB |
288 |
// |
289 |
dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data()); |
290 |
if( !dbc ) throw -2; |
291 |
// |
292 |
bool connect = dbc->IsConnected(); |
293 |
// |
294 |
if( !connect ){ |
295 |
printf(" Error, not connected to DB\n"); |
296 |
exit(-1); |
297 |
}; |
298 |
// |
299 |
myquery.str(""); |
300 |
myquery << "SET time_zone='+0:00';"; |
301 |
delete dbc->Query(myquery.str().c_str()); |
302 |
delete dbc->Query("SET sql_mode = 'NO_UNSIGNED_SUBTRACTION';"); |
303 |
// |
304 |
GL_ROOT *glroot = new GL_ROOT(); |
305 |
GL_RUN *glrun = new GL_RUN(); |
306 |
GL_TIMESYNC *dbtime = new GL_TIMESYNC(); |
307 |
UInt_t nr = 0; |
308 |
UInt_t rlist[500]; |
309 |
// |
310 |
// At which date correspond the DB time "dbti"? |
311 |
// |
312 |
if ( convert ){ |
313 |
printf("\n DB time %u is %s %s \n",dbti,dbtime->ConvertTime(tzone,dbti).Data(),tzone.Data()); |
314 |
}; |
315 |
// |
316 |
// convert OBT to date |
317 |
// |
318 |
if ( convobt && convobtts ){ |
319 |
UInt_t id = 0; |
320 |
myquery.str(""); |
321 |
myquery << "select "; |
322 |
myquery << " ID"; |
323 |
myquery << " from GL_ROOT where NAME=\"" << tsfile.Data() << "\";"; |
324 |
pResult = dbc->Query(myquery.str().c_str()); |
325 |
if ( pResult ){ |
326 |
Row = pResult->Next(); |
327 |
if ( Row ){ |
328 |
id = (UInt_t)atoll(Row->GetField(0)); |
329 |
delete pResult; |
330 |
GL_TIMESYNC *ctime = new GL_TIMESYNC(id,"ID",dbc); |
331 |
UInt_t abtime = ctime->DBabsTime(obt); |
332 |
TString thetime = dbtime->ConvertTime(tzone,abtime); |
333 |
printf("\n OBT %u in the file %s corresponds to DBtime %u and date %s %s \n",obt,tsfile.Data(),abtime,thetime.Data(),tzone.Data()); |
334 |
delete ctime; |
335 |
}; |
336 |
}; |
337 |
}; |
338 |
if ( (convobt && !convobtts) ){ |
339 |
printf("\n To convert a OBT to a date you must provide both OBT and file to be used for time sync \n"); |
340 |
}; |
341 |
// |
342 |
// convert OBT to Resurs seconds |
343 |
// |
344 |
if ( convres && convobtts ){ |
345 |
UInt_t id = 0; |
346 |
myquery.str(""); |
347 |
myquery << "select "; |
348 |
myquery << " ID"; |
349 |
myquery << " from GL_ROOT where NAME=\"" << tsfile.Data() << "\";"; |
350 |
pResult = dbc->Query(myquery.str().c_str()); |
351 |
if ( pResult ){ |
352 |
Row = pResult->Next(); |
353 |
if ( Row ){ |
354 |
id = (UInt_t)atoll(Row->GetField(0)); |
355 |
delete pResult; |
356 |
GL_TIMESYNC *ctime = new GL_TIMESYNC(id,"ID",dbc); |
357 |
UInt_t restime = ctime->ResursTime(res); |
358 |
//UInt_t abtime = ctime->DBabsTime(obt); |
359 |
// |
360 |
// TString thetime = dbtime->ConvertTime(tzone,abtime); |
361 |
printf("\n OBT %u in the file %s corresponds to Resurs time %u \n",res,tsfile.Data(),restime); |
362 |
delete ctime; |
363 |
}; |
364 |
}; |
365 |
}; |
366 |
if ( (convres && !convobtts) ){ |
367 |
printf("\n To convert a OBT to the Resurs time you must provide both OBT and file to be used for time sync \n"); |
368 |
}; |
369 |
// |
370 |
// Which run contains the date "runtime"? |
371 |
// |
372 |
if ( ruti ){ |
373 |
// |
374 |
TDatime ti = TDatime(runtime.Data()); |
375 |
// |
376 |
TTimeStamp *time = new TTimeStamp((UInt_t)ti.GetYear(),(UInt_t)ti.GetMonth(),(UInt_t)ti.GetDay(),(UInt_t)ti.GetHour(),(UInt_t)ti.GetMinute(),(UInt_t)ti.GetSecond(),0,true,0); |
377 |
// |
378 |
UInt_t dbti = time->GetSec(); |
379 |
// |
380 |
TString thetime = dbtime->UnConvertTime(tzone,dbti); |
381 |
// |
382 |
ti = TDatime(thetime.Data()); |
383 |
TTimeStamp *time2 = new TTimeStamp((UInt_t)ti.GetYear(),(UInt_t)ti.GetMonth(),(UInt_t)ti.GetDay(),(UInt_t)ti.GetHour(),(UInt_t)ti.GetMinute(),(UInt_t)ti.GetSecond(),0,true,0); |
384 |
// |
385 |
UInt_t mytime = time2->GetSec(); |
386 |
// |
387 |
Bool_t found = false; |
388 |
// |
389 |
myquery.str(""); |
390 |
myquery << "select "; |
391 |
myquery << " ID "; |
392 |
myquery << " from GL_RUN where RUNHEADER_TIME<=" << mytime << " AND " |
393 |
<< " RUNTRAILER_TIME>=" << mytime << " ;"; |
394 |
// printf("myquery is %s \n",myquery.str().c_str()); |
395 |
pResult = dbc->Query(myquery.str().c_str()); |
396 |
for( r=0; r < 1000; r++){ |
397 |
Row = pResult->Next(); |
398 |
if( Row == NULL ) break; |
399 |
found = true; |
400 |
printf("\n Date %s %s (DB time %u ) is contained in run %u \n",runtime.Data(),tzone.Data(),mytime,(UInt_t)atoll(Row->GetField(0))); |
401 |
}; |
402 |
myquery.str(""); |
403 |
myquery << "select "; |
404 |
myquery << " ID "; |
405 |
myquery << " from GL_RUN_TRASH where BELONGED_TO='GL_RUN' AND RUNHEADER_TIME<=" << mytime << " AND " |
406 |
<< " RUNTRAILER_TIME>=" << mytime << " ;"; |
407 |
// printf("myquery is %s \n",myquery.str().c_str()); |
408 |
pResult = dbc->Query(myquery.str().c_str()); |
409 |
for( r=0; r < 1000; r++){ |
410 |
Row = pResult->Next(); |
411 |
if( Row == NULL ) break; |
412 |
printf("\n Date %s %s (DB time %u ) is contained in run %u in the GL_RUN_TRASH table \n",runtime.Data(),tzone.Data(),mytime,(UInt_t)atoll(Row->GetField(0))); |
413 |
found = true; |
414 |
}; |
415 |
// |
416 |
if ( !found ){ |
417 |
printf("\n No run contains date %s %s (DB time %u )\n",runtime.Data(),tzone.Data(),mytime); |
418 |
}; |
419 |
// |
420 |
}; |
421 |
// |
422 |
// Which tle must be used for the date "tletime"? |
423 |
// |
424 |
if ( dtle ){ |
425 |
// |
426 |
// |
427 |
TDatime ti = TDatime(tletime.Data()); |
428 |
// |
429 |
TTimeStamp *time = new TTimeStamp((UInt_t)ti.GetYear(),(UInt_t)ti.GetMonth(),(UInt_t)ti.GetDay(),(UInt_t)ti.GetHour(),(UInt_t)ti.GetMinute(),(UInt_t)ti.GetSecond(),0,true,0); |
430 |
// |
431 |
UInt_t dbti = time->GetSec(); |
432 |
// |
433 |
TString thetime = dbtime->UnConvertTime(tzone,dbti); |
434 |
// |
435 |
ti = TDatime(thetime.Data()); |
436 |
TTimeStamp *time2 = new TTimeStamp((UInt_t)ti.GetYear(),(UInt_t)ti.GetMonth(),(UInt_t)ti.GetDay(),(UInt_t)ti.GetHour(),(UInt_t)ti.GetMinute(),(UInt_t)ti.GetSecond(),0,true,0); |
437 |
// |
438 |
UInt_t mytime = time2->GetSec(); |
439 |
// |
440 |
myquery.str(""); |
441 |
myquery << " select ID,TLE1,TLE2,TLE3 from GL_TLE where FROM_TIME<='" << time2->AsString("s") << "' ORDER BY FROM_TIME DESC LIMIT 1;"; |
442 |
// myquery << " from GL_TLE where FROM_TIME>=" << mytime << " ORDER BY FROM_TIME ASC LIMIT 1;"; |
443 |
// printf("myquery is %s \n",myquery.str().c_str()); |
444 |
pResult = dbc->Query(myquery.str().c_str()); |
445 |
Row = pResult->Next(); |
446 |
if ( !Row ){ |
447 |
printf("\n No TLE in the DB for date %s %s (DB time %u )\n",tletime.Data(),tzone.Data(),mytime); |
448 |
}; |
449 |
printf("\n Date %s %s (DB time %u ) is contained in TLE %u \n",tletime.Data(),tzone.Data(),mytime,(UInt_t)atoll(Row->GetField(0))); |
450 |
printf("\n%s\n",Row->GetField(1)); |
451 |
printf("%s\n",Row->GetField(2)); |
452 |
printf("%s\n",Row->GetField(3)); |
453 |
// |
454 |
FILE *tlefile; |
455 |
tlefile = fopen("tle.txt","w"); |
456 |
fprintf(tlefile,"%s\n",Row->GetField(1)); |
457 |
fprintf(tlefile,"%s\n",Row->GetField(2)); |
458 |
fprintf(tlefile,"%s\n",Row->GetField(3)); |
459 |
fclose (tlefile); |
460 |
// |
461 |
printf("\n TLE has been dumped in file tle.txt \n"); |
462 |
}; |
463 |
// |
464 |
// Which run contains the dbtime "runtime2"? |
465 |
// |
466 |
if ( ruti2 ){ |
467 |
// |
468 |
UInt_t mytime = runtime2; |
469 |
// |
470 |
Bool_t found = false; |
471 |
myquery.str(""); |
472 |
myquery << "select "; |
473 |
myquery << " ID "; |
474 |
myquery << " from GL_RUN where RUNHEADER_TIME<=" << mytime << " AND " |
475 |
<< " RUNTRAILER_TIME>=" << mytime << " ;"; |
476 |
// printf("myquery is %s \n",myquery.str().c_str()); |
477 |
pResult = dbc->Query(myquery.str().c_str()); |
478 |
for( r=0; r < 1000; r++){ |
479 |
Row = pResult->Next(); |
480 |
if( Row == NULL ) break; |
481 |
printf("\n DB time %u is contained in run %u \n",mytime,(UInt_t)atoll(Row->GetField(0))); |
482 |
found = true; |
483 |
}; |
484 |
// |
485 |
myquery.str(""); |
486 |
myquery << "select "; |
487 |
myquery << " ID "; |
488 |
myquery << " from GL_RUN_TRASH where BELONGED_TO='GL_RUN' AND RUNHEADER_TIME<=" << mytime << " AND " |
489 |
<< " RUNTRAILER_TIME>=" << mytime << " ;"; |
490 |
// printf("myquery is %s \n",myquery.str().c_str()); |
491 |
pResult = dbc->Query(myquery.str().c_str()); |
492 |
for( r=0; r < 1000; r++){ |
493 |
Row = pResult->Next(); |
494 |
if( Row == NULL ) break; |
495 |
printf("\n DB time %u is contained in run %u in the GL_RUN_TRASH table \n",mytime,(UInt_t)atoll(Row->GetField(0))); |
496 |
found = true; |
497 |
}; |
498 |
// |
499 |
if ( !found ){ |
500 |
printf("\n No run contains DB time %u \n",mytime); |
501 |
}; |
502 |
// |
503 |
}; |
504 |
// |
505 |
// To which file the run "run" belongs? |
506 |
// |
507 |
if ( run != 0 ){ |
508 |
Bool_t found = false; |
509 |
glrun->Clear(); |
510 |
error = glrun->Query_GL_RUN(run,dbc); |
511 |
glroot->Clear(); |
512 |
error = glroot->Query_GL_ROOT(glrun->ID_ROOT_L0,dbc); |
513 |
if ( glrun->ID_ROOT_L0 ){ |
514 |
if ( error ){ |
515 |
printf(" Error querying the DB! \n"); |
516 |
exit(-4); |
517 |
}; |
518 |
printf("\n Run %u belongs to file %s \n",run,(glroot->PATH+glroot->NAME).Data()); |
519 |
//filename = glroot->NAME; |
520 |
found = true; |
521 |
}; |
522 |
// |
523 |
myquery.str(""); |
524 |
myquery << "select "; |
525 |
myquery << " ID,FILENAMEL0,FILENAMEL2 "; |
526 |
myquery << " from GL_RUN_TRASH where BELONGED_TO='GL_RUN' AND ID=" << run << "; "; |
527 |
pResult = dbc->Query(myquery.str().c_str()); |
528 |
for( r=0; r < 1000; r++){ |
529 |
Row = pResult->Next(); |
530 |
if( Row == NULL ) break; |
531 |
printf("\n RUN %u has been deleted and now is contained in the GL_RUN_TRASH table \n",run); |
532 |
printf("\n RUN %u belonged to L0 file %s and L2 file %s \n",run,Row->GetField(1),Row->GetField(2)); |
533 |
found = true; |
534 |
}; |
535 |
// |
536 |
if ( !found ) printf("\n No run with ID=%u in the DB!\n",run); |
537 |
}; |
538 |
// |
539 |
// Which runs are contained in the file "filename"? |
540 |
// |
541 |
if ( strcmp(filename.Data(),"") ){ |
542 |
// ---------------- |
543 |
Bool_t found = false; |
544 |
Int_t ID = 0; |
545 |
Int_t ID_RAW = 0; |
546 |
// |
547 |
const char *rawpath = ""; |
548 |
const char *rawname = ""; |
549 |
// |
550 |
myquery.str(""); |
551 |
myquery << "select "; |
552 |
myquery << " ID"; |
553 |
myquery << ",ID_RAW"; |
554 |
myquery << ",PATH"; |
555 |
myquery << ",NAME"; |
556 |
myquery << " from GL_ROOT where NAME=\"" << filename.Data() << "\";"; |
557 |
pResult = dbc->Query(myquery.str().c_str()); |
558 |
for( r=0; r < 1000; r++){ |
559 |
Row = pResult->Next(); |
560 |
if( Row == NULL ) break; |
561 |
for( t = 0; t < pResult->GetFieldCount(); t++){ |
562 |
if(t==0) ID = atoi(Row->GetField(t)); |
563 |
if(t==1) ID_RAW = atoi(Row->GetField(t)); |
564 |
}; |
565 |
}; |
566 |
delete pResult; |
567 |
if ( !ID && !ID_RAW ){ |
568 |
if ( !neat ) printf("\n No file with name %s in the DB!\n",filename.Data()); |
569 |
} else { |
570 |
myquery.str(""); |
571 |
myquery << "select "; |
572 |
myquery << " ID,RUNHEADER_TIME,RUNTRAILER_TIME,NEVENTS "; |
573 |
myquery << " from GL_RUN where ID_ROOT_L0=" << ID << ";"; |
574 |
pResult = dbc->Query(myquery.str().c_str()); |
575 |
for( r=0; r < 1000; r++){ |
576 |
Row = pResult->Next(); |
577 |
if( Row == NULL ) break; |
578 |
found = true; |
579 |
if ( !r && !neat ) printf("\n File %s contains the following runs: \n\n",filename.Data()); |
580 |
TString UTC=tzone.Data(); |
581 |
if ( !neat ) printf(" => ID = %u --> the run started at %s %s ended at %s %s NEV = %u \n\n",(UInt_t)atoll(Row->GetField(0)),dbtime->ConvertTime(UTC,(UInt_t)atoll(Row->GetField(1))).Data(),tzone.Data(),dbtime->ConvertTime(UTC,(UInt_t)atoll(Row->GetField(2))).Data(),tzone.Data(),(UInt_t)atoll(Row->GetField(3))); |
582 |
if ( neat ) printf("%u\n",(UInt_t)atoll(Row->GetField(0))); |
583 |
}; |
584 |
delete pResult; |
585 |
myquery.str(""); |
586 |
myquery << "select "; |
587 |
myquery << " ID,RUNHEADER_TIME,RUNTRAILER_TIME "; |
588 |
myquery << " from GL_RUN_TRASH where BELONGED_TO='GL_RUN' AND FILENAMEL0='" << filename.Data() << "' ;"; |
589 |
pResult = dbc->Query(myquery.str().c_str()); |
590 |
for( r=0; r < 1000; r++){ |
591 |
Row = pResult->Next(); |
592 |
if( Row == NULL ) break; |
593 |
if ( !r && !neat ) printf("\n File %s contains the following DELETED runs: \n\n",filename.Data()); |
594 |
TString UTC=tzone.Data(); |
595 |
if ( !neat ) printf(" => ID = %i --> the run started at %s %s ended at %s %s \n\n",(UInt_t)atoll(Row->GetField(0)),dbtime->ConvertTime(UTC,(UInt_t)atoll(Row->GetField(1))).Data(),tzone.Data(),dbtime->ConvertTime(UTC,(UInt_t)atoll(Row->GetField(2))).Data(),tzone.Data()); |
596 |
}; |
597 |
// |
598 |
if ( !found ){ |
599 |
if ( !neat ) printf("\n No run associated to the file %s \n",filename.Data()); |
600 |
}; |
601 |
myquery.str(""); |
602 |
myquery << "select "; |
603 |
myquery << " PATH,NAME"; |
604 |
myquery << " from GL_RAW where ID=" << ID_RAW << ";"; |
605 |
pResult = dbc->Query(myquery.str().c_str()); |
606 |
for( r=0; r < 1000; r++){ |
607 |
Row = pResult->Next(); |
608 |
if( Row == NULL ) break; |
609 |
for( t = 0; t < pResult->GetFieldCount(); t++){ |
610 |
if(t==0) rawpath = Row->GetField(t); |
611 |
if(t==1) rawname = Row->GetField(t); |
612 |
}; |
613 |
}; |
614 |
delete pResult; |
615 |
if ( !neat ) printf("\n File %s belongs to raw data file %s/%s \n",filename.Data(),rawpath,rawname); |
616 |
}; |
617 |
}; |
618 |
// |
619 |
// Which runs are contained in the file "l2filename"? |
620 |
// |
621 |
if ( strcmp(l2filename.Data(),"") ){ |
622 |
// ---------------- |
623 |
// Bool_t found = false; |
624 |
Int_t ID = 0; |
625 |
// Int_t ID_RAW = 0; |
626 |
// |
627 |
const char *rawpath = ""; |
628 |
const char *rawname = ""; |
629 |
// |
630 |
myquery.str(""); |
631 |
myquery << "select "; |
632 |
myquery << " ID"; |
633 |
myquery << ",ID_RAW"; |
634 |
myquery << ",PATH"; |
635 |
myquery << ",NAME"; |
636 |
myquery << " from GL_ROOT where NAME=\"" << l2filename.Data() << "\" order by INSERT_TIME desc limit 1;"; |
637 |
pResult = dbc->Query(myquery.str().c_str()); |
638 |
for( r=0; r < 1000; r++){ |
639 |
Row = pResult->Next(); |
640 |
if( Row == NULL ) break; |
641 |
for( t = 0; t < pResult->GetFieldCount(); t++){ |
642 |
if(t==0) ID = atoi(Row->GetField(t)); |
643 |
// if(t==1) ID_RAW = atoi(Row->GetField(t)); |
644 |
}; |
645 |
}; |
646 |
delete pResult; |
647 |
if ( !ID ){ |
648 |
printf("\n No file with name %s in the DB!\n",l2filename.Data()); |
649 |
} else { |
650 |
nr = 0; |
651 |
rlist[0] = 0; |
652 |
myquery.str(""); |
653 |
myquery << "select "; |
654 |
myquery << " ID,ID_ROOT_L0 "; |
655 |
myquery << " from GL_RUN where ID_ROOT_L2=" << ID << ";"; |
656 |
pResult = dbc->Query(myquery.str().c_str()); |
657 |
for( r=0; r < pResult->GetRowCount(); r++){ |
658 |
Row = pResult->Next(); |
659 |
if( Row == NULL ) break; |
660 |
// found = true; |
661 |
if ( !r ) printf("\n File %s contains the following runs: \n\n",l2filename.Data()); |
662 |
printf(" %u ",(UInt_t)atoll(Row->GetField(0))); |
663 |
if ( r < (pResult->GetRowCount()-1) ){ |
664 |
printf("-"); |
665 |
} else { |
666 |
printf("\n\n"); |
667 |
}; |
668 |
if ( (UInt_t)atoll(Row->GetField(1)) != rlist[nr] ){ |
669 |
nr++; |
670 |
rlist[nr] = (UInt_t)atoll(Row->GetField(1)); |
671 |
}; |
672 |
}; |
673 |
delete pResult; |
674 |
// |
675 |
printf("\n File %s contains runs from L0 files: \n\n",l2filename.Data()); |
676 |
for (UInt_t l=1; l<nr+1; l++){ |
677 |
myquery.str(""); |
678 |
myquery << "select "; |
679 |
myquery << " PATH,NAME"; |
680 |
myquery << " from GL_ROOT where ID=" << rlist[l] << ";"; |
681 |
pResult = dbc->Query(myquery.str().c_str()); |
682 |
for( r=0; r < pResult->GetRowCount(); r++){ |
683 |
Row = pResult->Next(); |
684 |
if( Row == NULL ) break; |
685 |
for( t = 0; t < pResult->GetFieldCount(); t++){ |
686 |
if(t==0) rawpath = Row->GetField(t); |
687 |
if(t==1) rawname = Row->GetField(t); |
688 |
}; |
689 |
printf(" %s/%s \n",rawpath,rawname); |
690 |
}; |
691 |
delete pResult; |
692 |
}; |
693 |
}; |
694 |
}; |
695 |
// |
696 |
// |
697 |
// Show relationship between files around file "splitat" |
698 |
// |
699 |
if ( strcmp(splitat.Data(),"") ){ |
700 |
// ---------------- |
701 |
Int_t ID = 0; |
702 |
Int_t IDR = 0; |
703 |
TString PATH; |
704 |
TString NAME; |
705 |
UInt_t atime1 = 0; |
706 |
UInt_t atime2 = 0; |
707 |
UInt_t minid = 0; |
708 |
UInt_t maxid = 0; |
709 |
// |
710 |
myquery.str(""); |
711 |
myquery << " select ID,PATH,NAME from GL_ROOT where NAME>=\"" << splitat.Data() << "\" order by NAME asc limit " << nsplit << ";"; |
712 |
// printf(" myquery is %s \n",myquery.str().c_str()); |
713 |
pResult = dbc->Query(myquery.str().c_str()); |
714 |
GL_TIMESYNC *dbtime; |
715 |
GL_S4_CALIB *glS4calib; |
716 |
GL_ROOT *glroot; |
717 |
for( r=0; r < 1000; r++){ |
718 |
Row = pResult->Next(); |
719 |
if( Row == NULL ) break; |
720 |
printf("\n\n#################################################\n"); |
721 |
Int_t s=0; |
722 |
Int_t t=0; |
723 |
Int_t c=0; |
724 |
TString *s4files[500]; |
725 |
TString *tfiles[500]; |
726 |
TString *cfiles[2000]; |
727 |
ID = atoi(Row->GetField(0)); |
728 |
PATH = Row->GetField(1); |
729 |
NAME = Row->GetField(2); |
730 |
printf("\n FILE: %s/%s \n",PATH.Data(),NAME.Data()); |
731 |
myquery.str(""); |
732 |
myquery << " select ID,RUNHEADER_OBT,RUNTRAILER_OBT from GL_RUN where ID_ROOT_L0=" << ID << " order by ID;"; |
733 |
pResult2 = dbc->Query(myquery.str().c_str()); |
734 |
for( Int_t run=0; run < pResult2->GetRowCount(); run++){ |
735 |
Row2 = pResult2->Next(); |
736 |
if( Row2 == NULL ) break; |
737 |
IDR = atoi(Row2->GetField(0)); |
738 |
if ( run == 0 ) minid = IDR; |
739 |
if ( run == pResult2->GetRowCount()-1 ) maxid = IDR; |
740 |
// |
741 |
// here we make queries to DB looking for needed files and print results |
742 |
// |
743 |
dbtime = new GL_TIMESYNC(ID,"ID",dbc); |
744 |
// |
745 |
atime1 = dbtime->DBabsTime((UInt_t)atoll(Row2->GetField(1))); |
746 |
atime2 = dbtime->DBabsTime((UInt_t)atoll(Row2->GetField(2))); |
747 |
// printf(" atime1 %u atime2 %u field 1 %s field2 %s \n",atime1,atime2,Row2->GetField(1),Row2->GetField(2)); |
748 |
// |
749 |
// S4 |
750 |
// |
751 |
glS4calib = new GL_S4_CALIB(); |
752 |
glS4calib->Query_GL_S4_CALIB(atime1, dbc); |
753 |
glroot = new GL_ROOT(); |
754 |
glroot->Query_GL_ROOT(glS4calib->ID_ROOT_L0,dbc); |
755 |
// |
756 |
if ( s > 0 ){ |
757 |
Bool_t found = false; |
758 |
for (Int_t g=0; g<s; g++){ |
759 |
if ( !strcmp(glroot->NAME.Data(),s4files[g]->Data()) ){ |
760 |
found = true; |
761 |
}; |
762 |
}; |
763 |
if ( !found && strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
764 |
s4files[s] = new TString(glroot->NAME.Data()); |
765 |
s++; |
766 |
}; |
767 |
} else { |
768 |
if ( strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
769 |
s4files[s] = new TString(glroot->NAME.Data()); |
770 |
s++; |
771 |
}; |
772 |
}; |
773 |
// |
774 |
delete glS4calib; |
775 |
delete dbtime; |
776 |
delete glroot; |
777 |
// |
778 |
// TRACKER |
779 |
// |
780 |
GL_TRK_CALIB q2; |
781 |
q2.Query_GL_TRK_CALIB(atime1,dbc); |
782 |
GL_ROOT q3; |
783 |
q3.Query_GL_ROOT(q2.ID_ROOT_L0,dbc); |
784 |
// |
785 |
if ( t > 0 ){ |
786 |
Bool_t found = false; |
787 |
for (Int_t gt=0; gt<t; gt++){ |
788 |
if ( !strcmp(q3.NAME.Data(),tfiles[gt]->Data()) ){ |
789 |
found = true; |
790 |
}; |
791 |
}; |
792 |
if ( !found && strcmp(q3.NAME.Data(),NAME.Data()) ){ |
793 |
tfiles[t] = new TString(q3.NAME.Data()); |
794 |
t++; |
795 |
}; |
796 |
} else { |
797 |
if ( strcmp(q3.NAME.Data(),NAME.Data()) ){ |
798 |
tfiles[t] = new TString(q3.NAME.Data()); |
799 |
t++; |
800 |
}; |
801 |
}; |
802 |
// |
803 |
// CALO |
804 |
// |
805 |
GL_CALO_CALIB *glcalo = new GL_CALO_CALIB(); |
806 |
GL_CALOPULSE_CALIB *glp = new GL_CALOPULSE_CALIB(); |
807 |
GL_ROOT *glroot = new GL_ROOT(); |
808 |
// |
809 |
for (Int_t sc=0; sc<4; sc++){ |
810 |
// |
811 |
UInt_t pampli = 0; |
812 |
glcalo->Query_GL_CALO_CALIB(atime1,pampli,sc,dbc); |
813 |
glroot->Query_GL_ROOT(glcalo->ID_ROOT_L0,dbc); |
814 |
if ( c > 0 ){ |
815 |
Bool_t found = false; |
816 |
for (Int_t gt=0; gt<c; gt++){ |
817 |
if ( !strcmp(glroot->NAME.Data(),cfiles[gt]->Data()) ){ |
818 |
found = true; |
819 |
}; |
820 |
}; |
821 |
if ( !found && strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
822 |
cfiles[c] = new TString(glroot->NAME.Data()); |
823 |
c++; |
824 |
}; |
825 |
} else { |
826 |
if ( strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
827 |
cfiles[c] = new TString(glroot->NAME.Data()); |
828 |
c++; |
829 |
}; |
830 |
}; |
831 |
// |
832 |
pampli = 0; |
833 |
glcalo->Query_GL_CALO_CALIB(atime2,pampli,sc,dbc); |
834 |
glroot->Query_GL_ROOT(glcalo->ID_ROOT_L0,dbc); |
835 |
if ( c > 0 ){ |
836 |
Bool_t found = false; |
837 |
for (Int_t gt=0; gt<c; gt++){ |
838 |
if ( !strcmp(glroot->NAME.Data(),cfiles[gt]->Data()) ){ |
839 |
found = true; |
840 |
}; |
841 |
}; |
842 |
if ( !found && strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
843 |
cfiles[c] = new TString(glroot->NAME.Data()); |
844 |
c++; |
845 |
}; |
846 |
} else { |
847 |
if ( strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
848 |
cfiles[c] = new TString(glroot->NAME.Data()); |
849 |
c++; |
850 |
}; |
851 |
}; |
852 |
// |
853 |
pampli = 2; |
854 |
glp->Query_GL_CALOPULSE_CALIB(atime1,sc,pampli,dbc); |
855 |
glroot->Query_GL_ROOT(glp->ID_ROOT_L0,dbc); |
856 |
if ( c > 0 ){ |
857 |
Bool_t found = false; |
858 |
for (Int_t gt=0; gt<c; gt++){ |
859 |
if ( !strcmp(glroot->NAME.Data(),cfiles[gt]->Data()) ){ |
860 |
found = true; |
861 |
}; |
862 |
}; |
863 |
if ( !found && strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
864 |
cfiles[c] = new TString(glroot->NAME.Data()); |
865 |
c++; |
866 |
}; |
867 |
} else { |
868 |
if ( strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
869 |
cfiles[c] = new TString(glroot->NAME.Data()); |
870 |
c++; |
871 |
}; |
872 |
}; |
873 |
// |
874 |
pampli = 0; |
875 |
glp->Query_GL_CALOPULSE_CALIB(atime1,sc,pampli,dbc); |
876 |
glroot->Query_GL_ROOT(glp->ID_ROOT_L0,dbc); |
877 |
if ( c > 0 ){ |
878 |
Bool_t found = false; |
879 |
for (Int_t gt=0; gt<c; gt++){ |
880 |
if ( !strcmp(glroot->NAME.Data(),cfiles[gt]->Data()) ){ |
881 |
found = true; |
882 |
}; |
883 |
}; |
884 |
if ( !found && strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
885 |
cfiles[c] = new TString(glroot->NAME.Data()); |
886 |
c++; |
887 |
}; |
888 |
} else { |
889 |
if ( strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
890 |
cfiles[c] = new TString(glroot->NAME.Data()); |
891 |
c++; |
892 |
}; |
893 |
}; |
894 |
// |
895 |
pampli = 2; |
896 |
glp->Query_GL_CALOPULSE_CALIB(atime2,sc,pampli,dbc); |
897 |
glroot->Query_GL_ROOT(glp->ID_ROOT_L0,dbc); |
898 |
if ( c > 0 ){ |
899 |
Bool_t found = false; |
900 |
for (Int_t gt=0; gt<c; gt++){ |
901 |
if ( !strcmp(glroot->NAME.Data(),cfiles[gt]->Data()) ){ |
902 |
found = true; |
903 |
}; |
904 |
}; |
905 |
if ( !found && strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
906 |
cfiles[c] = new TString(glroot->NAME.Data()); |
907 |
c++; |
908 |
}; |
909 |
} else { |
910 |
if ( strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
911 |
cfiles[c] = new TString(glroot->NAME.Data()); |
912 |
c++; |
913 |
}; |
914 |
}; |
915 |
// |
916 |
pampli = 0; |
917 |
glp->Query_GL_CALOPULSE_CALIB(atime2,sc,pampli,dbc); |
918 |
glroot->Query_GL_ROOT(glp->ID_ROOT_L0,dbc); |
919 |
if ( c > 0 ){ |
920 |
Bool_t found = false; |
921 |
for (Int_t gt=0; gt<c; gt++){ |
922 |
if ( !strcmp(glroot->NAME.Data(),cfiles[gt]->Data()) ){ |
923 |
found = true; |
924 |
}; |
925 |
}; |
926 |
if ( !found && strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
927 |
cfiles[c] = new TString(glroot->NAME.Data()); |
928 |
c++; |
929 |
}; |
930 |
} else { |
931 |
if ( strcmp(glroot->NAME.Data(),NAME.Data()) ){ |
932 |
cfiles[c] = new TString(glroot->NAME.Data()); |
933 |
c++; |
934 |
}; |
935 |
}; |
936 |
// |
937 |
}; |
938 |
}; |
939 |
printf(" S4 requires:"); |
940 |
for (Int_t f=0;f<s; f++){ |
941 |
printf(" %s",s4files[f]->Data()); |
942 |
delete s4files[f]; |
943 |
}; |
944 |
printf("\n"); |
945 |
printf(" Tracker requires:"); |
946 |
for (Int_t f=0;f<t; f++){ |
947 |
printf(" %s",tfiles[f]->Data()); |
948 |
delete tfiles[f]; |
949 |
}; |
950 |
printf("\n"); |
951 |
printf(" Calorimeter requires:"); |
952 |
for (Int_t f=0;f<c; f++){ |
953 |
printf(" %s",cfiles[f]->Data()); |
954 |
delete cfiles[f]; |
955 |
}; |
956 |
printf("\n\n"); |
957 |
printf(" Minimum ID number in the file: %u \n",minid); |
958 |
printf(" Maximum ID number in the file: %u \n",maxid); |
959 |
printf("\n#################################################\n"); |
960 |
// |
961 |
}; |
962 |
delete pResult; |
963 |
if ( pResult2 ) delete pResult2; |
964 |
}; |
965 |
// |
966 |
// Close the DB connection |
967 |
// |
968 |
if ( dbc ){ |
969 |
dbc->Close(); |
970 |
delete dbc; |
971 |
dbc = 0; |
972 |
} |
973 |
// |
974 |
if ( !neat ) printf("\n"); |
975 |
// |
976 |
exit(0); |
977 |
} |