EgammaSshapeCalibration::Builder Class Reference

Inheritance diagram for EgammaSshapeCalibration::Builder:

Collaboration diagram for EgammaSshapeCalibration::Builder:

Public Member Functions
	Builder (const EgammaSshapeCalibration &, double, double, int)
	Constructor.
virtual float	calculate (int energy_ndx, bool &good) const
	Calculate the correction for tabulated energy ENERGY_NDX.

Private Member Functions
float	calc0 (float, float, const CaloRec::Array< 2 > &) const
float	calc3 (float, float, const CaloRec::Array< 2 > &) const
float	calc4 (float, float, const CaloRec::Array< 2 > &) const
float	calc5 (float, float, const CaloRec::Array< 2 > &) const
float	calc10 (float, float, const CaloRec::Array< 2 > &) const
float	calc11 (float, float, const CaloRec::Array< 2 > &) const

Private Attributes
const EgammaSshapeCalibration &	m_corr
double	m_aeta
double	m_u
int	m_region_ndx
int	m_form

Detailed Description

Definition at line 50 of file EgammaSshapeCalibration.h.

Constructor & Destructor Documentation

Builder()

EgammaSshapeCalibration::Builder::Builder	(	const EgammaSshapeCalibration &	corr,
		double	aeta,
		double	u,
		int	region_ndx )

Constructor.

Gets the parent correction object, the abs(eta) at which the correction is being evaluated (in cal-local coordinates), the fractional cell offset u, and the index of the region for this correction.

Definition at line 234 of file EgammaSshapeCalibration.cxx.

  : m_corr(corr),
    m_aeta(aeta),
    m_u(u),
    m_region_ndx(region_ndx),
    m_form (m_corr.m_regions()[region_ndx][REG_FORM])
{ }

Member Function Documentation

calc0()

float EgammaSshapeCalibration::Builder::calc0 ( float aeta, float u, const CaloRec::Array< 2 > & coef ) const

private

Definition at line 301 of file EgammaSshapeCalibration.cxx.

{
  const unsigned int NPAR = 5;
  double pars[NPAR];
 
  unsigned int degree = coef.size(1);
  for (unsigned int i=0; i < NPAR; i++) {
    CaloRec::Array<1> xcoef = coef[i];      
      double out = xcoef[0];
      for (unsigned int j=1; j<degree; j++)
    out = out * aeta + xcoef[j];
      pars[i] = out;
  }
 
  double c = (u>=0) ? pars[2] : pars[3];
  return pars[0] * (std::atan (pars[1] * u) + c*u + pars[4]);
}

calc10()

float EgammaSshapeCalibration::Builder::calc10 ( float aeta, float u, const CaloRec::Array< 2 > & coef ) const

private

Definition at line 399 of file EgammaSshapeCalibration.cxx.

{
  const unsigned int NPAR = 7;
  double pars[NPAR];
  
  unsigned int degree = coef.size(1);
  for (unsigned int i=0; i < NPAR; i++) {
    CaloRec::Array<1> xcoef = coef[i];      
    double out = xcoef[0];
    for (unsigned int j=1; j<degree; j++)
      out = out * aeta + xcoef[j];
    pars[i] = out;
  }
 
  double c = pars[2];
  if (u < 0)
    c = -c - 2*std::atan(pars[1]);
 
  double offs = pars[0] * (atan (pars[1] * u) + c*u + pars[3]) +
    pars[4] * std::cos (u * pars[5] * M_PI + pars[6]);
  if (u > 0.9 || u < -0.9) {
    double xdiff = pars[4] * (std::cos (pars[5] * M_PI + pars[6]) -
                              std::cos (- pars[5] * M_PI + pars[6]));
    if (u > 0.9)
      offs -= (u-0.9)/0.1 * xdiff/2;
    else
      offs += (-0.9-u)/0.1 * xdiff/2;
  }
  return offs;
}

calc11()

float EgammaSshapeCalibration::Builder::calc11 ( float aeta, float u, const CaloRec::Array< 2 > & coef ) const

private

Definition at line 435 of file EgammaSshapeCalibration.cxx.

{
  const unsigned int NPAR = 8;
  double pars[NPAR];
 
 unsigned int degree = coef.size(1);
  for (unsigned int i=0; i < NPAR; i++) {
    CaloRec::Array<1> xcoef = coef[i];      
      double out = xcoef[0];
      for (unsigned int j=1; j<degree; j++)
    out = out * aeta + xcoef[j];
      pars[i] = out;
  }
 
  double c = (u>=0) ? pars[2] : pars[3];
  return pars[0] * (std::atan (pars[1] * u) + c*u + pars[4])
    + pars[5] * std::cos (u * pars[6] * M_PI + pars[7]);
}

calc3()

float EgammaSshapeCalibration::Builder::calc3 ( float aeta, float u, const CaloRec::Array< 2 > & coef ) const

