NagFitStrategy.h

#ifndef NAGFITSTRATEGY_H
#define NAGFITSTRATEGY_H
#include "../FitStrategy.h"
#include <nag.h>
#include <nage04.h>
#include <vector>

using namespace Sct;
using namespace std;

namespace SctFitter{
    class NagFitStrategy : public FitStrategy {
    public:
    NagFitStrategy(string opt=string("")) throw();
    ~NagFitStrategy() throw();
    virtual void fitTH1(const TH1& hist, TF1& fit) const throw(LogicError, MathsError);
    virtual void fitTGraph(const TGraph& graph, TF1& fit) const throw(LogicError, MathsError);
    virtual void fitTGraphErrors(const TGraphErrors& graph, TF1& fit) const throw(LogicError, MathsError);
    virtual void fitTGraphAsymmErrors(const TGraphAsymmErrors& graph, TF1& fit) const throw(LogicError, MathsError);
    
    virtual const string& getName() const throw();
    private:
    NagFitStrategy( const NagFitStrategy& ); 
    const string name; 
    static bool inMap; 
    
        const int getNumberFitPoints(TH1& hist, const TF1& fit, 
                        int& iFirstBin) const throw(LogicError);
    void getFitPoints(TGraph& graph, const TF1& fit, vector<bool>& active) const throw(LogicError);
    class Cache{
    public:
        Cache(const TH1& hist, const unsigned npoints, 
          const unsigned firstpoint, TF1& fit) throw(LogicError, MathsError) ;

        Cache(const TGraph& graph, TF1& fit, const vector<bool>& active) throw(LogicError, MathsError) ;

        Cache(const TGraphErrors& graph, TF1& fit, const vector<bool>& active) throw(LogicError, MathsError) ;

        Cache(const TGraphAsymmErrors& graph, TF1& fit, const vector<bool>& active) throw(LogicError, MathsError) ;

        void convertPars(double inPars[]);
                
        void setupPars(TF1& fit); 
        
        ~Cache() throw();                
        vector<double> m_y; 
        vector<double> m_x; 
        vector<double> m_ey; 
        TF1* function;           
        double* pars;       
        unsigned int nVarPars;  
        double* inPars;     
        
    private:        
        int* map;               
    };

    void nagFit(Cache* cache, TF1& fit, bool quiet) const throw(MathsError);

    bool ranged() const throw();

    bool quiet() const throw();

    static void chiSquaredFunctionForNag(Integer m, Integer n, double x[], double fvec[],
                         Nag_Comm *comm) throw();
    };

} // end of namespace SctFitter;


#endif // #ifndef NAGFITSTRATEGY_H

---