Changeset 4255 for trunk/MagicSoft/Cosy/main

Timestamp:

06/01/04 16:48:19 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Cosy/main

Files:

• MBending.cc (modified) (7 diffs)
• MCosy.cc (modified) (1 diff)
• MCosy.h (modified) (1 diff)
• MStarguider.cc (modified) (2 diffs)

trunk/MagicSoft/Cosy/main/MBending.cc

@@ -118 +118 @@
     fDescr[15] =  "Alt/Az Coude Displacement (N-S)";
     fDescr[16] =  "Alt/Az Coude Displacement (E-W)";
-    fDescr[17] =  "n/a";
+    fDescr[17] =  "n/a <ZA Hysteresis>";
     fDescr[18] =  "n/a";
 }
@@ -204 +204 @@
         val *= kDeg2Rad;
 
-        Double_t *dest=NULL;
-
-        if (str=="IA")   dest = &fIa;
-        if (str=="IE")   dest = &fIe;
-        if (str=="FLOP") dest = &fFlop;
-        if (str=="NPAE") dest = &fNpae;
-        if (str=="CA")   dest = &fCa;
-        if (str=="AN")   dest = &fAn;
-        if (str=="AW")   dest = &fAw;
-        if (str=="TF")   dest = &fTf;
-        if (str=="TX")   dest = &fTx;
-        if (str=="NRX")  dest = &fNrx;
-        if (str=="NRY")  dest = &fNry;
-        if (str=="CRX")  dest = &fCrx;
-        if (str=="CRY")  dest = &fCry;
-        if (str=="ECES") dest = &fEces;
-        if (str=="ACES") dest = &fAces;
-        if (str=="ECEC") dest = &fEcec;
-        if (str=="ACEC") dest = &fAcec;
-
-        if (dest)
-            *dest = val;
+        // Find parameter
+        Int_t n = -1;
+        for (int i=0; i<fNumPar; i++)
+            if (str==fName[i])
+            {
+                n = i;
+                *fCoeff[i] = val;
+                break;
+            }
 
         fin >> val;
         cout << setw(9) << val << "°" << endl;
 
-        // Find corresponding error
-        for (int i=0; i<MBending::GetNumPar(); i++)
-            if (dest==fCoeff[i])
-            {
-                fError[i] = val*kDeg2Rad;
-                break;
-            }
+        // corresponding error
+        fError[n] = val*kDeg2Rad;
     }
     cout << endl;
@@ -490 +472 @@
     p += FLOP;
 
-    //const AltAz MAGIC1(fMagic1*sin(p.Az()), 0);
-    //p += MAGIC1;
+    const AltAz MAGIC1(fMagic1*TMath::Sign(1., sin(p.Az())), 0);
+    p += MAGIC1;
 
     const AltAz I(fIe, fIa);
@@ -512 +494 @@
     const AltAz I(fIe, fIa);
     p -= I;
+
+    const AltAz MAGIC1(fMagic1*TMath::Sign(1., sin(p.Az())), 0);
+    p -= MAGIC1;
 
     //const AltAz MAGIC1(fMagic1*sin(p.Az()), 0);
@@ -622 +607 @@
 }
 
-void MBending::SetMinuitParameters(TMinuit &m, Int_t n=-1) const
+void MBending::SetMinuitParameters(TMinuit &m, Int_t n) const
 {
     if (n<0)
@@ -633 +618 @@
 }
 
-void MBending::GetMinuitParameters(TMinuit &m, Int_t n=-1)
+void MBending::GetMinuitParameters(TMinuit &m, Int_t n)
 {
     if (n<0 || n>m.GetNumPars())
@@ -665 +650 @@
 }
 */
-void MBending::PrintMinuitParameters(TMinuit &m, Int_t n=-1) const
+void MBending::PrintMinuitParameters(TMinuit &m, Int_t n) const
 {
     if (n<0)

trunk/MagicSoft/Cosy/main/MCosy.cc

@@ -1447 +1447 @@
     cout << "Open Repfile: " << name << endl;
     fOutRep = new MLog(name, kTRUE);
+    *fOutRep << "[Drive Report File]" << endl;
+    *fOutRep << "Version <cvs>" << endl;
+    *fOutRep << "Date " << MTime(-1) << endl;
+    *fOutRep << "[Reports]" << endl;
 
 /*

trunk/MagicSoft/Cosy/main/MCosy.h

@@ -201 +201 @@
     }
 
+    MLog *GetOutRep() { return fOutRep; }
+
     // static ZdAz CorrectTarget(const ZdAz &src, const ZdAz &dst);
     //    static ZdAz RaDec2ZdAz(const double mjd, const RaDec &pos, const RaDec &pm=RaDec(0,0));

trunk/MagicSoft/Cosy/main/MStarguider.cc

@@ -922 +922 @@
     // Get tracking coordinates
     const XY    xy = fCRaDec->GetCoordinates();
-    const RaDec rd(xy.X()*TMath::DegToRad()*360/24, xy.Y()*TMath::DegToRad());
+    const RaDec rd(xy.X()*TMath::DegToRad()*15, xy.Y()*TMath::DegToRad());
 
     ac.SetRaDec(rd.Ra(), rd.Dec());
 
-    // Adapt coordinate system (GUI and humans are counting Y in different directions)
+    // Adapt coordinate system (GUIs and humans are counting Y in different directions)
     Double_t x = star->GetX()-center.GetX();
     Double_t y = center.GetY()-star->GetY();
+
+    cout << "STAR-Offset: " << MTime(-1) << " dx=" << x << "pix  dy=" << dy << "pix" << endl;
 
     // Convert from Pixel to millimeter (1pix=2.9mm)
@@ -993 +995 @@
     *fOutTp << " " << setprecision(11) << t.GetMjd();
     *fOutTp << endl;
+
+    MLog &outrep = *fCosy->GetOutRep();
+    if (outrep.Lock("MStarguider::FindStar"))
+    {
+        outrep << "FINDSTAR-REPORT " << MTime(-1) << " " << setprecision(7);
+        outrep << 90-za0.Alt() << " " << za0.Az() << " ";
+        outrep << za1.Zd() << " " << za1.Az() << " ";
+        outrep << xy.X() << " " << xy.Y() << " ";
+        outrep << -dzd << " " << -daz << " ";
+        outrep << star->GetX() << " " << star->GetY() << " ";
+        outrep << center.GetX() << " " << center.GetY() << " ";
+        outrep << x*1.0713/2.58427 << " " << y*1.0713/2.58427 << " " << star->GetMag();
+        outrep << setprecision(11) << t.GetMjd() << endl;
+        outrep.UnLock("MStarguider::FindStar")
+    }
 }