HBHERecalibration Class Reference

#include <HBHERecalibration.h>

Public Member Functions
float	getCorr (int ieta, int depth) const
	HBHERecalibration (float intlumi, float cutoff, std::string meanenergies)
int	maxDepth () const
void	setup (const std::vector< std::vector< int >> &m_segmentation, const HBHEDarkening *darkening)
	~HBHERecalibration ()

Private Member Functions
void	initialize ()

Private Attributes
std::vector< std::vector< float > >	corr_
float	cutoff_
const HBHEDarkening *	darkening_
std::vector< std::vector< int > >	dsegm_
int	ieta_shift_
float	intlumi_
int	max_depth_
std::vector< std::vector< float > >	meanenergies_

Detailed Description

Definition at line 14 of file HBHERecalibration.h.

Constructor & Destructor Documentation

HBHERecalibration::HBHERecalibration	(	float	intlumi,
		float	cutoff,
		std::string	meanenergies
	)

Definition at line 7 of file HBHERecalibration.cc.

    : intlumi_(intlumi),
      cutoff_(cutoff),
      ieta_shift_(0),
      max_depth_(0),
      meanenergies_(HBHEDarkening::readDoseMap(meanenergies)),
      darkening_(nullptr) {}

HBHERecalibration::~HBHERecalibration

(

)

Definition at line 15 of file HBHERecalibration.cc.

{}

Member Function Documentation

float HBHERecalibration::getCorr	(	int	ieta,
		int	depth
	)		const

Definition at line 40 of file HBHERecalibration.cc.

References funct::abs(), corr_, cutoff_, HLT_FULL_cff::depth, ieta_shift_, and findQualityFiles::size.

                                                          {
  ieta = abs(ieta);
  //shift ieta tower index to act as array index
  ieta -= ieta_shift_;

  //shift depth index to act as array index (depth = 0 - not used!)
  depth -= 1;

  //bounds check
  if (ieta < 0 or ieta >= int(corr_.size()))
    return 1.0;
  if (depth < 0 or depth >= int(corr_[ieta].size()))
    return 1.0;

  if (cutoff_ > 1 and corr_[ieta][depth] > cutoff_)
    return cutoff_;
  else
    return corr_[ieta][depth];
}

void HBHERecalibration::initialize ( )

private

Definition at line 60 of file HBHERecalibration.cc.

References corr_, darkening_, HBHEDarkening::degradation(), HLT_FULL_cff::depth, dsegm_, ieta_shift_, intlumi_, max_depth_, meanenergies_, min(), and findQualityFiles::size.

Referenced by setup().

                                   {
  std::vector<std::vector<float>> vtmp(dsegm_.size(), std::vector<float>(max_depth_, 0.0));
  auto dval =
      vtmp;  //conversion of meanenergies into depths-averaged values - denominator (including degradation for intlumi)
  auto nval = vtmp;  // conversion of meanenergies into depths-averaged values - numerator (no degradation)
  corr_ = vtmp;

  //converting energy values from layers into depths
  for (unsigned int ieta = 0; ieta < dsegm_.size(); ++ieta) {
    //fill sum(means(layer,0)) and sum(means(layer,lumi)) for each depth
    for (unsigned int ilay = 0; ilay < std::min(meanenergies_[ieta].size(), dsegm_[ieta].size()); ++ilay) {
      int depth = dsegm_[ieta][ilay] - 1;  // depth = 0 - not used!
      nval[ieta][depth] += meanenergies_[ieta][ilay];
      dval[ieta][depth] +=
          meanenergies_[ieta][ilay] *
          darkening_->degradation(intlumi_, ieta + ieta_shift_, ilay + 1);  //be careful of eta and layer numbering
    }

    //compute factors, w/ safety checks
    for (int depth = 0; depth < max_depth_; ++depth) {
      if (dval[ieta][depth] > 0)
        corr_[ieta][depth] = nval[ieta][depth] / dval[ieta][depth];
      else
        corr_[ieta][depth] = 1.0;

      if (corr_[ieta][depth] < 1.0)
        corr_[ieta][depth] = 1.0;
    }
  }
}

int HBHERecalibration::maxDepth ( ) const

inline

Definition at line 22 of file HBHERecalibration.h.

References max_depth_.

{ return max_depth_; }

void HBHERecalibration::setup	(	const std::vector< std::vector< int >> &	m_segmentation,
		const HBHEDarkening *	darkening
	)

Definition at line 17 of file HBHERecalibration.cc.

References darkening_, HLT_FULL_cff::depth, dsegm_, HBHEDarkening::get_ieta_shift(), ieta_shift_, initialize(), max_depth_, meanenergies_, and findQualityFiles::size.

Referenced by o2olib.O2OTool::execute().

                                                                                                             {
  darkening_ = darkening;
  ieta_shift_ = darkening_->get_ieta_shift();

  //infer eta bounds
  int min_ieta = ieta_shift_ - 1;
  int max_ieta = min_ieta + meanenergies_.size();
  dsegm_.reserve(max_ieta - min_ieta);
  for (int ieta = min_ieta; ieta < max_ieta; ++ieta) {
    dsegm_.push_back(m_segmentation[ieta]);
    //find maximum
    for (unsigned lay = 0; lay < dsegm_.back().size(); ++lay) {
      if (lay >= meanenergies_[0].size())
        break;
      int depth = dsegm_.back()[lay];
      if (depth > max_depth_)
        max_depth_ = depth;
    }
  }

  initialize();
}

Member Data Documentation

std::vector<std::vector<float> > HBHERecalibration::corr_

private

Definition at line 36 of file HBHERecalibration.h.

Referenced by getCorr(), and initialize().

float HBHERecalibration::cutoff_

private

Definition at line 30 of file HBHERecalibration.h.

Referenced by getCorr().

const HBHEDarkening* HBHERecalibration::darkening_

private

Definition at line 34 of file HBHERecalibration.h.

Referenced by initialize(), and setup().

std::vector<std::vector<int> > HBHERecalibration::dsegm_

private

Definition at line 35 of file HBHERecalibration.h.

Referenced by initialize(), and setup().

int HBHERecalibration::ieta_shift_

private

Definition at line 31 of file HBHERecalibration.h.

Referenced by getCorr(), initialize(), and setup().

float HBHERecalibration::intlumi_

private

Definition at line 29 of file HBHERecalibration.h.

Referenced by initialize().

int HBHERecalibration::max_depth_

private

Definition at line 32 of file HBHERecalibration.h.

Referenced by initialize(), maxDepth(), and setup().

std::vector<std::vector<float> > HBHERecalibration::meanenergies_

private

Definition at line 33 of file HBHERecalibration.h.

Referenced by initialize(), and setup().