Changeset 4107 for trunk/MagicSoft


Ignore:
Timestamp:
05/20/04 05:40:59 (20 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Cosy
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Cosy/Makefile

    r4105 r4107  
    145145        @ln -sf ../../Mars/mbase/MGMenu.h  base/MGMenu.h
    146146        @ln -sf ../../Mars/mbase/MGMenu.cc base/MGMenu.cc
     147        @ln -sf ../../Mars/mbase/MString.h  base/MString.h
     148        @ln -sf ../../Mars/mbase/MString.cc base/MString.cc
    147149        @ln -sf ../../Mars/mgeom/MGeomMirror.h  base/MGeomMirror.h
    148150        @ln -sf ../../Mars/mgeom/MGeomMirror.cc  base/MGeomMirror.cc
    … …  
    217219        @cp stars.txt /home/stesy/Cosy/
    218220        @chmod a-w /home/stesy/Cosy/stars.txt
    219         @echo Copy SAO catalog
    220         @chmod u+w /home/stesy/Cosy/sao/sao-*
    221         @cp sao/sao-* /home/stesy/Cosy/sao/
    222         @chmod a-w /home/stesy/Cosy/sao/sao-*
     221        @echo Copy PPM catalog
     222        @chmod u+w /home/stesy/Cosy/ppm9.bin
     223        @cp ppm9.bin /home/stesy/Cosy/
     224        @chmod a-w /home/stesy/Cosy/ppm9.bin
    223225        @echo finished...
    224226

  • trunk/MagicSoft/Cosy/caos/Led.cc

    r2278 r4107  
    77
    88#include "Ring.h"
     9#include "MString.h"
    910
    1011ClassImp(Led);
    … …  
    1718void Led::Print(Option_t *o=NULL) const
    1819{
     20    MString str;
    1921    cout << "Led: ";
    20     cout << "x="   << Form("%5.1f", fX)   << "+-" << fDx   << ", ";
    21     cout << "y="   << Form("%5.1f", fY)   << "+-" << fDy   << ", ";
    22     cout << "phi=" << Form("%6.1f", fPhi) << "+-" << fDphi << ", ";
     22    cout << "x="   << str.Print("%5.1f", fX)   << "+-" << fDx   << ", ";
     23    cout << "y="   << str.Print("%5.1f", fY)   << "+-" << fDy   << ", ";
     24    cout << "phi=" << str.Print("%6.1f", fPhi) << "+-" << fDphi << ", ";
    2325    cout << "mag=" << fMag << endl;
    2426}

  • trunk/MagicSoft/Cosy/caos/Makefile

    r1799 r4107  
    2020# @endcode
    2121
    22 INCLUDES = -I. -I..
     22INCLUDES = -I. -I.. -I../base
     23#base: MString
    2324
    2425# @code

  • trunk/MagicSoft/Cosy/caos/Ring.cc

    r2278 r4107  
    77
    88#include "Rings.h"
     9#include "MString.h"
    910
    1011ClassImp(Ring);
    1112
    12 Ring::Ring() :
    13     fX(0), fY(0), fR(0), fPhi(0), fDx(-1), fDy(-1), fDr(-1), fDphi(-1)
     13Ring::Ring(Double_t x, Double_t y) :
     14    fX(x), fY(y), fR(0), fPhi(0), fDx(-1), fDy(-1), fDr(-1), fDphi(-1)
    1415{
    1516}
    … …  
    114115void Ring::Print(Option_t *o=NULL) const
    115116{
     117    MString str;
    116118    cout << "Ring: ";
    117     cout << "x="   << Form("%5.1f", fX) << "+-" << Form("%.1f", fDx) << ", ";
    118     cout << "y="   << Form("%5.1f", fY) << "+-" << Form("%.1f", fDy) << ", ";
    119     cout << "r="   << Form("%5.1f", fR) << "+-" << Form("%.1f", fDr) << ", ";
     119    cout << "x="   << str.Print("%5.1f", fX) << "+-" << Form("%.1f", fDx) << ", ";
     120    cout << "y="   << str.Print("%5.1f", fY) << "+-" << Form("%.1f", fDy) << ", ";
     121    cout << "r="   << str.Print("%5.1f", fR) << "+-" << Form("%.1f", fDr) << ", ";
    120122    cout << "phi=" << fPhi              << "+-" << fDphi << endl;
    121123}

  • trunk/MagicSoft/Cosy/caos/Ring.h

    r2278 r4107  
    3232
    3333public:
    34     Ring();
     34    Ring(Double_t x=0, Double_t y=0);
    3535
    3636    Double_t GetX() const   { return fX; }

  • trunk/MagicSoft/Cosy/main/MCosy.cc

    r4105 r4107  
    883883    cout << " *  Y: " << kGear.Y() << "*" << kResRE.Y()/4 << "/" << kResSE.Y() << "=4*" << kGearTot.Y() << "/" << kResSE.Y() << endl;
    884884}
    885 
     885/*
    886886void MCosy::InitSync()
    887887{
    … …  
    899899    fMac3->StartPosSync();
    900900}
    901 
     901*/
    902902void MCosy::TalkThreadSeTest()
    903903{
    … …  
    10421042    }
    10431043
    1044     InitSync();
     1044    //InitSync();
    10451045
    10461046    /*** FOR DEMO MODE ***/

  • trunk/MagicSoft/Cosy/main/MTracking.cc

    r4105 r4107  
    550550        //        time.Zd(fCosy->fMac2->GetMjd());
    551551
    552         Double_t mjd1 = fCosy->fZd1->GetMjd();
    553         Double_t mjd2 = fCosy->fZd2->GetMjd();
    554         Double_t mjd0 = fCosy->fAz->GetMjd();
     552        //Double_t mjd1 = fCosy->fZd1->GetMjd();
     553        //Double_t mjd2 = fCosy->fZd2->GetMjd();
     554        //Double_t mjd0 = fCosy->fAz->GetMjd();
     555
     556        Double_t mjdaz = fCosy->fMac1->GetPdoMjd();//mjd0;
     557        Double_t mjdzd = fCosy->fMac2->GetPdoMjd();//TMath::Max(mjd1, mjd2);
    555558
    556559        // get current position of shaftencoders (interpolated
    … …  
    565568        if (phca1 || phca2)
    566569        {
    567             ZdAz dummy = sla.CalcZdAz(fCosy->fRaDec, fCosy->fMac2->GetPdoMjd());//TMath::Max(mjd1, mjd2));
     570            ZdAz dummy = sla.CalcZdAz(fCosy->fRaDec, mjdzd);
    568571            dummy = fCosy->AlignTrackingPos(dummy);
    569572            dummy = fCosy->fBending(dummy);
    … …  
    572575        if (phcaz)
    573576        {
    574             const Double_t mjd = fCosy->fMac1->GetMjd();
    575             ZdAz dummy = sla.CalcZdAz(fCosy->fRaDec, fCosy->fMac1->GetPdoMjd());//mjd0);
     577            ZdAz dummy = sla.CalcZdAz(fCosy->fRaDec, mjdaz);
    576578            dummy = fCosy->AlignTrackingPos(dummy);
    577579            dummy = fCosy->fBending(dummy);

  • trunk/MagicSoft/Cosy/tpoint/gui.C

    r2568 r4107  
     1#include <fstream.h>
     2#include <fstream.h>
    13#include <fstream.h>
    24#include <iostream.h>
    35#include <iomanip.h>
     6
     7#include <TError.h>
    48
    59#include <TGFrame.h>
    … …  
    3842public:
    3943    Set(Double_t sel=0, Double_t saz=0, Double_t rel=0, Double_t raz=0) :
    40         fStarAz(saz*TMath::Pi()/180),
    41         fStarEl(sel*TMath::Pi()/180),
    42         fRawAz(raz*TMath::Pi()/180),
    43         fRawEl(rel*TMath::Pi()/180)
    44     {
    45     }
    46 
    47     Double_t GetResidual() const
    48     {
    49          Double_t del = fRawEl-fStarEl;
    50          Double_t daz = fRawAz-fStarAz;
    51          Double_t dphi2 = daz/2.;
    52          Double_t cos2  = cos(dphi2)*cos(dphi2);
    53          Double_t sin2  = sin(dphi2)*sin(dphi2);
    54          Double_t d = cos(del)*cos2 - cos(fRawEl+fStarEl)*sin2;
    55 
    56          Double_t dist = acos(d);
    57 
    58          return dist * 180 / TMath::Pi();
     44        fStarAz(saz*TMath::DegToRad()),
     45        fStarEl(sel*TMath::DegToRad()),
     46        fRawAz(raz*TMath::DegToRad()),
     47        fRawEl(rel*TMath::DegToRad())
     48    {
     49    }
     50
     51    Double_t GetResidual(Double_t *err=0) const
     52    {
     53      /*
     54       TVector3 v1, v2;
     55       v1.SetMagThetaPhi(1, TMath::Pi()/2-fRawEl, fRawAz);
     56       v2.SetMagThetaPhi(1, TMath::Pi()/2-fStarEl, fStarAz);
     57
     58       return v1.Angle(v2)*TMath::RadToDeg();
     59       */
     60
     61        const Double_t del = fRawEl-fStarEl;
     62        const Double_t daz = fRawAz-fStarAz;
     63        /*
     64        const Double_t dphi2 = daz/2.;
     65        const Double_t cos2  = cos(dphi2)*cos(dphi2);
     66        const Double_t sin2  = sin(dphi2)*sin(dphi2);
     67        const Double_t d = cos(del)*cos2 - cos(fRawEl+fStarEl)*sin2;
     68        */
     69
     70        const Double_t d  = cos(del) - cos(fRawEl)*cos(fStarEl)*(1.-cos(daz));
     71
     72        if (err)
     73        {
     74            // Error of one pixel in the CCD
     75            const Double_t e1 = 32./3600*TMath::DegToRad()   * 0.5;
     76
     77            // Error of one SE unit
     78            const Double_t e2 = 360./16384*TMath::DegToRad() * 0.5;
     79
     80            const Double_t e11 =  sin(del)+cos(fRawEl)*sin(fStarEl)*(1-cos(daz));
     81            const Double_t e12 =  cos(fRawEl)*cos(fStarEl)*sin(daz);
     82
     83            const Double_t e21 = -sin(del)+sin(fRawEl)*cos(fStarEl)*(1-cos(daz));
     84            const Double_t e22 = -cos(fRawEl)*cos(fStarEl)*sin(daz);
     85
     86            const Double_t err1  = sqrt(1-d*d);
     87            const Double_t err2  = (e11*e11 + e12*e12)*e1*e1;
     88            const Double_t err3  = (e21*e21 + e22*e22)*e2*e2;
     89
     90            *err = sqrt(err2+err3)/err1 * TMath::RadToDeg();
     91        }
     92
     93        const Double_t dist = acos(d);
     94        return dist * TMath::RadToDeg();
    5995    }
    6096
    … …  
    67103    }
    68104
    69     Double_t GetDEl() const     { return (fRawEl-fStarEl)*180/TMath::Pi(); }
     105    Double_t GetDEl() const     { return (fRawEl-fStarEl)*TMath::RadToDeg(); }
    70106    Double_t GetDZd() const     { return -GetDEl(); }
    71     Double_t GetDAz() const     { return (fRawAz-fStarAz)*180/TMath::Pi(); }
    72     Double_t GetStarEl() const  { return fStarEl*180/TMath::Pi(); }
    73     Double_t GetStarZd() const  { return 90.-fStarEl*180/TMath::Pi(); }
    74     Double_t GetStarAz() const  { return fStarAz*180/TMath::Pi(); }
    75     Double_t GetRawEl() const   { return fRawEl*180/TMath::Pi(); }
    76     Double_t GetRawAz() const   { return fRawAz*180/TMath::Pi(); }
    77     Double_t GetRawZd() const   { return 90.-fRawEl*180/TMath::Pi(); }
     107    Double_t GetDAz() const     { return (fRawAz-fStarAz)*TMath::RadToDeg(); }
     108    Double_t GetStarEl() const  { return fStarEl*TMath::RadToDeg(); }
     109    Double_t GetStarZd() const  { return 90.-fStarEl*TMath::RadToDeg(); }
     110    Double_t GetStarAz() const  { return fStarAz*TMath::RadToDeg(); }
     111    Double_t GetRawEl() const   { return fRawEl*TMath::RadToDeg(); }
     112    Double_t GetRawAz() const   { return fRawAz*TMath::RadToDeg(); }
     113    Double_t GetRawZd() const   { return 90.-fRawEl*TMath::RadToDeg(); }
    78114
    79115    ZdAz  GetStarZdAz() const   { return ZdAz(TMath::Pi()/2-fStarEl, fStarAz); }
    … …  
    83119    AltAz GetRawAltAz() const   { return AltAz(fRawEl, fRawAz); }
    84120
    85     void AdjustEl(Double_t del) { fStarEl += del*TMath::Pi()/180; }
    86     void AdjustAz(Double_t daz) { fStarAz += daz*TMath::Pi()/180; }
     121    void AdjustEl(Double_t del) { fStarEl += del*TMath::DegToRad(); }
     122    void AdjustAz(Double_t daz) { fStarAz += daz*TMath::DegToRad(); }
    87123
    88124    void Adjust(const MBending &bend)
    … …  
    105141ifstream &operator>>(ifstream &fin, Set &set)
    106142{
    107     Double_t v[4];
    108     fin >> v[0];
    109     fin >> v[1];
    110     fin >> v[2];
    111     fin >> v[3];
    112 
    113     Double_t dummy;
    114     fin>>dummy;
    115     fin>>dummy;
    116     fin>>dummy;
    117 
    118     set.fStarAz = v[0]*TMath::Pi()/180;
    119     set.fStarEl = v[1]*TMath::Pi()/180;
    120 
    121     set.fRawAz  = v[2]*TMath::Pi()/180;
    122     set.fRawEl  = v[3]*TMath::Pi()/180;
     143    TString str;
     144    str.ReadLine(fin);
     145
     146    Float_t v[4];
     147    sscanf(str.Data(), "%f %f %f %f", v, v+1, v+2, v+3);
     148
     149    set.fStarAz = v[0]*TMath::DegToRad();
     150    set.fStarEl = v[1]*TMath::DegToRad();
     151
     152    set.fRawAz  = v[2]*TMath::DegToRad();
     153    set.fRawEl  = v[3]*TMath::DegToRad();
     154
     155    if (fin)
     156    {
     157        Double_t res, err;
     158        res = set.GetResidual(&err);
     159        cout << "Read: " << v[0] << " " << v[1] << "  :  " << v[2] << " " << v[3] << "  :  " << v[2]-v[0] << " " << v[3]-v[1] << "  :  " << res << " " << err << " " << err/res << endl;
     160    }
    123161
    124162    return fin;
    … …  
    160198            set.Adjust(bend);
    161199
    162             Double_t err = 0.02; // [deg] = 1SE
    163             Double_t res = set.GetResidual()/err;
     200            Double_t err;// = 0.005; // [deg] = 0.25SE
     201            Double_t res = set.GetResidual(&err);
     202            res /= err;
    164203
    165204            f += res*res;
    … …  
    168207        //f /= (fCoordinates.GetSize()-7)*(fCoordinates.GetSize()-7);
    169208        //f /= fCoordinates.GetSize()*fCoordinates.GetSize();
     209        //cout << f << ": " << fCoordinates.GetSize() << endl;
    170210        f /= fCoordinates.GetSize();
    171211    }
    … …  
    187227
    188228        TMarker mark0;
    189         TMarker mark1;
     229        //TMarker mark1;
    190230        mark0.SetMarkerStyle(kStar);
    191         mark1.SetMarkerStyle(kStar);
    192         mark1.SetMarkerColor(kRed);
     231        //mark1.SetMarkerStyle(kStar);
     232        //mark1.SetMarkerColor(kRed);
    193233   
    194234        r0 /= 90;
    195         r1 /= 90;
    196         phi0 *= TMath::Pi()/180;
    197         phi1 *= TMath::Pi()/180;
    198    
     235        //r1 /= 90;
     236        phi0 *= TMath::DegToRad();
     237        //phi1 *= TMath::DegToRad();
     238
    199239        Double_t x0[3] = { r0*cos(phi0), r0*sin(phi0), 0};
    200         Double_t x1[3] = { r1*cos(phi1), r1*sin(phi1), 0};
    201    
     240        //Double_t x1[3] = { r1*cos(phi1), r1*sin(phi1), 0};
     241
    202242        Double_t y0[3], y1[3];
    203243   
    204244        view->WCtoNDC(x0, y0);
    205         view->WCtoNDC(x1, y1);
     245        //view->WCtoNDC(x1, y1);
    206246   
    207247        mark0.DrawMarker(y0[0], y0[1]);
    208         mark1.DrawMarker(y1[0], y1[1]);
     248        //mark1.DrawMarker(y1[0], y1[1]);
    209249    }
    210250   
    … …  
    256296        Double_t dp = p1-p0;
    257297   
    258         Double_t x0[3] = { r0*cos(p0*TMath::Pi()/180), r0*sin(p0*TMath::Pi()/180), 0};
     298        Double_t x0[3] = { r0*cos(p0*TMath::DegToRad()), r0*sin(p0*TMath::DegToRad()), 0};
    259299   
    260300        for (double i=p0+10; i<p1+10; i+=10)
    … …  
    264304   
    265305            Double_t r = dr/dp*(i-p0)+r0;
    266             Double_t p = TMath::Pi()/180*i;
     306            Double_t p = TMath::DegToRad()*i;
    267307   
    268308            Double_t x1[3] = { r*cos(p), r*sin(p), 0};
    … …  
    280320    }
    281321   
    282     void DrawSet(TVirtualPad *pad, Set &set, Float_t scale=-1)
     322    void DrawSet(TVirtualPad *pad, Set &set, Float_t scale=-1, Float_t angle=0)
    283323    {
    284324        Double_t r0   = set.GetRawZd();
    285         Double_t phi0 = set.GetRawAz();
     325        Double_t phi0 = set.GetRawAz()-angle;
    286326        Double_t r1   = set.GetStarZd();
    287         Double_t phi1 = set.GetStarAz();
     327        Double_t phi1 = set.GetStarAz()-angle;
    288328   
    289329        if (r0<0)
    … …  
    323363    }
    324364
    325     void Fit()
     365    void Fit(Double_t &before, Double_t &after, Double_t &backw)
    326366    {
    327367        if (fOriginal.GetSize()==0)
    … …  
    336376
    337377        cout << "-----------------------------------------------------------------------" << endl;
    338    
    339         TH1F hres1("Res1", " Residuals before correction ", fOriginal.GetSize()/2, 0, 3);
    340         TH1F hres2("Res2", " Residuals after correction ",  fOriginal.GetSize(), 0, 3);
    341         TH1F hres3("Res3", " Residuals after backward correction ",  fOriginal.GetSize(), 0, 3);
     378
     379        gStyle->SetOptStat("emro");
     380
     381        TH1F hres1("Res1", " Residuals before correction ", fOriginal.GetSize()/3, 0, 0.3);
     382        TH1F hres2("Res2", " Residuals after correction ",  fOriginal.GetSize()/3, 0, 0.3);
     383        TH1F hres3("Res3", " Residuals after backward correction ",  fOriginal.GetSize()/3, 0, 0.3);
    342384   
    343385        hres1.SetXTitle("\\Delta [\\circ]");
    … …  
    505547        cout << endl;
    506548   
    507 
    508 
    509549
    510550        TCanvas *c1;
    … …  
    576616        cout << "before: " << hres1.GetMean() << " \xb1 " << hres1.GetRMS() << " deg " << endl;
    577617        cout << "after:  " << hres2.GetMean() << " \xb1 " << hres2.GetRMS() << " deg " << endl;
    578         cout << "before: " << (int)(hres1.GetMean()*60+.5) << " \xb1 " << (int)(hres1.GetRMS()*60+.5) << " arcmin" << endl;
    579         cout << "after:  " << (int)(hres2.GetMean()*60+.5) << " \xb1 " << (int)(hres2.GetRMS()*60+.5) << " arcmin" << endl;
    580         cout << "before: " << (int)(hres1.GetMean()*60*60/23.4+.5) << " \xb1 " << (int)(hres1.GetRMS()*60*60/23.4+.5) << " pix" << endl;
    581         cout << "after:  " << (int)(hres2.GetMean()*60*60/23.4+.5) << " \xb1 " << (int)(hres2.GetRMS()*60*60/23.4+.5) << " pix" << endl;
    582         cout << "after:  " << (int)(hres2.GetMean()*16384/360+.5) << " \xb1 " << (int)(hres2.GetRMS()*16384/360+.5) << " SE" << endl;
    583618        cout << endl;
    584         cout << "backw:  " << (int)(hres3.GetMean()*60+.5) << " \xb1 " << (int)(hres3.GetRMS()*60+.5) << " arcmin" << endl;
     619        cout << "before: " << Form("%.1f", hres1.GetMean()*60) << " \xb1 " << Form("%.1f", hres1.GetRMS()*60) << " arcmin" << endl;
     620        cout << "after:  " << Form("%.1f", hres2.GetMean()*60) << " \xb1 " << Form("%.1f", hres2.GetRMS()*60) << " arcmin" << endl;
     621        cout << endl;
     622        cout << "before: " << Form("%.1f", hres1.GetMean()*60*60/23.4) << " \xb1 " << Form("%.1f", hres1.GetRMS()*60*60/23.4) << " pix" << endl;
     623        cout << "after:  " << Form("%.1f", hres2.GetMean()*60*60/23.4) << " \xb1 " << Form("%.1f", hres2.GetRMS()*60*60/23.4) << " pix" << endl;
     624        cout << "after:  " << Form("%.1f", hres2.GetMean()*16384/360) << " \xb1 " << Form("%.1f", hres2.GetRMS()*16384/360) << " SE" << endl;
     625        cout << endl;
     626        cout << "backw:  " << Form("%.1f", hres3.GetMean()*60) << " \xb1 " << Form("%.1f", hres3.GetRMS()*60) << " arcmin" << endl;
    585627        cout << endl;                                            // ±
    586628   
    587    
     629
     630        before = hres1.GetMean()*16384/360;
     631        after  = hres2.GetMean()*16384/360;
     632        backw  = hres3.GetMean()*16384/360;
     633
    588634   
    589635        gStyle->SetOptStat(1110);
    … …  
    615661        cout << "Gaus-Fit  Sigma: " << f.GetParameter(2) << "\xb0" << endl;
    616662        cout << "Fit-Probability: " << f.GetProb()*100 << "%" << endl;
    617         cout << "      Chi^2/NDF: " << f.GetChisquare()/f.GetNDF() << endl;
     663        cout << "      Chi^2/NDF: " << f.GetChisquare() << "/" << f.GetNDF() << " = " << f.GetChisquare()/f.GetNDF() << endl;
    618664
    619665        c1->cd(1);
    … …  
    626672        gPad->Modified();
    627673        gPad->Update();
    628         for (int i=0; i<fCoordinates.GetSize(); i++)
    629             DrawSet(gPad, *(Set*)list.At(i));//, 10./hres1.GetMean());
     674        for (int i=0; i<fOriginal.GetSize(); i++)
     675            DrawSet(gPad, *(Set*)fOriginal.At(i));//, 10./hres1.GetMean());
    630676   
    631677        c1->cd(3);
    … …  
    638684        gPad->Update();
    639685        for (int i=0; i<fCoordinates.GetSize(); i++)
    640             DrawSet(gPad, *(Set*)fCoordinates.At(i), 10./hres2.GetMean());
     686            DrawSet(gPad, *(Set*)fCoordinates.At(i), 10./hres2.GetMean(), par[0]);
    641687
    642688        RaiseWindow();
    … …  
    663709
    664710            fin >> set;  // Read data from file [deg], it is stored in [rad]
    665 
    666711            if (!fin)
    667712                break;
    … …  
    687732                {
    688733                case kTbFit:
    689                     Fit();
    690                     DisplayBending();
     734                    {
     735                        Double_t before=0;
     736                        Double_t after=0;
     737                        Double_t backw=0;
     738                        Fit(before, after, backw);
     739                        DisplayBending();
     740                        DisplayResult(before, after, backw);
     741                    }
    691742                    return kTRUE;
    692743                case kTbLoad:
    693                     fBending.Load("bending_magic.txt");
     744                    fBending.Load("bending_new.txt");
    694745                    DisplayBending();
    695746                    return kTRUE;
    696747                case kTbSave:
    697                     fBending.Save("bending_magic.txt");
     748                    fBending.Save("bending_new.txt");
    698749                    return kTRUE;
    699750                case kTbLoadStars:
    … …  
    767818    }
    768819
     820    void DisplayResult(Double_t before, Double_t after, Double_t backw)
     821    {
     822        TGLabel *l1 = (TGLabel*)fLabel.At(3*MBending::GetNumPar()+1);
     823        l1->SetText(Form("Before: %.1f +- %.1f SE", before));
     824
     825        TGLabel *l2 = (TGLabel*)fLabel.At(3*MBending::GetNumPar()+2);
     826        l2->SetText(Form("After:  %.1f +- %.1f SE", after));
     827
     828        TGLabel *l3 = (TGLabel*)fLabel.At(3*MBending::GetNumPar()+3);
     829        l3->SetText(Form("Backw:  %.1f +- %.1f SE", backw));
     830    }
     831
    769832public:
    770833    ~MFit()
    … …  
    866929        fLabel.Add(l);
    867930
     931        l = new TGLabel(grp1, "");
     932        l->SetTextJustify(kTextLeft);
     933        grp1->AddFrame(l, hints5);
     934        fList->Add(l);
     935        fLabel.Add(l);
     936
     937        l = new TGLabel(grp1, "");
     938        l->SetTextJustify(kTextLeft);
     939        grp1->AddFrame(l, hints5);
     940        fList->Add(l);
     941        fLabel.Add(l);
     942
     943        l = new TGLabel(grp1, "");
     944        l->SetTextJustify(kTextLeft);
     945        grp1->AddFrame(l, hints5);
     946        fList->Add(l);
     947        fLabel.Add(l);
     948
    868949
    869950        ((TGCheckButton*)fList->FindWidget(0))->SetState(kButtonDown);
    … …  
    893974void gui()
    894975{
     976    gErrorIgnoreLevel = kError;
    895977    new MFit;
    896978}
Note: See TracChangeset for help on using the changeset viewer.