57 |
// CORE ROUTINE |
// CORE ROUTINE |
58 |
// |
// |
59 |
// |
// |
60 |
int NDCore(UInt_t run, TFile *file, TSQLServer *dbc, Int_t NDargc, char *NDargv[]){ |
int NDCore(UInt_t run, TFile *file, GL_TABLES *glt, Int_t NDargc, char *NDargv[]){ |
61 |
// |
// |
62 |
// Set these to true to have a verbose output. |
// Set these to true to have a verbose output. |
63 |
// |
// |
380 |
// |
// |
381 |
// prepare the timesync for the db |
// prepare the timesync for the db |
382 |
// |
// |
383 |
|
TString host = glt->CGetHost(); |
384 |
|
TString user = glt->CGetUser(); |
385 |
|
TString psw = glt->CGetPsw(); |
386 |
|
TSQLServer *dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data()); |
387 |
if ( !dbc->IsConnected() ) throw -604; |
if ( !dbc->IsConnected() ) throw -604; |
388 |
|
stringstream myquery; |
389 |
|
myquery.str(""); |
390 |
|
myquery << "SET time_zone='+0:00'"; |
391 |
|
dbc->Query(myquery.str().c_str()); |
392 |
dbtime = new GL_TIMESYNC(runinfo->ID_ROOT_L0,"ID",dbc); |
dbtime = new GL_TIMESYNC(runinfo->ID_ROOT_L0,"ID",dbc); |
393 |
// |
// |
394 |
// Search in the DB the path and name of the LEVEL0 file to be processed. |
// Search in the DB the path and name of the LEVEL0 file to be processed. |
395 |
// |
// |
396 |
if ( !dbc->IsConnected() ) throw -604; |
// if ( !dbc->IsConnected() ) throw -604; |
397 |
glroot->Query_GL_ROOT(runinfo->ID_ROOT_L0,dbc); |
glroot->Query_GL_ROOT(runinfo->ID_ROOT_L0,dbc); |
398 |
// |
// |
399 |
ftmpname.str(""); |
ftmpname.str(""); |
411 |
printf(" %i events to be processed for run %u: from %i to %i \n\n",totevent,idRun,runinfo->EV_FROM,runinfo->EV_FROM+totevent); |
printf(" %i events to be processed for run %u: from %i to %i \n\n",totevent,idRun,runinfo->EV_FROM,runinfo->EV_FROM+totevent); |
412 |
} |
} |
413 |
// |
// |
414 |
|
// if ( !totevent ) goto closeandexit; |
415 |
|
// |
416 |
// Open Level0 file |
// Open Level0 file |
417 |
// |
// |
418 |
l0File = new TFile(fname.Data()); |
l0File = new TFile(fname.Data()); |
449 |
// |
// |
450 |
nevents = l0ND->GetEntries(); |
nevents = l0ND->GetEntries(); |
451 |
// |
// |
452 |
if ( nevents < 1 ) { |
if ( nevents < 1 && totevent ) { |
453 |
if ( debug ) printf(" ND - ERROR: Level0 file is empty\n\n"); |
if ( debug ) printf(" ND - ERROR: Level0 file is empty\n\n"); |
454 |
l0File->Close(); |
l0File->Close(); |
455 |
code = -11; |
code = -11; |
456 |
goto closeandexit; |
goto closeandexit; |
457 |
}; |
}; |
458 |
// |
// |
459 |
if ( runinfo->EV_TO > nevents-1 ) { |
if ( runinfo->EV_TO > nevents-1 && totevent ) { |
460 |
if ( debug ) printf(" ND - ERROR: too few entries in the tree\n"); |
if ( debug ) printf(" ND - ERROR: too few entries in the tree\n"); |
461 |
l0File->Close(); |
l0File->Close(); |
462 |
code = -12; |
code = -12; |
467 |
// |
// |
468 |
if (verbose) printf("\n Ready to start! \n\n Processed events: \n\n"); |
if (verbose) printf("\n Ready to start! \n\n Processed events: \n\n"); |
469 |
// |
// |
470 |
|
if ( dbc ){ |
471 |
|
dbc->Close(); |
472 |
|
// delete dbc; |
473 |
|
}; |
474 |
|
// |
475 |
for ( re = runinfo->EV_FROM; re < (runinfo->EV_FROM+runinfo->NEVENTS); re++){ |
for ( re = runinfo->EV_FROM; re < (runinfo->EV_FROM+runinfo->NEVENTS); re++){ |
476 |
// |
// |
477 |
if ( procev%1000 == 0 && procev > 0 && verbose ) printf(" %iK \n",procev/1000); |
if ( procev%1000 == 0 && procev > 0 && verbose ) printf(" %iK \n",procev/1000); |
485 |
// |
// |
486 |
// paranoid check |
// paranoid check |
487 |
// |
// |
488 |
if ( (atime > runinfo->RUNTRAILER_TIME) || (atime < runinfo->RUNHEADER_TIME) ) { |
if ( (atime > (runinfo->RUNTRAILER_TIME+1)) || (atime < (runinfo->RUNHEADER_TIME-1)) ) { |
489 |
if (verbose) printf(" ND - WARNING: event at time outside the run time window, skipping it\n"); |
if (verbose) printf(" ND - WARNING: event at time outside the run time window, skipping it\n"); |
490 |
goto jumpev; |
goto jumpev; |
491 |
}; |
}; |