CaloSwPhimod_g3 Class Reference

#include <CaloSwPhimod_g3.h>

Inheritance diagram for CaloSwPhimod_g3:

Collaboration diagram for CaloSwPhimod_g3:

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 Member Functions
	CaloSwPhimod_g3 ()=delete
std::vector< float >	qphmop (const Context &myctx, float eta) const

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

Detailed Description

Definition at line 41 of file CaloSwPhimod_g3.h.

Constructor & Destructor Documentation

CaloSwPhimod_g3()

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

overridevirtual

Implements CaloClusterCorrection.

Definition at line 38 of file CaloSwPhimod_g3.cxx.

{
 
  float eta2 = cluster->etaBE(2); // use second sampling
  float phi2 = cluster->phiBE(2); // use the average over Bar and EC.
  if (eta2 < -900 || phi2 < -900)
    return;
 
  float aeta2 = fabs(eta2);
  float nabs;
  if (aeta2 < 1.5) 
      nabs = 1024; 
  else if (aeta2 < 2.5)
      nabs = 768;
  else
      nabs = 256.;
 
  std::vector<float> tab = qphmop (myctx, aeta2);
  assert (tab.size() == 5);
  float qphimod = 
    tab[0] + 
    tab[1] * cos(nabs*phi2) +
    tab[2] * cos(2.*nabs*phi2)+
    tab[3] * cos(3.*nabs*phi2)+ 
    tab[4] * sin(nabs*phi2);
  
  // set energy, and rescale each sampling
  setenergy(cluster, qphimod * cluster->e());
  
  
}

qphmop()

std::vector< float > CaloSwPhimod_g3::qphmop ( const Context & myctx, float eta ) const

private

Definition at line 71 of file CaloSwPhimod_g3.cxx.

{
  std::vector<float> coef;
 
  const CxxUtils::Array<2> correction = m_correction (myctx);
  if (correction.size() == 0) {
    ATH_MSG_ERROR("Empty correction array");
    coef.resize(5);
    return coef;
  }
 
  coef.resize (correction.size(1)-1);
 
  int i1 = 0;
  int i2 = 0;
  if (aeta < correction[0][0]) {
    for (unsigned int j = 1; j < correction.size(1); j++)
      coef[j-1] = correction[0][j];
    return coef;
  }
    
  else if (aeta >= correction[0][0]) {
    for (unsigned int i = 0; i < correction.size()-1; i++) {
      if (aeta >= correction[i][0] && aeta < correction[i+1][0]) {
        i1 = i;
        i2 = i+1;
        break;
      }
    }
  }
  if (i1 == i2) {
    i2 = correction.size()-1;
    i1 = i2-1;
  }
 
  float deta = correction[i2][0] - correction[i1][0];
  if (deta == 0)
    deta = 1;
 
  float m = (aeta - correction[i1][0]) / deta;
  for (unsigned int j = 1; j < correction.size(1); j++) {
    float dy = correction[i2][j] - correction[i1][j];
    coef[j-1] = (correction[i1][j] + dy * m);
  }
 
  return coef;
}

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_correction

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

private

Definition at line 56 of file CaloSwPhimod_g3.h.

{ this, "correction", "" };

---

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

• CaloSwPhimod_g3.h
• CaloSwPhimod_g3.cxx