CaloSwEta2b_g3 Class Reference

#include <CaloSwEta2b_g3.h>

Inheritance diagram for CaloSwEta2b_g3:

Collaboration diagram for CaloSwEta2b_g3:

Public Member Functions
virtual void	makeCorrection (const Context &myctx, xAOD::CaloCluster *cluster) const override
virtual const std::string &	toolType () const override
virtual StatusCode	initialize () override
virtual void	setsample (xAOD::CaloCluster *cluster, CaloSampling::CaloSample sampling, float em, float etam, float phim, float emax, float etamax, float phimax, float etas, float phis) const
virtual void	setenergy (xAOD::CaloCluster *cluster, float energy) const
StatusCode	execute (const EventContext &ctx, xAOD::CaloCluster *cluster) const override
virtual StatusCode	execute (const EventContext &ctx, xAOD::CaloClusterContainer *collection) const
	Execute on an entire collection of clusters.

Protected Attributes
SG::ReadCondHandleKey< CaloDetDescrManager >	m_caloMgrKey {this,"CaloDetDescrManager", "CaloDetDescrManager"}

Private Member Functions
	CaloSwEta2b_g3 ()=delete

Private Attributes
Constant< CxxUtils::Array< 2 > >	m_correction { this, "correction", "" }
Constant< CxxUtils::Array< 2 > >	m_residuals { this, "residuals", "" }
Constant< CxxUtils::Array< 1 > >	m_residual_bins { this, "residual_bins", "" }
Constant< int >	m_correction_degree { this, "correction_degree", "" }
Constant< int >	m_residual_eval_degree { this, "residual_eval_degree", ""}
Constant< int >	m_residual_degree { this, "residual_degree", "" }
Constant< float >	m_correction_coef { this, "correction_coef", "" }
Constant< float >	m_residual_coef { this, "residual_coef", "" }
Constant< int >	m_region { this, "region", "" }

Static Private Attributes
static const float	s_middle_layer_granularity = 0.025

Detailed Description

Definition at line 40 of file CaloSwEta2b_g3.h.

Constructor & Destructor Documentation

CaloSwEta2b_g3()

CaloSwEta2b_g3::CaloSwEta2b_g3 ( )

privatedelete

Member Function Documentation

execute() [1/2]

StatusCode CaloClusterCorrection::execute ( const EventContext & ctx, xAOD::CaloCluster * cluster ) const

overrideinherited

Definition at line 52 of file CaloClusterCorrection.cxx.

{
  this->makeCorrection (context(ctx), cluster);
 
#if 0
  ATH_MSG_DEBUG( " ...... e, et " << cluster->e() << " " << cluster->et() << endmsg);
  ATH_MSG_DEBUG( " ...... eta, etaBE, etaSmp " << cluster->eta() << " " << cluster->etaBE(2) 
      << " " << cluster->etaSample(CaloSampling::EMB1) 
      << " " << cluster->etaSample(CaloSampling::EMB2) 
      << " " << cluster->etaSample(CaloSampling::EMB3) << endmsg);
  ATH_MSG_DEBUG( " ...... phi, phiBE, phiSmp " << cluster->phi() << " " << cluster->phiBE(2) 
      << " " << cluster->phiSample(CaloSampling::EMB1) 
      << " " << cluster->phiSample(CaloSampling::EMB2) 
      << " " << cluster->phiSample(CaloSampling::EMB3) << endmsg);
#endif
  
  return StatusCode::SUCCESS;
}

execute() [2/2]

StatusCode CaloClusterProcessor::execute ( const EventContext & ctx, xAOD::CaloClusterContainer * collection ) const

inherited

Execute on an entire collection of clusters.

Parameters

collection

The container of clusters.

This will iterate over all the clusters in collection and call execute on each one individually.

Parameters

collection	The container of clusters.
ctx	The event context.

This will iterate over all the clusters in collection and call execute on each one individually.

