/[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.51 by pam-mep, Tue Nov 15 09:58:30 2011 UTC revision 1.52 by mocchiut, Wed Nov 23 21:19:34 2011 UTC
# Line 231  int OrbitalInfoCore(UInt_t run, TFile *f Line 231  int OrbitalInfoCore(UInt_t run, TFile *f
231    cEci eCi;    cEci eCi;
232        
233    // Initialize fortran routines!!!    // Initialize fortran routines!!!
234      Int_t ltp1 = 0;
235    Int_t ltp2 = 0;    Int_t ltp2 = 0;
236    Int_t ltp3 = 0;    Int_t ltp3 = 0;
237    Int_t uno = 1;    Int_t uno = 1;
238    const char *niente = " ";    const char *niente = " ";
239    GL_PARAM *glparam = new GL_PARAM();    GL_PARAM *glparam = new GL_PARAM();
240    GL_PARAM *glparam2 = new GL_PARAM();    GL_PARAM *glparam2 = new GL_PARAM();
241      GL_PARAM *glparam3 = new GL_PARAM();
242    
243    //    //
244    // Orientation variables. Vitaly    // Orientation variables. Vitaly
# Line 287  int OrbitalInfoCore(UInt_t run, TFile *f Line 289  int OrbitalInfoCore(UInt_t run, TFile *f
289    }    }
290    if ( verbose ) cout<<"We have read recovered data"<<endl;    if ( verbose ) cout<<"We have read recovered data"<<endl;
291    
292      UInt_t anni5 = 60. * 60. * 24. * 365. * 5 ;//1576800
293    parerror=glparam->Query_GL_PARAM(1,301,dbc); // parameters stored in DB in GL_PRAM table      parerror=glparam->Query_GL_PARAM(1,301,(dbc-anni5)); // parameters stored in DB in GL_PRAM table  
294    if ( parerror<0 ) {    if ( parerror<0 ) {
295      code = parerror;      code = parerror;
296      goto closeandexit;      goto closeandexit;
297    };    };
298    ltp2 = (Int_t)(glparam->PATH+glparam->NAME).Length();    ltp1 = (Int_t)(glparam->PATH+glparam->NAME).Length();
299    if ( verbose ) printf(" Reading Earth's Magnetic Field parameter file: %s \n",(glparam->PATH+glparam->NAME).Data());    if ( verbose ) printf(" Reading Earth's Magnetic Field parameter file: %s \n",(glparam->PATH+glparam->NAME).Data());
300    //    //
301    parerror=glparam2->Query_GL_PARAM(1,302,dbc); // parameters stored in DB in GL_PRAM table    parerror=glparam2->Query_GL_PARAM(1,301,dbc); // parameters stored in DB in GL_PRAM table  
302    if ( parerror<0 ) {    if ( parerror<0 ) {
303      code = parerror;      code = parerror;
304      goto closeandexit;      goto closeandexit;
305    };    };
306    ltp3 = (Int_t)(glparam2->PATH+glparam2->NAME).Length();    ltp2 = (Int_t)(glparam2->PATH+glparam->NAME).Length();
307    if ( verbose ) printf(" Reading Earth's Magnetic Field parameter file: %s \n",(glparam2->PATH+glparam2->NAME).Data());    if ( verbose ) printf(" Reading Earth's Magnetic Field parameter file: %s \n",(glparam2->PATH+glparam2->NAME).Data());
308    //    //
309    initize_((char *)niente,&uno,(char *)(glparam->PATH+glparam->NAME).Data(),&ltp2,(char *)(glparam2->PATH+glparam2->NAME).Data(),&ltp3);    parerror=glparam3->Query_GL_PARAM(1,302,dbc); // parameters stored in DB in GL_PRAM table
310      if ( parerror<0 ) {
311        code = parerror;
312        goto closeandexit;
313      };
314      ltp3 = (Int_t)(glparam3->PATH+glparam2->NAME).Length();
315      if ( verbose ) printf(" Reading Earth's Magnetic Field parameter file: %s \n",(glparam3->PATH+glparam3->NAME).Data());
316      //
317      initize_((char *)(glparam->PATH+glparam->NAME).Data(),&ltp1,(char *)(glparam2->PATH+glparam2->NAME).Data(),&ltp2,(char *)(glparam3->PATH+glparam3->NAME).Data(),&ltp3);
318    //    //
319    // End IGRF stuff//    // End IGRF stuff//
320    //    //
Legend:
Removed from v.1.51  
changed lines
  Added in v.1.52
  ViewVC Help
Powered by ViewVC 1.1.23