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

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.56 by mocchiut, Thu Apr 12 12:27:02 2012 UTC revision 1.57 by mocchiut, Tue May 15 14:31:13 2012 UTC
# Line 224  int OrbitalInfoCore(UInt_t run, TFile *f Line 224  int OrbitalInfoCore(UInt_t run, TFile *f
224    //    //
225    //Quaternions classes    //Quaternions classes
226    //    //
227    Quaternions *L_QQ_Q_l_lower = new Quaternions();    Quaternions *L_QQ_Q_l_lower = 0;
228    InclinationInfo *RYPang_lower = new InclinationInfo();    InclinationInfo *RYPang_lower = 0;
229    Quaternions *L_QQ_Q_l_upper = new Quaternions();    Quaternions *L_QQ_Q_l_upper = 0;
230    InclinationInfo *RYPang_upper = new InclinationInfo();    InclinationInfo *RYPang_upper = 0;
231        
232    cEci eCi;    cEci eCi;
233        
# Line 398  int OrbitalInfoCore(UInt_t run, TFile *f Line 398  int OrbitalInfoCore(UInt_t run, TFile *f
398        //        //
399        reprocall = true;        reprocall = true;
400        //        //
401        if (verbose) printf("\n OrbitalInfo - WARNING: Reprocessing all runs\n");        if (verbose) printf("\n OrbitalInfo - WARNING: Reprocessing all runs\n Deleting old tree...\n");
402        //        //
403      } else {      } else {
404        //        //
# Line 416  int OrbitalInfoCore(UInt_t run, TFile *f Line 416  int OrbitalInfoCore(UInt_t run, TFile *f
416        tempOrbitalInfo = OrbitalInfotrclone->CloneTree(-1,"fast");        tempOrbitalInfo = OrbitalInfotrclone->CloneTree(-1,"fast");
417        tempOrbitalInfo->SetName("OrbitalInfo-old");        tempOrbitalInfo->SetName("OrbitalInfo-old");
418        tempfile->Write();        tempfile->Write();
419          tempOrbitalInfo->Delete();
420        tempfile->Close();          tempfile->Close();  
421      }      }
422      //      //
423      // Delete the old tree from old file and memory      // Delete the old tree from old file and memory
424      //      //
425        OrbitalInfotrclone->Clear();
426      OrbitalInfotrclone->Delete("all");      OrbitalInfotrclone->Delete("all");
427      //      //
428      if (verbose) printf(" ...done!\n");      if (verbose) printf(" ...done!\n");
# Line 466  int OrbitalInfoCore(UInt_t run, TFile *f Line 468  int OrbitalInfoCore(UInt_t run, TFile *f
468          //          //
469        };        };
470        if (verbose) printf(" Finished successful copying!\n");        if (verbose) printf(" Finished successful copying!\n");
471      };                };
472    };    };
473    //    //
474    //    //
# Line 477  int OrbitalInfoCore(UInt_t run, TFile *f Line 479  int OrbitalInfoCore(UInt_t run, TFile *f
479    // Loop over the run to be processed    // Loop over the run to be processed
480    //    //
481    for (UInt_t irun=0; irun < numbofrun; irun++){    for (UInt_t irun=0; irun < numbofrun; irun++){
482    
483        L_QQ_Q_l_lower = new Quaternions();
484        RYPang_lower = new InclinationInfo();
485        L_QQ_Q_l_upper = new Quaternions();
486        RYPang_upper = new InclinationInfo();
487    
488      //      //
489      // retrieve the first run ID to be processed using the RunInfo list      // retrieve the first run ID to be processed using the RunInfo list
490      //      //
# Line 878  int OrbitalInfoCore(UInt_t run, TFile *f Line 886  int OrbitalInfoCore(UInt_t run, TFile *f
886              + (hour*3600.+min*60.+(float)sec)/(24.*3600.*365.);              + (hour*3600.+min*60.+(float)sec)/(24.*3600.*365.);
887            //            //
888            if ( debug ) printf(" %i compute magnetic dipole moment get dipole moment for year\n",procev);                  if ( debug ) printf(" %i compute magnetic dipole moment get dipole moment for year\n",procev);      
889              if ( debug ) printf(" %i jyear %f dimo %f \n",procev,jyear,dimo);      
890            feldcof_(&jyear, &dimo); // get dipole moment for year            feldcof_(&jyear, &dimo); // get dipole moment for year
891            if ( debug ) printf(" %i compute magnetic dipole moment end\n",procev);                  if ( debug ) printf(" %i compute magnetic dipole moment end\n",procev);      
892          } else {          } else {
# Line 1221  int OrbitalInfoCore(UInt_t run, TFile *f Line 1230  int OrbitalInfoCore(UInt_t run, TFile *f
1230          orbitalinfo->BB0 = babs/bequ;          orbitalinfo->BB0 = babs/bequ;
1231          orbitalinfo->L = xl;                orbitalinfo->L = xl;      
1232          // Set Stormer vertical cutoff using L shell.          // Set Stormer vertical cutoff using L shell.
1233          orbitalinfo->cutoffsvl = 14.9/(xl*xl);          orbitalinfo->cutoffsvl = 14.295 / (xl*xl); //
1234            /*
1235    ---------- Forwarded message ----------
1236    Date: Wed, 09 May 2012 12:16:47 +0200
1237    From: Alessandro Bruno <alessandro.bruno@ba.infn.it>
1238    To: Mirko Boezio <mirko.boezio@ts.infn.it>
1239    Cc: Francesco S. Cafagna <Francesco.Cafagna@ba.infn.it>
1240    Subject: Störmer vertical cutoff
1241    
1242    Ciao Mirko,
1243    volevo segnalarti che il valore dello Störmer vertical cutoff nel Level2 č
1244    sovrastimato di circa il 4%.
1245    Dopo un'approfondita analisi con l'IGRF-05 abbiamo ricavano un valore pari
1246    a: 14.295 / L^2 anzichč 14.9 / L^2, valore obsoleto in quanto riferito agli
1247    anni '50.
1248    */
1249            //14.9/(xl*xl);
1250          orbitalinfo->igrf_icode = icode;          orbitalinfo->igrf_icode = icode;
1251          //          //
1252        };              };      
# Line 1366  int OrbitalInfoCore(UInt_t run, TFile *f Line 1391  int OrbitalInfoCore(UInt_t run, TFile *f
1391      c1->SaveAs("9.Rollhyst.png");      c1->SaveAs("9.Rollhyst.png");
1392      delete c1;*/      delete c1;*/
1393    
1394        if ( verbose ) printf(" Clear before new run \n");
1395      delete dbtime;      delete dbtime;
1396      if ( L_QQ_Q_l_upper ) delete L_QQ_Q_l_upper;  
1397        mcmdrc->Clear();
1398        mcmdrc = 0;
1399        
1400        if ( verbose ) printf(" Clear before new run1 \n");
1401      if ( L_QQ_Q_l_lower ) delete L_QQ_Q_l_lower;      if ( L_QQ_Q_l_lower ) delete L_QQ_Q_l_lower;
1402        if ( verbose ) printf(" Clear before new run2 \n");
1403        if ( L_QQ_Q_l_upper ) delete L_QQ_Q_l_upper;
1404        if ( verbose ) printf(" Clear before new run3 \n");
1405      if ( RYPang_upper ) delete RYPang_upper;      if ( RYPang_upper ) delete RYPang_upper;
1406        if ( verbose ) printf(" Clear before new run4 \n");
1407      if ( RYPang_lower ) delete RYPang_lower;      if ( RYPang_lower ) delete RYPang_lower;
1408    
1409        if ( l0tr ) l0tr->Delete();
1410        
1411        if ( verbose ) printf(" End run \n");
1412    
1413    }; // process all the runs    }; // process all the runs
1414        
1415    if (verbose) printf("\n Finished processing data \n");    if (verbose) printf("\n Finished processing data \n");
# Line 1404  int OrbitalInfoCore(UInt_t run, TFile *f Line 1443  int OrbitalInfoCore(UInt_t run, TFile *f
1443        };        };
1444        if (verbose) printf(" Finished successful copying!\n");        if (verbose) printf(" Finished successful copying!\n");
1445      };      };
1446        //if ( OrbitalInfotrclone )    OrbitalInfotrclone->Clear();        
1447        //if ( OrbitalInfotrclone )    OrbitalInfotrclone->Delete();        
1448    };    };
1449    //    //
1450    // Close files, delete old tree(s), write and close level2 file    // Close files, delete old tree(s), write and close level2 file
1451    //    //
1452    if ( l0File ) l0File->Close();    if ( l0File ) l0File->Close();
   if ( tempfile ) tempfile->Close();              
