#include <src/CondFormats/interface/EcalShowerContainmentCorrections.h>

Classes
struct	Coefficients
	Structure defining the container for correction coefficients. More...

Public Member Functions
const double	correction3x3 (const EBDetId &xtal, const math::XYZPoint &pos) const
	The correction factor for 3x3 matrix. More...

const double	correction5x5 (const EBDetId &xtal, const math::XYZPoint &pos) const
	The correction factor for 5x5 matrix. More...

const Coefficients	correctionCoefficients (const EBDetId &centerxtal) const
	Get the correction coefficients for the given xtal. More...

void	fillCorrectionCoefficients (const EBDetId &xtal, int group, const Coefficients &coefficients)
	Fill the correction coefficients for a given xtal, part of group @group. More...

void	fillCorrectionCoefficients (const int supermodule, const int module, const Coefficients &coefficients)
	Fill the correction coefficients for a given Ecal module. More...

Private Types
enum	Direction { eX, eY }

typedef std::map< EBDetId, int >	GroupMap

enum	Type { e3x3, e5x5 }

Private Member Functions
const double	correctionXY (const EBDetId &xtal, double position, Direction dir, Type type) const

template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Private Attributes
std::vector< Coefficients >	coefficients_
	Holds the coeffiecients. The index corresponds to the group. More...

GroupMap	groupmap_
	Maps in which group a particular xtal has been placed. More...

Friends
class	boost::serialization::access

template<typename CondSerializationT , typename Enabled >
struct	cond::serialization::access

Detailed Description

Description: Holds the coefficients of a polynomial that describes the shower containment.

Usage example(for real code see CalibCalorimetry/EcalCorrectionModules/test) :

ESHandle<EcalShowerContainmentCorrections> pGapCorr;
iESetup.get<EcalShowerContainmentCorrectionsRcd>().get(pGapCorr);
 
PositionCalc pos(...);
 
Hep3Vector clusterPos= pos.CalculatePos(...);         
math::XYZPoint mathpoint(clusterPos.x(),clusterPos.y(),clusterPos.z());
 
double correction3x3 = pGapCorr->correction3x3(centerXtal,mathpoint);
double correction5x5 = pGapCorr->correction5x5(centerXtal,mathpoint);

Author: Stefano Argiro'

Date: Fri Mar 2 16:50:49 CET 2007 \id

Id: EcalShowerContainmentCorrections.h,v 1.1 2007/05/15 20:37:22 argiro Exp

Definition at line 45 of file EcalShowerContainmentCorrections.h.

Member Typedef Documentation

◆ GroupMap

typedef std::map<EBDetId, int> EcalShowerContainmentCorrections::GroupMap

private

Definition at line 121 of file EcalShowerContainmentCorrections.h.

Member Enumeration Documentation

◆ Direction

enum EcalShowerContainmentCorrections::Direction

private

Enumerator
eX
eY

Definition at line 115 of file EcalShowerContainmentCorrections.h.

115 { eX, eY };

◆ Type

enum EcalShowerContainmentCorrections::Type

private

Enumerator
e3x3
e5x5

Definition at line 116 of file EcalShowerContainmentCorrections.h.

116 { e3x3, e5x5 };

Member Function Documentation

◆ correction3x3()

const double EcalShowerContainmentCorrections::correction3x3	(	const EBDetId &	xtal,
		const math::XYZPoint &	pos
	)		const

The correction factor for 3x3 matrix.

Parameters

pos	is the distance in cm from the center of the xtal as calculated in RecoEcal/EgammaCoreTools/interface/PositionCalc.h The valid return value is in the range (0,1] (divide by this value to apply the correction) Returns -1 if correction is not avaiable for that xtal

Definition at line 90 of file EcalShowerContainmentCorrections.cc.

                                                                                                                {
   double x = pos.X() * 10;  // correction functions use mm
   double y = pos.Y() * 10;
  
   double corrx = correctionXY(xtal, x, eX, e3x3);
   double corry = correctionXY(xtal, y, eY, e3x3);
  
   return corrx * corry;
 }

References correctionXY(), e3x3, eX, eY, x, and y.

◆ correction5x5()

const double EcalShowerContainmentCorrections::correction5x5	(	const EBDetId &	xtal,
		const math::XYZPoint &	pos
	)		const

The correction factor for 5x5 matrix.

Parameters

pos	is the distance in cm from the center of the xtal as calculated in RecoEcal/EgammaCoreTools/interface/PositionCalc.h The return value is in the range (0,1] (divide by this value to apply the correction) Returns -1 if correction is not avaiable for that xtal

