TGPCostFunction.h

Go to the documentation of this file.

// Copyright (C) 2013-2024 CEA/DES
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as published
// by the Free Software Foundation, either version 3 of the License, or any
// later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with this program. If not, see <http://www.gnu.org/licenses/>.
// $Id$
// $Author$
// $Date$
// $Revision$
// $State$

#ifndef TGPCOSTFUNCTION_H
#define TGPCOSTFUNCTION_H    

#include "Modeler.h"
#include "TGPBuilder.h"
#include "TShareEval.h"

namespace URANIE
{
namespace Modeler
{
//class TGPBuilder;
class TGPCostFunction: public URANIE::Relauncher::TSimpleEval
{

    // Attributes
protected:
    Double_t* _H;    
    Double_t* _K;    
    Double_t* _C;    
    Double_t* _yObs; 
    Int_t _n;    
    Int_t _p;    

public:
    Bool_t _blog;          
    TGPBuilder* gpBuilder; 

    // Operations
public:
    //---------------------------------------------

    TGPCostFunction(TGPBuilder* gpb);

    virtual ~TGPCostFunction();

    //---------------------------------------------
    TGPBuilder* getBuilder()
    {
        return gpBuilder;
    }
    
    Double_t* getCovarianceMatrix()
    {
        return _K;
    }
    
    Double_t* getCorrelationMatrix()
    {
        return _C;
    }
    
    Double_t* getObservationVector()
    {
        return _yObs;
    }
    
    Int_t getNObs()
    {
        return _n;
    }

    void setBuilder(TGPBuilder* gpb)
    {
        gpBuilder = gpb;
    }

    //---------------------------------------------
    virtual Int_t eval(Double_t* theta, Double_t* res, int=0) = 0;

    //---------------------------------------------
    void setLog()
    {
        _blog = kTRUE;
    }
    void unsetLog()
    {
        _blog = kFALSE;
    }
    void changeLog()
    {
        _blog = !_blog;
    }
    Bool_t getLog()
    {
        return _blog;
    }
    virtual void printLog(Option_t *option = "");

    ClassDef(URANIE::Modeler::TGPCostFunction, ID_MODELER)
    //Classe de

protected:
    void duplicateGPB();
};

} // Fin du namespace Modeler
} // Fin du namespace URANIE
#endif