Ignore:
Timestamp:
05/30/05 10:31:23 (19 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Cosy/tpoint/gui.C

    r4617 r7104  
    1 #include <fstream.h>
    2 #include <fstream.h>
    3 #include <fstream.h>
    4 #include <iostream.h>
    5 #include <iomanip.h>
     1#include <fstream>
     2#include <iostream>
     3#include <iomanip>
    64
    75#include <TError.h>
     
    2624#include <TCanvas.h>
    2725
    28 #include "coord.h"
     26//#include "coord.h"
    2927
    3028#include "MGList.h"
    31 #include "MBending.h"
     29#include "MPointing.h"
     30
     31using namespace std;
    3232
    3333class Set : public TObject
     
    126126    void AdjustAz(Double_t daz) { fStarAz += daz*TMath::DegToRad(); }
    127127
    128     void Adjust(const MBending &bend)
     128    void Adjust(const MPointing &bend)
    129129    {
    130130        AltAz p = bend(GetStarAltAz());
     
    132132        fStarAz = p.Az();
    133133    }
    134     void AdjustBack(const MBending &bend)
     134    void AdjustBack(const MPointing &bend)
    135135    {
    136136        AltAz p = bend.CorrectBack(GetRawAltAz());
     
    182182private:
    183183    enum {
    184         kTbFit = 19, //MBending::GetNumPar(), // FIXME!!!
     184        kTbFit = 19, //MPointing::GetNumPar(), // FIXME!!!
    185185        kTbLoad,
    186186        kTbSave,
     
    196196    TList fLabel;
    197197
    198     MBending fBending;
     198    MPointing fBending;
    199199
    200200    FontStruct_t fFont;
     
    204204        f = 0;
    205205
    206         MBending bend;
    207         bend.SetParameters(par); // Set Parameters [deg] to MBending
     206        MPointing bend;
     207        bend.SetParameters(par); // Set Parameters [deg] to MPointing
    208208
    209209        for (int i=0; i<fCoordinates.GetSize(); i++)
     
    430430        grmag.SetTitle(" \\Delta vs. Mag ");
    431431   
    432         TMinuit minuit(MBending::GetNumPar());  //initialize TMinuit with a maximum of 5 params
     432        TMinuit minuit(MPointing::GetNumPar());  //initialize TMinuit with a maximum of 5 params
    433433        minuit.SetObjectFit(this);
    434434        minuit.SetPrintLevel(-1);
    435435        minuit.SetFCN(fcn);
    436436   
    437         fBending.SetMinuitParameters(minuit, MBending::GetNumPar()); // Init Parameters [deg]
    438    
    439         for (int i=0; i<MBending::GetNumPar(); i++)
     437        fBending.SetMinuitParameters(minuit, MPointing::GetNumPar()); // Init Parameters [deg]
     438   
     439        for (int i=0; i<MPointing::GetNumPar(); i++)
    440440        {
    441441            TGButton *l = (TGButton*)fList->FindWidget(i);
     
    482482        TArrayD par;
    483483        fBending.GetParameters(par);
    484         for (int i=2; i<MBending::GetNumPar(); i++)
     484        for (int i=2; i<MPointing::GetNumPar(); i++)
    485485            par[i]=0;
    486486   
    487         MBending b2;
     487        MPointing b2;
    488488        b2.SetParameters(par);
    489489   
     
    495495            Set &set0 = *(Set*)fCoordinates.At(i);
    496496   
    497             ZdAz za = set0.GetStarZdAz()*kRad2Deg;
     497            ZdAz za(set0.GetStarZdAz());
     498            za *=kRad2Deg;
    498499   
    499500            //
     
    615616
    616617        c1=new TCanvas("CanvGraphs", "Graphs");
    617         c1->Divide(3,3,0,0);
     618        c1->Divide(3,3,1e-10,1e-10);
    618619
    619620        TLine line;
     
    755756
    756757        c1=new TCanvas("CanvResiduals", "Residuals", 800, 800);
    757         c1->Divide(2, 2, 0, 0);
     758        c1->Divide(2, 2, 1e-10, 1e-10);
    758759   
    759760        c1->cd(2);
     
    952953        TGLabel *l;
    953954
    954         for (int i=0; i<MBending::GetNumPar(); i++)
     955        for (int i=0; i<MPointing::GetNumPar(); i++)
    955956        {
    956957            l = (TGLabel*)fLabel.At(i);
    957958            l->SetText(Form("%.4f\xb0", par[i]));
    958959
    959             l = (TGLabel*)fLabel.At(MBending::GetNumPar()+i);
     960            l = (TGLabel*)fLabel.At(MPointing::GetNumPar()+i);
    960961            l->SetText(Form("\xb1 %8.4f\xb0", err[i]));
    961962        }
     
    964965    void DisplayData()
    965966    {
    966         TGLabel *l = (TGLabel*)fLabel.At(3*MBending::GetNumPar());
     967        TGLabel *l = (TGLabel*)fLabel.At(3*MPointing::GetNumPar());
    967968        l->SetText(Form("%d data sets loaded.", fOriginal.GetSize()));
    968969    }
     
    970971    void DisplayResult(Double_t before, Double_t after, Double_t backw)
    971972    {
    972         TGLabel *l1 = (TGLabel*)fLabel.At(3*MBending::GetNumPar()+1);
     973        TGLabel *l1 = (TGLabel*)fLabel.At(3*MPointing::GetNumPar()+1);
    973974        l1->SetText(Form("Before: %.1f +- %.1f SE", before, 0));
    974975
    975         TGLabel *l2 = (TGLabel*)fLabel.At(3*MBending::GetNumPar()+2);
     976        TGLabel *l2 = (TGLabel*)fLabel.At(3*MPointing::GetNumPar()+2);
    976977        l2->SetText(Form("After:  %.1f +- %.1f SE", after, 0));
    977978
    978         TGLabel *l3 = (TGLabel*)fLabel.At(3*MBending::GetNumPar()+3);
     979        TGLabel *l3 = (TGLabel*)fLabel.At(3*MPointing::GetNumPar()+3);
    979980        l3->SetText(Form("Backw:  %.1f +- %.1f SE", backw, 0));
    980981    }
     
    10381039        fList->Add(vframe);
    10391040
    1040         for (int i=0; i<MBending::GetNumPar(); i++)
    1041             AddCheckButton(vframe, fBending.GetName(i), i);
     1041        for (int i=0; i<MPointing::GetNumPar(); i++)
     1042            AddCheckButton(vframe, fBending.GetVarName(i), i);
    10421043
    10431044        vframe = new TGVerticalFrame(comp, 1, 1);
     
    10571058        vframe->AddFrame(l,h);
    10581059
    1059         for (int i=1; i<MBending::GetNumPar(); i++)
     1060        for (int i=1; i<MPointing::GetNumPar(); i++)
    10601061            AddLabel(vframe, "+000.0000", h)->SetTextJustify(kTextRight);
    10611062
     
    10641065        fList->Add(vframe);
    10651066
    1066         for (int i=0; i<MBending::GetNumPar(); i++)
     1067        for (int i=0; i<MPointing::GetNumPar(); i++)
    10671068            AddLabel(vframe, "\xb1 00.0000\xb0", h)->SetTextJustify(kTextRight);
    10681069
     
    10711072        fList->Add(vframe);
    10721073
    1073         for (int i=0; i<MBending::GetNumPar(); i++)
     1074        for (int i=0; i<MPointing::GetNumPar(); i++)
    10741075            AddLabel(vframe, fBending.GetDescription(i), h);
    10751076
Note: See TracChangeset for help on using the changeset viewer.