genie::supernova::OrderedSampler Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/GenieSNova/src/Sampler.h"

Public Member Functions
	OrderedSampler ()
	~OrderedSampler ()
void	Init (TH1 *hist, bool integral=false)
bool	End ()
double	GenerateNext ()
int	GetBin ()
void	Reload ()
void	SetProbScale (double scale)
const size_t	Nbins ()
const double *	GetBinEdges ()
const double	GetBinValue (size_t n)

Protected Member Functions
double	BinStart ()
double	BinEnd ()
void	Step ()
void	StartNextBin ()

Protected Attributes
size_t	fBin
double	fPos
std::vector< double >	fBinEdges
std::vector< double >	fBinValues
double	fFreqCurrent
double	fProbScale

Detailed Description

Definition at line 22 of file Sampler.h.

Constructor & Destructor Documentation

OrderedSampler::OrderedSampler

(

)

Definition at line 44 of file Sampler.cxx.

: fProbScale(1) {}

OrderedSampler::~OrderedSampler

(

)

Definition at line 48 of file Sampler.cxx.

{}

Member Function Documentation

double OrderedSampler::BinEnd ( )

protected

Definition at line 155 of file Sampler.cxx.

References fBin, and fBinEdges.

Referenced by GenerateNext(), StartNextBin(), and Step().

{
    return fBinEdges[fBin + 1];
}

double OrderedSampler::BinStart ( )

protected

Definition at line 148 of file Sampler.cxx.

References fBin, and fBinEdges.

Referenced by GenerateNext(), and StartNextBin().

{
    return fBinEdges[fBin];
}

bool OrderedSampler::End

(

void

)

Definition at line 77 of file Sampler.cxx.

References fBin, and fBinEdges.

Referenced by genie::supernova::GFluxGenerator::End(), GenerateNext(), and StartNextBin().

{
    return fBin >= fBinEdges.size() - 1;
}

double OrderedSampler::GenerateNext

(

void

)

Definition at line 84 of file Sampler.cxx.

References BinEnd(), BinStart(), End(), fPos, LOG, pDEBUG, StartNextBin(), and Step().

Referenced by genie::supernova::GFluxGenerator::GenerateNext().

{
    LOG("Sampler", pDEBUG) << "from pos=" << fPos;

    if (this->End()) return 0;
    do {
        if (fPos >= this->BinEnd()) this->StartNextBin();
        if (End()) return 0;

        // Generate random step, using the average distance in current bin
        this->Step();
        LOG("Sampler", pDEBUG) << "now  pos=" << fPos
                               << " (bin =["  << this->BinStart()
                               << "," << this->BinEnd() << "])";
   } while (fPos >= this->BinEnd());

   LOG("Sampler", pDEBUG) << "return pos=" << fPos;
   return fPos;
}

int OrderedSampler::GetBin

(

)

Definition at line 106 of file Sampler.cxx.

References fBin.

Referenced by genie::supernova::GFluxGenerator::GenerateNext().

{
    return fBin;
}

const double * OrderedSampler::GetBinEdges

(

)

Definition at line 135 of file Sampler.cxx.

References fBinEdges.

{
    return fBinEdges.data();
}

const double OrderedSampler::GetBinValue

(

size_t

n

)

Definition at line 142 of file Sampler.cxx.

References fBinValues, fProbScale, and getGoodRuns4SAM::n.

Referenced by Step().

                                                 {
    return fBinValues[n] * fProbScale;
}

void OrderedSampler::Init	(	TH1 *	hist,
		bool	integral = false
	)

Definition at line 52 of file Sampler.cxx.

References fBin, fBinEdges, fBinValues, MECModelEnuComparisons::i, LOG, Nbins(), pDEBUG, pNOTICE, and StartNextBin().

Referenced by genie::supernova::GFluxGenerator::GFluxGenerator().

