#include <TFCS2DFunctionHistogram.h>

Inheritance diagram for TFCS2DFunctionHistogram:

Collaboration diagram for TFCS2DFunctionHistogram:

Public Member Functions
	TFCS2DFunctionHistogram (TH2 *hist=nullptr)
	~TFCS2DFunctionHistogram ()
void	Initialize (TH2 *hist)
virtual void	rnd_to_fct (float &valuex, float &valuey, float rnd0, float rnd1) const
const std::vector< float > &	get_HistoBordersx () const
std::vector< float > &	get_HistoBordersx ()
const std::vector< float > &	get_HistoBordersy () const
std::vector< float > &	get_HistoBordersy ()
const std::vector< float > &	get_HistoContents () const
std::vector< float > &	get_HistoContents ()
virtual void	rnd_to_fct (float value[], const float rnd[]) const
	Function gets array of random numbers rnd[] in the range [0,1) as arguments and returns function value in array value.
virtual int	ndim () const
	Return the number of dimensions for the function.
virtual std::size_t	MemorySize () const
	Gives the total memory size, including the size of additional memory allocated inside the class.
bool	msgLvl (const MSG::Level lvl) const
	Check whether the logging system is active at the provided verbosity level.
MsgStream &	msg () const
	Return a stream for sending messages directly (no decoration)
MsgStream &	msg (const MSG::Level lvl) const
	Return a decorated starting stream for sending messages.
MSG::Level	level () const
	Retrieve output level.
virtual void	setLevel (MSG::Level lvl)
	Update outputlevel.

Static Public Member Functions
static void unit_test	ATLAS_NOT_THREAD_SAFE (TH2 *hist=nullptr)
static double	CheckAndIntegrate2DHistogram (const TH2 *hist, std::vector< double > &integral_vec, int &first, int &last)
static void unit_test	ATLAS_NOT_THREAD_SAFE (TH2 hist=nullptr, TFCS2DFunction rtof=nullptr, const char *outfilename="TFCS2DFunction_unit_test.root", int nrnd=10000000)
static void unit_tests	ATLAS_NOT_THREAD_SAFE (TH2 hist=nullptr, const char outfilename="TFCS2DFunction_unit_test.root", int nrnd=10000000)
static std::string	startMsg (MSG::Level lvl, const std::string &file, int line)
	Make a message to decorate the start of logging.

Protected Attributes
std::vector< float >	m_HistoBorders
std::vector< float >	m_HistoBordersy
std::vector< float >	m_HistoContents

Private Attributes
std::string	m_nm
	Message source name.

Static Private Attributes
static boost::thread_specific_ptr< MsgStream > m_msg_tls	ATLAS_THREAD_SAFE
	Do not persistify!

Detailed Description

Definition at line 15 of file TFCS2DFunctionHistogram.h.

Constructor & Destructor Documentation

◆ TFCS2DFunctionHistogram()

TFCS2DFunctionHistogram::TFCS2DFunctionHistogram ( TH2 * hist = nullptr )

inline

Definition at line 17 of file TFCS2DFunctionHistogram.h.

                                               {
    if (hist)
      Initialize(hist);
  };

◆ ~TFCS2DFunctionHistogram()

TFCS2DFunctionHistogram::~TFCS2DFunctionHistogram ( )

inline

Definition at line 21 of file TFCS2DFunctionHistogram.h.

21{};

Member Function Documentation

◆ ATLAS_NOT_THREAD_SAFE() [1/3]

void unit_tests TFCS2DFunction::ATLAS_NOT_THREAD_SAFE	(	TH2 *	hist = nullptr,
		const char *	outfilename = "TFCS2DFunction_unit_test.root",
		int	nrnd = 10000000 )

staticinherited

◆ ATLAS_NOT_THREAD_SAFE() [2/3]

void unit_test TFCS2DFunction::ATLAS_NOT_THREAD_SAFE	(	TH2 *	hist = nullptr,
		TFCS2DFunction *	rtof = nullptr,
		const char *	outfilename = "TFCS2DFunction_unit_test.root",
		int	nrnd = 10000000 )

staticinherited

◆ ATLAS_NOT_THREAD_SAFE() [3/3]

void unit_test TFCS2DFunctionHistogram::ATLAS_NOT_THREAD_SAFE ( TH2 * hist = nullptr )

static

◆ CheckAndIntegrate2DHistogram()

double TFCS2DFunction::CheckAndIntegrate2DHistogram	(	const TH2 *	hist,
		std::vector< double > &	integral_vec,
		int &	first,
		int &	last )

staticinherited

