ATLAS Offline Software
Loading...
Searching...
No Matches
CaloSwEta2e_g3 Class Reference

#include <CaloSwEta2e_g3.h>

Inheritance diagram for CaloSwEta2e_g3:
Collaboration diagram for CaloSwEta2e_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< CaloDetDescrManagerm_caloMgrKey {this,"CaloDetDescrManager", "CaloDetDescrManager"}

Private Member Functions

 CaloSwEta2e_g3 ()=delete

Private Attributes

Constant< CxxUtils::Array< 2 > > m_correction { this, "correction", "" }
Constant< int > m_correction_degree { this, "correction_degree", "" }
Constant< int > m_region { this, "region", 1, "" }

Static Private Attributes

static const float s_middle_layer_granularity = 0.025

Detailed Description

Definition at line 36 of file CaloSwEta2e_g3.h.

Constructor & Destructor Documentation

◆ CaloSwEta2e_g3()

CaloSwEta2e_g3::CaloSwEta2e_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.

54{
55 this->makeCorrection (context(ctx), cluster);
56
57#if 0
58 ATH_MSG_DEBUG( " ...... e, et " << cluster->e() << " " << cluster->et() << endmsg);
59 ATH_MSG_DEBUG( " ...... eta, etaBE, etaSmp " << cluster->eta() << " " << cluster->etaBE(2)
60 << " " << cluster->etaSample(CaloSampling::EMB1)
61 << " " << cluster->etaSample(CaloSampling::EMB2)
62 << " " << cluster->etaSample(CaloSampling::EMB3) << endmsg);
63 ATH_MSG_DEBUG( " ...... phi, phiBE, phiSmp " << cluster->phi() << " " << cluster->phiBE(2)
64 << " " << cluster->phiSample(CaloSampling::EMB1)
65 << " " << cluster->phiSample(CaloSampling::EMB2)
66 << " " << cluster->phiSample(CaloSampling::EMB3) << endmsg);
67#endif
68
69 return StatusCode::SUCCESS;
70}
#define endmsg
#define ATH_MSG_DEBUG(x)
virtual void makeCorrection(const Context &myctx, xAOD::CaloCluster *) const =0
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.

46{
47 for (xAOD::CaloCluster* clu : *collection) {
48 ATH_CHECK( execute (ctx, clu) );
49 }
50 return StatusCode::SUCCESS;
51}
#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 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.

46 {
47 ATH_CHECK(m_caloMgrKey.initialize());
48 ATH_CHECK(base_class::initialize());
49 return StatusCode::SUCCESS;
50}
SG::ReadCondHandleKey< CaloDetDescrManager > m_caloMgrKey

◆ makeCorrection()

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

Implements CaloClusterCorrection.

Definition at line 44 of file CaloSwEta2e_g3.cxx.

46{
47 // Only for endcap
48 if (!cluster->inEndcap())
49 return;
50
51 float eta = cluster->etaSample(CaloSampling::EME2);
52 if (eta == -999.) return;
53 float etamax = cluster->etamax(CaloSampling::EME2);
54 float aeta = fabs(eta);
55
56 // peak around center of cell (+=0.0125)
57 float u2 = fmod(eta-etamax,s_middle_layer_granularity/2.);
59 // be careful with negative eta values, position correction is not symmetric.
60 if (eta < 0.) u2 = s_middle_layer_granularity-u2;
61
62 // First order (average) S-shape correction
63 float etc2 = aeta - interpolate (m_correction(myctx),
64 u2,
65 m_correction_degree(myctx));
66
67 if (eta < 0.)
68 etc2 = -etc2;
69
70 // make correction to second sampling
71 cluster->setEta( CaloSampling::EME2, etc2);
72}
Scalar eta() const
pseudorapidity method
Constant< int > m_correction_degree
Constant< CxxUtils::Array< 2 > > m_correction
static const float s_middle_layer_granularity
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 inEndcap() const
Returns true if at least one clustered cell in the endcap.
float etamax(const CaloSample sampling) const
Retrieve of cell with maximum energy in given sampling.
float interpolate(const CaloRec::Array< 2 > &a, float x, unsigned int degree, unsigned int ycol=1, const CaloRec::Array< 1 > &regions=CaloRec::Array< 1 >(), int n_points=-1, bool fixZero=false)
Polynomial interpolation in a table.

◆ setenergy()

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

Definition at line 93 of file CaloClusterCorrection.cxx.

94{
95 if (cluster->e() == 0) {
96 if (energy != 0)
97 REPORT_MESSAGE (MSG::WARNING)
98 << "Attempt to rescale zero-energy cluster to energy " << energy
99 << " ignored.";
100 return;
101 }
102
103 float correction = energy/cluster->e();
104 cluster->setE(energy);
105
106 // also correct individual sampling energies:
107
108 for (int iSample=CaloSampling::PreSamplerB;
109 iSample < CaloSampling::Unknown;
110 ++iSample)
111 {
112 CaloSampling::CaloSample sampling=static_cast<CaloSampling::CaloSample>(iSample);
113 if (cluster->hasSampling (sampling)) {
114 double e = cluster->eSample(sampling);
115 cluster->setEnergy(sampling,e*correction) ;
116 }
117 }
118}
#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.
correction(mu, runmode, campaign, run=None)
Definition zlumi_mc_cf.py:4

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

78{
79 cluster->setEnergy(sampling, em);
80 cluster->setEta(sampling, etam);
81 cluster->setPhi(sampling, phim);
82
83 cluster->setEmax(sampling,emax);
84 cluster->setEtamax(sampling,etamax);
85 cluster->setPhimax(sampling,phimax);
86
87 cluster->setEtasize(sampling, etas);
88 cluster->setPhisize(sampling, phis);
89}
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 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()

const std::string & CaloSwEta2e_g3::toolType ( ) const
overridevirtual

Definition at line 75 of file CaloSwEta2e_g3.cxx.

76{
77 static const std::string typeName = "CaloSwEta_g3";
78 return typeName;
79}

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_correction

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

Definition at line 54 of file CaloSwEta2e_g3.h.

54{ this, "correction", "" };

◆ m_correction_degree

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

Definition at line 55 of file CaloSwEta2e_g3.h.

55{ this, "correction_degree", "" };

◆ m_region

Constant<int> CaloSwEta2e_g3::m_region { this, "region", 1, "" }
private

Definition at line 56 of file CaloSwEta2e_g3.h.

56{ this, "region", 1, "" };

◆ s_middle_layer_granularity

const float CaloSwEta2e_g3::s_middle_layer_granularity = 0.025
staticprivate

Definition at line 58 of file CaloSwEta2e_g3.h.


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