TSamplerField.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 TSamplerField_H
#define TSamplerField_H

#include "TMatrixD.h"
#include "TVectorD.h"
#include "TVector3.h"
#include "TGraph2D.h"

#include "TSamplerStochastic.h"

namespace URANIE
{
namespace Sampler
{
class TSamplerField: public TSamplerStochastic
{
    enum SimulationMethod
    {
        kSphericalVariogram,
        kExponentialVariogram,
        kStableVariogram,
        kHyperbolicVariogram,
        kGaussianVariogram,
        kSinusCardinalVariogram,
        kSpectralMethod,
        kDilutionMethod,
        kTessellationMethod,
        kTurningBandsMethod
    };

protected:
    TMatrixD _vdc;                 
    TMatrixD _ap;                  
    TVectorD _scale;               
    TVectorD _seed;                
    TGraph2D* _graph2D;             
    vector<URANIE::DataServer::UEntry*> in;
    
public:
    Int_t _nXSize;
    Int_t _nYSize;
    Int_t _nZSize;
    Int_t _nRandomFunction;        
    Int_t _nMethod;                
    Int_t _nVariogram;             
    Bool_t _bIsotropy;             
    Double_t _dVariance;           
    Double_t _dSigma;              
    TVector3 _vScaleFactor;        
    Double_t _valcrt;
protected:
//       virtual void createTuple();

    // The spectral method of simulation
    void SpectralInit();
    void SpectralPrepare();
    void SpectralSimulate(unsigned int ievt);

    void VDC();

public:
    //---------------------------------------------
    TSamplerField(URANIE::DataServer::TDataServer *tds, Option_t *option =
                  "Gauss", int nS=1);
    virtual ~TSamplerField();

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

    void generateSample(Option_t *option = "");

    //---------------------------------------------
    /*
     In this case, the space is represented by TAttribute attributes.
     */
    void createListOfAttributes();

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

    void setMethod(TString str);

    TString getMethodName();


    void setVariogram(TString str);

    TString getVariogramName();


    // Setter and Getter
    void SetRandomFunction(Int_t nind)
    {
        _nRandomFunction = nind;
    }
    Int_t GetRandomFunction()
    {
        return _nRandomFunction;
    }
    void SetIsotropy(Bool_t bval)
    {
        _bIsotropy = bval;
    }
    Bool_t IsIsotropy()
    {
        return _bIsotropy;
    }
    void SetVariance(Double_t dval);
 
    Double_t GetVariance()
    {
        return _dVariance;
    }
    void SetScaleFactor(Double_t dval);

    void Draw2D(TString sx, TString sy, Option_t* opt = "TRI2");

//       void AddAttribute(Option_t *option, URANIE::DataServer::TAttribute *att, Int_t ind = 100);

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

    ClassDef(URANIE::Sampler::TSamplerField, ID_SAMPLER)
    //Abstract class for generation of random design.
};

} // Fin du namespace Sampler
} // Fin du namespace URANIE
#endif
// fin du fichier $RCSfile$.