/[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.85 by mocchiut, Mon Mar 30 13:37:22 2015 UTC revision 1.86 by pam-ts, Thu Apr 2 11:10:55 2015 UTC
# Line 1667  if(recqtime[mu]>=1160987921.75 && recqti Line 1667  if(recqtime[mu]>=1160987921.75 && recqti
1667          // If initial orientation data have reason to be inaccurate          // If initial orientation data have reason to be inaccurate
1668         Float_t tg = 0.00;         Float_t tg = 0.00;
1669         Float_t tmptg;         Float_t tmptg;
1670            Bool_t tgpar=false;
1671            Bool_t tgpar0=false;
1672            if (orbitalinfo->TimeGap>10.0 || ((modf(orbitalinfo->TimeGap,&tmptg)*1000>10 || modf(orbitalinfo->TimeGap,&tmptg)*1000==0.0) && orbitalinfo->TimeGap>2.0)) tgpar=true;
1673            if (orbitalinfo->TimeGap>180.0) tgpar0=true;
1674         if(MU!=0){         if(MU!=0){
1675  //      if(orbitalinfo->TimeGap>0 && errq==0 && azim==0){               // 10RED CHECK  (comparison between three metod of recovering orientation)  //      if(orbitalinfo->TimeGap>0 && errq==0 && azim==0){               // 10RED CHECK  (comparison between three metod of recovering orientation)
1676         if((atime>=RTstart[MU] && atime<RTstart[MU+1] && RTbank1[MU]==0 && RTbank2[MU]==0 && TMath::Abs(orbitalinfo->etha)>0.1) || ((RTbank1[MU]!=0 || RTbank2[MU]!=0) && atime>=RTstart[MU] && atime<RTstart[MU+1] && azim==0 && (errq!=0 || orbitalinfo->TimeGap>10.0 || ((modf(orbitalinfo->TimeGap,&tmptg)*1000>10 || modf(orbitalinfo->TimeGap,&tmptg)*1000==0.0) && orbitalinfo->TimeGap>2.0)))){         if((atime>=RTstart[MU] && atime<RTstart[MU+1] && RTbank1[MU]==0 && RTbank2[MU]==0 && {TMath::Abs(orbitalinfo->etha)>0.1 || tgpar0)) || ((RTbank1[MU]!=0 || RTbank2[MU]!=0) && atime>=RTstart[MU] && atime<RTstart[MU+1] && azim==0 && (errq!=0 || tgpar))){
1677          //found in Rotation Table this data for this time interval          //found in Rotation Table this data for this time interval
1678          if(atime<RTtime1[0])          if(atime<RTtime1[0])
1679            orbitalinfo->azim = 5;        //means that RotationTable no started yet            orbitalinfo->azim = 5;        //means that RotationTable no started yet
# Line 1760  if(recqtime[mu]>=1160987921.75 && recqti Line 1764  if(recqtime[mu]>=1160987921.75 && recqti
1764    
1765          orbitalinfo->Iij.ResizeTo(Iij);          orbitalinfo->Iij.ResizeTo(Iij);
1766          orbitalinfo->Iij = Iij;          orbitalinfo->Iij = Iij;
1767    
1768            Bool_t saso=true;
1769            if (orbitalinfo->qkind==1) saso=true;
1770            if (orbitalinfo->qkind==0) && orbitalinfo->azim>=0 && orbitalinfo->azim!=5 && tgpar) saso=false;
1771            if (orbitalinfo->qkind==0 && orbitalinfo->axim==5 && TMath::Abs(orbitalinfo->etha>0.1) && tgpar) saso=false;
1772            if (orbitalinfo->qkind==0 && orbitalinfo->azim==5 && TMath::Abs(orbitalinfo->etha<=0.1) && tgpar0) saso=false;
1773            if (saso) orbitalinfo->mode=orbitalinfo->rtqual; else orbitalinfo->mode=2;
1774    
1775          //          //
1776          //      A1 = Iij(0,2);          //      A1 = Iij(0,2);
1777          //      A2 = Iij(1,2);          //      A2 = Iij(1,2);

Legend:
Removed from v.1.85  
changed lines
  Added in v.1.86

  ViewVC Help
Powered by ViewVC 1.1.23