1453    if ( myfold ) gSystem->Unlink(tempname.str().c_str());    if ( myfold ) gSystem->Unlink(tempname.str().c_str());
1454    //    //
   if ( runinfo ) runinfo->Close();      
1455    if ( OrbitalInfotr ) OrbitalInfotr->SetName("OrbitalInfo");        if ( OrbitalInfotr ) OrbitalInfotr->SetName("OrbitalInfo");    
   if ( tof ) tof->Delete();  
   if ( ttof ) ttof->Delete();  
1456    //    //
1457    if ( file ){    if ( file ){
1458      file->cd();      file->cd();
1459      file->Write();      OrbitalInfotr->Write("OrbitalInfo", TObject::kOverwrite);
1460    };    };
1461    //    //
1462      if (verbose) printf("\n Exiting...\n");
1463    
1464    if ( myfold ) gSystem->Unlink(OrbitalInfofolder.str().c_str());    if ( myfold ) gSystem->Unlink(OrbitalInfofolder.str().c_str());
1465    //    //
1466    // the end    // the end
# Line 1430  int OrbitalInfoCore(UInt_t run, TFile *f Line 1469  int OrbitalInfoCore(UInt_t run, TFile *f
1469      dbc->Close();      dbc->Close();
1470      delete dbc;      delete dbc;
1471    };    };
   if (verbose) printf("\n Exiting...\n");  
   if(OrbitalInfotr)OrbitalInfotr->Delete();  
