xAOD::TrigHisto2D_v1 Class Reference

#include <TrigHisto2D_v1.h>

Inheritance diagram for xAOD::TrigHisto2D_v1:

Collaboration diagram for xAOD::TrigHisto2D_v1:

Public Member Functions
	TrigHisto2D_v1 ()
	~TrigHisto2D_v1 ()
	Destructor.
	TrigHisto2D_v1 (unsigned int nbins_x, float min_x, float max_x, unsigned int nbins_y, float min_y, float max_y)
	Standard constructor used by FEX algorithms.
const std::vector< float > &	contents () const
	Fill a 2D histogram.
unsigned int	nbinsX () const
float	minX () const
	Return the minimum along the x-axis.
float	maxX () const
	Return the maximum along the x-axis.
unsigned int	nbinsY () const
	Return the number of bins along the y-axis, not including the under and overflow.
float	minY () const
	Return the minimum along the y-axis.
float	maxY () const
	Return the maximum along the y-axis.
void	initialize (unsigned int nbins_x, float min_x, float max_x, unsigned int nbins_y, float min_y, float max_y)
	creates empty histogram
void	fill (float value_x, float value_y, float weight)
	fill histogram
std::vector< float >	profileX () const
	Sum the number of entries within the cut range.
std::vector< float >	profileY () const
	Collapse the x-axis and return a profile from the y-axis.
double	sumEntries (float value_x, float value_y, int cutType) const
	Sum the number of entries within the cut range.
void	clear ()
	clear m_contents vector
void	dump ()
	dump() function, for testing
void	setContents (const std::vector< float > &cont)
	contents setter
void	setMinX (float val)
void	setMaxX (float val)
void	setNbinsX (unsigned int nx)
	NbinsX setter.
void	setMinY (float val)
void	setMaxY (float val)
void	setNbinsY (unsigned int ny)
	NbinsY setter.

Private Member Functions
unsigned int	findBinX (float val) const
	returns x bin index
unsigned int	findBinY (float val) const
	returns y bin index

Private Attributes
float	m_binWidthX = 0.0F
float	m_binWidthY = 0.0F

Static Private Attributes
static const int	EXTRA_BINS =2
	additional bins for underflow and overflow bins

Detailed Description

Definition at line 23 of file TrigHisto2D_v1.h.

Constructor & Destructor Documentation

TrigHisto2D_v1() [1/2]

xAOD::TrigHisto2D_v1::TrigHisto2D_v1

(

)

Definition at line 16 of file TrigHisto2D_v1.cxx.

  : SG::AuxElement(),m_binWidthX(0),m_binWidthY(0) {
  }

~TrigHisto2D_v1()

xAOD::TrigHisto2D_v1::~TrigHisto2D_v1

(

)

Destructor.

Definition at line 25 of file TrigHisto2D_v1.cxx.

{}

TrigHisto2D_v1() [2/2]

xAOD::TrigHisto2D_v1::TrigHisto2D_v1	(	unsigned int	nbins_x,
		float	min_x,
		float	max_x,
		unsigned int	nbins_y,
		float	min_y,
		float	max_y )

Standard constructor used by FEX algorithms.

Definition at line 20 of file TrigHisto2D_v1.cxx.

  : SG::AuxElement() {
    initialize(nbins_x, min_x, max_x, nbins_y, min_y, max_y);
  }

Member Function Documentation

clear()

void xAOD::TrigHisto2D_v1::clear

(

)

clear m_contents vector

Definition at line 207 of file TrigHisto2D_v1.cxx.

                             {
    static const Accessor< std::vector<float> > acc_contents( "contents" );
    
    for(std::vector<float>::iterator contents_iter = acc_contents(*this).begin(); contents_iter !=acc_contents(*this).end(); ++contents_iter)
      *contents_iter = 0;
  }

contents()

const std::vector< float > & xAOD::TrigHisto2D_v1::contents

(

)

const

Fill a 2D histogram.

Return the bin contents of the histogram, including the under and overflow bins.

dump()

