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 |
// |
// |
165 |
TTree *tempND = 0; |
TTree *tempND = 0; |
166 |
stringstream tempname; |
stringstream tempname; |
167 |
stringstream NDfolder; |
stringstream NDfolder; |
168 |
|
Bool_t myfold = false; |
169 |
tempname.str(""); |
tempname.str(""); |
170 |
tempname << outDir; |
tempname << outDir; |
171 |
tempname << "/" << processFolder.Data(); |
tempname << "/" << processFolder.Data(); |
172 |
NDfolder.str(""); |
NDfolder.str(""); |
173 |
NDfolder << tempname.str().c_str(); |
NDfolder << tempname.str().c_str(); |
|
gSystem->MakeDirectory(NDfolder.str().c_str()); |
|
174 |
tempname << "/NDtree_run"; |
tempname << "/NDtree_run"; |
175 |
tempname << run << ".root"; |
tempname << run << ".root"; |
176 |
// |
// |
281 |
// |
// |
282 |
// copying old tree to a new file |
// copying old tree to a new file |
283 |
// |
// |
284 |
|
gSystem->MakeDirectory(NDfolder.str().c_str()); |
285 |
|
myfold = true; |
286 |
tempfile = new TFile(tempname.str().c_str(),"RECREATE"); |
tempfile = new TFile(tempname.str().c_str(),"RECREATE"); |
287 |
tempND = NDtrclone->CloneTree(-1,"fast"); |
tempND = NDtrclone->CloneTree(-1,"fast"); |
288 |
tempND->SetName("NeutronD-old"); |
tempND->SetName("NeutronD-old"); |
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()); |
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); |
565 |
// |
// |
566 |
if ( l0File ) l0File->Close(); |
if ( l0File ) l0File->Close(); |
567 |
if ( tempfile ) tempfile->Close(); |
if ( tempfile ) tempfile->Close(); |
568 |
gSystem->Unlink(tempname.str().c_str()); |
if ( myfold ) gSystem->Unlink(tempname.str().c_str()); |
569 |
// |
// |
570 |
if ( runinfo ) runinfo->Close(); |
if ( runinfo ) runinfo->Close(); |
571 |
if ( NDtr ) NDtr->SetName("NeutronD"); |
if ( NDtr ) NDtr->SetName("NeutronD"); |
574 |
file->Write(); |
file->Write(); |
575 |
}; |
}; |
576 |
// |
// |
577 |
gSystem->Unlink(NDfolder.str().c_str()); |
if ( myfold ) gSystem->Unlink(NDfolder.str().c_str()); |
578 |
// |
// |
579 |
// the end |
// the end |
580 |
// |
// |