CaloClusterUpdate Class Reference

#include <CaloClusterUpdate.h>

Inheritance diagram for CaloClusterUpdate:

Collaboration diagram for CaloClusterUpdate:

Public Member Functions
virtual void	makeCorrection (const Context &myctx, xAOD::CaloCluster *cluster) 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 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?" }

Detailed Description

Definition at line 46 of file CaloClusterUpdate.h.

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 CaloClusterUpdate::makeCorrection ( const Context & myctx, xAOD::CaloCluster * cluster ) const

overridevirtual

Implements CaloClusterCorrection.

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);
}

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);
}

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_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