void xAOD::TrigHisto2D_v1::dump

(

)

dump() function, for testing

Definition at line 214 of file TrigHisto2D_v1.cxx.

                           {
    static const Accessor< std::vector<float> > acc_contents( "contents" );
    
    std::cout<<"Dump contets vector of size:: "<<acc_contents(*this).size()<<std::endl;
    for( unsigned int i = 0 ; i < acc_contents(*this).size(); i++)
      std::cout<<acc_contents(*this).at(i)<<" ";
    std::cout<<std::endl;
    
    std::cout<<"NbinX:: "<< nbinsX()<<"\tRangeX:: ["<<minX()<<","<<maxX()<<"]"<<std::endl;
    std::cout<<"NbinY:: "<< nbinsY()<<"\tRangeY:: ["<<minY()<<","<<maxY()<<"]"<<std::endl;
    
    
    for(unsigned int iy = 0 ; iy < nbinsY() + TrigHisto2D_v1::EXTRA_BINS; iy++){
      std::cout<<std::endl;
      for(unsigned int ix = 0; ix < nbinsX() + TrigHisto2D_v1::EXTRA_BINS; ix++){
    std::cout<< acc_contents(*this).at(ix + iy*(nbinsX() + TrigHisto2D_v1::EXTRA_BINS))<<"\t";
      }
    }
    std::cout<<std::endl;
  }

fill()

void xAOD::TrigHisto2D_v1::fill	(	float	value_x,
		float	value_y,
		float	weight = 1. )

fill histogram

Definition at line 73 of file TrigHisto2D_v1.cxx.

                                                                      {
    
    unsigned int ibin_x = findBinX(value_x);
    unsigned int ibin_y = findBinY(value_y);
 
    unsigned int ibin = ibin_y*(nbinsX()+TrigHisto2D_v1::EXTRA_BINS) + ibin_x;
   
    static const Accessor< std::vector<float> > acc_contents( "contents" );
 
    acc_contents(*this).at(ibin)+= weight;
  }

findBinX()

unsigned int xAOD::TrigHisto2D_v1::findBinX ( float val ) const

private

returns x bin index

Definition at line 171 of file TrigHisto2D_v1.cxx.

                                                         {
    
    unsigned int ibin = 0;
    
    if(value < minX()) { // Underflow
      ibin = 0;
    }
    else if( !(value < maxX()) ) { // Overflow (catches NaN)
      ibin = nbinsX()+1;
    }
    else {
      while(value > (ibin*m_binWidthX+minX()) && ibin <= nbinsX()) { // None under/overflow from 1 to nbins
    ibin++;
      }
    }
    return ibin;
  }

findBinY()

unsigned int xAOD::TrigHisto2D_v1::findBinY ( float val ) const

private

returns y bin index

Definition at line 189 of file TrigHisto2D_v1.cxx.

                                                        {
    
    unsigned int ibin = 0;
    
    if(value < minY()) { // Underflow
      ibin = 0;
    }
    else if( !(value < maxY()) ) { // Overflow (catches NaN)
      ibin = nbinsY()+1;
    }
    else {
      while(value > (ibin*m_binWidthY+minY()) && ibin <= nbinsY()) { // None under/overflow from 1 to nbins
    ibin++;
      }
    }
    return ibin;
  }

initialize()

void xAOD::TrigHisto2D_v1::initialize	(	unsigned int	nbins_x,
		float	min_x,
		float	max_x,
		unsigned int	nbins_y,
		float	min_y,
		float	max_y )

creates empty histogram

compute contents_size

create empty contents vector

compute binwdiths'

