CaloClusterUpdate Class Reference

#include <CaloClusterUpdate.h>

Inheritance diagram for CaloClusterUpdate:

Collaboration diagram for CaloClusterUpdate:

Public Types
typedef ToolWithConstants	base_class
	Shorthand for derived classes.
using	Context
	Convenient alias for the Context type.
using	Constant
	Alias for the Constant type.

Public Member Functions
virtual void	makeCorrection (const Context &myctx, xAOD::CaloCluster *cluster) const override
virtual StatusCode	initialize () override
	Initialize method.
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
virtual void	makeCorrection (const Context &myctx, xAOD::CaloCluster *) const =0
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.
Context	context (const EventContext &ctx) const
	Create a Context object.
virtual void	writeConstants (std::ostream &stream, const std::string &name, const EventContext &ctx) const
	Dump method (for debugging)
virtual StatusCode	mergeConstants (CaloRec::ToolConstants &out, const EventContext &ctx) const override
	Merge our constants into out with the proper prefix.
virtual int	toolVersion () const
	Return the version number for this tool.
virtual const std::string &	toolType () const
	Return the name of the type of this tool.

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

Private Attributes
Gaudi::Property< float >	m_w1 { this, "layer1_weight", 3, "layer 1 additional weight in eta average" }
Constant< bool >	m_update_energy { this, "update_energy", "Should cluster total energy be updated?" }
SG::ReadCondHandleKey< CaloRec::ToolConstants >	m_DBHandle
	Handle to a ToolConstants conditions object.
StringProperty	m_prefix
	Prefix for finding our constants within the ToolConstants object.
ToolWithConstantsImpl	m_impl
	Internal implementation object.
Constant< int >	m_order
	Used to fix the ordering of tools when we're initializing from COOL based on a hierarchical tag.
Constant< bool >	m_isdummy
	If true, then this is a dummy tool that should not be executed.

Detailed Description

Definition at line 46 of file CaloClusterUpdate.h.

Member Typedef Documentation

base_class

typedef ToolWithConstants CaloUtils::ToolWithConstants< BASE >::base_class

inherited

Shorthand for derived classes.

Definition at line 450 of file ToolWithConstants.h.

Constant

using CaloUtils::ToolWithConstants< BASE >::Constant

inherited

Alias for the Constant type.

Definition at line 463 of file ToolWithConstants.h.

Context

using CaloUtils::ToolWithConstants< BASE >::Context

inherited

Convenient alias for the Context type.

Definition at line 458 of file ToolWithConstants.h.

Member Function Documentation

context()

Context CaloUtils::ToolWithConstants< BASE >::context ( const EventContext & ctx ) const

inherited

Create a Context object.

This can then be passed to Constant::operator().

execute() [1/2]

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

overrideinherited

Definition at line 53 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

Initialize method.

Derived classes must call this.

Reimplemented from CaloUtils::ToolWithConstants< CaloClusterProcessor >.

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

Definition at line 47 of file CaloClusterCorrection.cxx.

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

makeCorrection() [1/2]

virtual void CaloClusterCorrection::makeCorrection ( const Context & myctx, xAOD::CaloCluster *  ) const

pure virtualinherited

Implemented in CaloClusterBadChannelList, and CaloFillRectangularCluster.

makeCorrection() [2/2]

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

overridevirtual

Definition at line 48 of file CaloClusterUpdate.cxx.

{
  float energy=0; 
  float eta=0; 
  float weta=0;
 
  // set eta to be weighted average of eta1 and eta2
  ATH_MSG_DEBUG(" inBarrel "<<cluster->inBarrel()
        << " inEndcap "<<cluster->inEndcap()) ;
 
  for(int i=0; i<5; i=i+4 )
  { 
  
   if (i==0 && !cluster->inBarrel()) continue;
   if (i==4 && !cluster->inEndcap()) continue;
 
   CaloSampling::CaloSample sam0 = (CaloSampling::CaloSample)(CaloSampling::PreSamplerB+i);   
   CaloSampling::CaloSample sam1 = (CaloSampling::CaloSample)(CaloSampling::PreSamplerB+i+1);   
   CaloSampling::CaloSample sam2 = (CaloSampling::CaloSample)(CaloSampling::PreSamplerB+i+2);   
   CaloSampling::CaloSample sam3 = (CaloSampling::CaloSample)(CaloSampling::PreSamplerB+i+3);
   
   float e0 = cluster->eSample(sam0);
   float e1 = cluster->eSample(sam1);
   float e2 = cluster->eSample(sam2);
   float e3 = cluster->eSample(sam3);
 
   // total energy is the sum of each sampling, which had all corrections
   energy += (e0 + e1 + e2 + e3);
 
   // do not consider including in the average if both energies are negative
   if (e1 <= 0 && e2 <= 0) continue;
 
   // reject weird clusters. Not even sure it can ever happen
   float eta2 = cluster->etaSample(sam2);
   if (cluster->hasSampling(sam2) && eta2 == -999.) continue;
 
   float eta1 = cluster->etaSample(sam1);
   if (cluster->hasSampling(sam1) && eta1 == -999.) continue;
 
   // eta1 has better resolution, so weight it differently
   float w1 = e1*m_w1;
   float w2 = e2;
 
   // do not include layer if energy is negative
   if (e1 <= 0) {
     w1 = 0;
   }
   else if (e2 <= 0) {
     w2 = 0;
   }
 
   eta  += (eta1*w1+eta2*w2) ;
   weta += w1+w2; 
  }
 
  //
  // set them in the cluster. 
  //
  if (eta  != -999. && weta != 0)
    eta  = eta/weta; 
  else {
    ATH_MSG_DEBUG("Weird cluster "
          " EB1 = " << cluster->eSample(CaloSampling::EMB1)
          << " etaB1 = " << cluster->etaSample(CaloSampling::EMB1)
          << " EE1 = " << cluster->eSample(CaloSampling::EME1)
          << " etaE1 = " << cluster->etaSample(CaloSampling::EME1)
          << " EB2 = " << cluster->eSample(CaloSampling::EMB2)
          << " etaB2 = " << cluster->etaSample(CaloSampling::EMB2)
          << " EE2 = " << cluster->eSample(CaloSampling::EME2)
          << " etaE2 = " << cluster->etaSample(CaloSampling::EME2));
    if (cluster->inBarrel() && !cluster->inEndcap())
      eta = cluster->etaSample(CaloSampling::EMB2);
    else if (cluster->inEndcap() && !cluster->inBarrel())
      eta = cluster->etaSample(CaloSampling::EME2);
    else {
      if (cluster->eSample(CaloSampling::EMB2) >
      cluster->eSample(CaloSampling::EME2))
    eta = cluster->etaSample(CaloSampling::EMB2);
      else
    eta = cluster->etaSample(CaloSampling::EME2);
    }
  }
 
  cluster->setEta(eta);
  cluster->setPhi(cluster->phiBE(2));
 
  if (m_update_energy (myctx))
    cluster->setE(energy);
}

