ATLAS Offline Software
Loading...
Searching...
No Matches
CaloSwGap_g3 Class Referenceabstract

#include <CaloSwGap_g3.h>

Inheritance diagram for CaloSwGap_g3:
Collaboration diagram for CaloSwGap_g3:

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 StatusCode initialize () override
 Standard Gaudi initialize method.
virtual void makeCorrection (const Context &myctx, xAOD::CaloCluster *cluster) const override
virtual StatusCode setCaloCellContainerName (const std::string &name) override
 Change the name of the CaloCellContainer used by this tool.
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< CaloDetDescrManagerm_caloMgrKey {this,"CaloDetDescrManager", "CaloDetDescrManager"}

Private Member Functions

 CaloSwGap_g3 ()=delete

Private Attributes

Constant< float > m_etamin_crack { this, "etamin_crack", "" }
Constant< float > m_etamax_crack { this, "etamax_crack", "" }
Constant< CxxUtils::Array< 1 > > m_scint_weight { this, "scint_weight", "" }
Constant< CxxUtils::Array< 1 > > m_correction { this, "correction", "" }
SG::ReadHandleKey< CaloCellContainerm_cells_name { this, "cells_name", "AllCalo", "" }
 Property: The name of the container in which to look to find tile cells.
SG::ReadCondHandleKey< CaloRec::ToolConstantsm_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 34 of file CaloSwGap_g3.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.

Constructor & Destructor Documentation

◆ CaloSwGap_g3()

CaloSwGap_g3::CaloSwGap_g3 ( )
privatedelete

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.

55{
56 this->makeCorrection (context(ctx), cluster);
57
58#if 0
59 ATH_MSG_DEBUG( " ...... e, et " << cluster->e() << " " << cluster->et() << endmsg);
60 ATH_MSG_DEBUG( " ...... eta, etaBE, etaSmp " << cluster->eta() << " " << cluster->etaBE(2)
61 << " " << cluster->etaSample(CaloSampling::EMB1)
62 << " " << cluster->etaSample(CaloSampling::EMB2)
63 << " " << cluster->etaSample(CaloSampling::EMB3) << endmsg);
64 ATH_MSG_DEBUG( " ...... phi, phiBE, phiSmp " << cluster->phi() << " " << cluster->phiBE(2)
65 << " " << cluster->phiSample(CaloSampling::EMB1)
66 << " " << cluster->phiSample(CaloSampling::EMB2)
67 << " " << cluster->phiSample(CaloSampling::EMB3) << endmsg);
68#endif
69
70 return StatusCode::SUCCESS;
71}
#define endmsg
#define ATH_MSG_DEBUG(x)
virtual void makeCorrection(const Context &myctx, xAOD::CaloCluster *) const =0
Context context(const EventContext &ctx) const
float phiBE(const unsigned layer) const
Get the phi in one layer of the EM Calo.
float phiSample(const CaloSample sampling) const
Retrieve barycenter in a given sample.
virtual double eta() const
The pseudorapidity ( ) of the particle.
virtual double e() const
The total energy of the particle.
virtual double phi() const
The azimuthal angle ( ) of the particle.
float etaSample(const CaloSample sampling) const
Retrieve barycenter in a given sample.
float etaBE(const unsigned layer) const
Get the eta in one layer of the EM Calo.

◆ execute() [2/2]

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

Execute on an entire collection of clusters.

Parameters
collectionThe container of clusters.

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

Parameters
collectionThe container of clusters.
ctxThe 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.

47{
48 for (xAOD::CaloCluster* clu : *collection) {
49 ATH_CHECK( execute (ctx, clu) );
50 }
51 return StatusCode::SUCCESS;
52}
#define ATH_CHECK
Evaluate an expression and check for errors.
StatusCode execute(const EventContext &ctx, xAOD::CaloCluster *cluster) const override
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.

◆ initialize()

StatusCode CaloSwGap_g3::initialize ( )
overridevirtual

Standard Gaudi initialize method.

Reimplemented from CaloClusterCorrection.

Definition at line 54 of file CaloSwGap_g3.cxx.

