Ignore:
Timestamp:
03/03/03 10:59:22 (22 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Cosy/main/MBending.cc

    r1805 r1806  
    178178    p += CRY;
    179179
    180     const AltAz NRX(fNrx*kDeg2Rad*sin(p.Alt()), -fNrx*kDeg2Rad);
    181     const AltAz NRY(fNry*kDeg2Rad*cos(p.Alt()), -fNry*kDeg2Rad*tan(p.Alt()));
     180    const AltAz NRX(fNrx*sin(p.Alt()), -fNrx);
     181    const AltAz NRY(fNry*cos(p.Alt()), -fNry*tan(p.Alt()));
    182182    p += NRX;
    183183    p += NRY;
     
    186186    p += MAGIC;
    187187
    188     const AltAz AW( fAw*sin(p.Az()), -fAw*kDeg2Rad*cos(p.Az())*tan(p.Alt()));
     188    const AltAz AW( fAw*sin(p.Az()), -fAw*cos(p.Az())*tan(p.Alt()));
    189189    const AltAz AN(-fAn*cos(p.Az()), -fAn*sin(p.Az())*tan(p.Alt()));
    190190    p += AW;
    191191    p += AN;
    192192
    193     //zd1 += 11.22;
    194     //zd1 += 1.15*sin((az0+0.05)*3.1415/180);
     193//    const AltAz MAGIC(-fMagic1*sin(p.Az()-fMagic2), 0);
     194//    p += MAGIC;
     195
    195196
    196197    const AltAz CA(0, -fCa/cos(p.Alt()));
     
    212213    AltAz p = aa;
    213214
    214     const AltAz I(fIe/**kDeg2Rad*/, fIa/**kDeg2Rad*/);
     215    const AltAz I(fIe, fIa);
    215216    p -= I;
    216217
    217     const AltAz NPAE(0, -fNpae*kDeg2Rad*tan(p.Alt()));
     218    const AltAz NPAE(0, -fNpae*tan(p.Alt()));
    218219    p -= NPAE;
    219220
    220     const AltAz CA(0, -fCa*kDeg2Rad/cos(p.Alt()));
     221    const AltAz CA(0, -fCa/cos(p.Alt()));
    221222    p -= CA;
    222223
    223     const AltAz AN(-fAn*kDeg2Rad*cos(p.Az()), -fAn*kDeg2Rad*sin(p.Az())*tan(p.Alt()));
    224     const AltAz AW( fAw*kDeg2Rad*sin(p.Az()), 0/*-fAw*kDeg2Rad*cos(p.Az())*tan(p.Alt())*/);
     224    const AltAz AN(-fAn*cos(p.Az()), -fAn*sin(p.Az())*tan(p.Alt()));
     225    const AltAz AW( fAw*sin(p.Az()), -fAw*cos(p.Az())*tan(p.Alt()));
    225226    p -= AN;
    226227    p -= AW;
    227228
    228     const AltAz NRY(fNry*kDeg2Rad*cos(p.Alt()), -fNry*kDeg2Rad*tan(p.Alt()));
    229     const AltAz NRX(fNrx*kDeg2Rad*sin(p.Alt()), -fNrx*kDeg2Rad);
     229    const AltAz MAGIC(-fMagic1*sin(p.Az()-fMagic2), 0);
     230    p -= MAGIC;
     231
     232    const AltAz NRY(fNry*cos(p.Alt()), -fNry*tan(p.Alt()));
     233    const AltAz NRX(fNrx*sin(p.Alt()), -fNrx);
    230234    p -= NRY;
    231235    p -= NRX;
    232236
    233     const AltAz CRY(-fCry*kDeg2Rad*cos(p.Az()-p.Alt()), -fCry*kDeg2Rad*sin(p.Az()-p.Alt())/cos(p.Alt()));
    234     const AltAz CRX(-fCrx*kDeg2Rad*sin(p.Az()-p.Alt()),  fCrx*kDeg2Rad*cos(p.Az()-p.Alt())/cos(p.Alt()));
     237    const AltAz CRY(-fCry*cos(p.Az()-p.Alt()), -fCry*sin(p.Az()-p.Alt())/cos(p.Alt()));
     238    const AltAz CRX(-fCrx*sin(p.Az()-p.Alt()),  fCrx*cos(p.Az()-p.Alt())/cos(p.Alt()));
    235239    p -= CRY;
    236240    p -= CRX;
    237241
    238     const AltAz CEC(-fEcec*kDeg2Rad*cos(p.Alt()), -fAcec*kDeg2Rad*cos(p.Az()));
    239     const AltAz CES(-fEces*kDeg2Rad*sin(p.Alt()), -fAces*kDeg2Rad*sin(p.Az()));
     242    const AltAz CEC(-fEcec*cos(p.Alt()), -fAcec*cos(p.Az()));
     243    const AltAz CES(-fEces*sin(p.Alt()), -fAces*sin(p.Az()));
    240244    p -= CEC;
    241245    p -= CES;
Note: See TracChangeset for help on using the changeset viewer.