private

Definition at line 324 of file EgammaSshapeCalibration.cxx.

{
  const unsigned int NPAR = 4;
  double pars[NPAR];
 
  unsigned int degree = coef.size(1);
  for (unsigned int i=0; i < NPAR; i++) {
    CaloRec::Array<1> xcoef = coef[i];      
    double out = xcoef[0];
    for (unsigned int j=1; j<degree; j++)
      out = out * aeta + xcoef[j];
    pars[i] = out;
  }
  
  double c = pars[2];
  if (u < 0)
    c = -c - 2*std::atan(pars[1]);
  return pars[0] * (std::atan (pars[1] * u) + c*u + pars[3]);
}

calc4()

float EgammaSshapeCalibration::Builder::calc4 ( float aeta, float u, const CaloRec::Array< 2 > & coef ) const

private

Definition at line 349 of file EgammaSshapeCalibration.cxx.

{
  const unsigned int NPAR = 3;
  double pars[NPAR];
 
  unsigned int degree = coef.size(1);
  for (unsigned int i=0; i < NPAR; i++) {
    CaloRec::Array<1> xcoef = coef[i];      
    double out = xcoef[0];
    for (unsigned int j=1; j<degree; j++)
      out = out * aeta + xcoef[j];
    pars[i] = out;
  }
  
  double b = std::max ((double)pars[1], 1e-5);
  double atanb = std::atan(b);
  double sq = std::sqrt (b/atanb - 1);
  double den = (sq/b*atanb - std::atan(sq));
  return  pars[0]* ((- std::atan (b*u) + u*atanb) / den +
                    pars[2]*(1-std::abs(u)));
}

calc5()

float EgammaSshapeCalibration::Builder::calc5 ( float aeta, float , const CaloRec::Array< 2 > & coef ) const

private

Definition at line 376 of file EgammaSshapeCalibration.cxx.

{
  const unsigned int NPAR = 1;
  double pars[NPAR];
 
  //  poly_eval (aeta, coef, NPAR, pars);
 unsigned int degree = coef.size(1);
  for (unsigned int i=0; i < NPAR; i++) {
    CaloRec::Array<1> xcoef = coef[i];      
      double out = xcoef[0];
      for (unsigned int j=1; j<degree; j++)
    out = out * aeta + xcoef[j];
      pars[i] = out;
  }
 
  return pars[0];
}

calculate()

float EgammaSshapeCalibration::Builder::calculate ( int energy_ndx, bool & good ) const

virtual

Calculate the correction for tabulated energy ENERGY_NDX.

Implements CaloClusterCorrectionCommon::TableBuilder.

Definition at line 248 of file EgammaSshapeCalibration.cxx.

                                            {
 
  // Find the proper array of coefficients.
  CaloRec::Array<2> coef = m_corr.m_correction()[energy_ndx][m_region_ndx];
  
  // If we don't have coefficients for this energy/region, skip it.
  if(coef[0].end()[-1]==0) {
    good = false;
    return 0;
  }
 
  // Which functional form to use?
  int form;
  CxxUtils::Array<2> forms = m_corr.m_forms();
  if (forms.size(0) != 0 && forms.size(1) != 0)
    form = forms[m_region_ndx][energy_ndx];
  else
    form = m_form;
 
 
 // Evaluate the correction!
  good = true;
  float ret = 0;
  switch (form) {
  case 0:
    ret = calc0 (m_aeta, m_u, coef);
    break;
  case 3:
    ret = calc3 (m_aeta, m_u, coef);
    break;
  case 4:
    ret = calc4 (m_aeta, m_u, coef);
    break;
  case 5:
    ret = calc5 (m_aeta, m_u, coef);
    break;
  case 10:
    ret = calc10 (m_aeta, m_u, coef);
    break;
  case 11:
    ret = calc11 (m_aeta, m_u, coef);
    break;
  default:
    abort();
  }
 
  return ret;
}

Member Data Documentation

m_aeta

double EgammaSshapeCalibration::Builder::m_aeta

private

Definition at line 72 of file EgammaSshapeCalibration.h.

m_corr

const EgammaSshapeCalibration& EgammaSshapeCalibration::Builder::m_corr

private

Definition at line 71 of file EgammaSshapeCalibration.h.

m_form

int EgammaSshapeCalibration::Builder::m_form

private

Definition at line 75 of file EgammaSshapeCalibration.h.

m_region_ndx

int EgammaSshapeCalibration::Builder::m_region_ndx

private

Definition at line 74 of file EgammaSshapeCalibration.h.

m_u

double EgammaSshapeCalibration::Builder::m_u

private

Definition at line 73 of file EgammaSshapeCalibration.h.

---

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

• EgammaSshapeCalibration.h
• EgammaSshapeCalibration.cxx