TABC.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/>.
// TABC
// $Id$
// $Author$
// $Date$
// $Revision 1.2 $
// $State$

#ifndef TABC_H
#define TABC_H
 
#include "TTree.h"
#include "TMatrixD.h"
#include "TRandom3.h"
#include "TMath.h"
 
// Uranie
#include "TCalibration.h"
#include "TDataServer.h"
 
namespace URANIE
{
namespace Calibration
{
 
class TABC : public URANIE::Calibration::TCalibration
{
public :
    vector<double>                                     _values;
    URANIE::DataServer::TDataServer *_tdsPosterior;     
    
protected:    

    double _threshold; 
    Int_t _kPosterior;
    Int_t _kEval;
    
    TRandom3 *_rand;

    Bool_t _boolEpsilon;
    Double_t _Epsilon;


public :

    //---------------------------------------------
    
    TABC(URANIE::DataServer::TDataServer *tds, URANIE::Relauncher::TRun *run, int ns = 100, Option_t *option = "");
    
    TABC(URANIE::DataServer::TDataServer *tds, void (*fcn)(Double_t*,Double_t*),  const char *varexpinput, const char *varexpoutput, int ns = 100, Option_t *option = "");
    
    TABC(URANIE::DataServer::TDataServer *tds, const char *fcn, const char *varexpinput, const char *varexpoutput, int ns = 100, Option_t *option = "");
    
    TABC(URANIE::DataServer::TDataServer *tds, URANIE::Launcher::TCode *fcode, int ns = 1, Option_t *option = "");
    
    virtual ~TABC();

    //---------------------------------------------
    virtual void setDistance(const char *distName, URANIE::DataServer::TDataServer *tdsRef, const char *input, const char *reference, const char *weight="");

    virtual void setDistance(URANIE::Calibration::TDistanceLikelihoodFunction *distFunc, URANIE::DataServer::TDataServer *tdsRef, const char *input, const char *reference, const char *weight="");
    

    void initABC();
    
    
    //---------------------------------------------
    //---------------------------------------------
    void setUpABC();
//     void evalABC();
    void computeABC(int nEval);
    void computeABCscore(Double_t newEpsilon);
    
    void posteriorToPar();

    void checktdsParContent();
    
    //---------------------------------------------
    void setGaussianNoise(const char *stdname);
       
    void setPercentile(Double_t eps)
    {
        try
        {
            if ((eps>1.0) || (eps<0.0))
            {
                throw URANIE::Exceptions::UErrorExceptions(__FILE__, __LINE__,
                                                           Form( "TABC::setPercentile: the percentile must be set between 0 and 1: however it is fix at %4.2f ! ", eps));
            }
        }
        catch (URANIE::Exceptions::UErrorExceptions& ue)
        {
            ue.printMessage();
            throw ue;
        }
        
        _boolEpsilon = true;
        _Epsilon = eps;
    }

    URANIE::DataServer::TDataServer *getTdsPosterior()
    {
        return _tdsPosterior;
    }
    double getThreshold()
    {
        try
        {
            if (_threshold==12345678.9)
            {
                throw URANIE::Exceptions::UErrorExceptions(__FILE__, __LINE__,
                         "TABC::getThreshold: the threshold is not yet computed");
            }
        }
        catch (URANIE::Exceptions::UErrorExceptions& ue)
        {
            ue.printMessage();
            throw ue;
        } 
        return _threshold;
    }
    int getNEval()
    {
        return _kEval;
    }
    
    
    
    

    ClassDef(URANIE::Calibration::TABC, ID_CALIBRATION)
 
};
 
} // End of namespace ABC
} // End of namespace URANIE
#endif