Definition at line 100 of file EcalShowerContainmentCorrections.cc.

                                                                                                                {
   double x = pos.X() * 10;  // correction functions use mm
   double y = pos.Y() * 10;
  
   double corrx = correctionXY(xtal, x, eX, e5x5);
   double corry = correctionXY(xtal, y, eY, e5x5);
  
   return corrx * corry;
 }

References correctionXY(), e5x5, eX, eY, x, and y.

◆ correctionCoefficients()

const EcalShowerContainmentCorrections::Coefficients EcalShowerContainmentCorrections::correctionCoefficients ( const EBDetId & centerxtal ) const

Get the correction coefficients for the given xtal.

Return zero coefficients in case the correction is not available for that xtal

Definition at line 13 of file EcalShowerContainmentCorrections.cc.

                                      {
   GroupMap::const_iterator iter = groupmap_.find(centerxtal.rawId());
  
   if (iter != groupmap_.end()) {
     int group = iter->second;
     return coefficients_[group - 1];
   }
  
   edm::LogError("ShowerContaiment Correction not found");
   return Coefficients();
 }

References coefficients_, watchdog::group, groupmap_, and DetId::rawId().

◆ correctionXY()

const double EcalShowerContainmentCorrections::correctionXY	(	const EBDetId &	xtal,
		double	position,
		EcalShowerContainmentCorrections::Direction	dir,
		EcalShowerContainmentCorrections::Type	type
	)		const

private

Calculate the correction for the given direction and type

Definition at line 61 of file EcalShowerContainmentCorrections.cc.

                                                                                                              {
   GroupMap::const_iterator iter = groupmap_.find(xtal);
   if (iter == groupmap_.end())
     return -1;
  
   int group = iter->second;
   EcalShowerContainmentCorrections::Coefficients coeff = coefficients_[group - 1];
  
   int offset = 0;
  
   if (dir == eY)
     offset += 2 * Coefficients::kPolynomialDegree;
   if (position < 0)
     offset += Coefficients::kPolynomialDegree;
   if (type == e5x5)
     offset += 4 * Coefficients::kPolynomialDegree;
  
   double corr = 0;
  
   for (int i = offset; i < offset + Coefficients::kPolynomialDegree; ++i) {
     corr += coeff.data[i] * pow(position, i - offset);
   }
  
   return corr;
 }

References coefficients_, alignCSCRings::corr, EcalShowerContainmentCorrections::Coefficients::data, DeadROC_duringRun::dir, e5x5, eY, watchdog::group, groupmap_, mps_fire::i, EcalShowerContainmentCorrections::Coefficients::kPolynomialDegree, hltrates_dqm_sourceclient-live_cfg::offset, position, and funct::pow().

Referenced by correction3x3(), and correction5x5().

◆ fillCorrectionCoefficients() [1/2]

void EcalShowerContainmentCorrections::fillCorrectionCoefficients	(	const EBDetId &	xtal,
		int	group,
		const Coefficients &	coefficients
	)

Fill the correction coefficients for a given xtal, part of group @group.

Do not replace if xtal is already there

Definition at line 26 of file EcalShowerContainmentCorrections.cc.

                                                                                                     {
   // do not replace if we already have the xtal
   if (groupmap_.find(xtal) != groupmap_.end())
     return;
   groupmap_[xtal] = group;
  
   if (coefficients_.size() < (unsigned int)(group)) {
     coefficients_.resize(group);
     coefficients_[group - 1] = coefficients;
   }
  
   // we don't need to fill coefficients if the group has already been inserted
 }

References coefficients_, watchdog::group, groupmap_, and createfilelist::int.

Referenced by fillCorrectionCoefficients().

◆ fillCorrectionCoefficients() [2/2]

void EcalShowerContainmentCorrections::fillCorrectionCoefficients	(	const int	supermodule,
		const int	module,
		const Coefficients &	coefficients
	)

Fill the correction coefficients for a given Ecal module.

Assume that corresponding modules in different supermodules use the same coefficients

Definition at line 42 of file EcalShowerContainmentCorrections.cc.

                                                                                                     {
   if (module > EBDetId::kModulesPerSM) {
     edm::LogError("Invalid Module Number");
     return;
   }
  
   // what is EBDetID::kModuleBoundaries ? we better redefine them here ...
   const int kModuleLow[] = {1, 501, 901, 1301};
   const int kModuleHigh[] = {500, 900, 1300, 1700};
  
   for (int xtal = kModuleLow[module - 1]; xtal <= kModuleHigh[module - 1]; ++xtal) {
     EBDetId detid(supermodule, xtal, EBDetId::SMCRYSTALMODE);
     fillCorrectionCoefficients(detid, module, coefficients);
   }
 }