Changeset 18627 for trunk/FACT++/drive


Ignore:
Timestamp:
09/18/16 14:53:16 (8 years ago)
Author:
tbretz
Message:
Disentengle from root and Mars were possible, clean up of gui from all obsolete stuff
Location:
trunk/FACT++/drive
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/FACT++/drive/MStarguider.cc

    r18618 r18627  
    44#include "MStarguider.h"
    55
    6 #include <fstream>    // ifstream
    7 
    8 #include <TH2F.h>
    9 #include <TGraph.h>
    10 #include <TTimer.h>
     6#ifdef HAVE_PNG
     7#include <png.h>
     8#endif
     9#include <iostream>
     10#include <fstream>
     11
     12#include <math.h>
     13
    1114#include <TEnv.h>
    1215#include <TSystem.h>
    13 #include <TFile.h> // temp writeout of histos
    14 //#include <TSocket.h>
    1516
    1617#include <TGMenu.h>
     
    2021#include <TGTextEntry.h>
    2122#include <TGLayout.h>
    22 #include "MString.h"
    23 
    24 #include "MLog.h"
    25 #include "MLogManip.h"
    26 
    27 #include "MAstro.h"
    28 #include "MString.h"
    29 
    30 
    31 #include "MCosy.h"
     23
    3224#include "MCaos.h"
    33 // #include "MStargLeds.h"
    34 #include "MStargHistograms.h"
    35 #include "MDriveCom.h"
    36 
    37 #include "MGStarg.h"
    38 #include "MGNumStars.h"
    39 
    4025#include "MGImage.h"
    41 #include "MGCoordinates.h"
    42 
    4326#include "Camera.h"
    44 #include "PngReader.h"
    45 
     27#include "../src/DimSetup.h"
    4628#include "Led.h"
    47 #include "Writer.h"
     29//#include "Writer.h"
    4830#include "FilterLed.h"
    49 #include "MStarList.h"
    50 #include "CaosFilter.h"
    51 #include "StarCatalog.h"
    52 #include "MGeomCamMagic.h"
    53 #include "MAstroCamera.h"
    54 
    55 #include "MGCosy.h"
    56 
    57 ClassImp(MStarguider);
     31//#include "CaosFilter.h"
    5832
    5933using namespace std;
     
    6337    IDM_kFindStar,
    6438    IDM_kCaosFilter,
    65     IDM_kCatalog,
    6639    IDM_kStarguider,
    67     IDM_kStart,
    68     IDM_kStop,
    69     IDM_kFileType,
    70     IDM_kPPM,
    71     IDM_kPNG,
    72     IDM_kOnce,
    7340    IDM_kStretch,
    7441    IDM_kInput,
    75     IDM_kChannel1,
    76     IDM_kChannel2,
    77     IDM_kChannel3,
    78     IDM_kContinous,
    79     IDM_kRate25ps,
    80     IDM_kRate5ps,
    81     IDM_kRate1s,
    82     IDM_kRate5s,
    83     IDM_kRate30s,
    84     IDM_kRate1m,
    85     IDM_kRate5m,
    8642    IDM_kSetup,
    87     IDM_kLimMag3,
    88     IDM_kLimMag4,
    89     IDM_kLimMag5,
    90     IDM_kLimMag6,
    91     IDM_kLimMag7,
    92     IDM_kLimMag8,
    93     IDM_kLimMag9,
    94     IDM_kPixSize,
    95     IDM_kAngle,
    9643    IDM_kCut,
    9744    IDM_kInterpol250,
     
    10350    IDM_kInterpol2,
    10451    IDM_kInterpol1,
    105     IDM_kCaosPrintRings,
    106     IDM_kCaosPrintLeds,
    107     IDM_kCaosAnalStart,
    108     IDM_kCaosAnalStop,
     52    //IDM_kCaosPrintRings,
     53    //IDM_kCaosPrintLeds,
    10954    IDM_kCaosWriteStart,
    11055    IDM_kCaosWriteStop,
    11156    IDM_kResetHistograms,
    112     IDM_kStargAnalysis,
    113     IDM_kStargCaosFilter,
    114     IDM_kStarguiderMode,
    115     IDM_kTpointMode
    116 
    11757};
    11858
    119 Bool_t MStarguider::HandleTimer(TTimer *t)
     59Bool_t MStarguider::HandleTimer(TTimer *)
    12060{
    12161    if (IsMapped())
    122     {
    12362        fImage->DoRedraw();
    124         fZoomImage->DoRedraw();
    125     }
    126 
    127     if (fCosy && fCosy->GetWin()->IsMapped())
    128         fCosy->GetWin()->GetImage()->DoRedraw();
    129 
    130     fGStarg->Update(fPos, fD);   
    131     if (fDisplay->IsEntryChecked(IDM_kStarguider))
    132         fGNumStars->Update(fNumStarsDetected, fNumStarsCorrelated);
    13363 
    13464    return kTRUE;
     
    13666
    13767#define kZOOM 96
    138 
     68/*
    13969XY MStarguider::GetCoordinates() const
    14070{
    14171    return fPZdAz->GetCoordinates();
    14272}
    143 
    144 void MStarguider::InitGui(Int_t channel)
    145 {
    146     fList = new MGList;
    147     fList->SetOwner();
    148 
    149     const TGWindow *p=gClient->GetRoot();
     73*/
     74void MStarguider::InitGui(Int_t)
     75{
     76    //fList = new MGList;
     77    //fList->SetOwner();
     78
     79/*
     80 const TGWindow *p=gClient->GetRoot();
    15081
    15182    fChannel = new TGPopupMenu(p);
    15283    fChannel->AddEntry("Starfield Camera", IDM_kChannel1);
    15384    fChannel->AddEntry("TPoint Camera",    IDM_kChannel2);
    154     fChannel->AddEntry("Read from File",   IDM_kChannel3);
    155     if (channel<0)
    156         fChannel->CheckEntry(IDM_kChannel3);
    157     else
     85    //fChannel->AddEntry("Read from File",   IDM_kChannel3);
     86    //if (channel<0)
     87    //    fChannel->CheckEntry(IDM_kChannel3);
     88    //else
    15889        fChannel->CheckEntry(channel==0?IDM_kChannel1:IDM_kChannel2);
    15990    fChannel->Associate(this);
     
    199130    fLimMag->Associate(this);
    200131    fList->Add(fLimMag);
    201 
    202     fSao->SetLimitMag(9.0);
     132    */
     133    //fSao->SetLimitMag(9.0);
     134
     135    const TGWindow *p=gClient->GetRoot();
    203136
    204137    fInterpol = new TGPopupMenu(p);
     
    212145    fInterpol->AddEntry("Off", IDM_kInterpol1);
    213146    fInterpol->Associate(this);
    214     fList->Add(fInterpol);
     147    //fList->Add(fInterpol);
    215148
    216149    TString disp=gVirtualX->DisplayName();
    217     gLog << all << "Display: " << disp << endl;
     150    cout << "Display: " << disp << endl;
    218151    if (disp.First(':')>=0)
    219152        disp=disp(0, disp.First(':'));
     
    221154    if (disp.IsNull() || disp==(TString)"localhost")
    222155    {
    223         fInterpol->CheckEntry(IDM_kInterpol5);
    224         fIntRate = 50;
     156        fInterpol->CheckEntry(IDM_kInterpol25);
     157        fIntRate = 25;
    225158    }
    226159    else
     
    230163    }
    231164
     165/*
    232166    fCaosPrint = new TGPopupMenu(p);
    233167    fCaosPrint->AddEntry("&Leds",  IDM_kCaosPrintLeds);
     
    249183    fCaosAnalyse->Associate(this);
    250184    fList->Add(fCaosAnalyse);
    251 
     185*/
    252186    fMenu = new TGMenuBar(this, 0, 0, kHorizontalFrame);
    253187    fDisplay       = fMenu->AddPopup("&Display");
    254     fMode          = fMenu->AddPopup("&Mode");
    255     fWritePictures = fMenu->AddPopup("&WritePics");
     188    //fMode          = fMenu->AddPopup("&Mode");
     189    //fWritePictures = fMenu->AddPopup("&WritePics");
    256190    fSetup         = fMenu->AddPopup("&Setup");
    257     fOperations    = fMenu->AddPopup("&Operations");
     191    //fOperations    = fMenu->AddPopup("&Operations");
    258192    fMenu->Resize(fMenu->GetDefaultSize());
    259     AddFrame(fMenu); 
     193    AddFrame(fMenu);
    260194
    261195    //
     
    268202    fDisplay->AddEntry("Find &Star",            IDM_kFindStar);
    269203    fDisplay->AddEntry("C&aos Filter",          IDM_kCaosFilter);
    270     fDisplay->AddSeparator();
    271     fDisplay->AddEntry("SAO &Catalog",          IDM_kCatalog);
    272     fDisplay->AddEntry("Starguider",            IDM_kStarguider);
    273     fDisplay->AddEntry("Starguider LED Filter", IDM_kStargCaosFilter);
    274     fDisplay->AddSeparator();   
    275     if (channel>=0)
    276         fDisplay->AddPopup("&Input",   fChannel);
    277     fDisplay->CheckEntry(IDM_kStretch);
     204    //fDisplay->AddSeparator();
     205    //if (channel>=0)
     206    //    fDisplay->AddPopup("&Input",   fChannel);
     207    // fDisplay->CheckEntry(IDM_kStretch);
     208    fDisplay->CheckEntry(IDM_kFindStar);
     209    fDisplay->CheckEntry(IDM_kCaosFilter);
    278210    fDisplay->Associate(this);
    279211
    280     fMode->AddEntry("Starguider", IDM_kStarguiderMode);
    281     fMode->AddEntry("Tpoint",     IDM_kTpointMode);
    282     fMode->Associate(this);
    283 
     212    //fMode->AddEntry("Tpoint",     IDM_kTpointMode);
     213    //fMode->Associate(this);
     214/*
    284215    fWritePictures->AddEntry("&Start",      IDM_kStart);
    285216    fWritePictures->AddEntry("Sto&p",       IDM_kStop);
    286217    fWritePictures->AddSeparator();
    287     fWritePictures->AddPopup("File &Type",  fFileType);
     218    //fWritePictures->AddPopup("File &Type",  fFileType);
    288219    fWritePictures->AddPopup("&Write Type", fWriteType);
    289220    fWritePictures->AddPopup("Write &Rate", fWriteRate);
    290221    fWritePictures->DisableEntry(IDM_kStop);
    291222    fWritePictures->Associate(this);
    292 
    293     fSetup->AddPopup("Lim. &Magnitude",      fLimMag);
     223    */
    294224    fSetup->AddPopup("Disp. &Interpolation", fInterpol);
    295225    fSetup->Associate(this);
    296226
    297     fOperations->AddEntry("Starguider Analysis", IDM_kStargAnalysis);
    298     fOperations->DisableEntry(IDM_kStargAnalysis);
    299     fOperations->Associate(this);
    300 
     227/*
    301228    fCaOs = new TGPopupMenu(p);
    302     fCaOs->AddPopup("&Write",   fCaosWrite);
     229    //fCaOs->AddPopup("&Write",   fCaosWrite);
    303230    fCaOs->AddPopup("&Print",   fCaosPrint);
    304     fCaOs->AddPopup("&Analyse", fCaosAnalyse);
     231    //fCaOs->AddPopup("&Analyse", fCaosAnalyse);
    305232    fCaOs->Associate(this);
    306233    fList->Add(fCaOs);
    307 
    308     fGStarg = new MGStarg(this, 235);
    309     fGStarg->Move(530,596+5);
    310     fList->Add(fGStarg);
    311 
     234*/
    312235    //RA,Dec for catalog
     236    /*
    313237    fCRaDec = new MGCoordinates(this, kETypeRaDec);
    314238    fCRaDec->Move(4, fMenu->GetDefaultHeight()+584);
     
    344268    fTPoint->AllowStayDown(kTRUE);
    345269    AddFrame(fTPoint);
    346     /*
    347      fStargTPoint = new TGTextButton(this, "StargTPoint");
     270
     271    fStargTPoint = new TGTextButton(this, "StargTPoint");
    348272     fStargTPoint->Move(170, fMenu->GetDefaultHeight()+785);
    349273     fStargTPoint->AllowStayDown(kTRUE);
    350274     AddFrame(fStargTPoint);
    351     */
    352 
    353     fFps = new TGLabel(this, "---fps");
     275
     276     fFps = new TGLabel(this, "---fps");
    354277    fFps->SetTextJustify(kTextRight);
    355278    fFps->Move(650-495, fMenu->GetDefaultHeight()+714+23);
     
    395318    fDZdAzText->Move(240+12+20+7, fMenu->GetDefaultHeight()+676+2*20-5-46);
    396319    AddFrame(fDZdAzText);
    397 
     320*/
    398321    // Set input box for rotation angle
     322    /*
    399323    fAngle = new TGTextEntry(this, "           ", IDM_kAngle);
    400324    fAngle->SetAlignment(kTextCenterX);
     
    402326    AddFrame(fAngle);
    403327
    404     SetRotationAngle(-0.2);
     328    //SetRotationAngle(-0.2);
    405329
    406330    // Set input box for pixel size
     
    409333    fPixSize->Move(547-410, fMenu->GetDefaultHeight()+690);
    410334    AddFrame(fPixSize);
    411 
    412     SetPixSize(48.9);
     335    */
     336    //SetPixSize(48.9);
    413337
    414338    // Set input box for cleaning cut
    415     fCut = new TGTextEntry(this, "           ", IDM_kCut);
    416     fCut->SetAlignment(kTextCenterX);
    417     fCut->Move(547-410, fMenu->GetDefaultHeight()+713);
    418     AddFrame(fCut);
    419 
    420     SetCut(3.0);
     339    //fCut = new TGTextEntry(this, "           ", IDM_kCut);
     340    //fCut->SetAlignment(kTextCenterX);
     341    //fCut->Move(547-410, fMenu->GetDefaultHeight()+713);
     342    //AddFrame(fCut);
     343
     344    //SetCut(3.0);
    421345
    422346    // TGHorizontal3DLine *fLineSep = new TGHorizontal3DLine(this);
     
    426350    //
    427351    // Create Image Display
    428     //
     352    /*
    429353    fZoomImage = new MGImage(this, kZOOM, kZOOM);
    430354    // fZoomImage->Move(768-kZOOM-2, 700-kZOOM-2);
    431355    fZoomImage->Move(4, 700-kZOOM-2+85);
    432356    AddFrame(fZoomImage);
    433 
     357    */
    434358    fImage = new MGImage(this, 768, 576);
    435359    fImage->Move(0, fMenu->GetDefaultHeight());
     
    437361
    438362    const Int_t w = 768;
    439     const Int_t h = 840;
     363    const Int_t h = 576;
    440364    SetWMSizeHints(w, h, w, h, 1, 1);  // set the smallest and biggest size of the Main frame
    441365
     
    443367    // Make everything visible
    444368    //
    445     SetWindowName("MStarguider Main Window");
    446     SetIconName("MStarguider");
     369    SetWindowName("TPoint Main Window");
     370    SetIconName("TPoint");
    447371
    448372    MapSubwindows();
    449     fTPoint->UnmapWindow();
     373    //fTPoint->UnmapWindow();
    450374    //fStargTPoint->UnmapWindow();
    451     fGStarg->UnmapWindow();
    452     fGNumStars->UnmapWindow();
    453     fCRaDec->UnmapWindow();
    454     fCZdAz->UnmapWindow();
    455     fCZdAzText->UnmapWindow();
    456     fPZdAz->UnmapWindow();
    457     fPZdAzText->UnmapWindow();
    458     fDZdAz->UnmapWindow();
    459     fDZdAzText->UnmapWindow();
    460     fSZdAz->UnmapWindow();
    461     fSkyBright->UnmapWindow();
     375    //fGStarg->UnmapWindow();
     376    //fGNumStars->UnmapWindow();
     377    //fCRaDec->UnmapWindow();
     378    //fCZdAz->UnmapWindow();
     379    //fCZdAzText->UnmapWindow();
     380    //fPZdAz->UnmapWindow();
     381    //fPZdAzText->UnmapWindow();
     382    //fDZdAz->UnmapWindow();
     383    //fDZdAzText->UnmapWindow();
     384    //fSZdAz->UnmapWindow();
     385    //fSkyBright->UnmapWindow();
    462386    MapWindow();
    463387
     
    472396}
    473397
    474 MStarguider::MStarguider(MObservatory::LocationName_t obs, Int_t channel)
    475     : TGMainFrame(gClient->GetRoot(), 768, 840),
    476       fNumStarsDetected(0),
    477       fNumStarsCorrelated(0),
    478       fCosy(NULL),
    479       fOutTp(0),
    480       fOutStargTp(0),
    481 //      fOutRq(0),
    482       fDx((768-kZOOM)/2),
    483       fDy((512-kZOOM)/2),
    484       fStatus(MDriveCom::kStandby),
    485       fRadius(200),
    486       fTPointFromCC(-1),
    487       fFindStarCut(3.0), fFindStarBox(70)
    488 {
    489     gLog << warn << " #### FIXME: Make MCaos Thread safe!" << endl;
     398MStarguider::MStarguider(Int_t channel) : TGMainFrame(gClient->GetRoot(), 768, 840),
     399fDimData("TPOINT/DATA", "D:11", (void*)NULL, 0),
     400fDimTPoint("TPOINT/EXECUTE", "", this),
     401fDimScreenshot("TPOINT/SCREENSHOT", "B:1;C", this),
     402fRadius(200), fFindStarCut(2.), fFindStarBox(30), fTPointMode(0)
     403{
     404    cout << " #### FIXME: Make MCaos Thread safe!" << endl;
    490405
    491406    // This means that all objects added with AddFrame are deleted
     
    493408    SetCleanup();
    494409
    495     fSao = new StarCatalog(obs);
    496     fRaDec = new RaDec(180, 40);
    497 
    498410    fCaos = new MCaos;
    499     fStargCaos = new MCaos;
    500 
    501     fStargHistograms = new MStargHistograms();
     411    fCaos->ReadResources("leds_fact.txt");
    502412
    503413    InitGui(channel);
     
    506416    fTimer->TurnOn();
    507417
    508     fTime.Now();
    509 
    510     fTimeFromTp.Clear();
    511     fAltAzOffsetFromTp = AltAz(-1000,-1000);
    512                                              
    513418    gVirtualX->GrabButton(fId, kButton2, 0, 0, 0, 0, kTRUE);
    514419
    515     if (channel<0)
    516         fGetter = new PngReader(*this);
    517     else
    518         fGetter = new Camera(*this, channel);
    519 }
    520 
    521 void MStarguider::SetRotationAngle(Double_t angle)
    522 {
    523     fSao->SetRotationAngle(angle);
    524 
    525     TString txt;
    526     txt += angle;
    527 
    528     fAngle->SetText(txt);
    529 }
    530 
    531 void MStarguider::SetPixSize(Double_t size)
    532 {
    533     fSao->SetPixSize(size);
    534 
    535     TString txt;
    536     txt += size;
    537 
    538     fPixSize->SetText(txt);
    539 }
    540 
    541 void MStarguider::SetCut(Double_t cut)
    542 {
    543     TString txt;
    544     txt += cut;
    545 
    546     fCut->SetText(txt);
    547 }
    548 
     420    fGetter = new Camera(*this, channel);
     421
     422    DimClient::setNoDataCopy();
     423    DimServer::start("TPOINT");
     424}
     425
     426MStarguider::~MStarguider()
     427{
     428    DimServer::stop();
     429
     430    fGetter->ExitLoop();
     431    delete fGetter;
     432
     433    gVirtualX->GrabButton(fId, kButton2, 0, 0, 0, 0, kFALSE);
     434
     435    fTimer->TurnOff();
     436    delete fTimer;
     437
     438    delete fInterpol;
     439
     440    //delete fList;
     441    delete fCaos;
     442
     443    cout << "Camera Display destroyed." << endl;
     444}
     445
     446/*
    549447void MStarguider::SetupEnv(TEnv &env)
    550448{
    551449    fCaos->ReadEnv(env, "TPointLeds", kTRUE);
    552     fStargCaos->ReadEnv(env, "StarguiderLeds", kTRUE);
    553 
    554     SetRotationAngle(env.GetValue("Starguider.RotationAngle", fSao->GetRotationAngle()));
    555     SetCut(env.GetValue("Starguider.CleaningLevel", atof(fCut->GetText())));
    556 
    557     SetPixSize(env.GetValue("StarguiderLeds.ArcsecPerPixel", fSao->GetPixSize()));
     450    //fStargCaos->ReadEnv(env, "StarguiderLeds", kTRUE);
     451
     452    //SetRotationAngle(env.GetValue("Starguider.RotationAngle", fSao->GetRotationAngle()));
     453    //SetCut(env.GetValue("Starguider.CleaningLevel", atof(fCut->GetText())));
     454
     455    //SetPixSize(env.GetValue("StarguiderLeds.ArcsecPerPixel", fSao->GetPixSize()));
    558456
    559457    fRadius = env.GetValue("Leds.Radius", fRadius);
     
    570468    fFindStarCut = env.GetValue("FindStar.CleaningLevel", fFindStarCut);
    571469}
    572 
    573 MStarguider::~MStarguider()
    574 {
    575     fGetter->ExitLoop();
    576     delete fGetter;
    577 
    578     gVirtualX->GrabButton(fId, kButton2, 0, 0, 0, 0, kFALSE);
    579 
    580     fTimer->TurnOff();
    581     delete fTimer;
    582 
    583     delete fList;
    584 
    585     delete fCaos;
    586     delete fStargCaos;
    587     delete fStargHistograms;
    588     delete fSao;
    589     delete fRaDec;
    590 
    591     if (fOutTp)
    592         delete fOutTp;
    593 
    594     if (fOutStargTp)
    595         delete fOutStargTp;
    596 
    597 //    if (fOutRq)
    598 //        delete fOutRq;
    599 
    600     gLog << inf2 << "Camera Display destroyed." << endl;
    601 }
     470*/
    602471
    603472void MStarguider::Layout()
     
    608477void MStarguider::CloseWindow()
    609478{
    610     gLog << inf2 << "EventDisplay::CloseWindow: Exit Application Loop." << endl;
     479    cout << "EventDisplay::CloseWindow: Exit Application Loop." << endl;
    611480
    612481    //fClient.ExitLoop();
     
    622491}
    623492
     493/*
    624494void MStarguider::SetChannel()
    625495{
     
    643513            fGetter->SetChannel(ch);
    644514    }
    645 }
     515}*/
    646516
    647517void MStarguider::Toggle(TGPopupMenu *p, UInt_t id)
     
    653523}
    654524
    655 void MStarguider::ToggleStargAnalysis()
    656 {
    657     if (fOperations->IsEntryChecked(IDM_kStargAnalysis))
    658         fStargHistograms->OpenFile();
    659     else
    660         fStargHistograms->CloseFile();
    661 }
    662 
    663 /*
    664 void MStarguider::ToggleFindStar()
    665 {
    666     if (fDisplay->IsEntryChecked(IDM_kFindStar) && fCosy)
    667     fTPoint->MapWindow();
    668     else
    669     {
    670         fTPoint->UnmapWindow();
    671         fTPoint->SetDown(kFALSE);
    672     }
    673 }
    674 */
    675 void MStarguider::ToggleStarguider()
    676 {
    677     if (fDisplay->IsEntryChecked(IDM_kStarguider))
    678     {
    679         fLastBright = 0xff;
    680 
    681         fChannel->EnableEntry(IDM_kChannel1);
    682         fDisplay->DisableEntry(IDM_kFindStar);
    683         fDisplay->DisableEntry(IDM_kCaosFilter);
    684         fOperations->EnableEntry(IDM_kStargAnalysis);
    685 
    686         //fStargTPoint->MapWindow();
    687 
    688         fPZdAz->MapWindow();
    689         fDZdAz->MapWindow();
    690         fSkyBright->MapWindow();
    691 
    692         for (int i=IDM_kLimMag3; i<=IDM_kLimMag9; i++)
    693             fLimMag->UnCheckEntry(i);
    694         fLimMag->CheckEntry(IDM_kLimMag9);
    695         fSao->SetLimitMag(9.0);
    696 
    697         for (int i=IDM_kInterpol250; i<=IDM_kInterpol1; i++)
    698             fInterpol->UnCheckEntry(i);
    699         fInterpol->CheckEntry(IDM_kInterpol125);
    700         fIntRate = 125;
    701 
    702         //catalog
    703         fDisplay->CheckEntry(IDM_kCatalog);
    704         fCRaDec->MapWindow();
    705         fCZdAz->MapWindow();
    706         fCZdAzText->MapWindow();
    707 
    708         fPZdAz->MapWindow();
    709         fPZdAzText->MapWindow();
    710         fDZdAz->MapWindow();
    711         fDZdAzText->MapWindow();
    712 
    713         fGStarg->MapWindow();
    714         fGNumStars->MapWindow();
    715 
    716         SwitchOff(fChannel, IDM_kChannel2);
    717         fChannel->CheckEntry(IDM_kChannel1);
    718     }
    719     else
    720     {
    721         fStatus = MDriveCom::kStandby;
    722 
    723         //fStargTPoint->UnmapWindow();
    724         //fStargTPoint->SetDown(kFALSE);
    725 
    726         fPZdAz->UnmapWindow();
    727         fPZdAzText->UnmapWindow();
    728         fDZdAz->UnmapWindow();
    729         fDZdAzText->UnmapWindow();
    730         fSkyBright->UnmapWindow();
    731         //catalog
    732         fCRaDec->UnmapWindow();
    733         fCZdAz->UnmapWindow();
    734         fCZdAzText->UnmapWindow();
    735         fDisplay->UnCheckEntry(IDM_kCatalog);
    736 
    737         fGStarg->UnmapWindow();
    738         fGNumStars->UnmapWindow();
    739 
    740         fChannel->EnableEntry(IDM_kChannel2);
    741         fDisplay->EnableEntry(IDM_kFindStar);
    742         fDisplay->EnableEntry(IDM_kCaosFilter);
    743         fOperations->DisableEntry(IDM_kStargAnalysis);
    744     }
    745 
    746     SetChannel();
    747     //gSystem->Unlink("tracking_error.txt");
    748 }
    749 
    750 void MStarguider::ToggleCaosFilter()
    751 {
    752     if (fDisplay->IsEntryChecked(IDM_kCaosFilter))
    753     {
    754         if (!fMode->IsEntryChecked(IDM_kStarguiderMode)
    755             && !fMode->IsEntryChecked(IDM_kTpointMode))
    756             fMenu->AddPopup("&CaOs", fCaOs, NULL);
    757 
    758         SwitchOff(fDisplay, IDM_kStarguider);
    759         SwitchOff(fDisplay, IDM_kStargCaosFilter);
    760     }
    761     else
    762     {
    763         fCaosWrite->UnCheckEntry(IDM_kCaosPrintLeds);
    764         fCaosWrite->UnCheckEntry(IDM_kCaosPrintRings);
    765 
    766         if (fCaosAnalyse->IsEntryEnabled(IDM_kCaosAnalStop))
    767         {
    768             fCaosAnalyse->DisableEntry(IDM_kCaosAnalStop);
    769             fCaosAnalyse->EnableEntry(IDM_kCaosAnalStart);
    770             fCaos->DeleteHistograms();
    771         }
    772         if (fCaosWrite->IsEntryEnabled(IDM_kCaosWriteStop))
    773         {
    774             fCaosWrite->DisableEntry(IDM_kCaosWriteStop);
    775             fCaosWrite->EnableEntry(IDM_kCaosWriteStart);
    776             fCaos->CloseFile();
    777         }
    778         fMenu->RemovePopup("CaOs");
    779 
    780         fDisplay->EnableEntry(IDM_kStargCaosFilter);
    781         fDisplay->EnableEntry(IDM_kStarguider);
    782     }
    783     fMenu->Resize(fMenu->GetDefaultSize());
    784     MapSubwindows(); // maps everything, but we don't want that
    785     fTPoint->UnmapWindow();
    786     //fStargTPoint->UnmapWindow();
    787     fGStarg->UnmapWindow();
    788     fGNumStars->UnmapWindow();
    789     fCRaDec->UnmapWindow();
    790     fCZdAz->UnmapWindow();
    791     fCZdAzText->UnmapWindow();
    792     fPZdAz->UnmapWindow();
    793     fPZdAzText->UnmapWindow();
    794     fDZdAz->UnmapWindow();
    795     fDZdAzText->UnmapWindow();
    796     fSZdAz->UnmapWindow();
    797     fSkyBright->UnmapWindow();
    798     MapWindow();
    799 }
    800 
    801 Bool_t MStarguider::ProcessMessage(Long_t msg, Long_t mp1, Long_t mp2)
     525
     526Bool_t MStarguider::ProcessMessage(Long_t msg, Long_t mp1, Long_t)
    802527{
    803528    switch (GET_MSG(msg))
     
    807532            switch (mp1)
    808533            {
     534                /*
    809535            case IDM_kPixSize:
    810536                {
     
    826552                    gLog << all << "Starguider cleaning level changed to " << cut << " sigma." << endl;
    827553                    return kTRUE;
    828                 }
     554                }*/
    829555            }
    830556        return kTRUE;
     
    836562            switch (mp1)
    837563            {
    838             case IDM_kCatalog:
    839                 Toggle(fDisplay, IDM_kCatalog);
    840                 if (fDisplay->IsEntryChecked(IDM_kCatalog))
    841                 {
    842                     fCRaDec->MapWindow();
    843                     fCZdAz->MapWindow();
    844                     fCZdAzText->MapWindow();
    845                 }
    846                 else
    847                 {
    848                     fCRaDec->UnmapWindow();
    849                     fCZdAz->UnmapWindow();
    850                     fCZdAzText->UnmapWindow();
    851                 }
    852                 return kTRUE;
    853 
    854             case IDM_kStarguider:
    855                 Toggle(fDisplay, IDM_kStarguider);
    856                 ToggleStarguider();
    857                 return kTRUE;
    858 
    859             case IDM_kStargAnalysis:
    860                 Toggle(fOperations, IDM_kStargAnalysis);
    861                 ToggleStargAnalysis();
    862                 return kTRUE;
    863 
    864             case IDM_kStarguiderMode:
    865                 Toggle(fMode, IDM_kStarguiderMode);
    866 
    867                 if (fMode->IsEntryChecked(IDM_kStarguiderMode))
    868                 {
    869                     //uncheck not needed items
    870                     //general
    871                     SwitchOff(fDisplay, IDM_kFilter);
    872                     SwitchOff(fChannel, IDM_kChannel3);
    873                     SwitchOff(fOperations, IDM_kStargAnalysis);
    874 
    875                     //tpoint
    876                     SwitchOff(fDisplay, IDM_kFindStar);
    877                     SwitchOff(fDisplay, IDM_kCaosFilter);
    878                     ToggleCaosFilter();
    879                     fMode->UnCheckEntry(IDM_kTpointMode);
    880 
    881                     //enable starguider items
    882                     fDisplay->EnableEntry(IDM_kStargCaosFilter);
    883                     fDisplay->EnableEntry(IDM_kStarguider);
    884                     fDisplay->EnableEntry(IDM_kCatalog);
    885 
    886                     //check needed items
    887                     fDisplay->UnCheckEntry(IDM_kStretch);
    888                     fDisplay->CheckEntry(IDM_kStargCaosFilter);
    889                     fDisplay->CheckEntry(IDM_kStarguider);
    890                     ToggleStarguider();
    891 
    892                     fTPoint->MapWindow();
    893                 }
    894                 else
    895                 {
    896                     //uncheck
    897                     //starguider items
    898                     fDisplay->UnCheckEntry(IDM_kStargCaosFilter);
    899                     fDisplay->UnCheckEntry(IDM_kStarguider);
    900                     fTPoint->UnmapWindow();
    901                     ToggleStarguider();
    902 
    903                     //general
    904                     fDisplay->UnCheckEntry(IDM_kCatalog);
    905                     fDisplay->EnableEntry(IDM_kFilter);
    906                     fDisplay->EnableEntry(IDM_kCaosFilter);
    907                     fDisplay->EnableEntry(IDM_kFindStar);
    908                     fChannel->EnableEntry(IDM_kChannel3);
    909                 }
    910                 return kTRUE;
    911 
     564/*
    912565            case IDM_kTpointMode:
    913566                Toggle(fMode, IDM_kTpointMode);
     
    921574
    922575                    //from starguider
    923                     SwitchOff(fDisplay, IDM_kStargCaosFilter);
    924                     SwitchOff(fDisplay, IDM_kCatalog);
    925                     SwitchOff(fDisplay, IDM_kStarguider);
    926                     ToggleStarguider();
    927                     fMode->UnCheckEntry(IDM_kStarguiderMode);
    928                     SwitchOff(fOperations, IDM_kStargAnalysis);
    929                     ToggleStargAnalysis();
     576                    //SwitchOff(fDisplay, IDM_kStargCaosFilter);
     577                    //SwitchOff(fDisplay, IDM_kCatalog);
     578                    //SwitchOff(fDisplay, IDM_kStarguider);
     579                    //ToggleStarguider();
     580                    //fMode->UnCheckEntry(IDM_kStarguiderMode);
     581                    //SwitchOff(fOperations, IDM_kStargAnalysis);
     582                    //ToggleStargAnalysis();
    930583
    931584                    //switch camera
     
    946599                    //enable
    947600                    //starguider items
    948                     fDisplay->EnableEntry(IDM_kStargCaosFilter);
    949                     fDisplay->EnableEntry(IDM_kCatalog);
    950                     fDisplay->EnableEntry(IDM_kStarguider);
    951                     fOperations->EnableEntry(IDM_kStargAnalysis);
     601                    //fDisplay->EnableEntry(IDM_kStargCaosFilter);
     602                    //fDisplay->EnableEntry(IDM_kCatalog);
     603                    //fDisplay->EnableEntry(IDM_kStarguider);
     604                    //fOperations->EnableEntry(IDM_kStargAnalysis);
    952605
    953606                    //general
     
    963616                }
    964617                return kTRUE;
    965 
     618*/
    966619            case IDM_kFilter:
    967620                Toggle(fDisplay, IDM_kFilter);           
     
    979632            case IDM_kCaosFilter:
    980633                Toggle(fDisplay, IDM_kCaosFilter);
    981                 ToggleCaosFilter();
    982                 return kTRUE;
    983 
    984             case IDM_kStargCaosFilter:
    985                 Toggle(fDisplay, IDM_kStargCaosFilter);               
    986                 if (fDisplay->IsEntryChecked(IDM_kStargCaosFilter))
    987                 {
    988                     SwitchOff(fDisplay, IDM_kCaosFilter);
    989                     SwitchOff(fDisplay, IDM_kFindStar);
    990                     SwitchOff(fDisplay, IDM_kFilter);
    991                 }
    992                 else
    993                 {
    994                     fDisplay->EnableEntry(IDM_kFilter);
    995                     fDisplay->EnableEntry(IDM_kFindStar);
    996                     fDisplay->EnableEntry(IDM_kCaosFilter);
    997                 }
    998                 return kTRUE;
    999 
     634                //ToggleCaosFilter();
     635                return kTRUE;
     636/*
    1000637            case IDM_kCaosPrintLeds:
    1001638            case IDM_kCaosPrintRings:
     
    1006643                fCaosAnalyse->DisableEntry(IDM_kCaosAnalStart);
    1007644                fCaosAnalyse->EnableEntry(IDM_kCaosAnalStop);
    1008                 fCaos->InitHistograms();
     645                //fCaos->InitHistograms();
    1009646                return kTRUE;
    1010647
     
    1111748                }
    1112749                return kTRUE;
    1113 
     750             */
    1114751            case IDM_kInterpol250:
    1115752            case IDM_kInterpol125:
     
    1153790                }
    1154791                return kTRUE;
    1155 
    1156             case IDM_kLimMag3:
    1157             case IDM_kLimMag4:
    1158             case IDM_kLimMag5:
    1159             case IDM_kLimMag6:
    1160             case IDM_kLimMag7:
    1161             case IDM_kLimMag8:
    1162             case IDM_kLimMag9:
    1163                 for (int i=IDM_kLimMag3; i<=IDM_kLimMag9; i++)
    1164                     if (mp1==i)
    1165                         fLimMag->CheckEntry(i);
    1166                     else
    1167                         fLimMag->UnCheckEntry(i);
    1168 
    1169                 fSao->SetLimitMag(mp1-IDM_kLimMag3+3);
    1170                 return kTRUE;
    1171792            }
    1172793            break;
     
    1177798    return kTRUE;
    1178799}
    1179 
    1180 void MStarguider::SetPointingPosition(RaDec rd)
    1181 {
    1182     rd.Ra(rd.Ra()*24/360);
    1183     fCRaDec->SetCoordinates(rd);
    1184 }
    1185 
    1186 ZdAz MStarguider::TrackingError(TArrayF &x, TArrayF &y, TArrayF &mag,
    1187                                 Int_t &num) const
    1188 {
    1189     num = -1;
    1190 
    1191     // Width of the final 2D gaussian
    1192     const Double_t width = 0.8;
    1193 
    1194     // The binning is made 1.7 sigma (which, in the 1D case, gives the
    1195     // highest significance of a gaussian signal)  (1bin equiv 3x3 sigma)
    1196     const Double_t max =  50;
    1197     const Int_t    n   = TMath::Nint(2*max/(1.77*1.7*width));
    1198     //1.77=sqrt(pi) comes from pir=bin
    1199 
    1200     // Fill a histogram with the dx/dy values from the file
    1201     TH2F hist("Hist", "dX/dY", n, -max, max, n, -max, max);
    1202 //    hist.SetDirectory(0);
    1203 
    1204     //
    1205     // Search for matching Magnitudes
    1206     //
    1207     //for (int i=0; i<mag.GetSize(); i++)
    1208     //{
    1209     //   if (mag[i]>48-15 && mag[i]<48+15)
    1210     //        h.Fill(x[i], y[i]);
    1211     //}
    1212 
    1213     // fill dx and dy into histogram
    1214     for (int i=0; i<mag.GetSize(); i++)
    1215         hist.Fill(x[i], y[i]);
    1216 
    1217     // Remove all bins which have content of only a single event
    1218     // including under- and overflow bins
    1219     for (int i=0; i<hist.GetSize(); i++)
    1220         if (hist.GetBinContent(i)<1.5)
    1221             hist.SetBinContent(i, 0);
    1222 
    1223     // Find the bin containing the maximum
    1224     Int_t bx, by, bz;
    1225     hist.GetMaximumBin(bx, by, bz);
    1226 
    1227     // In the worst case the events are spread through the
    1228     // bins which are the neighbors of the maximum
    1229     // Normally neighbors which do not belong to the signal are empty!
    1230     hist.GetXaxis()->SetRange(bx-1, bx+1); // 3x3bins  ~8x8pix  ~9x9sigma
    1231     hist.GetYaxis()->SetRange(by-1, by+1); // 3x3bins  ~8x8pix  ~9x9sigma
    1232 
    1233     // Check whether this region contains events
    1234     num = TMath::Nint(hist.Integral());
    1235     if (num<1)
    1236         return ZdAz();
    1237 
    1238     // Get Mean for the given axis range
    1239     Double_t mx = hist.GetMean(1);
    1240     Double_t my = hist.GetMean(2);
    1241 
    1242     const Int_t max1 = TMath::Nint(hist.GetBinContent(bx, by));
    1243 
    1244     // Check if the maximum is unique
    1245     Int_t bx2, by2, bz2;
    1246     hist.GetXaxis()->SetRange(-1, 9999);
    1247     hist.GetYaxis()->SetRange(-1, 9999);
    1248     hist.SetBinContent(bx, by, 0);
    1249     hist.GetMaximumBin(bx2, by2, bz2);
    1250 
    1251     const Int_t max2 = TMath::Nint(hist.GetBinContent(bx2, by2));
    1252     if (max1==max2 && TMath::Abs(bx2-bx)+TMath::Abs(by2-by)>1)
    1253         return ZdAz();
    1254 
    1255     // loop again over the file and fill the histogram again.
    1256     // to get closer to the correct value
    1257     Double_t sumx = 0;
    1258     Double_t sumy = 0;
    1259              num  = 0;
    1260     for (int i=0; i<mag.GetSize(); i++)
    1261     {
    1262         // only fill values into the histogram which
    1263         // are inside 2*1.7=3.4 sigma (99.7%)
    1264         if (TMath::Hypot(x[i]-mx, y[i]-my)>width*3.4)
    1265             continue;
    1266 
    1267         sumx += x[i];
    1268         sumy += y[i];
    1269         num++;
    1270     }
    1271 
    1272     if (num<1)
    1273         return ZdAz();
    1274 
    1275     // calc MeanX and MeanY
    1276     mx = sumx/num;
    1277     my = sumy/num;
    1278 
    1279     // loop again to fill the mispointing corrected histograms
    1280     // and fill another histogram to check the quality of the calculated
    1281     // mispointing.
    1282     sumx = 0;
    1283     sumy = 0;
    1284     num  = 0;
    1285     for (int i=0; i<mag.GetSize(); i++)
    1286     {
    1287         // only fill values into the histogram which
    1288         // are inside 1.7=3.4 sigma (92%)
    1289         // Cut at 1.7 sigma which gives the best background supression
    1290         if (TMath::Hypot(x[i]-mx, y[i]-my)>width*1.7)
    1291             continue;
    1292 
    1293         sumx += x[i];
    1294         sumy += y[i];
    1295         num++;
    1296     }
    1297 
    1298     if (num<1)
    1299         return ZdAz(); // FIXME!!!!!!!!!!!!
    1300 
    1301     // Mispointing
    1302     mx = sumx/num;
    1303     my = sumy/num;
    1304 
    1305 #ifdef EXPERT
    1306     cout << "Offset-XY:    " << mx << " +- " << my << endl;
    1307 #endif
    1308 
    1309     return fSao->CalcDeltaZdAzFromPix(mx, my)*TMath::RadToDeg();
    1310800/*
    1311     AltAz pos0 = fSao->CalcAltAzFromPix(768/2,    576/2)*kRad2Deg;
    1312     AltAz pos1 = fSao->CalcAltAzFromPix(768/2+mx, 576/2+my)*kRad2Deg;
    1313 
    1314     pos1 -= pos0;
    1315 
    1316     ZdAz res2 = ZdAz(-pos1.Alt(), pos1.Az());
    1317 
    1318     return res1;
    1319 */
    1320 }
    1321 
    1322 Int_t MStarguider::CalcTrackingError(const Leds &leds, MStarList &stars,
    1323                                      ZdAz &d, const MTime &t, const double &bright,
    1324                                      Int_t &num, const Ring &center, Int_t numleds, Int_t numrings)
     801void MStarguider::DrawZoomImage(const byte *img)
     802{
     803    byte zimg[kZOOM*kZOOM];
     804    for (int y=0; y<kZOOM; y++)
     805        for (int x=0; x<kZOOM; x++)
     806            zimg[x+y*kZOOM] = img[(fDx+(x-kZOOM/2)/2)+(fDy+(y-kZOOM/2)/2)*768];
     807
     808    fZoomImage->DrawImg(zimg);
     809}
     810
     811void MStarguider::DrawCosyImage(const byte *img)
     812{
     813    if (!fCosy)
     814        return;
     815
     816    byte simg[(768/2-1)*(576/2-1)];
     817    for (int y=0; y<576/2-1; y++)
     818        for (int x=0; x<768/2-1; x++)
     819            simg[x+y*(768/2-1)] = ((unsigned int)img[2*x+2*y*768]+img[2*x+2*y*768+1]+img[2*x+2*(y+1)*768]+img[2*x+2*(y+1)*768+1])/4;
     820
     821    fCosy->GetWin()->GetImage()->DrawImg(simg);
     822}
     823
     824Led MStarguider::FindStar(const FilterLed &f, const FilterLed &f2, const Ring &center, Int_t numleds, Int_t numrings)
    1325825{
    1326826    // Get tracking coordinates
    1327     const XY xy = fCRaDec->GetCoordinates();
    1328 
    1329     num = leds.GetEntries();
    1330     if (num < 3) //was 1
    1331     {
    1332         gLog << warn << "Sorry, less than 3 detected spot in FOV!" << endl;
    1333         if (fTPoint->IsDown() && fCosy && fCosy->GetDriveCom())
    1334             fCosy->GetDriveCom()->SendTPoint(false, 'S', fTPointStarMag, fTPointStarName, AltAz(), ZdAz(), xy, 0, 0, t, center, Led(), numleds, numrings, num);    // Report
    1335         fTPoint->SetDown(kFALSE);
    1336         return 0;
    1337     }
    1338 
    1339     if (stars.GetRealEntries() < 3)
    1340     {
    1341         gLog << warn << "Sorry, less than 3 stars in FOV!" << endl;
    1342         if (fTPoint->IsDown() && fCosy && fCosy->GetDriveCom())
    1343             fCosy->GetDriveCom()->SendTPoint(false, 'S', fTPointStarMag, fTPointStarName, AltAz(), ZdAz(), xy, 0, 0, t, center, Led(), numleds, numrings, num);    // Report
    1344         fTPoint->SetDown(kFALSE);
    1345         return 0;
    1346     }
    1347 
    1348     stars.Sort(); // Sort by magnitude
    1349 
    1350 #ifdef EXPERT
    1351     TString str = "data/tracking_";
    1352     str += fSao->GetMjd();
    1353     str += ".txt";
    1354 
    1355     ofstream fout(str);
    1356 #endif
    1357 
    1358     TArrayF x, y, mag;
    1359 
    1360     // FIXME: Is predefined value 3 a good idea?
    1361     Int_t cnt=0;
    1362 
    1363     MStar *star=0;
    1364     MStarListIter NextStar(&stars);
    1365     while ((star=NextStar()) && cnt<num+5)
    1366     {
    1367         TIter NextSp(&leds);
    1368         Led *spot=NULL;
    1369         while ((spot=(Led*)NextSp()))
    1370         {
    1371             const XY dpos(spot->GetX()-(768-star->GetX()),
    1372                           spot->GetY()-star->GetY());
    1373            
    1374             const Int_t idx = x.GetSize();
    1375 
    1376             x.Set(idx+1);
    1377             y.Set(idx+1);
    1378             mag.Set(idx+1);
    1379 
    1380             x.AddAt(dpos.X(), idx);
    1381             y.AddAt(dpos.Y(), idx);
    1382             mag.AddAt(1, idx);
    1383 //            mag.AddAt(log(spot->GetMag())/star->GetMag(), idx);
    1384 /*
    1385 #ifdef EXPERT
    1386             if (fout)
    1387             {
    1388                 fout << spot->GetX() << " "
    1389                     << spot->GetY()  << " "
    1390                     << -2.5*log(spot->GetMag())  << " "
    1391                     << star->GetX() << " "
    1392                     << star->GetY()  << " "
    1393                     << star->GetMag()  << " ";
    1394                 fout << x[idx] << " " << y[idx] << " " << mag[idx] << endl;
    1395             }
    1396 #endif
    1397 */
    1398         }
    1399         cnt++;
    1400     }
    1401 
    1402     Int_t numcor = 0;
    1403     d = TrackingError(x, y, mag, numcor);
    1404     if (numcor<1)
    1405     {
    1406         fTPoint->SetDown(kFALSE);
    1407         if (fTPoint->IsDown() && fCosy && fCosy->GetDriveCom())
    1408             fCosy->GetDriveCom()->SendTPoint(false, 'S', fTPointStarMag, fTPointStarName, AltAz(), ZdAz(), xy, 0, 0, t, center, Led(), numleds, numrings, num, numcor);    // Report
    1409         return 0;
    1410     }
    1411 
    1412     fDZdAz->SetCoordinates(d);
    1413 
    1414     //
    1415     // Current Pointing position
    1416     //
    1417     ZdAz cpos = fSao->GetZdAz()-d;
    1418     fPZdAz->SetCoordinates(cpos);
    1419 
    1420     // Check StargTPoint data set request
    1421     if (!fMode->IsEntryChecked(IDM_kStarguiderMode) || !fTPoint->IsDown())
    1422         return numcor;
    1423     fTPoint->SetDown(kFALSE);
    1424 
    1425     // If no file open: open new file
    1426     if (!fOutStargTp)
    1427     {
    1428         // open tpoint file
    1429         const TString name = MCosy::GetFileName("tpoint", "starg", "txt");
    1430         gLog << all << "TPoint-Starg File ********* " << name << " ********** " << endl;
    1431 
    1432         fOutStargTp = new ofstream(name);
    1433         *fOutStargTp << "Magic Model  STARGTPOINT data file" << endl;
    1434         *fOutStargTp << ": ALTAZ" << endl;
    1435         *fOutStargTp << "49 48 0 ";
    1436         *fOutStargTp << t << endl;
    1437     }
    1438 
    1439     const RaDec rd(xy.X()*MAstro::HorToRad(), xy.Y()*TMath::DegToRad());
    1440 
    1441     // From the star position in the camera we calculate the Alt/Az
    1442     // position we are currently tracking (real pointing position)
    1443     fSao->SetMjd(t.GetMjd());
    1444     AltAz za0 = fSao->CalcAltAz(rd)*kRad2Deg;
    1445 
    1446     //correction with offset from last tpoint measurement
    1447     za0 -= fAltAzOffsetFromTp;
    1448 
    1449     //if the difference between the tpoint and the starguider tpoint
    1450     //is too big, the starguider tpoint is not stored
    1451     if ((t.GetMjd()-fTimeFromTp.GetMjd())>0.001) //1min20sec
    1452     {
    1453         cout << "     time difference between tpoint and starguider-tpoint > 1m20s" << endl;
    1454         cout << "     => starguider tpoint hasn't been stored. " << endl;
    1455         cout << "        Please repeat whole procedure. " << endl;
    1456         return numcor;
    1457     }
    1458 
    1459 
    1460     // Write real pointing position
    1461     *fOutStargTp << setprecision(7) << za0.Az() << " " << za0.Alt() << " ";
    1462 
    1463     // Write system pointing position
    1464     *fOutStargTp << fmod(cpos.Az()+360, 360) << " " << 90-cpos.Zd();
    1465 
    1466     *fOutStargTp << " " << xy.X() << " " << xy.Y();
    1467     *fOutStargTp << " " << d.Zd() << " " << d.Az();
    1468     *fOutStargTp << " " << setprecision(11) << t.GetMjd();
    1469     *fOutStargTp << " " << setprecision(4) << center.GetMag();
    1470     *fOutStargTp << " 0";// << star->GetMag();
    1471     *fOutStargTp << " " << center.GetX() << " " << center.GetY();
    1472     *fOutStargTp << " 0 0";
    1473     *fOutStargTp << " " << numleds << " " << numrings;
    1474     *fOutStargTp << " " << num << " " << numcor << " " << bright;
    1475     *fOutStargTp << " " << fTPointStarMag << " " << fTPointStarName;
    1476     *fOutStargTp << endl;
    1477 
    1478     gLog << all << "Starguider TPoint successfully taken." << endl;
    1479 
    1480     fTimeFromTp.Clear();
    1481 
    1482     if (!fCosy)
    1483         return numcor;
    1484 
    1485     MDriveCom *com = fCosy->GetDriveCom();
    1486     if (!com)
    1487         return numcor;
    1488 
    1489     com->SendTPoint(true, 'S', fTPointStarMag, fTPointStarName, za0, cpos, xy, d.Zd(), d.Az(), t, center, Led(), numleds, numrings, num, numcor);    // Report
    1490 
    1491     return numcor;
    1492 }
    1493 
    1494 void MStarguider::FindStar(const FilterLed &f, const FilterLed &f2, const Ring &center, const MTime &t, Int_t numleds, Int_t numrings)
    1495 {
    1496     // Get tracking coordinates
    1497     const XY xy = fCRaDec->GetCoordinates();  // [h, deg]
     827    //const XY xy = fCRaDec->GetCoordinates();  // [h, deg]
    1498828
    1499829    if (center.GetX()<=0 && center.GetY()<=0)
    1500830    {
    1501         gLog << warn << "Couldn't determine center of the camera." << endl;
    1502         if (fTPoint->IsDown() && fCosy && fCosy->GetDriveCom())
    1503             fCosy->GetDriveCom()->SendTPoint(false, 'T', fTPointStarMag, fTPointStarName, AltAz(), ZdAz(), xy, 0, 0, t, center, Led(), numleds, numrings);    // Report
     831        cout << "Couldn't determine center of the camera." << endl;
     832        //if (fTPoint->IsDown() && fCosy && fCosy->GetDriveCom())
     833        //    fCosy->GetDriveCom()->SendTPoint(false, 'T', fTPointStarMag, fTPointStarName, AltAz(), ZdAz(), xy, 0, 0, t, center, Led(), numleds, numrings);    // Report
    1504834        return;
    1505835    }
     
    1510840
    1511841    // Check whether star found
    1512     Led *star = (Led*)leds.At(0);
    1513     if (!star)
    1514     {
    1515         gLog << warn << "No star found." << endl;
    1516         if (fTPoint->IsDown() && fCosy && fCosy->GetDriveCom())
    1517             fCosy->GetDriveCom()->SendTPoint(false, 'T', fTPointStarMag, fTPointStarName, AltAz(), ZdAz(), xy, 0, 0, t, center, Led(), numleds, numrings);    // Report
     842    if (leds.size()==0)
     843    {
     844        cout << "No star found." << endl;
     845        //if (fTPoint->IsDown() && fCosy && fCosy->GetDriveCom())
     846        //    fCosy->GetDriveCom()->SendTPoint(false, 'T', fTPointStarMag, fTPointStarName, AltAz(), ZdAz(), xy, 0, 0, t, center, Led(), numleds, numrings);    // Report
    1518847        return;
    1519848    }
    1520849
    1521850    cout << "Found star @ " << flush;
    1522     star->Print();
    1523     f2.MarkPoint(star->GetX(), star->GetY(), 2<<2);
     851
     852    const Led &star = leds.front();
     853
     854    star.Print();
     855    f2.MarkPoint(star.GetX(), star.GetY(), 2<<2);
     856
     857    return star;
    1524858
    1525859    const RaDec rd(xy.X()*MAstro::HorToRad(), xy.Y()*TMath::DegToRad());
     
    1530864    ac.SetGeom(geom);
    1531865    ac.SetRadiusFOV(3);
    1532     ac.SetObservatory(*fSao);
     866    //ac.SetObservatory(*fSao);
    1533867    ac.SetTime(t);
    1534868    ac.SetRaDec(rd.Ra(), rd.Dec());
     
    1536870    // Convert from Pixel to millimeter (1pix=2.6mm) [deg/pix / deg/mm = mm/pix]
    1537871    // Correct for abberation.
    1538     const Double_t conv = fCaos->GetArcsecPerPixel()/3600/geom.GetConvMm2Deg()/ 1.0713;
     872    const double sec_per_pix = 45.311;
     873    const double mm_to_deg   = 0.011693;
     874    const Double_t conv = sec_per_pix/3600/mm_to_deg / 1.0713;
    1539875
    1540876    // Adapt coordinate system (GUIs and humans are counting Y in different directions)
     
    1558894        // open tpoint file
    1559895        //
     896
    1560897        const TString name = MCosy::GetFileName("tpoint", "tpoint", "txt");
    1561898        cout << "TPoint File ********* " << name << " ********** " << endl;
     
    1566903        *fOutTp << "49 48 0 ";
    1567904        *fOutTp << t << endl;
    1568         // temp(°C) pressure(mB) height(m) humidity(1) wavelength(microm) troplapserate(K/m)
     905        // temp(°C) pressure(mB) height(m) humidity(1) wavelength(microm) troplapserate(K/m
    1569906    }
    1570907
     
    1608945
    1609946    gLog << all << "TPoint successfully taken." << endl;
    1610 /*
     947
    1611948    MLog &outrep = *fCosy->GetOutRep();
    1612949    if (outrep.Lock("MStarguider::FindStar"))
     
    1622959        outrep << setprecision(11) << t.GetMjd() << endl;
    1623960        outrep.UnLock("MStarguider::FindStar");
    1624     }*/
     961    }
    1625962   
    1626963//    return zdaz;
     
    1635972    com->SendTPoint(true, 'T', fTPointStarMag, fTPointStarName, za0, za1, xy, -dzd, -daz, t, center, *star, numleds, numrings);    // Report
    1636973}
     974*/
     975
     976void MStarguider::WritePNG(const char *name, const byte *gbuf, const byte *cbuf)
     977{
     978    string fname(name);
     979    if (fname.size()<4 || fname.compare(fname.size()-4,4, ".png")==0)
     980        fname += ".png";
     981
     982    cout << "Writing PNG '" << fname << "'" << endl;
     983
     984#ifdef HAVE_PNG
     985    //
     986    // open file
     987    //
     988    FILE *fd = fopen(fname.c_str(), "w");
     989    if (!fd)
     990    {
     991        cout << "Warning: Cannot open file for writing." << endl;
     992        return;
     993    }
     994
     995    //
     996    // allocate memory
     997    //
     998    png_structp fPng = png_create_write_struct(PNG_LIBPNG_VER_STRING,
     999                                               NULL, NULL, NULL);
     1000
     1001    if (!fPng)
     1002    {
     1003        cout << "Warning: Unable to create PNG structure" << endl;
     1004        fclose(fd);
     1005        return;
     1006    }
     1007
     1008
     1009    png_infop fInfo = png_create_info_struct(fPng);
     1010
     1011    if (!fInfo)
     1012    {
     1013        cout << "Warning: Unable to create PNG info structure" << endl;
     1014        png_destroy_write_struct (&fPng, NULL);
     1015        fclose(fd);
     1016        return;
     1017    }
     1018
     1019    fInfo->width      = 768;
     1020    fInfo->height     = 576;
     1021    fInfo->bit_depth  = 8;
     1022    fInfo->color_type = PNG_COLOR_TYPE_RGB;
     1023//    fInfo->color_type = PNG_COLOR_TYPE_GRAY;
     1024
     1025    //
     1026    // set jump-back point in case of errors
     1027    //
     1028    if (setjmp(fPng->jmpbuf))
     1029    {
     1030        cout << "longjmp Warning: PNG encounterd an error!" << endl;
     1031        png_destroy_write_struct (&fPng, &fInfo);
     1032        fclose(fd);
     1033        return;
     1034    }
     1035
     1036    //
     1037    // connect file to PNG-Structure
     1038    //
     1039    png_init_io(fPng, fd);
     1040
     1041    // png_set_compression_level (fPng, Z_BEST_COMPRESSION);
     1042
     1043    //
     1044    // Write header
     1045    //
     1046    png_write_info(fPng, fInfo);
     1047
     1048    png_byte buf[768*576*3];
     1049
     1050    png_byte *d = buf;
     1051    const byte *g = gbuf;
     1052    const byte *c = cbuf;
     1053
     1054    // d=destination, s1=source1, s2=source2, e=end
     1055    while (d<buf+768*576*3)
     1056    {
     1057        if (fScreenshotColor && *c)
     1058        {
     1059            *d++ = ((*c>>4)&0x3)*85;
     1060            *d++ = ((*c>>2)&0x3)*85;
     1061            *d++ = ((*c++ )&0x3)*85;
     1062            g++;
     1063        }
     1064        else
     1065        {
     1066            *d++ = *g;
     1067            *d++ = *g;
     1068            *d++ = *g++;
     1069            c++;
     1070        }
     1071    }
     1072
     1073    //
     1074    // Write bitmap data
     1075    //
     1076    for (unsigned int y=0; y<768*576*3; y+=768*3)
     1077        png_write_row (fPng, buf+y);
     1078
     1079    //
     1080    // Write footer
     1081    //
     1082    png_write_end (fPng, fInfo);
     1083
     1084    //
     1085    // free memory
     1086    //
     1087    png_destroy_write_struct (&fPng, &fInfo);
     1088
     1089    fclose(fd);
     1090#else
     1091    cout << "Sorry, no png support compiled into tpoint." << endl;
     1092#endif
     1093}
    16371094
    16381095bool MStarguider::Interpolate(const unsigned long n, byte *img) const
     
    16411098        return true;
    16421099
     1100    const int rate = fTPointMode>0 ? 5*25 : fIntRate;
     1101
    16431102    static unsigned short myimg[768*576];
    16441103
     
    16471106    const byte *end = img+768*576;
    16481107
    1649     if (n%fIntRate)
     1108    if (n%rate)
    16501109    {
    16511110        while (img<end)
     
    16571116        while (img<end)
    16581117        {
    1659             *img = (*img + *f)/fIntRate;
     1118            *img = (*img + *f)/rate;
    16601119            ++img;
    16611120            *f++ = 0;
     
    16661125}
    16671126
    1668 void MStarguider::DrawZoomImage(const byte *img)
    1669 {
    1670     byte zimg[kZOOM*kZOOM];
    1671     for (int y=0; y<kZOOM; y++)
    1672         for (int x=0; x<kZOOM; x++)
    1673             zimg[x+y*kZOOM] = img[(fDx+(x-kZOOM/2)/2)+(fDy+(y-kZOOM/2)/2)*768];
    1674 
    1675     fZoomImage->DrawImg(zimg);
    1676 }
    1677 
    1678 void MStarguider::DrawCosyImage(const byte *img)
    1679 {
    1680     if (!fCosy)
    1681         return;
    1682 
    1683     byte simg[(768/2-1)*(576/2-1)];
    1684     for (int y=0; y<576/2-1; y++)
    1685         for (int x=0; x<768/2-1; x++)
    1686             simg[x+y*(768/2-1)] = ((unsigned int)img[2*x+2*y*768]+img[2*x+2*y*768+1]+img[2*x+2*(y+1)*768]+img[2*x+2*(y+1)*768+1])/4;
    1687 
    1688     fCosy->GetWin()->GetImage()->DrawImg(simg);
    1689 }
    1690 
    1691 void MStarguider::StartTPoint(char *cmd)
    1692 {
    1693     fTPointStarName = "Dummy";
    1694     fTPointStarMag  = 0;
    1695 
    1696     TString str(cmd);
    1697 
    1698     TObjArray *arr = str.Tokenize(' ');
    1699 
    1700     if (arr->GetEntries()==2)
    1701     {
    1702         fTPointStarName = (*arr)[0]->GetName();
    1703         fTPointStarMag  = atof((*arr)[1]->GetName());
    1704     }
    1705 
    1706     delete arr;
    1707 
    1708     if (fTPointFromCC<0)
    1709         fTPointFromCC=0;
    1710 }
    1711 
    1712 void MStarguider::StartStarguider(bool on)
    1713 {
    1714     // Switch to starguider mode
    1715     if (on)
    1716     {
    1717         cout << " * Switching Starguider mode on" << endl;
    1718         fMode->UnCheckEntry(IDM_kStarguiderMode);
    1719         ProcessMessage(MK_MSG(kC_COMMAND, kCM_MENU), IDM_kStarguiderMode, 0);
    1720     }
    1721     else
    1722     {
    1723         cout << " * Switching Starguider mode off" << endl;
    1724         fMode->CheckEntry(IDM_kStarguiderMode);
    1725         ProcessMessage(MK_MSG(kC_COMMAND, kCM_MENU), IDM_kStarguiderMode, 0);
    1726     }
    1727 }
    1728 
    1729 Int_t MStarguider::GetStarguiderMode() const
    1730 {
    1731     return fMode->IsEntryChecked(IDM_kStarguiderMode) ? 1 : 2;
    1732 
    1733 }
    1734 
    1735 Bool_t MStarguider::DoTPoint()
    1736 {
    1737     if (fTPoint->IsDown() && fTPointFromCC<0)
    1738     {
    1739         fTPointFromCC = 0;
    1740         fTPoint->SetDown(kFALSE);
    1741     }
    1742 
    1743     if (fTPointFromCC<0)
    1744         return kTRUE;
    1745 
    1746     switch (++fTPointFromCC)
    1747     {
    1748     case 1:
    1749         fTimeFromTp.Clear();
    1750         fNumStarsDetected = 0;
    1751 
    1752         if (!fMode->IsEntryChecked(IDM_kTpointMode))
    1753         {
    1754             cout << " * Switching to TPoint mode" << endl;
    1755             // Switch to tpoint mode
    1756             fMode->UnCheckEntry(IDM_kTpointMode);
    1757             ProcessMessage(MK_MSG(kC_COMMAND, kCM_MENU), IDM_kTpointMode, 0);
    1758             return kFALSE;
    1759         }
    1760         fTPointFromCC++;
    1761 
    1762     case 2:
    1763     case 3:
    1764         cout << " * Waiting one frame" << endl;
    1765         // Wait one frame
    1766         return kFALSE;
    1767 
    1768     case 4:
    1769         cout << " * Taking TPoint" << endl;
    1770         fTPoint->SetDown(); // kTRUE
    1771         return kTRUE;
    1772 
    1773     case 5:
    1774         if (!fTimeFromTp) // TPoint failed
    1775             break;
    1776 #ifdef FACT
    1777         fTPointFromCC = -1;
    1778         return kFALSE;
    1779 #else
    1780         // Switch to starguider mode
    1781         cout << " * Switching to Starguider mode" << endl;
    1782         fMode->UnCheckEntry(IDM_kStarguiderMode);
    1783         ProcessMessage(MK_MSG(kC_COMMAND, kCM_MENU), IDM_kStarguiderMode, 0);
    1784         return kFALSE;
    1785 #endif
    1786     case 6:
    1787     case 7:
    1788         cout << " * Waiting one frame" << endl;
    1789         // Wait one frame
    1790         return kFALSE;
    1791 
    1792     case 8:
    1793         cout << " * Taking Starguider TPoint" << endl;
    1794         fTPoint->SetDown(); // kTRUE
    1795         return kTRUE;
    1796 
    1797 //    case 9:
    1798 //        cout << " * Send Report" << endl;
    1799 //        return kTRUE;
    1800 
    1801     case 9:
    1802         cout << " * Switching to TPoint mode" << endl;
    1803         // Switch to tpoint mode
    1804         fMode->UnCheckEntry(IDM_kTpointMode);
    1805         ProcessMessage(MK_MSG(kC_COMMAND, kCM_MENU), IDM_kTpointMode, 0);
    1806         fTPointFromCC = -1;
    1807         return kFALSE;
    1808     }
    1809 
    1810     // Send report
    1811 /*
    1812     if (!fCosy)
    1813         return kTRUE;
    1814 
    1815     MDriveCom *com = fCosy->GetDriveCom();
    1816     if (!com)
    1817         return kTRUE;
    1818 
    1819     // nominalaz, nominalel, realaz, realel, nomra, nomdec,
    1820     // diffaz, diffel, mjd, numleds, artmag
    1821 
    1822     //fTimeFromTp.Clear();
    1823     //fStatus==MDriveCom::kMonitoring
    1824     //fNumStarsDetected   = numstars;
    1825     //fNumStarsCorrelated = rc;
    1826     com->SendTPoint(fNumStarsCorrelated>0);
    1827 */
    1828     return kTRUE;
    1829 }
    1830 
    1831 void MStarguider::ProcessFrame(const unsigned long n, byte *img,
     1127void MStarguider::ProcessFrame(const unsigned long n, byte *img,
    18321128                               struct timeval *tm)
    18331129{
    1834     static unsigned long n0 = n;
    1835 
    1836     MTime t(*tm);
    1837 
    1838     const Double_t d = t-fTime;
    1839     if (d>1)
    1840     {
    1841         const TString txt = MString::Format("%dfps", (int)((n-n0)/d+.5));
    1842         fFps->SetText(txt);
    1843         fTime = t;
    1844         n0 = n;
    1845     }
    1846 
    18471130    if (!Interpolate(n, img))
    18481131        return;
     1132
     1133    if (fTPointMode==2)
     1134    {
     1135        fTPointMode=1;
     1136        return;
     1137    }
    18491138
    18501139    byte cimg[768*576];
    18511140    memset(cimg, 0, 768*576);
    18521141
    1853     FilterLed f(img,   768, 576, 2.5); // 2.5
    1854     FilterLed f2(cimg, 768, 576); // former color 0xb0
    1855 
    1856     if (fDisplay->IsEntryChecked(IDM_kStretch))
     1142    FilterLed f (img,  768, 576, 2.5); // 2.5
     1143    FilterLed f2(cimg, 768, 576);      // former color 0xb0
     1144
     1145    if (!fTPointMode && fScreenshotName.empty() && fDisplay->IsEntryChecked(IDM_kStretch))
    18571146        f.Stretch();
    18581147
    1859     if (!fWritePictures->IsEntryEnabled(IDM_kStart) &&
    1860         (!(n%fWrtRate) || fWriteType->IsEntryChecked(IDM_kOnce)))
    1861     {
    1862 
    1863         if (fFileType->IsEntryChecked(IDM_kPNG))
    1864             Writer::Png("pix/file", img, tm, fCRaDec->GetCoordinates());
    1865 
    1866         if (fFileType->IsEntryChecked(IDM_kPPM))
    1867             Writer::Ppm("pix/file", img, tm, fCRaDec->GetCoordinates());
    1868        
    1869         if (fWriteType->IsEntryChecked(IDM_kOnce))
    1870             ProcessMessage(MK_MSG(kC_COMMAND, kCM_MENU), IDM_kStop, 0);
    1871     }
    1872 
    18731148    // Visual Filter, whole FOV
    1874     if (fDisplay->IsEntryChecked(IDM_kFilter))
    1875         f.Execute();
    1876 
    1877     if (!DoTPoint())
    1878         return;
    1879 
    1880     Int_t numleds = 0;
    1881     Int_t numrings = 0;
     1149    if (!fTPointMode && fDisplay->IsEntryChecked(IDM_kFilter))
     1150    {
     1151        vector<Led> leds;
     1152        f.Execute(leds, 768/2, 576/2);
     1153        for (auto it=leds.begin(); it!=leds.end(); it++)
     1154            f.MarkPoint(*it);
     1155    }
    18821156
    18831157    // Find Center of Camera for Caos and Tpoints
    1884 //    Ring center;//(-1, -1);
    1885     Ring center(5, 5);//(-1, -1);
    1886     if (fDisplay->IsEntryChecked(IDM_kCaosFilter))
    1887     {
    1888         const bool printl = fCaosPrint->IsEntryChecked(IDM_kCaosPrintLeds);
    1889         const bool printr = fCaosPrint->IsEntryChecked(IDM_kCaosPrintRings);
    1890         ZdAz pos;
    1891         if (fCosy)
    1892             pos = fCosy->GetPointingPos();
    1893 
    1894         center = fCaos->Run(img, printl, printr, pos, t);
    1895 
     1158    int numleds  = 0;
     1159    int numrings = 0;
     1160    Ring center(-1, -1);//(5, 5);
     1161
     1162    if (fTPointMode || fDisplay->IsEntryChecked(IDM_kCaosFilter))
     1163    {
     1164        center   = fCaos->Run(img);
    18961165        numleds  = fCaos->GetNumDetectedLEDs();
    18971166        numrings = fCaos->GetNumDetectedRings();
    18981167    }
    18991168
     1169    //cout << "cx=" << center.GetX() << "   cy=" << center.GetY() << "   Nled=" << numleds << "   Nrings=" << numrings << endl;
     1170
    19001171    // Find Star at Center---for Tpoint Procedure
    1901     if (fDisplay->IsEntryChecked(IDM_kFindStar))
    1902     {
    1903         // Set search Paremeters (FIXME: Get them from user input!)
    1904         f.SetCut(fFindStarCut);
    1905         f.SetBox(fFindStarBox);
    1906 
    1907         FindStar(f, f2, center, t, numleds, numrings);
    1908     }
    1909 
    1910     DrawZoomImage(img);
    1911     DrawCosyImage(img);
    1912 
    1913     // Position corresponding to the camera center (53.2, 293.6)
    1914     Ring sgcenter;//(-1, -1); // Center of camera in SG picture [px]
    1915 
    1916     // Find Center of Camera in Starfield Camera picture
    1917     if (fDisplay->IsEntryChecked(IDM_kStargCaosFilter))
    1918     {
    1919         ZdAz pos;
    1920         if (fCosy)
    1921             pos = fCosy->GetPointingPos();
    1922 
    1923         sgcenter = fStargCaos->Run(img, kFALSE, kFALSE, pos, t); // [px]
    1924 
    1925         numleds  = fStargCaos->GetNumDetectedLEDs();
    1926         numrings = fStargCaos->GetNumDetectedRings();
    1927     }
    1928 
    1929     // we obtain a list of stars in the FOV from the SAO catalog
    1930     if (fDisplay->IsEntryChecked(IDM_kCatalog))
    1931     {
    1932         MTime time(*tm);
    1933 
    1934         XY xy = fCRaDec->GetCoordinates();  //[h,  deg]
    1935         fRaDec->Set(xy.X()*15, xy.Y());     //[deg,deg]
    1936 
    1937         UpdatePosZoom();
    1938 
    1939         // Always call SetMjd first!
    1940         fSao->SetPointing(time.GetMjd(), *fRaDec);
    1941         fCZdAz->SetCoordinates(fSao->GetZdAz());
    1942 
    1943         fSao->SetBox(fStarguiderW, fStarguiderH); // 280 Region of interest around center
    1944 
    1945         // If center of camera cannot be determined sgcenter is (0,0)
    1946         const Bool_t centerisvalid = sgcenter.GetX()>0 && sgcenter.GetY()>0;
    1947 
    1948         // We determine the ideal starfield using camera sagging info
    1949         // from the LEDs
    1950         //Please never change this offsets!!!
    1951         // 53.2 and 293.6 are the "preliminary" camera center
    1952         // -9 and 28.5 are the offsets of the pointing position in the sky
    1953         //const XY off(sgcenter.GetX()- 53.2-9,
    1954         //             sgcenter.GetY()-293.6+28.5);
    1955 
    1956         // Skycenter
    1957         const XY off(sgcenter.GetX() + fSkyOffsetX,
    1958                      sgcenter.GetY() + fSkyOffsetY);
    1959 
    1960         MStarList stars;
    1961         if (centerisvalid)
     1172    Led star(-1, -1);
     1173    if (center.GetX()>0 && center.GetY()>0)
     1174    {
     1175        if (fTPointMode || fDisplay->IsEntryChecked(IDM_kFindStar))
    19621176        {
    1963             // we obtain stars in the effective star FOV and draw them.
    1964             // coordinates are video frame coords.
    1965             fSao->CalcStars(stars, fStarguiderX, fStarguiderY, off.X(), off.Y());
    1966             fSao->DrawStars(stars, cimg);
     1177            // Set search Paremeters (FIXME: Get them from user input!)
     1178            f.SetCut(fFindStarCut+1.5);  // FindStar.CleaningLevel
     1179            f.SetBox(fFindStarBox+12/*+80*/);  // FindStar.SizeBox
     1180
     1181            // Try to find the star
     1182            vector<Led> leds;
     1183            f.FindStar(leds, (Int_t)center.GetX(), (Int_t)center.GetY(), true);
     1184
     1185            // Check whether star found
     1186            if (leds.size()>0)
     1187            {
     1188                //cout << "Found star @ " << flush;
     1189                //leds[0].Print();
     1190                f2.MarkPoint(leds[0].GetX(), leds[0].GetY(), 2<<2);
     1191                star = leds[0];
     1192            }
    19671193        }
    19681194
    1969         // There are two corrections to the mispointing
    1970         // - Sagging of the camera as measured with the LEDs
    1971         // - Star not ideally centered on MAGIC Camera
    1972 
    1973         // Next we evaluate the offset given by the LEDs. This we obtain
    1974         // in Zd/Az and add it to the tracking error.
    1975         if (fDisplay->IsEntryChecked(IDM_kStarguider))
     1195        // DrawZoomImage(img);
     1196        // DrawCosyImage(img);
     1197
     1198        // Position corresponding to the camera center (53.2, 293.6)
     1199        // Draw Circles around center of Camera
     1200        if (fTPointMode || fDisplay->IsEntryChecked(IDM_kCaosFilter))
    19761201        {
    1977             const Float_t cut = atof(fCut->GetText());
    1978 
    1979             Leds spots;
    1980             f.SetBox(fStarguiderW, fStarguiderH); // 280
    1981             f.SetCut(cut);
    1982 
    1983             double bright;
    1984             f.ExecuteAndMark(spots, fStarguiderX, fStarguiderY, bright);
    1985 
    1986             ULong_t color;
    1987             gClient->GetColorByName("Green", color);
    1988             if (bright> 60)
    1989                 gClient->GetColorByName("Yellow", color);
    1990             if (bright> 85)
    1991                 gClient->GetColorByName("Orange", color);
    1992             if (bright> 95)
    1993                 gClient->GetColorByName("Red", color);
    1994             fSkyBright->SetBackgroundColor(color);
    1995 
    1996             const MString txt = MString::Format("Sky Brightness: %.1f", bright);
    1997             fSkyBright->SetText(txt);
    1998 
    1999             const Bool_t brightnessisvalid = bright<1.75*fLastBright &&
    2000                 bright>30 && bright<110;
    2001 
    2002             fLastBright = bright;
    2003 
    2004             Int_t numstars = 0;
    2005             const Int_t rc = CalcTrackingError(spots, stars, fD, t, bright, numstars, sgcenter, numleds, numrings);
    2006 
    2007             const Bool_t monitoring = brightnessisvalid && centerisvalid && fNumStarsCorrelated>2;
    2008 
    2009             fStatus = monitoring ? MDriveCom::kMonitoring : MDriveCom::kError;
    2010 
    2011             if (fCosy)
    2012                 fPos = fCosy->GetPointingPos();
    2013 
    2014             if (fOperations->IsEntryChecked(IDM_kStargAnalysis))
    2015                 fStargHistograms->Fill(spots, stars, fD, fSao->GetZdAz(),
    2016                                        sgcenter, bright, fPos, t);
    2017 
    2018             fNumStarsDetected   = numstars;
    2019             fNumStarsCorrelated = rc;
    2020 
    2021             if (fCosy)
    2022             {
    2023                 MDriveCom *com = fCosy->GetDriveCom();
    2024                 if (com)
    2025                     com->SendStargReport(fStatus, fD, fSao->GetZdAz(),
    2026                                          sgcenter, numstars, rc, bright,
    2027                                          time.GetMjd(), numleds, numrings);    // Report
    2028             }
    2029 
    2030         } //kStarguider
    2031 
    2032         if (centerisvalid && fNumStarsCorrelated>2)
    2033         {
    2034             // Position around which the circles are drawn.
    2035             const Ring skycenter(off.X(), off.Y());
    2036 
    2037             f2.DrawCircle(skycenter,   2.0,     0x0a);
    2038 
    2039             const Double_t ap = fSao->GetPixSize()/3600; //[deg/pix]
    2040 
    2041             f2.DrawCircle(skycenter,  0.10/ap, 0x0a); //0.1deg
    2042             f2.DrawHexagon(skycenter, 2.06/ap, 0x0a);
    2043             f2.DrawHexagon(skycenter, 3.50/ap, 0x0a);
     1202            f2.DrawCircle(center, 0x0a);
     1203            f2.DrawCircle(center,   7.0,
     1204                          fDisplay->IsEntryChecked(IDM_kFindStar)?3:0xb0);
     1205            //f2.DrawCircle(center, 115.0, 0x0a);
     1206            //f2.DrawCircle(center, 230.0, 0x0a);
     1207            //f2.DrawCircle(center, 245.0, 0x0a);
    20441208        }
    2045 
    2046     } //CalcStars
    2047 
    2048     if (fTPoint->IsDown())
    2049         fTPoint->SetDown(kFALSE);
    2050 
    2051     // Draw Circles around center of Camera
    2052     if (fDisplay->IsEntryChecked(IDM_kCaosFilter) &&
    2053         center.GetX()>0 && center.GetY()>0)
    2054     {
    2055         f2.DrawCircle(center, 0x0a);
    2056         f2.DrawCircle(center,   7.0,
    2057                       fDisplay->IsEntryChecked(IDM_kFindStar)?3:0xb0);
    2058         f2.DrawCircle(center, 115.0, 0x0a);
    2059         f2.DrawCircle(center, 230.0, 0x0a);
    2060         f2.DrawCircle(center, 245.0, 0x0a);
    2061     }
    2062 
    2063     if (fDisplay->IsEntryChecked(IDM_kStargCaosFilter) &&
    2064         sgcenter.GetX()>0 && sgcenter.GetY()>0)
    2065     {
    2066         f2.DrawCircle(sgcenter, 0x0a);
    2067         f2.DrawCircle(sgcenter,   5.0,
    2068                       fDisplay->IsEntryChecked(IDM_kFindStar)?3:0xb0);
    2069     }
    2070 
    2071     if (fDisplay->IsEntryChecked(IDM_kCaosFilter) ||
    2072         fDisplay->IsEntryChecked(IDM_kCatalog)    ||
     1209    }
     1210
     1211    if (fTPointMode ||
     1212        fDisplay->IsEntryChecked(IDM_kCaosFilter) ||
    20731213        fDisplay->IsEntryChecked(IDM_kFindStar))
    20741214        fImage->DrawColImg(img, cimg);
    20751215    else
    20761216        fImage->DrawImg(img);
    2077 }
    2078 
     1217
     1218    if (fTPointMode && !fScreenshotName.empty())
     1219    {
     1220        WritePNG(fScreenshotName.c_str(), img, cimg);
     1221        fScreenshotName = "";
     1222        fTPointMode = 0;
     1223        return;
     1224    }
     1225
     1226
     1227    if (star.GetX()<0 || star.GetY()<0 || fTPointMode==0)
     1228        return;
     1229
     1230    if (fTPointMode==1)
     1231        fTPointMode=0;
     1232
     1233    // Convert from Pixel to millimeter (1pix=2.6mm) [deg/pix / deg/mm = mm/pix]
     1234    // Correct for abberation.
     1235
     1236    //const float Dleds =     510; // 5.96344 deg
     1237    //const float Dpix  = 2*237.58;
     1238
     1239    // The DC reflector elongates light from off axis sources
     1240    // This is a correction. It is 7% for MAGIC (1:1) and
     1241    // less for FACT (1:1.4). This is an estimate from a
     1242    // Orbit mode observation at 0.17deg distance to the
     1243    // camera center [4.90m might also not be very accurate
     1244    // depending on the position of the CCD camera]
     1245    const double abberation  = 1.0638; //1.0713;
     1246    const double sec_per_pix = 45.14;  //45.311;  (atan(510mm/2 / 4.90m) / 237.58)           // FACT LEDs
     1247
     1248    const double conv = sec_per_pix/abberation;
     1249
     1250    const double dx = star.GetX()-center.GetX();
     1251    const double dy = star.GetY()-center.GetY();
     1252
     1253    //const double dxx = - conv * (star.GetX()-center.GetX());
     1254    //const double dyy =   conv * (star.GetY()-center.GetY());
     1255
     1256    const double dphi = - center.GetPhi() * M_PI/180;
     1257
     1258    // The sign is because the pixels are not counted in
     1259    // both directions in the same direction as Zd/Az
     1260    const double dxx = - conv * (dx*cos(dphi) - dy*sin(dphi));
     1261    const double dyy =   conv * (dx*sin(dphi) + dy*cos(dphi));
     1262
     1263    double arr[11] = {
     1264        dxx, dyy,
     1265        double(numleds), double(numrings),
     1266        center.GetX(), center.GetY(), center.GetMag(),
     1267        star.GetX(),   star.GetY(),   star.GetMag(),
     1268        center.GetPhi()
     1269    };
     1270
     1271    fDimData.setData(arr, 11*sizeof(double));
     1272    fDimData.setQuality(0);
     1273    fDimData.setTimestamp(tm->tv_sec, tm->tv_usec/1000);
     1274    fDimData.updateService();
     1275}
     1276
     1277/*
    20791278void MStarguider::UpdatePosZoom()
    20801279{
    2081     MString txt;/*
     1280    MString txt;
    20821281    if (fDisplay->IsEntryChecked(IDM_kCatalog))
    20831282    {
     
    20891288        txt.Form("(%d, %d) %.1fd/%.1fd", fDx, fDy, -aa.Alt(), aa.Az()-180);
    20901289    }
    2091     else*/
     1290    else
    20921291        txt.Form("(%d, %d)", fDx, fDy);
    20931292    fPosZoom->SetText(txt);
     
    21151314    fDy = dy;
    21161315
    2117     UpdatePosZoom();
     1316    //UpdatePosZoom();
    21181317    return kTRUE;
    21191318}
    21201319
    2121 void MStarguider::Print(TString &str, Double_t deg) const
    2122 {
    2123     Char_t sgn;
    2124     UShort_t d, m, s;
    2125 
    2126     MAstro::Deg2Dms(deg, sgn, d, m, s);
    2127 
    2128     str += MString::Format("%c %03d %02d %03d ", sgn, d, m, s);
    2129 }
     1320*/
     1321
     1322/// Overwritten DimCommand::commandHandler
     1323void MStarguider::commandHandler()
     1324{
     1325    DimCommand *cmd = getCommand();
     1326    if (!cmd)
     1327        return;
     1328
     1329    if (cmd==&fDimTPoint)
     1330    {
     1331        fTPointMode = 2;
     1332        fScreenshotName = "";
     1333        cout << "DimCommand[TPOINT]: " << cmd->itsSize << " " << string((char*)cmd->itsData, cmd->itsSize) << endl;
     1334        return;
     1335    }
     1336
     1337    if (cmd==&fDimScreenshot && fTPointMode==0)
     1338    {
     1339        if (cmd->itsSize<2)
     1340            return;
     1341
     1342        fTPointMode = 2;
     1343        fScreenshotColor = ((uint8_t*)cmd->itsData)[0];
     1344        fScreenshotName = string((char*)cmd->itsData+1, cmd->itsSize-1);
     1345        cout << "DimCommand[SCREENSHOT]: " << fScreenshotName << endl;
     1346        return;
     1347    }
     1348
     1349    cout << "DimCommand[UNKNOWN]: " << cmd->itsSize << " " << string((char*)cmd->itsData, cmd->itsSize) << endl;
     1350}
  • trunk/FACT++/drive/MStarguider.h

    r18618 r18627  
    44#include "PixClient.h"
    55#include "MGImage.h"
     6#include "Led.h"
     7#include "Camera.h"
    68
    7 #include "MPointing.h"
     9#include "dic.hxx"
     10#include "dis.hxx"
    811
    9 #ifndef MARS_MObservatory
    10 #include "MObservatory.h"
    11 #endif
    12 #ifndef MARS_MTime
    13 #include "MTime.h"
    14 #endif
    15 #ifndef MARS_MGList
    16 #include "MGList.h"
    17 #endif
    18 
    19 class TArrayF;
    20 class TGLabel;
    21 class TGButton;
    22 class TString;
    23 //class TSocket;
    24 
    25 class TTimer;
     12#include <TGButton.h>
    2613
    2714class TGMenuBar;
     
    3219
    3320class MGImage;
    34 class MGCoordinates;
    35 
    36 class StarCatalog;
    37 class MStarList;
    38 class MCosy;
    3921class MCaos;
    40 // class MStargLeds;
    41 class MStargHistograms;
    42 class MGNumStars;
    43 class MGStarg;
    4422class FilterLed;
    4523class Ring;
    4624
    47 class Leds;
    48 #include <TGButton.h>
    49 class MStarguider : public PixClient, public TGMainFrame
     25//class Leds;
     26class MStarguider : public PixClient, public TGMainFrame, public DimCommandHandler
    5027{
    5128private:
    52     PixGetter     *fGetter;
     29    DimService fDimData;
     30    DimCommand fDimTPoint;
     31    DimCommand fDimScreenshot;
    5332
    54     MGList        *fList;
     33    Camera        *fGetter;
    5534
    5635    TGMenuBar     *fMenu;
    5736    MGImage       *fImage;
    58     MGImage       *fZoomImage;
    59 
    60     MGStarg       *fGStarg;
    61     MGNumStars    *fGNumStars;
    6237
    6338    TGPopupMenu   *fDisplay;
    64     TGPopupMenu   *fMode;
    65     TGPopupMenu   *fWritePictures;
    6639    TGPopupMenu   *fSetup;
    67     TGPopupMenu   *fOperations;
    68 
    69     TGPopupMenu   *fFileType;
    70     TGPopupMenu   *fWriteType;
    71     TGPopupMenu   *fWriteRate;
    7240    TGPopupMenu   *fInterpol;
    73     TGPopupMenu   *fLimMag;
    74     TGPopupMenu   *fChannel;
    7541
    7642    TGPopupMenu   *fCaosWrite;
     
    7945    TGPopupMenu   *fCaOs;
    8046
    81     MGCoordinates *fCRaDec;
    82     MGCoordinates *fCZdAz;
    83 
    84     MGCoordinates *fPZdAz;
    85     MGCoordinates *fDZdAz;
    86     MGCoordinates *fSZdAz;
    87 
    88     TGTextEntry   *fPixSize;
    89     TGTextEntry   *fAngle;
    90     TGTextEntry   *fCut;
    91 
    92     TGLabel       *fFps;
    93     TGLabel       *fPosZoom;
    94     TGLabel       *fSkyBright;
    95     TGLabel       *fCZdAzText;
    96     TGLabel       *fPZdAzText;
    97     TGLabel       *fDZdAzText;
    98 public:
    99     TGButton      *fTPoint;
    100 //    TGButton      *fStargTPoint;
    10147private:
    102     ZdAz          fPos;
    103     ZdAz          fD;
    104     AltAz         fAltAzOffsetFromTp;
    105 
    106     Int_t          fNumStarsDetected;
    107     Int_t          fNumStarsCorrelated;
    108 
    109     StarCatalog   *fSao;
    110     MCosy         *fCosy;
    11148    MCaos         *fCaos;
    112     MCaos         *fStargCaos;
    113     MStargHistograms *fStargHistograms;
    114    
    115     RaDec *fRaDec;
    116 
    117     TTimer *fTimer;
    118     MTime   fTime;
    119     MTime   fTimeFromTp;
    120 
    121     std::ofstream *fOutTp;
    122     std::ofstream *fOutStargTp;
    123 //    ofstream *fOutRq;
     49    TTimer        *fTimer;
    12450
    12551    Int_t fDx;
     
    12753
    12854    byte fIntRate;
    129     int  fWrtRate;
     55//    int  fWrtRate;
    13056
    131     UInt_t fStatus;
    13257    Double_t fLastBright;
    133 
    13458    Double_t fRadius; // LED radius [cm]
    135 
    136     Int_t fTPointFromCC;
    137 
    138     Int_t fStarguiderW;
    139     Int_t fStarguiderH;
    140     Int_t fStarguiderX;
    141     Int_t fStarguiderY;
    142 
    143     Float_t fSkyOffsetX;  // Offset between camera center and sky position
    144     Float_t fSkyOffsetY;  // Offset between camera center and sky position
    14559
    14660    Float_t fFindStarCut;
    14761    Int_t   fFindStarBox;
    14862
    149     TString fTPointStarName;
    150     Float_t fTPointStarMag;
     63    int fTPointMode;
     64
     65    bool fScreenshotColor;
     66    std::string fScreenshotName;
     67
     68    void WritePNG(const char *name, const byte *gbuf, const byte *cbuf);
    15169
    15270    void Toggle(TGPopupMenu *p, UInt_t id);
    15371    void SwitchOff(TGPopupMenu *p, UInt_t id);
    154     void ToggleStargAnalysis();
    155 //    void ToggleFindStar();
    156     void ToggleStarguider();
    157     void ToggleCaosFilter();
    158     void SetChannel();
    159     Int_t CalcTrackingError(const Leds &, MStarList &, ZdAz &, const MTime &, const double &bright, Int_t &num, const Ring &center, Int_t numleds, Int_t numrings);
    160     ZdAz TrackingError(TArrayF &alt, TArrayF &az, TArrayF &mag, Int_t &num) const;
    16172    bool Interpolate(const unsigned long n, byte *img) const;
    162     void FindStar(const FilterLed &f, const FilterLed &f2, const Ring &center, const MTime &t, Int_t numleds, Int_t numrings);
    16373
    16474    void InitGui(Int_t channel);
    16575
    166     void UpdatePosZoom();
    167     void DrawZoomImage(const byte *img);
    168     void DrawCosyImage(const byte *img);
     76    //void DrawZoomImage(const byte *img);
     77    //void DrawCosyImage(const byte *img);
    16978
    17079    Bool_t HandleTimer(TTimer *t);
    17180
    172     void SetRotationAngle(Double_t angle);
    173     void SetPixSize(Double_t size);
    17481    void SetCut(Double_t cut);
    17582
    176     Bool_t DoTPoint();
    177 
    17883public:
    179     MStarguider(MObservatory::LocationName_t obs, Int_t channel);
     84    MStarguider(Int_t channel);
    18085    virtual ~MStarguider();
    18186
    182     void SetupEnv(TEnv &env);
     87    //void SetupEnv(TEnv &env);
    18388
    18489    void Layout();
     
    18792    Bool_t ProcessMessage(Long_t msg, Long_t parm1, Long_t parm2);
    18893
    189     Bool_t HandleDoubleClick(Event_t *event);
    190 
    191     void SetPointingPosition(RaDec rd);
    192     void SetCosy(MCosy *c) { fCosy = c; }
    193 
    194     XY GetCoordinates() const;
     94    //Bool_t HandleDoubleClick(Event_t *event);
    19595
    19696    //
     
    19999    void ProcessFrame(const unsigned long n, byte *img, struct timeval *tm);
    200100
    201     void Print(TString &str, Double_t deg) const;
    202     void Print(Option_t *o) const { }
    203 
    204     void StartTPoint(char *cmd=0);// { if (fTPointFromCC<0) fTPointFromCC=0; }
    205     void StartStarguider(bool on);
    206 
    207     Int_t GetStarguiderMode() const;
    208 
    209     ClassDef(MStarguider, 0)
     101    void commandHandler();       /// Overwritten DimCommand::commandHandler
    210102};
    211103
Note: See TracChangeset for help on using the changeset viewer.