Definition at line 37 of file TrigHisto2D_v1.cxx.

                                                                                                                               {
    
    //the aux store has to exist if we wont to write anything there
    if(!this->hasStore())
      this->makePrivateStore();
    
    unsigned int contents_size = (nbins_x+TrigHisto2D_v1::EXTRA_BINS)*(nbins_y+TrigHisto2D_v1::EXTRA_BINS);
#ifdef DBG
    std::cout<<"Creating Thist2D with :"<<contents_size<< " elements"<<std::endl<<"nbins_x+2="<< nbins_x+TrigHisto2D_v1::EXTRA_BINS<<std::endl<<"nbins_y+2="<<nbins_y+TrigHisto2D_v1::EXTRA_BINS<<std::endl;
#endif
    setContents(std::vector<float>(contents_size,0.));
    setMinX(min_x);
    setMaxX(max_x);
    setNbinsX(nbins_x);
    setMinY(min_y);
    setMaxY(max_y);
    setNbinsY(nbins_y);
    
    if(nbins_x != 0) {
      m_binWidthX = (max_x - min_x)/((float)nbins_x); // Calculate bin size.
    }
    else {
      m_binWidthX = 0.;
    }
    
    if(nbins_y != 0) {
      m_binWidthY = (max_y - min_y)/((float)nbins_y); // Calculate bin size.
    }
    else {
      m_binWidthY = 0.;
    }   
  }

maxX()

float xAOD::TrigHisto2D_v1::maxX

(

)

const

Return the maximum along the x-axis.

maxY()

float xAOD::TrigHisto2D_v1::maxY

(

)

const

Return the maximum along the y-axis.

minX()

float xAOD::TrigHisto2D_v1::minX

(

)

const

Return the minimum along the x-axis.

minY()

float xAOD::TrigHisto2D_v1::minY

(

)

const

Return the minimum along the y-axis.

nbinsX()

unsigned int xAOD::TrigHisto2D_v1::nbinsX

(

)

const

nbinsY()

unsigned int xAOD::TrigHisto2D_v1::nbinsY

(

)

const

Return the number of bins along the y-axis, not including the under and overflow.

profileX()

std::vector< float > xAOD::TrigHisto2D_v1::profileX

(

)

const

Sum the number of entries within the cut range.

Collapse the y-axis and return a profile from the x-axis

Definition at line 85 of file TrigHisto2D_v1.cxx.

                                                  {
    std::vector<float> contentsX(nbinsX()+TrigHisto2D_v1::EXTRA_BINS,0.);
    
    static const Accessor< std::vector<float> > acc_contents( "contents" );
    
    for(unsigned int ix = 0 ; ix < nbinsX() + TrigHisto2D_v1::EXTRA_BINS; ix++){
      for(unsigned int iy = 0; iy < nbinsY() + TrigHisto2D_v1::EXTRA_BINS; iy++){
    contentsX.at(ix) += acc_contents(*this).at(ix + iy*(nbinsX() + TrigHisto2D_v1::EXTRA_BINS));
      }
    }
    return contentsX;
  }

profileY()

std::vector< float > xAOD::TrigHisto2D_v1::profileY

(

)

const

Collapse the x-axis and return a profile from the y-axis.

Definition at line 98 of file TrigHisto2D_v1.cxx.

                                                  {
    std::vector<float> contentsY(nbinsY()+TrigHisto2D_v1::EXTRA_BINS);
    
    static const Accessor< std::vector<float> > acc_contents( "contents" );
    
    for(unsigned int iy = 0 ; iy < nbinsY() + TrigHisto2D_v1::EXTRA_BINS; iy++){
      for(unsigned int ix = 0; ix < nbinsX() + TrigHisto2D_v1::EXTRA_BINS; ix++){
    contentsY.at(iy) += acc_contents(*this).at(ix + iy*(nbinsX() + TrigHisto2D_v1::EXTRA_BINS));
      }
    }
    return contentsY;
  }

setContents()

void xAOD::TrigHisto2D_v1::setContents

(

const std::vector< float > &

cont

)

contents setter

setMaxX()

void xAOD::TrigHisto2D_v1::setMaxX

(

float

val

)

setMaxY()

void xAOD::TrigHisto2D_v1::setMaxY

(

float

val

)

setMinX()

void xAOD::TrigHisto2D_v1::setMinX

(

float

val

)

setMinY()

void xAOD::TrigHisto2D_v1::setMinY

(

float

val

)

setNbinsX()

