TSequence.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 TSequence_H
#define TSequence_H

#include "Rtypes.h"
#include "TNamed.h"
#include "TSampler.h"

namespace URANIE
{
namespace Sampler
{
class TSequence: public TNamed
{

    // Attributes
public:
    ULong_t _nDim;             
    Double_t * _dValue;     //[_nDim]   ///< Pointeur
    static Int_t _seqCounter;  
    static Bool_t _blog;       

    // Operations
public:
    //---------------------------------------------
    TSequence();
    TSequence(ULong_t ndim);
    virtual ~TSequence();

    //---------------------------------------------
    virtual void setInstance(const ULong_t n) = 0;
    virtual void init() = 0;

    inline Double_t operator[](const ULong_t) const;
    inline Double_t getComponent(const ULong_t) const;

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

    ClassDef(URANIE::Sampler::TSequence, ID_SAMPLER)
    //Classe de g�n�ration par Quasi Monte-Carlo
};


inline Double_t TSequence::operator[](const ULong_t j) const
{
    if (j >= _nDim)
    {
        throw URANIE::Exceptions::UErrorExceptions(__FILE__, __LINE__,
                Form(
                        "TSequence::operator[] ERROR : the indice %lu is greater than the dimension[%lu]",
                        j, _nDim));
    }
    return (_dValue[j]);
}


inline Double_t TSequence::getComponent(const ULong_t j) const
{
    if (j >= _nDim)
    {
        throw URANIE::Exceptions::UErrorExceptions(__FILE__, __LINE__,
                Form(
                        "TSequence::getComponent ERROR : the indice %lu is greater than the dimension[%lu]",
                        j, _nDim));
    }
    return (_dValue[j]);
}

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