CaloSwGap_g3 Class Reference

#include <CaloSwGap_g3.h>

Inheritance diagram for CaloSwGap_g3:

Collaboration diagram for CaloSwGap_g3:

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
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
	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< CaloCellContainer >	m_cells_name { this, "cells_name", "AllCalo", "" }
	Property: The name of the container in which to look to find tile cells.

Detailed Description

Definition at line 34 of file CaloSwGap_g3.h.

Constructor & Destructor Documentation

CaloSwGap_g3()

CaloSwGap_g3::CaloSwGap_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 CaloSwGap_g3::initialize ( )

overridevirtual

Standard Gaudi initialize method.

Reimplemented from CaloClusterCorrection.

Definition at line 54 of file CaloSwGap_g3.cxx.

{
  ATH_CHECK( CaloClusterCorrection::initialize() );
  ATH_CHECK( m_cells_name.initialize() );
  return StatusCode::SUCCESS;
}

makeCorrection()

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

overridevirtual

Implements CaloClusterCorrection.

Definition at line 62 of file CaloSwGap_g3.cxx.

{
  float eta = cluster->eta();
  float phi = cluster->phi();
 
  float aeta = fabs(eta);
 
  const float etamin_crack = m_etamin_crack (myctx);
  const float etamax_crack = m_etamax_crack (myctx);
  if (aeta < etamin_crack || aeta > etamax_crack) 
    return; // no correction required
 
  const CxxUtils::Array<1> scint_weight = m_scint_weight (myctx);
  const CxxUtils::Array<1> correction = m_correction (myctx);
  assert (scint_weight.size() == correction.size());
  if (correction.size() == 0) return;
 
  SG::ReadHandle<CaloCellContainer> cc (m_cells_name, myctx.ctx());
 
  double eh_scint = 0;
  if(cc.isValid())
  {
    CaloCellContainer::const_iterator f_cell =
      cc->beginConstCalo(CaloCell_ID::TILE);
    CaloCellContainer::const_iterator l_cell =
      cc->endConstCalo(CaloCell_ID::TILE);
 
    for ( ; f_cell!=l_cell; ++f_cell)
    {
      const CaloCell* cell = (*f_cell) ; 
 
      int inum = cell->caloDDE()->getSampling();
 
      if (CaloSampling::TileGap3 == inum  ) {
        double phic = cell->phi();
        double etac = cell->eta(); 
        if (phic-phi < -pi) phic = phic + twopi;
        if (phic-phi > pi)  phic = phic - twopi;    
        if (etac >= eta-deta && etac <= eta+deta) {
          if (phic >= phi-dphi && phic <= phi+dphi) {
            eh_scint += cell->e();
          }
        }
      }
    }
  }
 
  // make the correction
  double granularity = (etamax_crack - etamin_crack)/correction.size();  
  int ind = static_cast<int> ((aeta - etamin_crack)/granularity);
 
// Scintillator energy calibration
  eh_scint = eh_scint * scint_weight[ind];
 
//Correct for normalization factor total cluster energy and samplings energies 
  double energy = (cluster->e())*correction[ind];
  setenergy(cluster,energy);
 
// Add the scintillator energy to the cluster and correct the total energy
  energy += eh_scint * correction[ind]; 
  cluster->setE(energy);  
}

setCaloCellContainerName()

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

overridevirtual

Change the name of the CaloCellContainer used by this tool.

Parameters

name	The new container name.

Implements ISetCaloCellContainerName.

Definition at line 132 of file CaloSwGap_g3.cxx.

{
  return this->setProperty (StringProperty ("cells_name", name));
}

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

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

{ this, "correction",   "" };

m_etamax_crack

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

private

Definition at line 61 of file CaloSwGap_g3.h.

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

{ this, "etamin_crack", "" };

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.

{ this, "scint_weight", "" };

---

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

• CaloSwGap_g3.h
• CaloSwGap_g3.cxx