Changeset 1699 for trunk/MagicSoft/Cosy/main/MBending.cc

Timestamp:

01/08/03 17:21:33 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Cosy/main/MBending.cc (modified) (3 diffs)

trunk/MagicSoft/Cosy/main/MBending.cc

@@ -7 +7 @@
 
 ClassImp(MBending);
+
+void MBending::Reset()
+{
+    fIa   = 0;
+    fIe   = 0;
+    fCa   = 0;
+    fAn   = 0;
+    fAw   = 0;
+    fNrx  = 0;
+    fNry  = 0;
+    fCrx  = 0;
+    fCry  = 0;
+    fNpae = 0;
+    fEces = 0;
+    fAces = 0;
+    fEcec = 0;
+    fAcec = 0;
+}
 
 void MBending::Load(const char *name)
@@ -141 +159 @@
     p += NPAE;
 
-    const AltAz I(fIe*kDeg2Rad, fIa*kDeg2Rad);
+    const AltAz I(fIe/**kDeg2Rad*/, fIa/**kDeg2Rad*/);
     p += I;
+
+    return p;
+}
+
+AltAz MBending::CorrectBack(const AltAz &aa) const
+{
+    // Correct [rad]
+    // zdaz    [rad]
+    AltAz p = aa;
+
+    const AltAz CES(-fEces*kDeg2Rad*sin(p.Alt()), -fAces*kDeg2Rad*sin(p.Az()));
+    const AltAz CEC(-fEcec*kDeg2Rad*cos(p.Alt()), -fAcec*kDeg2Rad*cos(p.Az()));
+    p -= CES;
+    p -= CEC;
+
+    const AltAz CRX(-fCrx*kDeg2Rad*sin(p.Az()-p.Alt()),  fCrx*kDeg2Rad*cos(p.Az()-p.Alt())/cos(p.Alt()));
+    const AltAz CRY(-fCry*kDeg2Rad*cos(p.Az()-p.Alt()), -fCry*kDeg2Rad*sin(p.Az()-p.Alt())/cos(p.Alt()));
+    p -= CRX;
+    p -= CRY;
+
+    const AltAz NRX(fNrx*kDeg2Rad*sin(p.Alt()), -fNrx*kDeg2Rad);
+    const AltAz NRY(fNry*kDeg2Rad*cos(p.Alt()), -fNry*kDeg2Rad*tan(p.Alt()));
+    p -= NRX;
+    p -= NRY;
+
+    const AltAz AW( fAw*kDeg2Rad*sin(p.Az()), -fAw*kDeg2Rad*cos(p.Az())*tan(p.Alt()));
+    const AltAz AN(-fAn*kDeg2Rad*cos(p.Az()), -fAn*kDeg2Rad*sin(p.Az())*tan(p.Alt()));
+    p -= AW;
+    p -= AN;
+
+    const AltAz CA(0, -fCa*kDeg2Rad/cos(p.Alt()));
+    p -= CA;
+
+    const AltAz NPAE(0, -fNpae*kDeg2Rad*tan(p.Alt()));
+    p -= NPAE;
+
+    const AltAz I(fIe/**kDeg2Rad*/, fIa/**kDeg2Rad*/);
+    p -= I;
 
     return p;
@@ -156 +212 @@
 }
 
+ZdAz MBending::CorrectBack(const ZdAz &zdaz) const
+{
+    // Correct [rad]
+    // zdaz    [rad]
+    AltAz p(TMath::Pi()/2-zdaz.Zd(), zdaz.Az());
+    AltAz c = CorrectBack(p);
+    return ZdAz(TMath::Pi()/2-c.Alt(), c.Az());
+}
+
 void MBending::SetParameters(const Double_t *par, Int_t n=14)
 {