Definition at line 26 of file TFCS2DFunction.cxx.

                                                                             {
  ISF_FCS::MLogging logger;
  Int_t nbinsx = hist->GetNbinsX();
  Int_t nbinsy = hist->GetNbinsY();
  Int_t nbins = nbinsx * nbinsy;
 
  float integral = 0.0;
  float hint = hist->Integral();
  integral_vec.resize(nbins);
 
  for (int ix = 1; ix <= nbinsx; ix++) {
    for (int iy = 1; iy <= nbinsy; iy++) {
      int globalbin = (ix - 1) * nbinsy + iy - 1;
      float binval = hist->GetBinContent(ix, iy);
      if (binval < 0) {
        // Can't work if a bin is negative, forcing bins to 0 in this case
        double fraction = binval / hint;
        if (std::abs(fraction) > 1e-5) {
          ATH_MSG_NOCLASS(logger,
                          "Warning: bin content is negative in histogram "
                              << hist->GetName() << " : " << hist->GetTitle()
                              << " binval=" << binval << " " << fraction * 100
                              << "% of integral=" << hist->Integral()
                              << ". Forcing bin to 0.");
        }
        binval = 0;
      }
      integral += binval;
      integral_vec[globalbin] = integral;
    }
  }
 
  for (first = 0; first < nbins; first++)
    if (integral_vec[first] != 0)
      break;
  for (last = nbins - 1; last > 0; last--)
    if (integral_vec[last] != integral)
      break;
  last++;
 
  if (integral <= 0) {
    ATH_MSG_NOCLASS(logger, "Error: histogram "
                                << hist->GetName() << " : " << hist->GetTitle()
                                << " integral=" << integral << " is <=0");
  }
 
  return integral;
}

◆ get_HistoBordersx() [1/2]

std::vector< float > & TFCS2DFunctionHistogram::get_HistoBordersx ( )

inline

Definition at line 32 of file TFCS2DFunctionHistogram.h.

32{ return m_HistoBorders; };

TFCS2DFunctionHistogram::m_HistoBorders

std::vector< float > m_HistoBorders

Definition TFCS2DFunctionHistogram.h:45

◆ get_HistoBordersx() [2/2]

const std::vector< float > & TFCS2DFunctionHistogram::get_HistoBordersx ( ) const

inline

Definition at line 29 of file TFCS2DFunctionHistogram.h.

                                                    {
    return m_HistoBorders;
  };

◆ get_HistoBordersy() [1/2]

std::vector< float > & TFCS2DFunctionHistogram::get_HistoBordersy ( )

inline

Definition at line 36 of file TFCS2DFunctionHistogram.h.

36{ return m_HistoBordersy; };

TFCS2DFunctionHistogram::m_HistoBordersy

std::vector< float > m_HistoBordersy

Definition TFCS2DFunctionHistogram.h:46

◆ get_HistoBordersy() [2/2]

const std::vector< float > & TFCS2DFunctionHistogram::get_HistoBordersy ( ) const

inline

Definition at line 33 of file TFCS2DFunctionHistogram.h.

                                                    {
    return m_HistoBordersy;
  };

◆ get_HistoContents() [1/2]

std::vector< float > & TFCS2DFunctionHistogram::get_HistoContents ( )

inline

Definition at line 40 of file TFCS2DFunctionHistogram.h.

40{ return m_HistoContents; };

TFCS2DFunctionHistogram::m_HistoContents

std::vector< float > m_HistoContents

Definition TFCS2DFunctionHistogram.h:47

◆ get_HistoContents() [2/2]

const std::vector< float > & TFCS2DFunctionHistogram::get_HistoContents ( ) const

inline

Definition at line 37 of file TFCS2DFunctionHistogram.h.

                                                    {
    return m_HistoContents;
  };

◆ Initialize()

void TFCS2DFunctionHistogram::Initialize ( TH2 * hist )