Definition at line 65 of file CaloClusterProcessor.cxx.

{
  for (xAOD::CaloCluster* clu : *collection) {
    ATH_CHECK( execute (ctx, clu) );
  }
  return StatusCode::SUCCESS;
}

initialize()

StatusCode CaloClusterCorrection::initialize ( )

overridevirtualinherited

Reimplemented in CaloClusterBadChannelList, CaloDummyCorrection, CaloFillRectangularCluster, CaloSwDeadOTX_back, CaloSwDeadOTX_ps, CaloSwGap_g3, CaloSwGap_v2, CaloSwGap_v3, CaloTopoEMGap, and CaloTopoEMlayers.

Definition at line 46 of file CaloClusterCorrection.cxx.

                                             {
  ATH_CHECK(m_caloMgrKey.initialize());
  ATH_CHECK(base_class::initialize());
  return StatusCode::SUCCESS;
}

makeCorrection()

void CaloSwEta2b_g3::makeCorrection ( const Context & myctx, xAOD::CaloCluster * cluster ) const

overridevirtual

Implements CaloClusterCorrection.

Definition at line 49 of file CaloSwEta2b_g3.cxx.

{
  // Only for barrel
  if (!cluster->inBarrel())
    return;
 
  const CxxUtils::Array<2> correction = m_correction (myctx);
  const CxxUtils::Array<2> residuals = m_residuals (myctx);
  const CxxUtils::Array<1> residual_bins = m_residual_bins (myctx);
  const int correction_degree = m_correction_degree(myctx);
  const int residual_eval_degree = m_residual_eval_degree(myctx);
  const int residual_degree = m_residual_degree(myctx);
  const float correction_coef = m_correction_coef(myctx);
  const float residual_coef = m_residual_coef(myctx);
 
  assert (residuals.size(1)-1 == residual_bins.size());
 
  // eta of second sampling
  float eta = cluster->etaSample(CaloSampling::EMB2); 
  if (eta == -999.) return;
  float aeta = fabs(eta);
 
  //   u2 is the distance to the inner edge of the cell (granularity 0.025)
  float u2 = fmod(aeta,s_middle_layer_granularity) ;
 
  //   First order (average) S-shape correction
  float deta = correction_coef * interpolate (correction,
                                              u2,
                                              correction_degree);
 
  //    Residual S-shape correction: six ranges in eta from 0 to 1.2
  //    Evaluated by 2nd order interpolation between 10 tabulated values
 
  // first we evaluate our residual in each eta range, given u2
 
  unsigned int shape[] = {residual_bins.size(), 2};
  CaloRec::WritableArrayData<2> w (shape);
  
  for (unsigned int i=0; i < residual_bins.size(); i++) {
    w[i][0] = residual_bins[i];
    w[i][1] = interpolate (residuals,
               u2, 
               residual_eval_degree,
               i+1);
  }
 
  // finally we interpolate (2nd order) for the actual eta
 
  deta += residual_coef * interpolate (w,
                                       aeta,
                                       residual_degree);
  if (eta < 0)
    deta = -deta;
 
  // make correction to second sampling
  cluster->setEta(CaloSampling::EMB2, eta - deta);
}

setenergy()

void CaloClusterCorrection::setenergy ( xAOD::CaloCluster * cluster, float energy ) const

virtualinherited

Definition at line 93 of file CaloClusterCorrection.cxx.

{
  if (cluster->e() == 0) {
    if (energy != 0)
      REPORT_MESSAGE (MSG::WARNING)
        << "Attempt to rescale zero-energy cluster to energy " << energy
        << " ignored.";
    return;
  }
 
  float correction = energy/cluster->e();
  cluster->setE(energy);
 
  // also correct individual sampling energies:
 
  for (int iSample=CaloSampling::PreSamplerB; 
        iSample < CaloSampling::Unknown; 
       ++iSample) 
    {
      CaloSampling::CaloSample sampling=static_cast<CaloSampling::CaloSample>(iSample);
      if (cluster->hasSampling (sampling)) {
        double e = cluster->eSample(sampling);
        cluster->setEnergy(sampling,e*correction) ;
      }
  }
}

