/[PAMELA software]/DarthVader/OrbitalInfo/src/OrbitalInfoCore.cpp
ViewVC logotype

Contents of /DarthVader/OrbitalInfo/src/OrbitalInfoCore.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.2 - (show annotations) (download)
Wed Jun 14 09:52:55 2006 UTC (18 years, 5 months ago) by mocchiut
Branch: MAIN
CVS Tags: v0r02, v1r00
Changes since 1.1: +21 -0 lines
Added OBT and pkt_num variables

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

  ViewVC Help
Powered by ViewVC 1.1.23