Definition at line 18 of file TFCS2DFunctionHistogram.cxx.

                                                  {
  Int_t nbinsx = hist->GetNbinsX();
  Int_t nbinsy = hist->GetNbinsY();
  Int_t nbins = nbinsx * nbinsy;
 
  float integral = 0;
  m_HistoBorders.resize(nbinsx + 1);
  m_HistoBordersy.resize(nbinsy + 1);
  m_HistoContents.resize(nbins);
  int ibin = 0;
  for (int iy = 1; iy <= nbinsy; iy++) {
    for (int ix = 1; ix <= nbinsx; ix++) {
      float binval = hist->GetBinContent(ix, iy);
      if (binval < 0) {
        // Can't work if a bin is negative, forcing bins to 0 in this case
        double fraction = binval / hist->Integral();
        if (TMath::Abs(fraction) > 1e-5) {
          ATH_MSG_WARNING("bin content is negative in histogram "
                          << hist->GetName() << " : " << hist->GetTitle()
                          << " binval=" << binval << " " << fraction * 100
                          << "% of integral=" << hist->Integral()
                          << ". Forcing bin to 0.");
        }
        binval = 0;
      }
      integral += binval;
      m_HistoContents[ibin] = integral;
      ++ibin;
    }
  }
  if (integral <= 0) {
    ATH_MSG_ERROR("histogram " << hist->GetName() << " : " << hist->GetTitle()
                               << " integral=" << integral << " is <=0");
    m_HistoBorders.resize(0);
    m_HistoBordersy.resize(0);
    m_HistoContents.resize(0);
    return;
  }
 
  for (int ix = 1; ix <= nbinsx; ix++)
    m_HistoBorders[ix - 1] = hist->GetXaxis()->GetBinLowEdge(ix);
  m_HistoBorders[nbinsx] = hist->GetXaxis()->GetXmax();
 
  for (int iy = 1; iy <= nbinsy; iy++)
    m_HistoBordersy[iy - 1] = hist->GetYaxis()->GetBinLowEdge(iy);
  m_HistoBordersy[nbinsy] = hist->GetYaxis()->GetXmax();
 
  for (ibin = 0; ibin < nbins; ++ibin)
    m_HistoContents[ibin] /= integral;
}

◆ level()

MSG::Level ISF_FCS::MLogging::level ( ) const

inlineinherited

Retrieve output level.

Definition at line 201 of file MLogging.h.

201{ return msg().level(); }

ISF_FCS::MLogging::msg

MsgStream & msg() const

Return a stream for sending messages directly (no decoration)

Definition MLogging.h:231

◆ MemorySize()

virtual std::size_t TFCSFunction::MemorySize ( ) const

inlinevirtualinherited

Gives the total memory size, including the size of additional memory allocated inside the class.

Definition at line 19 of file TFCSFunction.h.

19{ return sizeof(*this); };

◆ msg() [1/2]

MsgStream & ISF_FCS::MLogging::msg ( ) const

inlineinherited

Return a stream for sending messages directly (no decoration)

Definition at line 231 of file MLogging.h.

                                      {
  MsgStream *ms = m_msg_tls.get();
  if (!ms) {
    ms = new MsgStream(Athena::getMessageSvc(), m_nm);
    m_msg_tls.reset(ms);
  }
  return *ms;
}

◆ msg() [2/2]

MsgStream & ISF_FCS::MLogging::msg ( const MSG::Level lvl ) const

inlineinherited

Return a decorated starting stream for sending messages.

Definition at line 240 of file MLogging.h.

                                                        {
  return msg() << lvl;
}

◆ msgLvl()

bool ISF_FCS::MLogging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Check whether the logging system is active at the provided verbosity level.