setsample()

void CaloClusterCorrection::setsample ( xAOD::CaloCluster * cluster, CaloSampling::CaloSample sampling, float em, float etam, float phim, float emax, float etamax, float phimax, float etas, float phis ) const

virtualinherited

Definition at line 73 of file CaloClusterCorrection.cxx.

{
  cluster->setEnergy(sampling, em);
  cluster->setEta(sampling, etam);
  cluster->setPhi(sampling, phim);
 
  cluster->setEmax(sampling,emax);
  cluster->setEtamax(sampling,etamax);
  cluster->setPhimax(sampling,phimax);
 
  cluster->setEtasize(sampling, etas);
  cluster->setPhisize(sampling, phis);
}

toolType()

const std::string & CaloSwEta2b_g3::toolType ( ) const

overridevirtual

Definition at line 109 of file CaloSwEta2b_g3.cxx.

{
  static const std::string typeName = "CaloSwEta_g3";
  return typeName;
}

Member Data Documentation

m_caloMgrKey

SG::ReadCondHandleKey<CaloDetDescrManager> CaloClusterCorrection::m_caloMgrKey {this,"CaloDetDescrManager", "CaloDetDescrManager"}

protectedinherited

Definition at line 83 of file CaloClusterCorrection.h.

{this,"CaloDetDescrManager", "CaloDetDescrManager"};

m_correction

Constant<CxxUtils::Array<2> > CaloSwEta2b_g3::m_correction { this, "correction", "" }

private

Definition at line 59 of file CaloSwEta2b_g3.h.

{ this, "correction", "" };

m_correction_coef

Constant<float> CaloSwEta2b_g3::m_correction_coef { this, "correction_coef", "" }

private

Definition at line 65 of file CaloSwEta2b_g3.h.

{ this, "correction_coef", "" };

m_correction_degree

Constant<int> CaloSwEta2b_g3::m_correction_degree { this, "correction_degree", "" }

private

Definition at line 62 of file CaloSwEta2b_g3.h.

{ this, "correction_degree", "" };

m_region

Constant<int> CaloSwEta2b_g3::m_region { this, "region", "" }

private

Definition at line 67 of file CaloSwEta2b_g3.h.

{ this, "region",          "" };

m_residual_bins

Constant<CxxUtils::Array<1> > CaloSwEta2b_g3::m_residual_bins { this, "residual_bins", "" }

private

Definition at line 61 of file CaloSwEta2b_g3.h.

{ this, "residual_bins", "" };

m_residual_coef

Constant<float> CaloSwEta2b_g3::m_residual_coef { this, "residual_coef", "" }

private

Definition at line 66 of file CaloSwEta2b_g3.h.

{ this, "residual_coef",   "" };

m_residual_degree

Constant<int> CaloSwEta2b_g3::m_residual_degree { this, "residual_degree", "" }

private

Definition at line 64 of file CaloSwEta2b_g3.h.

{ this, "residual_degree", "" };

m_residual_eval_degree

Constant<int> CaloSwEta2b_g3::m_residual_eval_degree { this, "residual_eval_degree", ""}

private

Definition at line 63 of file CaloSwEta2b_g3.h.

{ this, "residual_eval_degree", ""};

m_residuals

Constant<CxxUtils::Array<2> > CaloSwEta2b_g3::m_residuals { this, "residuals", "" }

private

Definition at line 60 of file CaloSwEta2b_g3.h.

{ this, "residuals", "" };

s_middle_layer_granularity

const float CaloSwEta2b_g3::s_middle_layer_granularity = 0.025

staticprivate

Definition at line 69 of file CaloSwEta2b_g3.h.

---

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

• CaloSwEta2b_g3.h
• CaloSwEta2b_g3.cxx