Ignore:
Timestamp:
06/01/04 16:48:19 (20 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Cosy/main
Files:
4 edited

Legend:

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

    r4105 r4255  
    118118    fDescr[15] =  "Alt/Az Coude Displacement (N-S)";
    119119    fDescr[16] =  "Alt/Az Coude Displacement (E-W)";
    120     fDescr[17] =  "n/a";
     120    fDescr[17] =  "n/a <ZA Hysteresis>";
    121121    fDescr[18] =  "n/a";
    122122}
     
    204204        val *= kDeg2Rad;
    205205
    206         Double_t *dest=NULL;
    207 
    208         if (str=="IA")   dest = &fIa;
    209         if (str=="IE")   dest = &fIe;
    210         if (str=="FLOP") dest = &fFlop;
    211         if (str=="NPAE") dest = &fNpae;
    212         if (str=="CA")   dest = &fCa;
    213         if (str=="AN")   dest = &fAn;
    214         if (str=="AW")   dest = &fAw;
    215         if (str=="TF")   dest = &fTf;
    216         if (str=="TX")   dest = &fTx;
    217         if (str=="NRX")  dest = &fNrx;
    218         if (str=="NRY")  dest = &fNry;
    219         if (str=="CRX")  dest = &fCrx;
    220         if (str=="CRY")  dest = &fCry;
    221         if (str=="ECES") dest = &fEces;
    222         if (str=="ACES") dest = &fAces;
    223         if (str=="ECEC") dest = &fEcec;
    224         if (str=="ACEC") dest = &fAcec;
    225 
    226         if (dest)
    227             *dest = val;
     206        // Find parameter
     207        Int_t n = -1;
     208        for (int i=0; i<fNumPar; i++)
     209            if (str==fName[i])
     210            {
     211                n = i;
     212                *fCoeff[i] = val;
     213                break;
     214            }
    228215
    229216        fin >> val;
    230217        cout << setw(9) << val << "°" << endl;
    231218
    232         // Find corresponding error
    233         for (int i=0; i<MBending::GetNumPar(); i++)
    234             if (dest==fCoeff[i])
    235             {
    236                 fError[i] = val*kDeg2Rad;
    237                 break;
    238             }
     219        // corresponding error
     220        fError[n] = val*kDeg2Rad;
    239221    }
    240222    cout << endl;
     
    490472    p += FLOP;
    491473
    492     //const AltAz MAGIC1(fMagic1*sin(p.Az()), 0);
    493     //p += MAGIC1;
     474    const AltAz MAGIC1(fMagic1*TMath::Sign(1., sin(p.Az())), 0);
     475    p += MAGIC1;
    494476
    495477    const AltAz I(fIe, fIa);
     
    512494    const AltAz I(fIe, fIa);
    513495    p -= I;
     496
     497    const AltAz MAGIC1(fMagic1*TMath::Sign(1., sin(p.Az())), 0);
     498    p -= MAGIC1;
    514499
    515500    //const AltAz MAGIC1(fMagic1*sin(p.Az()), 0);
     
    622607}
    623608
    624 void MBending::SetMinuitParameters(TMinuit &m, Int_t n=-1) const
     609void MBending::SetMinuitParameters(TMinuit &m, Int_t n) const
    625610{
    626611    if (n<0)
     
    633618}
    634619
    635 void MBending::GetMinuitParameters(TMinuit &m, Int_t n=-1)
     620void MBending::GetMinuitParameters(TMinuit &m, Int_t n)
    636621{
    637622    if (n<0 || n>m.GetNumPars())
     
    665650}
    666651*/
    667 void MBending::PrintMinuitParameters(TMinuit &m, Int_t n=-1) const
     652void MBending::PrintMinuitParameters(TMinuit &m, Int_t n) const
    668653{
    669654    if (n<0)
  • trunk/MagicSoft/Cosy/main/MCosy.cc

    r4153 r4255  
    14471447    cout << "Open Repfile: " << name << endl;
    14481448    fOutRep = new MLog(name, kTRUE);
     1449    *fOutRep << "[Drive Report File]" << endl;
     1450    *fOutRep << "Version <cvs>" << endl;
     1451    *fOutRep << "Date " << MTime(-1) << endl;
     1452    *fOutRep << "[Reports]" << endl;
    14491453
    14501454/*
  • trunk/MagicSoft/Cosy/main/MCosy.h

    r4154 r4255  
    201201    }
    202202
     203    MLog *GetOutRep() { return fOutRep; }
     204
    203205    // static ZdAz CorrectTarget(const ZdAz &src, const ZdAz &dst);
    204206    //    static ZdAz RaDec2ZdAz(const double mjd, const RaDec &pos, const RaDec &pm=RaDec(0,0));
  • trunk/MagicSoft/Cosy/main/MStarguider.cc

    r4154 r4255  
    922922    // Get tracking coordinates
    923923    const XY    xy = fCRaDec->GetCoordinates();
    924     const RaDec rd(xy.X()*TMath::DegToRad()*360/24, xy.Y()*TMath::DegToRad());
     924    const RaDec rd(xy.X()*TMath::DegToRad()*15, xy.Y()*TMath::DegToRad());
    925925
    926926    ac.SetRaDec(rd.Ra(), rd.Dec());
    927927
    928     // Adapt coordinate system (GUI and humans are counting Y in different directions)
     928    // Adapt coordinate system (GUIs and humans are counting Y in different directions)
    929929    Double_t x = star->GetX()-center.GetX();
    930930    Double_t y = center.GetY()-star->GetY();
     931
     932    cout << "STAR-Offset: " << MTime(-1) << " dx=" << x << "pix  dy=" << dy << "pix" << endl;
    931933
    932934    // Convert from Pixel to millimeter (1pix=2.9mm)
     
    993995    *fOutTp << " " << setprecision(11) << t.GetMjd();
    994996    *fOutTp << endl;
     997
     998    MLog &outrep = *fCosy->GetOutRep();
     999    if (outrep.Lock("MStarguider::FindStar"))
     1000    {
     1001        outrep << "FINDSTAR-REPORT " << MTime(-1) << " " << setprecision(7);
     1002        outrep << 90-za0.Alt() << " " << za0.Az() << " ";
     1003        outrep << za1.Zd() << " " << za1.Az() << " ";
     1004        outrep << xy.X() << " " << xy.Y() << " ";
     1005        outrep << -dzd << " " << -daz << " ";
     1006        outrep << star->GetX() << " " << star->GetY() << " ";
     1007        outrep << center.GetX() << " " << center.GetY() << " ";
     1008        outrep << x*1.0713/2.58427 << " " << y*1.0713/2.58427 << " " << star->GetMag();
     1009        outrep << setprecision(11) << t.GetMjd() << endl;
     1010        outrep.UnLock("MStarguider::FindStar")
     1011    }
    9951012}
    9961013
Note: See TracChangeset for help on using the changeset viewer.