| 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 |  |