{
    LOG("Sampler", pDEBUG) << "init";

    fBin = 0;
    int Nbins = hist->GetNbinsX() + 1;

    fBinEdges.resize(Nbins);
    fBinValues.resize(Nbins);

    for (int i = 0; i < Nbins; ++i) {
        fBinEdges[i] = hist->GetBinLowEdge(i);
        fBinValues[i] = hist->GetBinContent(i);
        if(integral)
            fBinValues[i] /= hist->GetBinWidth(i);
    }

    // Start from bin#1, skipping underflow
    this->StartNextBin();
    LOG("Sampler", pNOTICE) << "Ordered sample with nbins=" << Nbins
                            << " [" << fBinEdges.front() << "," << fBinEdges.back() << "]";
}

const size_t OrderedSampler::Nbins

(

)

Definition at line 128 of file Sampler.cxx.

References fBinValues.

Referenced by Init().

{
    return fBinValues.size();
}

void OrderedSampler::Reload

(

)

Definition at line 113 of file Sampler.cxx.

References fBin, and StartNextBin().

Referenced by genie::supernova::GFluxGenerator::Reload().

{
    fBin = 0;
    this->StartNextBin();
}

void OrderedSampler::SetProbScale

(

double

scale

)

Definition at line 121 of file Sampler.cxx.

References fProbScale, and scale.

Referenced by genie::supernova::GFluxGenerator::SetProbScale().

{
    fProbScale = scale;
}

void OrderedSampler::StartNextBin ( )

protected

Definition at line 186 of file Sampler.cxx.

References BinEnd(), BinStart(), End(), fBin, fBinValues, fPos, LOG, and pDEBUG.

Referenced by GenerateNext(), Init(), and Reload().

 {
    // Go to next bin and update all the variables
    ++fBin;
    LOG("Sampler", pDEBUG) << "Starting bin#" << fBin << "/" << fBinValues.size()
                                              << "(" << this->End() << ")";

    // Start the steps from beginning of this bin
    fPos = this->BinStart();

    if (this->End()) return;

    LOG("Sampler", pDEBUG) << "Starting bin#" << fBin
                           << " [" << this->BinStart() << "," << this->BinEnd() << "]"
                           << " Freq=" << fBinValues[fBin];
}

void OrderedSampler::Step ( )

protected

Definition at line 162 of file Sampler.cxx.

References BinEnd(), e, fBin, fPos, GetBinValue(), genie::RandomGen::Instance(), LOG, pDEBUG, generate_hists::rnd, genie::RandomGen::RndFlux(), and fillBadChanDBTables::step.

Referenced by GenerateNext().

{
    double binValue = this->GetBinValue(fBin);

    // If the bin value is too small, move on to the next bin.
    if (binValue < 1e-15) {
        fPos = this->BinEnd();
        return;
    }

    genie::RandomGen * rnd = genie::RandomGen::Instance();

    float tau = 1. / binValue;
    double step = rnd->RndFlux().Exp(tau);

    LOG("Sampler", pDEBUG) << "step by " << step << ": " << fPos << "->" << fPos + step;

    fPos += step;

    return;
}

Member Data Documentation

size_t genie::supernova::OrderedSampler::fBin

protected

Definition at line 46 of file Sampler.h.

Referenced by BinEnd(), BinStart(), End(), GetBin(), Init(), Reload(), StartNextBin(), and Step().

std::vector<double> genie::supernova::OrderedSampler::fBinEdges

protected

Definition at line 48 of file Sampler.h.

Referenced by BinEnd(), BinStart(), End(), GetBinEdges(), and Init().

std::vector<double> genie::supernova::OrderedSampler::fBinValues

protected

Definition at line 49 of file Sampler.h.

Referenced by GetBinValue(), Init(), Nbins(), and StartNextBin().

double genie::supernova::OrderedSampler::fFreqCurrent

protected

Definition at line 51 of file Sampler.h.

double genie::supernova::OrderedSampler::fPos

protected

Definition at line 47 of file Sampler.h.

Referenced by GenerateNext(), StartNextBin(), and Step().

double genie::supernova::OrderedSampler::fProbScale

protected

Definition at line 54 of file Sampler.h.

Referenced by GetBinValue(), and SetProbScale().

---

The documentation for this class was generated from the following files:

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/GenieSNova/src/Sampler.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/GenieSNova/src/Sampler.cxx