1472    //    //
1473      if (verbose) printf("\n Exiting...\n");
1474      if ( tempfile ) tempfile->Close();            
1475      
1476    if ( PO ) delete PO;    if ( PO ) delete PO;
1477    if ( orbitalinfo ) delete orbitalinfo;    if ( gltle ) delete gltle;
1478    if ( orbitalinfoclone ) delete orbitalinfoclone;    if ( glparam ) delete glparam;
1479      if ( glparam2 ) delete glparam2;
1480      if ( glparam3 ) delete glparam3;
1481      if (verbose) printf("\n Exiting3...\n");
1482    if ( glroot ) delete glroot;    if ( glroot ) delete glroot;
1483      if (verbose) printf("\n Exiting4...\n");
1484      if ( runinfo ) runinfo->Close();    
1485    if ( runinfo ) delete runinfo;    if ( runinfo ) delete runinfo;
1486      
1487      cout << "1   0x" << OrbitalInfotr << endl;
1488      cout << "2   0x" << OrbitalInfotrclone << endl;
1489      cout << "3   0x" << l0tr << endl;
1490      cout << "4   0x" << tempOrbitalInfo << endl;
1491      cout << "5   0x" << ttof << endl;
1492    
1493      //
1494      file->ls();
1495    //    //
1496    if(code < 0)  throw code;    if(code < 0)  throw code;
1497    return(code);    return(code);
Legend:
Removed from v.1.56  
changed lines
  Added in v.1.57
  ViewVC Help
Powered by ViewVC 1.1.23