CaloLayerCalculator::Helper Struct Reference

Helper class for layer variable calculation. More...

Collaboration diagram for CaloLayerCalculator::Helper:

Public Member Functions
	Helper (CaloLayerCalculator &calc, double eta, double phi, double deta, double dphi, CaloSampling::CaloSample sampling, xAOD::CaloCluster *tofill, const CaloCellContainer *cell_container)
	Initialize for doing layer variable calculation.
void	cell (const CaloCell *cell, double weight)
	Accumulate results for another cell.
	~Helper ()
	Calculate the layer variables from the accumulated cells The variables in the parent CaloLayerCalculator instance (passed to the constructor) are updated.

Public Attributes
CaloLayerCalculator &	m_calc
	These members hold the values passed to the constructor.
double	m_eta
double	m_phi
double	m_deta
double	m_dphi
CaloSampling::CaloSample	m_sampling
xAOD::CaloCluster *	m_tofill
const CaloCellContainer *	m_cell_container
double	m_s00
	Temporary variables used to accumulate the results.
double	m_s10
double	m_s20
double	m_s01
double	m_s02
double	m_s10r
double	m_s01r

Detailed Description

Helper class for layer variable calculation.

This class is used to allow factoring out the layer variable calculation itself (which we only want to write once) from the loop over cells (of which there may be multiple varieties). To compute the layer variables, a Helper instance is created; it will contain the temporary variables used to accumulate the results. Then, for each cell in the list, call the cell method of Helper. At the end, the Helper destructor will update the layer variables.

Definition at line 285 of file CaloLayerCalculator.h.

Constructor & Destructor Documentation

Helper()

CaloLayerCalculator::Helper::Helper	(	CaloLayerCalculator &	calc,
		double	eta,
		double	phi,
		double	deta,
		double	dphi,
		CaloSampling::CaloSample	sampling,
		xAOD::CaloCluster *	tofill,
		const CaloCellContainer *	cell_container )

Initialize for doing layer variable calculation.

Parameters

calc	The parent CaloLayerCalculator instance.
eta	The \(\eta\) center of the cluster, for cell selection.
phi	The \(\phi\) center of the cluster, for cell selection.
deta	The \(\eta\) full width for cell selection.
dphi	The \(\phi\) full width for cell selection.
sampling	The sampling for cell selection.
tofill	If non-null, then selected cells will be added to this cluster.
cell_container	If tofill is non-null, then this is this container from which the cells came.

Definition at line 154 of file CaloLayerCalculator.cxx.

  : m_calc (calc),
    m_eta (eta),
    m_phi (CaloPhiRange::fix (phi)),
    m_deta (deta),
    m_dphi (dphi),
    m_sampling (sampling),
    m_tofill (tofill),
    m_cell_container (cell_container),
    m_s00 (0),
    m_s10 (0),
    m_s20 (0),
    m_s01 (0),
    m_s02 (0),
    // Added to cluster in the calorimeter frame
    m_s10r (0),
    m_s01r (0)
{
  m_calc.m_emax   = -999999;
}

~Helper()

CaloLayerCalculator::Helper::~Helper

(

)

Calculate the layer variables from the accumulated cells The variables in the parent CaloLayerCalculator instance (passed to the constructor) are updated.

Definition at line 244 of file CaloLayerCalculator.cxx.

{
  m_calc.m_em   = m_s00;
 
  if (m_s00 > 0.) {
    // When we computed m_s01 and m_s02 above, we had subtracted phi from each
    // angle to avoid wrapping problems.  We need to add it back in
    // when computing the mean, but not when we compute the width.
    // To be pedantic (but i got this wrong the first time, so maybe
    // someone else will also be confused):
    //  We want <phic^2>-<phic>^2.
    //    s2 = s02/s00 = <(phic-phi)^2>
    //    s1 = s01/s00 = <(phic-phi)>
    //  Expanding the square in s2 gives
    //    <phic^2>-<phic>^2 = s2 - <phic>^2 + phi*(2<phic>-phi)
    //  Expanding the square in s1^2 and substituting for <phic>^2 above gives
    //    <phic^2>-<phic>^2 = s2 - s1^2
    double s1 = m_s01/m_s00;
    m_calc.m_phim = CaloPhiRange::fix (s1 + m_phi);
    m_calc.m_phis = m_s02/m_s00 - s1*s1;
 
 
    m_calc.m_etam = m_s10/m_s00;
    m_calc.m_etas = m_s20/m_s00 - m_calc.m_etam*m_calc.m_etam;
 
    // In the calorimeter frame
    m_calc.m_etamr = m_s10r/m_s00;
    m_calc.m_phimr = CaloPhiRange::fix(m_s01r/m_s00 + m_phi);
 
    // Reset position variables if the calculated eta,phi is outside
    // of the window defined by deta,dphi. This usually happens
    // if the negative energy within the window nearly cancels
    // the positive energy.
    if (std::abs (m_calc.m_etamr - m_eta) > m_deta ||
        std::abs (m_calc.m_etam) > 10 ||
        std::abs (CaloPhiRange::fix (m_calc.m_phimr - m_phi)) > m_dphi)
    {
      m_calc.resetOnNegativeEnergy(m_eta,m_phi);
    }
  }
  else {
    //reset();//replaced by:
    m_calc.resetOnNegativeEnergy(m_eta,m_phi);
 
    //energy of a sampling CAN be negative (especially PS and back)
    //when there is the (gaussian) noise
    //TEMPORARY SOLUTION ( energy negative is not reset to 0 anymore,
    //but m_etam, m_phim, m_etas, m_phis are still null !! )
    // =>  << TO BE REDONE >>
  }
 
  if (m_calc.m_etas > 0.) m_calc.m_etas = sqrt(m_calc.m_etas);
  if (m_calc.m_phis > 0.) m_calc.m_phis = sqrt (m_calc.m_phis);
}

