1 |
// |
2 |
// C/C++ headers |
3 |
// |
4 |
#include <fstream> |
5 |
#include <string.h> |
6 |
#include <iostream> |
7 |
#include <cstring> |
8 |
#include <stdio.h> |
9 |
// |
10 |
// ROOT headers |
11 |
// |
12 |
#include <TTree.h> |
13 |
#include <TClassEdit.h> |
14 |
#include <TObject.h> |
15 |
#include <TList.h> |
16 |
#include <TArrayL.h> |
17 |
#include <TSystem.h> |
18 |
#include <TSystemDirectory.h> |
19 |
#include <TString.h> |
20 |
#include <TFile.h> |
21 |
#include <TClass.h> |
22 |
#include <TSQLServer.h> |
23 |
#include <TSQLRow.h> |
24 |
#include <TSQLResult.h> |
25 |
// |
26 |
// YODA headers |
27 |
// |
28 |
#include <PamelaRun.h> |
29 |
#include <RegistryEvent.h> |
30 |
#include <PscuHeader.h> |
31 |
#include <PscuEvent.h> |
32 |
#include <EventHeader.h> |
33 |
#include <RegistryEvent.h> |
34 |
// |
35 |
// RunInfo header |
36 |
// |
37 |
#include <RunInfo.h> |
38 |
#include <GLTables.h> |
39 |
// |
40 |
// This program headers |
41 |
// |
42 |
#include <OrbitalInfo.h> |
43 |
#include <OrbitalInfoCore.h> |
44 |
#include <OrbitalInfoVerl2.h> |
45 |
|
46 |
using namespace std; |
47 |
|
48 |
// |
49 |
// CORE ROUTINE |
50 |
// |
51 |
// |
52 |
|
53 |
int OrbitalInfoCore(ULong64_t run, TFile *file, TSQLServer *dbc, Int_t OrbitalInfoargc, char *OrbitalInfoargv[]){ |
54 |
Int_t i = 0; |
55 |
// |
56 |
TString processFolder = "OrbitalInfoFolder"; |
57 |
// |
58 |
// Set these to true to have a very verbose output. |
59 |
// |
60 |
Bool_t debug = false; |
61 |
// |
62 |
Bool_t verbose = false; |
63 |
|
64 |
if ( OrbitalInfoargc > 0 ){ |
65 |
i = 0; |
66 |
while ( i < OrbitalInfoargc ){ |
67 |
if ( !strcmp(OrbitalInfoargv[i],"-processFolder") ) { |
68 |
if ( OrbitalInfoargc < i+1 ){ |
69 |
throw -3; |
70 |
}; |
71 |
processFolder = (TString)OrbitalInfoargv[i+1]; |
72 |
i++; |
73 |
}; |
74 |
if ( (!strcmp(OrbitalInfoargv[i],"--debug")) || (!strcmp(OrbitalInfoargv[i],"-g")) ) { |
75 |
verbose = true; |
76 |
}; |
77 |
if ( (!strcmp(OrbitalInfoargv[i],"--verbose")) || (!strcmp(OrbitalInfoargv[i],"-v")) ) { |
78 |
verbose = true; |
79 |
}; |
80 |
i++; |
81 |
}; |
82 |
}; |
83 |
// |
84 |
const char* outDir = gSystem->DirName(gSystem->DirName(file->GetPath())); |
85 |
// |
86 |
TTree *OrbitalInfotr = 0; |
87 |
Long64_t nevents = 0LL; |
88 |
// |
89 |
// variables needed to reprocess data |
90 |
// |
91 |
TString OrbitalInfoversion; |
92 |
ItoRunInfo *runinfo = 0; |
93 |
TArrayL *runlist = 0; |
94 |
TTree *OrbitalInfotrclone = 0; |
95 |
Bool_t reproc = false; |
96 |
Bool_t reprocall = false; |
97 |
UInt_t nobefrun = 0; |
98 |
UInt_t noaftrun = 0; |
99 |
UInt_t numbofrun = 0; |
100 |
stringstream ftmpname; |
101 |
TString fname; |
102 |
Long64_t totfileentries = 0ULL; |
103 |
Long64_t idRun = 0LL; |
104 |
// |
105 |
// variables needed to handle error signals |
106 |
// |
107 |
Int_t code = 0; |
108 |
Int_t sgnl; |
109 |
// |
110 |
// OrbitalInfo classes |
111 |
// |
112 |
OrbitalInfo *orbitalinfo = new OrbitalInfo(); |
113 |
OrbitalInfo *orbitalinfoclone = new OrbitalInfo(); |
114 |
// |
115 |
// define variables for opening and reading level0 file |
116 |
// |
117 |
TFile *l0File = 0; |
118 |
TTree *l0tr = 0; |
119 |
TBranch *l0registry = 0; |
120 |
pamela::RegistryEvent *l0reg=0; |
121 |
// EM: open also header branch |
122 |
TBranch *l0head = 0; |
123 |
pamela::EventHeader *eh = 0; |
124 |
pamela::PscuHeader *ph = 0; |
125 |
// end EM |
126 |
// |
127 |
// Define other basic variables |
128 |
// |
129 |
UInt_t procev = 0; |
130 |
stringstream file2; |
131 |
stringstream file3; |
132 |
stringstream qy; |
133 |
Int_t totevent = 0; |
134 |
ULong64_t atime = 0ULL; |
135 |
Int_t ei = 0; |
136 |
Int_t re = 0; |
137 |
// |
138 |
// Working filename |
139 |
// |
140 |
TString outputfile; |
141 |
stringstream name; |
142 |
name.str(""); |
143 |
name << outDir << "/"; |
144 |
// |
145 |
// temporary file and folder |
146 |
// |
147 |
TFile *tempfile = 0; |
148 |
TTree *tempOrbitalInfo = 0; |
149 |
stringstream tempname; |
150 |
stringstream OrbitalInfofolder; |
151 |
tempname.str(""); |
152 |
tempname << outDir; |
153 |
tempname << "/" << processFolder.Data(); |
154 |
OrbitalInfofolder.str(""); |
155 |
OrbitalInfofolder << tempname.str().c_str(); |
156 |
gSystem->MakeDirectory(OrbitalInfofolder.str().c_str()); |
157 |
tempname << "/OrbitalInfotree_run"; |
158 |
tempname << run << ".root"; |
159 |
// |
160 |
// DB classes |
161 |
// |
162 |
GL_ROOT *glroot = new GL_ROOT(); |
163 |
// |
164 |
// Let's start! |
165 |
// |
166 |
// As a first thing we must check what we have to do: if run = 0 we must process all events in the file has been passed |
167 |
// if run != 0 we must process only that run but first we have to check if the tree MyDetector2 already exist in the file |
168 |
// if it exists we are reprocessing data and we must delete that entries, if not we must create it. |
169 |
// |
170 |
if ( run == 0ULL ) reproc = true; |
171 |
// |
172 |
// |
173 |
// Output file is "outputfile" |
174 |
// |
175 |
if ( !file->IsOpen() ){ |
176 |
//printf(" OrbitalInfo - ERROR: cannot open file for writing\n"); |
177 |
throw -901; |
178 |
}; |
179 |
// |
180 |
// Retrieve GL_RUN variables from the level2 file |
181 |
// |
182 |
OrbitalInfoversion = OrbitalInfoInfo(false); // we should decide how to handle versioning system |
183 |
// |
184 |
// create an interface to RunInfo called "runinfo" |
185 |
// |
186 |
runinfo = new ItoRunInfo(file); |
187 |
// |
188 |
// open "Run" tree in level2 file, if not existing return an error (sngl != 0) |
189 |
// |
190 |
sgnl = 0; |
191 |
sgnl = runinfo->Update(run, "ORB", OrbitalInfoversion); |
192 |
//sgnl = runinfo->Read(run); |
193 |
|
194 |
if ( sgnl ){ |
195 |
//printf("OrbitalInfo - ERROR: RunInfo exited with non-zero status\n"); |
196 |
code = sgnl; |
197 |
goto closeandexit; |
198 |
} else { |
199 |
sgnl = 0; |
200 |
}; |
201 |
// |
202 |
// number of events in the file BEFORE the first event of our run |
203 |
// |
204 |
nobefrun = runinfo->GetFirstEntry(); |
205 |
// |
206 |
// total number of events in the file |
207 |
// |
208 |
totfileentries = runinfo->GetFileEntries(); |
209 |
// |
210 |
// first file entry AFTER the last event of our run |
211 |
// |
212 |
noaftrun = runinfo->GetLastEntry() + 1; |
213 |
// |
214 |
// number of run to be processed |
215 |
// |
216 |
numbofrun = runinfo->GetNoRun(); |
217 |
// |
218 |
// Try to access the OrbitalInfo tree in the file, if it exists we are reprocessing data if not we are processing a new run |
219 |
// |
220 |
OrbitalInfotrclone = (TTree*)file->Get("OrbitalInfo"); |
221 |
// |
222 |
if ( !OrbitalInfotrclone ){ |
223 |
// |
224 |
// tree does not exist, we are not reprocessing |
225 |
// |
226 |
reproc = false; |
227 |
if ( run == 0ULL ){ |
228 |
if (verbose) printf(" OrbitalInfo - WARNING: you are reprocessing data but OrbitalInfo tree does not exist!\n"); |
229 |
} |
230 |
if ( runinfo->IsReprocessing() && run != 0ULL ) { |
231 |
if (verbose) printf(" OrbitalInfo - WARNING: it seems you are not reprocessing data but OrbitalInfo\n versioning information already exists in RunInfo.\n"); |
232 |
} |
233 |
} else { |
234 |
// |
235 |
// tree exists, we are reprocessing data. Are we reprocessing a single run or all the file? |
236 |
// |
237 |
reproc = true; |
238 |
// |
239 |
// |
240 |
if (verbose) printf("\n Preparing the pre-processing...\n"); |
241 |
// |
242 |
if ( run == 0ULL ){ |
243 |
// |
244 |
// we are reprocessing all the file |
245 |
// if we are reprocessing everything we don't need to copy any old event and we can just work with the new tree and delete the old one immediately |
246 |
// |
247 |
reprocall = true; |
248 |
// |
249 |
if (verbose) printf("\n OrbitalInfo - WARNING: Reprocessing all runs\n"); |
250 |
// |
251 |
} else { |
252 |
// |
253 |
// we are reprocessing a single run, we must copy to the new tree the events in the file which preceed the first event of the run |
254 |
// |
255 |
reprocall = false; |
256 |
// |
257 |
if (verbose) printf("\n OrbitalInfo - WARNING: Reprocessing run number %llu \n",run); |
258 |
// |
259 |
// copying old tree to a new file |
260 |
// |
261 |
tempfile = new TFile(tempname.str().c_str(),"RECREATE"); |
262 |
tempOrbitalInfo = OrbitalInfotrclone->CloneTree(-1,"fast"); |
263 |
tempOrbitalInfo->SetName("OrbitalInfo-old"); |
264 |
tempfile->Write(); |
265 |
tempfile->Close(); |
266 |
} |
267 |
// |
268 |
// Delete the old tree from old file and memory |
269 |
// |
270 |
OrbitalInfotrclone->Delete("all"); |
271 |
// |
272 |
if (verbose) printf(" ...done!\n"); |
273 |
// |
274 |
}; |
275 |
// |
276 |
// create mydetector tree mydect |
277 |
// |
278 |
file->cd(); |
279 |
OrbitalInfotr = new TTree("OrbitalInfo-new","PAMELA OrbitalInfo data"); |
280 |
OrbitalInfotr->Branch("OrbitalInfo","OrbitalInfo",&orbitalinfo); |
281 |
// |
282 |
if ( reproc && !reprocall ){ |
283 |
// |
284 |
// open new file and retrieve also tree informations |
285 |
// |
286 |
tempfile = new TFile(tempname.str().c_str(),"READ"); |
287 |
OrbitalInfotrclone = (TTree*)tempfile->Get("OrbitalInfo-old"); |
288 |
OrbitalInfotrclone->SetBranchAddress("OrbitalInfo",&orbitalinfoclone); |
289 |
// |
290 |
if ( nobefrun > 0 ){ |
291 |
if (verbose){ |
292 |
printf("\n Pre-processing: copying events from the old tree before the processed run\n"); |
293 |
printf(" Copying %u events in the file which are before the beginning of the run %llu \n",nobefrun,run); |
294 |
printf(" Start copying at event number 0, end copying at event number %u \n",nobefrun); |
295 |
} |
296 |
for (UInt_t j = 0; j < nobefrun; j++){ |
297 |
// |
298 |
OrbitalInfotrclone->GetEntry(j); |
299 |
// |
300 |
// copy orbitalinfoclone to mydec |
301 |
// |
302 |
orbitalinfo = new OrbitalInfo(); |
303 |
memcpy(&orbitalinfo,&orbitalinfoclone,sizeof(orbitalinfoclone)); |
304 |
// |
305 |
// Fill entry in the new tree |
306 |
// |
307 |
OrbitalInfotr->Fill(); |
308 |
// |
309 |
}; |
310 |
if (verbose) printf(" Finished successful copying!\n"); |
311 |
}; |
312 |
}; |
313 |
// |
314 |
// Get the list of run to be processed, if only one run has to be processed the list will contain one entry only. |
315 |
// |
316 |
runlist = runinfo->GetRunList(); |
317 |
// |
318 |
// Loop over the run to be processed |
319 |
// |
320 |
for (UInt_t irun=0; irun < numbofrun; irun++){ |
321 |
// |
322 |
// retrieve the first run ID to be processed using the RunInfo list |
323 |
// |
324 |
idRun = runlist->At(irun); |
325 |
if (verbose){ |
326 |
printf("\n\n\n ####################################################################### \n"); |
327 |
printf(" PROCESSING RUN NUMBER %i \n",(int)idRun); |
328 |
printf(" ####################################################################### \n\n\n"); |
329 |
} |
330 |
// |
331 |
runinfo->ID_REG_RUN = 0ULL; |
332 |
// |
333 |
// store in the runinfo class the GL_RUN variables for our run |
334 |
// |
335 |
sgnl = 0; |
336 |
sgnl = runinfo->GetRunInfo(idRun); |
337 |
if ( sgnl ){ |
338 |
//printf("\n OrbitalInfo - ERROR: RunInfo exited with non-zero status\n"); |
339 |
code = sgnl; |
340 |
goto closeandexit; |
341 |
} else { |
342 |
sgnl = 0; |
343 |
}; |
344 |
// |
345 |
// now you can access that variables using the RunInfo class this way runinfo->ID_REG_RUN |
346 |
// |
347 |
if ( runinfo->ID_REG_RUN == 0 ){ |
348 |
//printf("\n OrbitalInfo - ERROR: no run with ID_RUN = %i \n\n Exiting... \n\n",(int)idRun); |
349 |
code = -5; |
350 |
goto closeandexit; |
351 |
}; |
352 |
// |
353 |
// Search in the DB the path and name of the LEVEL0 file to be processed. |
354 |
// |
355 |
glroot->Query_GL_ROOT(runinfo->ID_REG_RUN,dbc); |
356 |
// |
357 |
ftmpname.str(""); |
358 |
ftmpname << glroot->PATH.Data() << "/"; |
359 |
ftmpname << glroot->NAME.Data(); |
360 |
fname = ftmpname.str().c_str(); |
361 |
// |
362 |
// print out informations |
363 |
// |
364 |
totevent = runinfo->EV_REG_PHYS_TO - runinfo->EV_REG_PHYS_FROM + 1; |
365 |
if (verbose){ |
366 |
printf("\n LEVEL0 data file: %s \n",fname.Data()); |
367 |
printf(" RUN HEADER absolute time is: %llu \n",runinfo->RUNHEADER_TIME); |
368 |
printf(" RUN TRAILER absolute time is: %llu \n",runinfo->RUNTRAILER_TIME); |
369 |
printf(" %i events to be processed for run %llu: from %i to %i (reg entries)\n\n",totevent,idRun,runinfo->EV_REG_PHYS_FROM,runinfo->EV_REG_PHYS_TO); |
370 |
}// |
371 |
// Open Level0 file |
372 |
l0File = new TFile(fname.Data()); |
373 |
if ( !l0File ) { |
374 |
//printf(" OrbitalInfo - ERROR: problems opening Level0 file\n"); |
375 |
code = -6; |
376 |
goto closeandexit; |
377 |
}; |
378 |
l0tr = (TTree*)l0File->Get("Physics"); |
379 |
if ( !l0tr ) { |
380 |
//printf(" OrbitalInfo - ERROR: no Physics tree in Level0 file\n"); |
381 |
l0File->Close(); |
382 |
code = -7; |
383 |
goto closeandexit; |
384 |
}; |
385 |
// EM: open header branch as well |
386 |
l0head = l0tr->GetBranch("Header"); |
387 |
if ( !l0head ) { |
388 |
//if ( verbose ) printf(" CALORIMETER - ERROR: no Header branch in Level0 tree\n"); |
389 |
l0File->Close(); |
390 |
code = -8; |
391 |
goto closeandexit; |
392 |
}; |
393 |
l0tr->SetBranchAddress("Header", &eh); |
394 |
// end EM |
395 |
l0registry = l0tr->GetBranch("Registry"); |
396 |
if ( !l0registry ) { |
397 |
//printf(" OrbitalInfo - ERROR: no Registry branch in Level0 tree\n"); |
398 |
l0File->Close(); |
399 |
code = -9; |
400 |
goto closeandexit; |
401 |
}; |
402 |
// |
403 |
l0tr->SetBranchAddress("Registry", &l0reg); |
404 |
// |
405 |
nevents = l0registry->GetEntries(); |
406 |
// |
407 |
if ( nevents < 1 ) { |
408 |
//printf(" OrbitalInfo - ERROR: Level0 file is empty\n\n"); |
409 |
l0File->Close(); |
410 |
code = -11; |
411 |
goto closeandexit; |
412 |
}; |
413 |
// |
414 |
if ( runinfo->EV_REG_PHYS_TO > nevents-1 ) { |
415 |
//printf(" OrbitalInfo - ERROR: too few entries in the registry tree\n"); |
416 |
l0File->Close(); |
417 |
code = -12; |
418 |
goto closeandexit; |
419 |
}; |
420 |
// |
421 |
// run over all the events of the run |
422 |
// |
423 |
if (verbose) printf("\n Ready to start! \n\n Processed events: \n\n"); |
424 |
// |
425 |
for ( re = runinfo->EV_REG_PHYS_FROM; re <= runinfo->EV_REG_PHYS_TO; re++){ |
426 |
// |
427 |
if ( procev%1000 == 0 && procev > 0 && verbose) printf(" %iK \n",procev/1000); |
428 |
// |
429 |
l0registry->GetEntry(re); |
430 |
// |
431 |
// absolute time of this event |
432 |
// |
433 |
atime = l0reg->absTime; |
434 |
// |
435 |
// physics events is at entry number ei where |
436 |
// |
437 |
ei = l0reg->event; |
438 |
// |
439 |
// paranoid check |
440 |
// |
441 |
if ( (atime > runinfo->RUNTRAILER_TIME) || (atime < runinfo->RUNHEADER_TIME) ) { |
442 |
if (verbose) printf(" OrbitalInfo - WARNING: event at time outside the run time window, skipping it\n"); |
443 |
goto jumpev; |
444 |
}; |
445 |
// |
446 |
procev++; |
447 |
// |
448 |
// start processing |
449 |
// |
450 |
orbitalinfo = new OrbitalInfo(); |
451 |
orbitalinfo->absTime = l0reg->absTime; |
452 |
// EM: add OBT and plt_num infos from the header |
453 |
l0head->GetEntry(ei); |
454 |
ph = eh->GetPscuHeader(); |
455 |
orbitalinfo->pkt_num = ph->GetCounter(); |
456 |
orbitalinfo->OBT = ph->GetOrbitalTime(); |
457 |
// end EM |
458 |
OrbitalInfotr->Fill(); |
459 |
// |
460 |
// |
461 |
jumpev: |
462 |
debug = false; |
463 |
// |
464 |
}; |
465 |
// |
466 |
// Here you may want to clear some variables before processing another run |
467 |
// |
468 |
ei = 0; |
469 |
}; // process all the runs |
470 |
// |
471 |
if (verbose) printf("\n Finished processing data \n"); |
472 |
// |
473 |
closeandexit: |
474 |
// |
475 |
// we have finished processing the run(s). If we processed a single run now we must copy all the events after our run from the old tree to the new one and delete the old tree. |
476 |
// |
477 |
if ( !reprocall && reproc && code >= 0 ){ |
478 |
if ( totfileentries > noaftrun ){ |
479 |
if (verbose){ |
480 |
printf("\n Post-processing: copying events from the old tree after the processed run\n"); |
481 |
printf(" Copying %i events in the file which are after the end of the run %i \n",(int)(totfileentries-noaftrun),(int)run); |
482 |
printf(" Start copying at event number %i end copying at event number %i \n",(int)noaftrun,(int)totfileentries); |
483 |
} |
484 |
for (UInt_t j = noaftrun; j < totfileentries; j++ ){ |
485 |
// |
486 |
// Get entry from old tree |
487 |
// |
488 |
OrbitalInfotrclone->GetEntry(j); |
489 |
// |
490 |
// copy orbitalinfoclone to OrbitalInfo |
491 |
// |
492 |
orbitalinfo = new OrbitalInfo(); |
493 |
memcpy(&orbitalinfo,&orbitalinfoclone,sizeof(orbitalinfoclone)); |
494 |
// |
495 |
// Fill entry in the new tree |
496 |
// |
497 |
OrbitalInfotr->Fill(); |
498 |
}; |
499 |
if (verbose) printf(" Finished successful copying!\n"); |
500 |
}; |
501 |
}; |
502 |
// |
503 |
// Close files, delete old tree(s), write and close level2 file |
504 |
// |
505 |
if ( l0File ) l0File->Close(); |
506 |
if ( tempfile ) tempfile->Close(); |
507 |
gSystem->Unlink(tempname.str().c_str()); |
508 |
|
509 |
//if ( code < 0 ) printf("\n OrbitalInfo - ERROR: an error occurred, try to save anyway...\n"); |
510 |
//printf("\n Writing and closing rootple\n"); |
511 |
if ( runinfo ) runinfo->Close(); |
512 |
if ( OrbitalInfotr ) OrbitalInfotr->SetName("OrbitalInfo"); |
513 |
if ( file ){ |
514 |
file->cd(); |
515 |
file->Write(); |
516 |
}; |
517 |
// |
518 |
gSystem->Unlink(OrbitalInfofolder.str().c_str()); |
519 |
// |
520 |
// the end |
521 |
// |
522 |
if (verbose) printf("\n Exiting...\n"); |
523 |
if(OrbitalInfotr)OrbitalInfotr->Delete(); |
524 |
if(code < 0) throw code; |
525 |
return(code); |
526 |
} |
527 |
|