mergeConstants()

virtual StatusCode CaloUtils::ToolWithConstants< BASE >::mergeConstants ( CaloRec::ToolConstants & out, const EventContext & ctx ) const

overridevirtualinherited

Merge our constants into out with the proper prefix.

Parameters

[out]	out	Object to receive our constants.
	ctx	Event context.

setenergy()

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

virtualinherited

Definition at line 94 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 74 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()

virtual const std::string & CaloUtils::ToolWithConstants< BASE >::toolType ( ) const

virtualinherited

Return the name of the type of this tool.

A saved set of constants includes both the C++ class name and a version number. Normally, the class name is taken from the Gaudi type() method, but that may be changed by overriding this method. This can be used, for example, when there are tools with distinct C++ classes but which are yet similar enough to combine together.

Reimplemented in CaloSwEta1b_g3, CaloSwEta1e_g3, CaloSwEta2b_g3, and CaloSwEta2e_g3.

toolVersion()

virtual int CaloUtils::ToolWithConstants< BASE >::toolVersion ( ) const

virtualinherited

Return the version number for this tool.

A saved set of constants includes both the C++ class name and a version number. The idea is that the version number can be bumped whenever there's a backwards-incompatible change; this gives some protection against trying to use an old version of a tool with an incompatible newer set of constants.

If you want a tool to have a version number, override this method. Otherwise, it will default to a version number of 0.

writeConstants()

virtual void CaloUtils::ToolWithConstants< BASE >::writeConstants ( std::ostream & stream, const std::string & name, const EventContext & ctx ) const

virtualinherited

Dump method (for debugging)

Parameters

stream	Ostream to which to write.
name	Name to go in output
ctx	Event context.

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_DBHandle

SG::ReadCondHandleKey<CaloRec::ToolConstants> CaloUtils::ToolWithConstants< BASE >::m_DBHandle

privateinherited

Handle to a ToolConstants conditions object.

Definition at line 534 of file ToolWithConstants.h.

{ this, "DBHandleKey", "", "" };

m_impl

ToolWithConstantsImpl CaloUtils::ToolWithConstants< BASE >::m_impl

privateinherited

Internal implementation object.

Definition at line 542 of file ToolWithConstants.h.

{ this->name(), m_prefix, m_DBHandle };

m_isdummy

Constant<bool> CaloUtils::ToolWithConstants< BASE >::m_isdummy

privateinherited

If true, then this is a dummy tool that should not be executed.

This is used for the case of reading from COOL using hierarchical tags: we need to have such tags associated with some object in each folder, regardless of whether or not the correction from that folder is actually used. [Every folder that IOVDbSvc knows about at configuration time needs to have a valid object for the configured tag, else IOVDbSvc will raise a fatal error. But we don't know at configuration time which folders we're actually going to need, so we gotta configure all of them.]

Definition at line 560 of file ToolWithConstants.h.

{this, "isDummy", false };

m_order

Constant<int> CaloUtils::ToolWithConstants< BASE >::m_order

privateinherited

Used to fix the ordering of tools when we're initializing from COOL based on a hierarchical tag.

Tools should be executed in order of increasing m_order.

Definition at line 547 of file ToolWithConstants.h.

{ this, "order", 0 };

m_prefix

StringProperty CaloUtils::ToolWithConstants< BASE >::m_prefix

privateinherited

Prefix for finding our constants within the ToolConstants object.

Definition at line 538 of file ToolWithConstants.h.

{ this, "prefix", "", "" };

m_update_energy

Constant<bool> CaloClusterUpdate::m_update_energy { this, "update_energy", "Should cluster total energy be updated?" }

private

Definition at line 64 of file CaloClusterUpdate.h.

{ this, "update_energy", "Should cluster total energy be updated?" };

m_w1

Gaudi::Property<float> CaloClusterUpdate::m_w1 { this, "layer1_weight", 3, "layer 1 additional weight in eta average" }

private

Definition at line 60 of file CaloClusterUpdate.h.

{ this, "layer1_weight", 3, "layer 1 additional weight in eta average" };

---

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

• CaloClusterUpdate.h
• CaloClusterUpdate.cxx