TTester.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 TTESTER_H
#define TTESTER_H

// ROOT
#include "TNamed.h"
#include "Math/ProbFuncMathCore.h"
#include "TList.h"
#include "TVector.h"

// Uranie
#include "UncertModeler.h"
#include "TDataServer.h"

namespace URANIE
{
namespace UncertModeler
{
class TTester: public TNamed
{
    // Attributes

protected:
    Long64_t _nS;                           
    Double_t _pValue;                       
    Double_t _dScore;                       
    Double_t _dModifiedScore; 
    TString _sAttributes; 
    TList* _listOfLaws;                   
    Bool_t _bDrawSame;               

public:
    enum ETestType
    {
        kKolmogorovSmirnov, kAndersonDarling, kCramerVonMises, kUnknown
    };
    ETestType _nTestType;                    
    enum ECase
    {
        kCase0, kCase1, kCase2, kCase3, kCaseUnknown
    };
    ECase _nCase;
    URANIE::DataServer::TDataServer *_tds;   
    Bool_t _blog;                            
    Bool_t _bquiet;                          

public:
    //---------------------------------------------
    TTester(URANIE::DataServer::TDataServer *tds, ETestType nTestType,
            const char *varexpinput, Option_t * option = "");

    virtual ~TTester();

    //---------------------------------------------
    Int_t getID()
    {
        return ID_UNCERTMODELER;
    }
    Double_t getScore()
    {
        return _dScore;
    }
    Double_t getModifiedScore()
    {
        return _dModifiedScore;
    }
    Double_t getPValue()
    {
        return _pValue;
    }

    ECase getCase()
    {
        return _nCase;
    }

    //---------------------------------------------

    Double_t computeScore(Option_t * option = "");

    virtual Double_t computeTheCurrentScore(Int_t i, Double_t dzi) = 0;
    virtual void computeModifiedScore(
            URANIE::DataServer::TStochasticAttribute::ELawType tlaw) = 0;

    virtual const char * getTitleStatistic() = 0;


    void addUniformLaw(Double_t min, Double_t max, Color_t color = kBlack, TString stitle = "");


    void addNormalLaw(Double_t mu, Double_t sigma, Color_t color = kBlack, TString stitle = "");


    void addLogNormalLaw(Double_t mu, Double_t sigma, Color_t color = kBlack, TString stitle = "");

    void addWeibullLaw(Double_t beta, Double_t m, Color_t color = kBlack, TString stitle = "");

    void addGammaLaw(Double_t alpha, Double_t beta, Color_t color = kBlack, TString stitle = "");


    void addBetaLaw(Double_t alpha, Double_t beta, Double_t min, Double_t max, Color_t color = kBlack, TString stitle = "");


    void treatOption(Option_t * option);


    //---------------------------------------------
    void computePValue(URANIE::DataServer::TStochasticAttribute::ELawType tlaw);


    Double_t computePValue_KS_CaseAll(Double_t dz);


    Double_t computePValue_AD_Case0(Double_t dz);

    Double_t computePValue_AD_Case3(Double_t dz, URANIE::DataServer::TStochasticAttribute::ELawType tlaw);

    Double_t computePValue_CvM_Case3(Double_t dz, URANIE::DataServer::TStochasticAttribute::ELawType tlaw);

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

protected:
    //---------------------------------------------

    void findParametersLaw(TString soption, TVectorD &vec);

    Int_t nTimesValueLTvalVector(Double_t dval, TVectorD vecValues);

    TGraph *createGraph(Color_t color = kBlack);


    ClassDef(URANIE::UncertModeler::TTester, ID_UNCERTMODELER)
    // Tester classe
};
} // Fin du namespace UncertModeler
} // Fin du namespace URANIE
#endif