Definition at line 222 of file MLogging.h.

                                                     {
  if (msg().level() <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

◆ ndim()

virtual int TFCS2DFunction::ndim ( ) const

inlinevirtualinherited

Return the number of dimensions for the function.

Reimplemented from TFCSFunction.

Definition at line 20 of file TFCS2DFunction.h.

20{ return 2; };

◆ rnd_to_fct() [1/2]

void TFCS2DFunctionHistogram::rnd_to_fct	(	float &	valuex,
		float &	valuey,
		float	rnd0,
		float	rnd1 ) const

virtual

Implements TFCS2DFunction.

Definition at line 69 of file TFCS2DFunctionHistogram.cxx.

                                                                       {
  if (m_HistoContents.empty()) {
    valuex = 0;
    valuey = 0;
    return;
  }
  auto it =
      std::upper_bound(m_HistoContents.begin(), m_HistoContents.end(), rnd0);
  int ibin = std::distance(m_HistoContents.begin(), it);
  if (ibin >= (int)m_HistoContents.size())
    ibin = m_HistoContents.size() - 1;
  Int_t nbinsx = m_HistoBorders.size() - 1;
  Int_t biny = ibin / nbinsx;
  Int_t binx = ibin - nbinsx * biny;
 
  float basecont = 0;
  if (ibin > 0)
    basecont = m_HistoContents[ibin - 1];
 
  float dcont = m_HistoContents[ibin] - basecont;
  if (dcont > 0) {
    valuex = m_HistoBorders[binx] +
             (m_HistoBorders[binx + 1] - m_HistoBorders[binx]) *
                 (rnd0 - basecont) / dcont;
  } else {
    valuex = m_HistoBorders[binx] +
             (m_HistoBorders[binx + 1] - m_HistoBorders[binx]) / 2;
  }
  valuey = m_HistoBordersy[biny] +
           (m_HistoBordersy[biny + 1] - m_HistoBordersy[biny]) * rnd1;
}

◆ rnd_to_fct() [2/2]

void TFCS2DFunction::rnd_to_fct	(	float	value[],
		const float	rnd[] ) const

virtual

Function gets array of random numbers rnd[] in the range [0,1) as arguments and returns function value in array value.

For a n-dimensional function, value and rnd should both have n elements.

Reimplemented from TFCS2DFunction.

Definition at line 24 of file TFCS2DFunction.cxx.

                                                                      {
  rnd_to_fct(value[0], value[1], rnd[0], rnd[1]);
}

◆ setLevel()

void ISF_FCS::MLogging::setLevel ( MSG::Level lvl )

virtualinherited

Update outputlevel.

Definition at line 105 of file MLogging.cxx.

                                    {
  lvl = (lvl >= MSG::NUM_LEVELS) ? MSG::ALWAYS
        : (lvl < MSG::NIL)       ? MSG::NIL
                                 : lvl;
  msg().setLevel(lvl);
}

◆ startMsg()

std::string ISF_FCS::MLogging::startMsg	(	MSG::Level	lvl,
		const std::string &	file,
		int	line )

staticinherited

Make a message to decorate the start of logging.

Print a message for the start of logging.

Definition at line 116 of file MLogging.cxx.

                                                                          {
  int col1_len = 20;
  int col2_len = 5;
  int col3_len = 10;
  auto last_slash = file.find_last_of('/');
  int path_len = last_slash == std::string::npos ? 0 : last_slash;
  int trim_point = path_len;
  int total_len = file.length();
  if (total_len - path_len > col1_len)
    trim_point = total_len - col1_len;
  std::string trimmed_name = file.substr(trim_point);
  const char *LevelNames[MSG::NUM_LEVELS] = {
      "NIL", "VERBOSE", "DEBUG", "INFO", "WARNING", "ERROR", "FATAL", "ALWAYS"};
  std::string level = LevelNames[lvl];
  std::string level_string = std::string("(") + level + ") ";
  std::stringstream output;
  output << std::setw(col1_len) << std::right << trimmed_name << ":"
         << std::setw(col2_len) << std::left << line << std::setw(col3_len)
         << std::right << level_string;
  return output.str();
}

Member Data Documentation

◆ ATLAS_THREAD_SAFE

boost::thread_specific_ptr<MsgStream> m_msg_tls ISF_FCS::MLogging::ATLAS_THREAD_SAFE

inlinestaticprivateinherited

Do not persistify!

MsgStream instance (a std::cout like with print-out levels)

Definition at line 215 of file MLogging.h.

◆ m_HistoBorders

std::vector<float> TFCS2DFunctionHistogram::m_HistoBorders

protected

Definition at line 45 of file TFCS2DFunctionHistogram.h.

◆ m_HistoBordersy

std::vector<float> TFCS2DFunctionHistogram::m_HistoBordersy

protected

Definition at line 46 of file TFCS2DFunctionHistogram.h.

◆ m_HistoContents

std::vector<float> TFCS2DFunctionHistogram::m_HistoContents

protected

Definition at line 47 of file TFCS2DFunctionHistogram.h.

◆ m_nm

std::string ISF_FCS::MLogging::m_nm

privateinherited

Message source name.

Definition at line 211 of file MLogging.h.

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

Public Member Functions

Static Public Member Functions

Protected Attributes

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TFCS2DFunctionHistogram()

◆ ~TFCS2DFunctionHistogram()

Member Function Documentation

◆ ATLAS_NOT_THREAD_SAFE() [1/3]

◆ ATLAS_NOT_THREAD_SAFE() [2/3]

◆ ATLAS_NOT_THREAD_SAFE() [3/3]

◆ CheckAndIntegrate2DHistogram()

◆ get_HistoBordersx() [1/2]

◆ get_HistoBordersx() [2/2]

◆ get_HistoBordersy() [1/2]

◆ get_HistoBordersy() [2/2]

◆ get_HistoContents() [1/2]

◆ get_HistoContents() [2/2]

◆ Initialize()

◆ level()

◆ MemorySize()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ ndim()

◆ rnd_to_fct() [1/2]

◆ rnd_to_fct() [2/2]

◆ setLevel()

◆ startMsg()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_HistoBorders

◆ m_HistoBordersy

◆ m_HistoContents

◆ m_nm