TMorris.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/>.
// TMorris mer sep 13 12:39:50 CEST 2006
// $Id$
// $Author$
// $Date$
// $Revision$
// $State$

#ifndef TMORRIS_H
#define TMORRIS_H

#ifndef ROOT_TMatrixD
#include "TMatrixD.h"
#endif
#include "TVectorD.h"
#include "TRandom2.h"

#include "TScreeming.h"

namespace URANIE
{
namespace Sensitivity
{
class TMorris: public TScreeming
{
private:
    Double_t _delta; 
    Int_t _level; 
    TMatrixD _matB;    
    TVectorD _vecDirection;    
    TNtupleD* _ntdatamorris;           
    URANIE::DataServer::TDSNtupleD* _ntmorrisresult;   
    TRandom2 * _rand;

    void stdInit();
public:
    //---------------------------------------------
    TMorris(URANIE::DataServer::TDataServer *tds, const char *fcn, Int_t nReplique,
            Int_t level, Double_t delta = 0, TString sout = "");
#ifdef ROOT5
    TMorris(URANIE::DataServer::TDataServer *tds, void *fcn, Int_t nReplique,
            Int_t level, Double_t delta = 0, TString sout = "");
#else
    TMorris(URANIE::DataServer::TDataServer *tds, void (*fcn)(Double_t*,Double_t*), 
            TString sin, TString sout,
            Int_t nReplique, Int_t level, Double_t delta = 0);
#endif
    TMorris(URANIE::DataServer::TDataServer *tds, const char *varexpinput,
            const char *varexpoutput, Option_t * option = "");
    TMorris(URANIE::DataServer::TDataServer *tds,
            URANIE::Launcher::TCode *fcode, Int_t nReplique, Int_t level,
            Double_t delta = 0);
    TMorris(URANIE::DataServer::TDataServer *tds,
            URANIE::Relauncher::TRun *run, Int_t nReplique, Int_t level,
            Double_t delta = 0);
    ~TMorris();
    //---------------------------------------------
    URANIE::DataServer::TDSNtupleD* getMorrisResults()
    {
        return _ntmorrisresult;
    }

    TTree *getResultTuple(bool commonresulttuple=false);

    void setLevel(Int_t val)
    {
        _level = val;
    }
    Int_t getLevel()
    {
        return _level;
    }
    void setDelta(Double_t val)
    {
        _delta = val;
    }
    Double_t getDelta()
    {
        return _delta;
    }
    Int_t getNReplica()
    {
        return _nS;
    }

    //---------------------------------------------
    virtual void generateSample(Option_t * option = "");


    void computeOrientationMatrixOld();
    void computeOrientationMatrix();

    //---------------------------------------------
    void evaluateIndexes(Option_t * option = "");


    void preTreatment();
    
    //---------------------------------------------
    virtual void createTuple();
    virtual void drawIndexes(TString sTitre, const char *select = "", Option_t * option = "", bool internalCall=kFALSE);
    void drawSample(TString svar = "", Int_t nreplique = 0, Option_t *option="");
    void drawIndexes(Option_t * option = "", bool internalCall=kFALSE);
    void drawIndexesAll(Option_t * option = "");
    void pyDrawIndexes(Option_t * option = "");

    //---------------------------------------------
    void printLog(Option_t *option = "");

    ClassDef(URANIE::Sensitivity::TMorris, ID_SENSITIVITY)
    //Classe de
};

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