void xAOD::TrigHisto2D_v1::setNbinsX

(

unsigned int

nx

)

NbinsX setter.

setNbinsY()

void xAOD::TrigHisto2D_v1::setNbinsY

(

unsigned int

ny

)

NbinsY setter.

sumEntries()

double xAOD::TrigHisto2D_v1::sumEntries	(	float	value_x,
		float	value_y,
		int	cutType ) const

Sum the number of entries within the cut range.

Definition at line 111 of file TrigHisto2D_v1.cxx.

                                                                                   {
    
    unsigned int ibin, ibin_x, ibin_y, ibin_x_selected, ibin_y_selected;
    double entries;
    
    static const Accessor< std::vector<float> > acc_contents( "contents" );
    
    // Find the x bin index that the cut corresponds to.
    ibin_x_selected = findBinX(value_x);
    
    // Find the y bin index that the cut corresponds to.
    ibin_y_selected = findBinY(value_y);
    entries = 0.;
    
    
    if( nbinsX()==0 || nbinsY()==0 ){
      return 0;
    }
    else{
      
      if(cutType == TrigHistoCutType::BELOW_X_BELOW_Y) {
    for(ibin_x = 0; ibin_x <= ibin_x_selected; ibin_x++) {
      for(ibin_y = 0; ibin_y <= ibin_y_selected; ibin_y++) {
        ibin = ibin_y*(nbinsX()+2) + ibin_x; // Find position in 1d data array
        entries += acc_contents(*this).at(ibin);
      }
    }
      }
      else if(cutType == TrigHistoCutType::ABOVE_X_BELOW_Y) {
    for(ibin_x = ibin_x_selected; ibin_x < nbinsX()+ TrigHisto2D_v1::EXTRA_BINS; ibin_x++) {
      for(ibin_y = 0; ibin_y <= ibin_y_selected; ibin_y++) {
        ibin = ibin_y*(nbinsX()+ TrigHisto2D_v1::EXTRA_BINS) + ibin_x; // Find position in 1d data array
        entries += acc_contents(*this).at(ibin);
      }
    }
      }
      else if(cutType == TrigHistoCutType::BELOW_X_ABOVE_Y) {
    for(ibin_x = 0; ibin_x <= ibin_x_selected; ibin_x++) {
      for(ibin_y = ibin_y_selected; ibin_y < nbinsY()+ TrigHisto2D_v1::EXTRA_BINS; ibin_y++) {
        ibin = ibin_y*(nbinsX()+ TrigHisto2D_v1::EXTRA_BINS) + ibin_x; // Find position in 1d data array
        entries += acc_contents(*this).at(ibin);
      }
    }
      }
      else if(cutType == TrigHistoCutType::ABOVE_X_ABOVE_Y) {   
    for(ibin_x = ibin_x_selected; ibin_x < nbinsX()+ TrigHisto2D_v1::EXTRA_BINS; ibin_x++) {
      for(ibin_y = ibin_y_selected; ibin_y < nbinsY()+ TrigHisto2D_v1::EXTRA_BINS; ibin_y++) {
        ibin = ibin_y*(nbinsX()+ TrigHisto2D_v1::EXTRA_BINS) + ibin_x; // Find position in 1d data array
        entries += acc_contents(*this).at(ibin);
      }
    }
      }
      else {
    return 0;
      }
    }// else of m_nbins!=0
    
    return entries;
  }

Member Data Documentation

EXTRA_BINS

const int xAOD::TrigHisto2D_v1::EXTRA_BINS =2

staticprivate

additional bins for underflow and overflow bins

Definition at line 108 of file TrigHisto2D_v1.h.

m_binWidthX

float xAOD::TrigHisto2D_v1::m_binWidthX = 0.0F

private

Definition at line 111 of file TrigHisto2D_v1.h.

m_binWidthY

float xAOD::TrigHisto2D_v1::m_binWidthY = 0.0F

private

Definition at line 114 of file TrigHisto2D_v1.h.

---

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

• TrigHisto2D_v1.h
• TrigHisto2D_v1.cxx