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