Member Function Documentation

cell()

void CaloLayerCalculator::Helper::cell	(	const CaloCell *	cell,
		double	weight )

Accumulate results for another cell.

Parameters

cell	The cell to accumulate.
weight	The weight for the cell.

Definition at line 190 of file CaloLayerCalculator.cxx.

{
  // Test to see if this cell should be selected.
  if (cell->caloDDE()->getSampling() == m_sampling) {
    double etac = cell->eta();
    double phic = cell->phi();
    double etar = cell->caloDDE()->eta_raw();
    double phir = cell->caloDDE()->phi_raw();
    // Ideally, phi is computed as a deltaphi + phi to be safe with range
    // the offset is taken to be the seed from the previous sampling
    // I take the same offset for the raw and real computations
    double dphic = CaloPhiRange::diff (phic, m_phi);
    double dphir = CaloPhiRange::diff (phir, m_phi);
 
    // The conditions are applied in the calorimeter frame, the biggest difference w.r.t. before...
    if (etar >= m_eta-m_deta/2. && etar <= m_eta+m_deta/2.) {
      if (dphir >= -m_dphi/2. && dphir <= m_dphi/2.) {
 
        // Is this the maximum-energy cell so far?
        double ene = cell->energy() * weight;
 
        if (ene > m_calc.m_emax) {
          m_calc.m_etamax = etac;
          m_calc.m_phimax = phic;
          m_calc.m_emax   = ene;
      // In the calorimeter frame
          m_calc.m_etarmax = etar;
          m_calc.m_phirmax = phir;
        }
 
        // Update accumulation.
        m_s00 += ene;
        m_s10 += etac*ene;
        m_s20 += ene*etac*etac;
        m_s01 += dphic*ene;
        m_s02 += ene*dphic*dphic;
    // In the calorimeter frame
    m_s10r += etar*ene;
    m_s01r += dphir*ene;
 
        // Add the cell to the cluster, if requested.
        if (m_tofill)
          m_tofill->addCell(m_cell_container->findIndex(cell->caloDDE()->calo_hash()), weight);
      }
    }
  }
}

Member Data Documentation

m_calc

CaloLayerCalculator& CaloLayerCalculator::Helper::m_calc

These members hold the values passed to the constructor.

Definition at line 327 of file CaloLayerCalculator.h.

m_cell_container

const CaloCellContainer* CaloLayerCalculator::Helper::m_cell_container

Definition at line 334 of file CaloLayerCalculator.h.

m_deta

double CaloLayerCalculator::Helper::m_deta

Definition at line 330 of file CaloLayerCalculator.h.

m_dphi

double CaloLayerCalculator::Helper::m_dphi

Definition at line 331 of file CaloLayerCalculator.h.

m_eta

double CaloLayerCalculator::Helper::m_eta

Definition at line 328 of file CaloLayerCalculator.h.

m_phi

double CaloLayerCalculator::Helper::m_phi

Definition at line 329 of file CaloLayerCalculator.h.

m_s00

double CaloLayerCalculator::Helper::m_s00

Temporary variables used to accumulate the results.

Definition at line 337 of file CaloLayerCalculator.h.

m_s01

double CaloLayerCalculator::Helper::m_s01

Definition at line 340 of file CaloLayerCalculator.h.

m_s01r

double CaloLayerCalculator::Helper::m_s01r

Definition at line 345 of file CaloLayerCalculator.h.

m_s02

double CaloLayerCalculator::Helper::m_s02

Definition at line 341 of file CaloLayerCalculator.h.

m_s10

double CaloLayerCalculator::Helper::m_s10

Definition at line 338 of file CaloLayerCalculator.h.

m_s10r

double CaloLayerCalculator::Helper::m_s10r

Definition at line 344 of file CaloLayerCalculator.h.

m_s20

double CaloLayerCalculator::Helper::m_s20

Definition at line 339 of file CaloLayerCalculator.h.

m_sampling

CaloSampling::CaloSample CaloLayerCalculator::Helper::m_sampling

Definition at line 332 of file CaloLayerCalculator.h.

m_tofill

xAOD::CaloCluster* CaloLayerCalculator::Helper::m_tofill

Definition at line 333 of file CaloLayerCalculator.h.

---

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

• CaloLayerCalculator.h
• CaloLayerCalculator.cxx