ISF::PDFcreator Class Reference

Creates random numbers with a distribution given as ROOT TF1. More...

#include <PDFcreator.h>

Collaboration diagram for ISF::PDFcreator:

Public Member Functions
	PDFcreator ()
	construct the class with a given TF1 and a random engine
	~PDFcreator ()
void	setName (const std::string &PDFname)
	all following is used to set up the class
void	addToEnergyEtaHist1DMap (int energy, int etaMin, TH1 *hist)
double	getRand (CLHEP::HepRandomEngine *rndmEngine, const std::vector< int > &inputParameters) const
	get the random value with this method, by providing the input parameters
const std::string &	getName () const

Private Attributes
std::string	m_name
	Give pdf a name for debug purposes.
std::map< int, std::map< int, TH1 * > >	m_energy_eta_hists1D
	map of energies to map of eta ranges to 1D histograms

Detailed Description

Creates random numbers with a distribution given as ROOT TF1.

The TF1 function parameters will be retrieved from a histogram given by addPar.

Author

Elmar Ritsch Elmar.nosp@m..Rit.nosp@m.sch@c.nosp@m.ern..nosp@m.ch @maintainer/updater Thomas Carter thoma.nosp@m.s.mi.nosp@m.chael.nosp@m..car.nosp@m.ter@c.nosp@m.ern..nosp@m.ch

Definition at line 36 of file PDFcreator.h.

Constructor & Destructor Documentation

PDFcreator()

ISF::PDFcreator::PDFcreator ( )

inline

construct the class with a given TF1 and a random engine

Definition at line 41 of file PDFcreator.h.

{};

~PDFcreator()

ISF::PDFcreator::~PDFcreator

(

)

Definition at line 32 of file PDFcreator.cxx.

                           {
  for(auto & each_outter : m_energy_eta_hists1D) {
    for(auto & each_inner : each_outter.second) {
      if(each_inner.second != nullptr)
        delete each_inner.second;
    }
  }
}

Member Function Documentation

addToEnergyEtaHist1DMap()

void ISF::PDFcreator::addToEnergyEtaHist1DMap	(	int	energy,
		int	etaMin,
		TH1 *	hist )

Definition at line 41 of file PDFcreator.cxx.

                                                                             {
  // Make a local copy
  TH1* localHist = (TH1*) hist->Clone();
  localHist->SetDirectory(0);
 
  if(m_energy_eta_hists1D.find(energy) != m_energy_eta_hists1D.end()){ //if energy entry exists, insert into inner eta map
    (m_energy_eta_hists1D.find(energy)->second).insert(std::make_pair(etaMin, localHist));
  }
  else{ //if energy entry does not exist create new full energy entry
    std::map< int, TH1*> inner;
    inner.insert(std::make_pair(etaMin, localHist));
    m_energy_eta_hists1D.insert(std::make_pair(energy, inner));
  }
}

getName()

const std::string & ISF::PDFcreator::getName ( ) const

inline

Definition at line 51 of file PDFcreator.h.

{return m_name;};

getRand()

double ISF::PDFcreator::getRand	(	CLHEP::HepRandomEngine *	rndmEngine,
		const std::vector< int > &	inputParameters ) const

get the random value with this method, by providing the input parameters

Definition at line 56 of file PDFcreator.cxx.

{
 
    //define variable to return from getRand call, should never return zero
    float randomHist = 0.;
 
    if( m_energy_eta_hists1D.find(inputParameters.at(0)) == m_energy_eta_hists1D.end()) {
        //this should never be reached
        return 0.;
    }
 
    const std::map< int, TH1*>& etaMin_hists = m_energy_eta_hists1D.at(inputParameters.at(0));
 
    if( etaMin_hists.find(inputParameters.at(1)) == etaMin_hists.end()) {
        //this should never be reached
        return 0.;
    }
 
    //get the chosen histogram from the map
    TH1*  hist = etaMin_hists.at(inputParameters.at(1));
 
    //Draw randomly from the histogram CDF distribution.
    double randomShoot = CLHEP::RandFlat::shoot(rndmEngine);
 
    //first select the matching CDF bin
    int iBinSelect = 0;
    for(int iBin = 0; iBin < hist->GetNbinsX(); iBin ++){
        iBinSelect = iBin;
        if(hist->GetBinContent(iBin) > randomShoot){
            break;
        }
    }
 
    //Select random value within bin
    randomHist = CLHEP::RandFlat::shoot(rndmEngine, hist->GetBinLowEdge(iBinSelect) , hist->GetBinLowEdge(iBinSelect+1));
 
 
    return randomHist;
 
}

setName()

void ISF::PDFcreator::setName ( const std::string & PDFname )

inline

all following is used to set up the class

Definition at line 46 of file PDFcreator.h.

{ m_name = PDFname; }; //get the pdf's name

Member Data Documentation

m_energy_eta_hists1D

std::map< int , std::map< int, TH1*> > ISF::PDFcreator::m_energy_eta_hists1D

private

map of energies to map of eta ranges to 1D histograms

Definition at line 55 of file PDFcreator.h.

m_name

std::string ISF::PDFcreator::m_name

private

Give pdf a name for debug purposes.

Definition at line 54 of file PDFcreator.h.

---

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

• PDFcreator.h
• PDFcreator.cxx