55{
57 ATH_CHECK( m_cells_name.initialize() );
58 return StatusCode::SUCCESS;
59}
virtual StatusCode initialize() override
Initialize method.
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.

◆ makeCorrection() [1/2]

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

◆ makeCorrection() [2/2]

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

Definition at line 62 of file CaloSwGap_g3.cxx.

64{
65 float eta = cluster->eta();
66 float phi = cluster->phi();
67
68 float aeta = fabs(eta);
69
70 const float etamin_crack = m_etamin_crack (myctx);
71 const float etamax_crack = m_etamax_crack (myctx);
72 if (aeta < etamin_crack || aeta > etamax_crack)
73 return; // no correction required
74
75 const CxxUtils::Array<1> scint_weight = m_scint_weight (myctx);
76 const CxxUtils::Array<1> correction = m_correction (myctx);
77 assert (scint_weight.size() == correction.size());
78 if (correction.size() == 0) return;
79
80 SG::ReadHandle<CaloCellContainer> cc (m_cells_name, myctx.ctx());
81
82 double eh_scint = 0;
83 if(cc.isValid())
84 {
86 cc->beginConstCalo(CaloCell_ID::TILE);
88 cc->endConstCalo(CaloCell_ID::TILE);
89
90 for ( ; f_cell!=l_cell; ++f_cell)
91 {
92 const CaloCell* cell = (*f_cell) ;
93
94 int inum = cell->caloDDE()->getSampling();
95
96 if (CaloSampling::TileGap3 == inum ) {
97 double phic = cell->phi();
98 double etac = cell->eta();
99 if (phic-phi < -pi) phic = phic + twopi;
100 if (phic-phi > pi) phic = phic - twopi;
101 if (etac >= eta-deta && etac <= eta+deta) {
102 if (phic >= phi-dphi && phic <= phi+dphi) {
103 eh_scint += cell->e();
104 }
105 }
106 }
107 }
108 }
109
110 // make the correction
111 double granularity = (etamax_crack - etamin_crack)/correction.size();
112 int ind = static_cast<int> ((aeta - etamin_crack)/granularity);
113
114// Scintillator energy calibration
115 eh_scint = eh_scint * scint_weight[ind];
116
117//Correct for normalization factor total cluster energy and samplings energies
118 double energy = (cluster->e())*correction[ind];
119 setenergy(cluster,energy);
120
121// Add the scintillator energy to the cluster and correct the total energy
122 energy += eh_scint * correction[ind];
123 cluster->setE(energy);
124}
Scalar eta() const
pseudorapidity method
Scalar phi() const
phi method
#define pi
constexpr double twopi
Data object for each calorimeter readout cell.
Definition CaloCell.h:57
virtual void setenergy(xAOD::CaloCluster *cluster, float energy) const
Read-only multidimensional array.
unsigned int size(unsigned int dim=0) const
Return the size of the array along one dimension.
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:838
correction(mu, runmode, campaign, run=None)
Definition zlumi_mc_cf.py:4
const EventContext & ctx() const

◆ 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]outObject to receive our constants.
ctxEvent context.

◆ setCaloCellContainerName()

StatusCode CaloSwGap_g3::setCaloCellContainerName ( const std::string & name)
overridevirtual

Change the name of the CaloCellContainer used by this tool.

Parameters
nameThe new container name.

Implements ISetCaloCellContainerName.

Definition at line 132 of file CaloSwGap_g3.cxx.

134{
135 return this->setProperty (StringProperty ("cells_name", name));
136}
void setProperty(columnar::PythonToolHandle &self, const std::string &key, nb::object value)

◆ setenergy()

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

Definition at line 94 of file CaloClusterCorrection.cxx.

