/[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.82 by malakhov, Tue Mar 3 10:58:13 2015 UTC revision 1.85 by mocchiut, Mon Mar 30 13:37:22 2015 UTC
# Line 198  int OrbitalInfoCore(UInt_t run, TFile *f Line 198  int OrbitalInfoCore(UInt_t run, TFile *f
198    Float_t dimo = 0.0; // dipole moment (computed from dat files) // EM GCC 4.7    Float_t dimo = 0.0; // dipole moment (computed from dat files) // EM GCC 4.7
199    Float_t bnorth, beast, bdown, babs;    Float_t bnorth, beast, bdown, babs;
200    Float_t xl; // L value    Float_t xl; // L value
201    Float_t icode; // code value for L accuracy (see fortran code)    Int_t icode; // code value for L accuracy (see fortran code)
202    Float_t bab1; // What's  the difference with babs?    Float_t bab1; // What's  the difference with babs?
203    Float_t stps = 0.005; // step size for field line tracing    Float_t stps = 0.005; // step size for field line tracing
204    Float_t bdel = 0.01; // required accuracy    Float_t bdel = 0.01; // required accuracy
# Line 1205  int OrbitalInfoCore(UInt_t run, TFile *f Line 1205  int OrbitalInfoCore(UInt_t run, TFile *f
1205        cCoordGeo coo;        cCoordGeo coo;
1206        Float_t jyear=0.;            Float_t jyear=0.;    
1207        //        //
1208        if(atime >= gltle->GetToTime()) {        if(atime >= gltle->GetToTime() || atime < gltle->GetFromTime() ) {  // AGH! bug when reprocessing??
1209          if ( !gltle->Query(atime, dbc) ){          if ( !gltle->Query(atime, dbc) ){
1210            //                  //      
1211            // Compute the magnetic dipole moment.            // Compute the magnetic dipole moment.
# Line 1300  int OrbitalInfoCore(UInt_t run, TFile *f Line 1300  int OrbitalInfoCore(UInt_t run, TFile *f
1300                              qRoll[sizeqmcmd]=RYPang_upper->Kren;                              qRoll[sizeqmcmd]=RYPang_upper->Kren;
1301                              qYaw[sizeqmcmd]=RYPang_upper->Ryskanie;                              qYaw[sizeqmcmd]=RYPang_upper->Ryskanie;
1302                              qPitch[sizeqmcmd]=RYPang_upper->Tangazh;                              qPitch[sizeqmcmd]=RYPang_upper->Tangazh;
 /* CHECK RECOVERED QUATERNIONS PROBLEM */  
 if(recqtime[mu]>=1160987921.75 && recqtime[mu]<=1160987932.00){  
   cout<<"We found it while checking all quaternions"<<"\t"<<recqtime[mu]<<endl;  
 }  
1303                            }                            }
1304                          }                          }
1305                          if(recqtime[mu]>=u_time){                          if(recqtime[mu]>=u_time){
# Line 1322  if(recqtime[mu]>=1160987921.75 && recqti Line 1318  if(recqtime[mu]>=1160987921.75 && recqti
1318                              qRoll[sizeqmcmd]=RYPang_upper->Kren;                              qRoll[sizeqmcmd]=RYPang_upper->Kren;
1319                              qYaw[sizeqmcmd]=RYPang_upper->Ryskanie;                              qYaw[sizeqmcmd]=RYPang_upper->Ryskanie;
1320                              qPitch[sizeqmcmd]=RYPang_upper->Tangazh;                              qPitch[sizeqmcmd]=RYPang_upper->Tangazh;
 /* CHECK RECOVERED QUATERNIONS PROBLEM */  
 if(recqtime[mu]>=1160987921.75 && recqtime[mu]<=1160987932.00){  
   cout<<"We found it while checking all quaternions"<<"\t"<<recqtime[mu]<<endl;  
 }  
1321                              break;                              break;
1322                            }                            }
1323                          }                          }
# Line 1394  if(recqtime[mu]>=1160987921.75 && recqti Line 1386  if(recqtime[mu]>=1160987921.75 && recqti
1386                             qYaw[sizeqmcmd]=RYPang_upper->Ryskanie;                             qYaw[sizeqmcmd]=RYPang_upper->Ryskanie;
1387                             qPitch[sizeqmcmd]=RYPang_upper->Tangazh;                             qPitch[sizeqmcmd]=RYPang_upper->Tangazh;
1388                             CopyQ(L_QQ_Q_l_lower,L_QQ_Q_l_upper);                             CopyQ(L_QQ_Q_l_lower,L_QQ_Q_l_upper);
 /* CHECK RECOVERED QUATERNIONS PROBLEM */  
 if(recqtime[mu]>=1160987921.75 && recqtime[mu]<=1160987932.00){  
   cout<<"We found it while checking all quaternions"<<"\t"<<recqtime[mu]<<endl;  
 }  
1389                             break;                             break;
1390                           }                           }
1391                        }                        }
# Line 1656  if(recqtime[mu]>=1160987921.75 && recqti Line 1644  if(recqtime[mu]>=1160987921.75 && recqti
1644  //           anni '50.  //           anni '50.
1645  //          //        
1646          //14.9/(xl*xl);          //14.9/(xl*xl);
1647          orbitalinfo->igrf_icode = icode;          orbitalinfo->igrf_icode = (Float_t)icode;
1648          //          //
1649        }              }      
1650        //        //
Legend:
Removed from v.1.82  
changed lines
  Added in v.1.85
  ViewVC Help
Powered by ViewVC 1.1.23