Ignore:
Timestamp:
01/09/08 11:45:59 (17 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

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

    r8766 r8802  
    485485        gmaz.SetMarkerStyle(kFullDotMedium);;
    486486        gmzd.SetMarkerStyle(kFullDotMedium);
    487    
     487
    488488        TMinuit minuit(MPointing::GetNumPar());  //initialize TMinuit with a maximum of 5 params
    489489        minuit.SetObjectFit(this);
    490490        minuit.SetPrintLevel(-1);
    491491        minuit.SetFCN(fcn);
    492    
     492
    493493        fBending.SetMinuitParameters(minuit, MPointing::GetNumPar()); // Init Parameters [deg]
    494    
     494
    495495        for (int i=0; i<MPointing::GetNumPar(); i++)
    496496        {
     
    500500                minuit.Release(i);
    501501        }
    502    
     502
    503503        //minuit.Command("SHOW PARAMETERS");
    504504        //minuit.Command("SHOW LIMITS");
    505    
     505
    506506        cout << endl;
    507507        cout << "Starting fit..." << endl;
     
    509509        cout << "0.02deg (=1SE) is assumed." << endl;
    510510        cout << endl;
    511    
     511
    512512        Int_t ierflg = 0;
    513513        ierflg = minuit.Migrad();
     
    516516        ierflg = minuit.Migrad();
    517517        cout << "Migrad returns " << ierflg << endl << endl;
    518    
     518
    519519        //
    520520        // Get Fit Results
     
    524524        cout << endl;
    525525        //fBending.Save("bending_magic.txt");
    526    
    527    
     526
     527
    528528        //
    529529        // Make a copy of all list entries
     
    533533        for (int i=0; i<fCoordinates.GetSize(); i++)
    534534            list.Add(new Set(*(Set*)fCoordinates.At(i)));
    535    
     535
    536536        //
    537537        // Correct for Offsets only
     
    541541        for (int i=2; i<MPointing::GetNumPar(); i++)
    542542            par[i]=0;
    543    
     543
    544544        MPointing b2;
    545545        b2.SetParameters(par);
     
    989989    }
    990990
    991     void LoadCollection(TString fname)
     991    Bool_t LoadCollection(TString fname)
    992992    {
    993993        ifstream fin(fname);
     
    995995        {
    996996            cout << "Collection '" << fname << "' not found!" << endl;
    997             return;
     997            return kFALSE;
    998998        }
    999999
     
    10131013            if (!line.EndsWith(".txt"))
    10141014            {
    1015                 cout << "WARNING: " << line << endl;
     1015                cout << "WARNING - Wring extension: " << line << endl;
    10161016                continue;
    10171017            }
    10181018
    1019             LoadStars(line);
    1020         }
    1021     }
    1022 
    1023 
    1024     void LoadStars(TString fname="tpoint.txt")
    1025     {
    1026         if (fname.EndsWith(".col"))
    1027         {
    1028             LoadCollection(fname);
    1029             fFileNameStars = fname;
    1030             return;
    1031         }
    1032 
     1019            if (!LoadStarFile(line))
     1020                return kFALSE;
     1021        }
     1022        return kTRUE;
     1023    }
     1024
     1025    Bool_t LoadStarFile(TString fname)
     1026    {
    10331027        const Int_t size = fOriginal.GetSize();
    10341028
     
    10411035        {
    10421036            cout << "File '" << fname << "' not found!" << endl;
    1043             return;
     1037            return kFALSE;
    10441038        }
    10451039
     
    10591053        cout << " (Total=" << fOriginal.GetSize() << ")" << endl;
    10601054
     1055        return kTRUE;
     1056    }
     1057
     1058    Bool_t LoadStars(TString fname="tpoint.txt")
     1059    {
     1060        const Bool_t iscol = fname.EndsWith(".col");
     1061
     1062        const Bool_t rc = iscol ? LoadCollection(fname) : LoadStarFile(fname);
     1063        if (!rc)
     1064            return kFALSE;
     1065
    10611066        fFileNameStars = fname;
     1067        SetWindowName(fname);
     1068
     1069        return kTRUE;
    10621070    }
    10631071
Note: See TracChangeset for help on using the changeset viewer.