95{
96 if (cluster->e() == 0) {
97 if (energy != 0)
98 REPORT_MESSAGE (MSG::WARNING)
99 << "Attempt to rescale zero-energy cluster to energy " << energy
100 << " ignored.";
101 return;
102 }
103
104 float correction = energy/cluster->e();
105 cluster->setE(energy);
106
107 // also correct individual sampling energies:
108
109 for (int iSample=CaloSampling::PreSamplerB;
110 iSample < CaloSampling::Unknown;
111 ++iSample)
112 {
113 CaloSampling::CaloSample sampling=static_cast<CaloSampling::CaloSample>(iSample);
114 if (cluster->hasSampling (sampling)) {
115 double e = cluster->eSample(sampling);
116 cluster->setEnergy(sampling,e*correction) ;
117 }
118 }
119}
#define REPORT_MESSAGE(LVL)
Report a message.
float eSample(const CaloSample sampling) const
bool setEnergy(const CaloSample sampling, const float e)
Set energy for a given sampling. Returns false if the sample isn't part of the cluster.
bool hasSampling(const CaloSample s) const
Checks if certain smapling contributes to cluster.

◆ 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.

79{
80 cluster->setEnergy(sampling, em);
81 cluster->setEta(sampling, etam);
82 cluster->setPhi(sampling, phim);
83
84 cluster->setEmax(sampling,emax);
85 cluster->setEtamax(sampling,etamax);
86 cluster->setPhimax(sampling,phimax);
87
88 cluster->setEtasize(sampling, etas);
89 cluster->setPhisize(sampling, phis);
90}
bool setPhi(const CaloSample sampling, const float phi)
Set in a given sampling. Returns false if the sample isn't part of the cluster.
bool setEtasize(const CaloSample sampling, const float etaSize)
Set the cluster size in for a given sampling.
bool setPhisize(const CaloSample sampling, const float phiSize)
Set the cluster size in for a given sampling.
bool setPhimax(const CaloSample sampling, const float phiMax)
Set the phi of the cell with the highest energy in a particular sampling.
bool setEta(const CaloSample sampling, const float eta)
Set in a given sampling. Returns false if the sample isn't part of the cluster.
bool setEmax(const CaloSample sampling, const float eMax)
Set the Energy of the cell with the highest energy in a particular sampling.
bool setEtamax(const CaloSample sampling, const float etaMax)
Set the eta of the cell with the highest energy in a particular sampling.

◆ 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
streamOstream to which to write.
nameName to go in output
ctxEvent context.

Member Data Documentation

◆ m_caloMgrKey

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

Definition at line 83 of file CaloClusterCorrection.h.

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

◆ m_cells_name

SG::ReadHandleKey<CaloCellContainer> CaloSwGap_g3::m_cells_name { this, "cells_name", "AllCalo", "" }
private

Property: The name of the container in which to look to find tile cells.

Definition at line 66 of file CaloSwGap_g3.h.

67{ this, "cells_name", "AllCalo", "" };

◆ m_correction

Constant<CxxUtils::Array<1> > CaloSwGap_g3::m_correction { this, "correction", "" }
private

Definition at line 63 of file CaloSwGap_g3.h.

63{ this, "correction", "" };

◆ m_DBHandle

Handle to a ToolConstants conditions object.

Definition at line 534 of file ToolWithConstants.h.

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

◆ m_etamax_crack

Constant<float> CaloSwGap_g3::m_etamax_crack { this, "etamax_crack", "" }
private

Definition at line 61 of file CaloSwGap_g3.h.

61{ this, "etamax_crack", "" };

◆ m_etamin_crack

Constant<float> CaloSwGap_g3::m_etamin_crack { this, "etamin_crack", "" }
private

Definition at line 60 of file CaloSwGap_g3.h.

60{ this, "etamin_crack", "" };

◆ m_impl

ToolWithConstantsImpl CaloUtils::ToolWithConstants< BASE >::m_impl
privateinherited

Internal implementation object.

Definition at line 542 of file ToolWithConstants.h.

542{ this->name(), m_prefix, m_DBHandle };
Hold constants for a tool.

◆ 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.

561{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.

548{ 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.

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

◆ m_scint_weight

Constant<CxxUtils::Array<1> > CaloSwGap_g3::m_scint_weight { this, "scint_weight", "" }
private

Definition at line 62 of file CaloSwGap_g3.h.

62{ this, "scint_weight", "" };

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