Class to hold the information ECL energy resolution derived from PERC. More...

#include <ECLPhotonEnergyResolution.h>

Inheritance diagram for ECLPhotonEnergyResolution:

Public Member Functions
	ECLPhotonEnergyResolution ()
	!

	~ECLPhotonEnergyResolution ()
	!

ParticleWeightingLookUpTable	getFullEnergyResolution () const
	Get ParticleWeightingLookUpTable with binned energy resolution.

void	addRelativeEnergyResolution (std::vector< double > energyResolution, Binning binning)
	Add energy bin [vector of three pairs] and energy resolution [vector of three values] Energy resolution structure: [value, uncertainty up, uncertainty down] Energy bin structure for bin limit pairs: [energy limits, theta limits, phi limits].

double	getRelativeEnergyResolution (double energy, double theta, double phi) const
	Get energy resolution for given energy, theta, phi.

void	addThetaPhiResolution (std::vector< double > thetaResolution, std::vector< double > phiResolution, Binning binning)
	For neutral hadrons Add theta resolution [vector of three values] and phi resolution [vector of three values] in bins of (clusterUncEnergy, theta, phi) [vector of three pairs] thetaResolution structure: [value, uncertainty up, uncertainty down] phiResolution structure: [value, uncertainty up, uncertainty down] bin structure for bin limit pairs: [energy limits, theta limits, phi limits].

double	getThetaPhiResolution (double energy, double theta, double phi, bool thetaOrPhi) const
	Get theta or phi resolution for given energy, theta, phi.

Private Member Functions
	ClassDef (ECLPhotonEnergyResolution, 2)
	ClassDef.

Private Attributes
ParticleWeightingLookUpTable	m_resolutionBinningTable
	Utilise ParticleWeighting modules as holders for binning of energy resolution.

Detailed Description

Class to hold the information ECL energy resolution derived from PERC.

Currently the energy resolution is given in energy bins and only valid for barrel region

Definition at line 32 of file ECLPhotonEnergyResolution.h.

Constructor & Destructor Documentation

◆ ECLPhotonEnergyResolution()

ECLPhotonEnergyResolution ( )

!

Default constructor

Define out of range resolution value as -1

Definition at line 15 of file ECLPhotonEnergyResolution.cc.

{
  WeightInfo energyResolutionInfo = {
    {"RelativeEnergyResolution", -1.},
    {"RelativeEnergyResolutionUncertUp", -1.},
    {"RelativeEnergyResolutionUncertDown", -1.}
  };
 
  m_resolutionBinningTable.defineOutOfRangeWeight(energyResolutionInfo);
}

◆ ~ECLPhotonEnergyResolution()

~ECLPhotonEnergyResolution ( )

inline

!

Destructor

Definition at line 42 of file ECLPhotonEnergyResolution.h.

42{};

Member Function Documentation

◆ addRelativeEnergyResolution()

void addRelativeEnergyResolution	(	std::vector< double >	energyResolution,
		Binning	binning )

Add energy bin [vector of three pairs] and energy resolution [vector of three values] Energy resolution structure: [value, uncertainty up, uncertainty down] Energy bin structure for bin limit pairs: [energy limits, theta limits, phi limits].

Parameters

energyResolution	energy resolution with uncertainties
binning	binning of energy resolution valid region in bins of [energy, theta phi]

Build energy resolution with uncertainties as info map

Extract energy, theta and phi bin values

Add resolution value to binning table

Definition at line 29 of file ECLPhotonEnergyResolution.cc.

{
  WeightInfo energyResolutionInfo = {
    {"RelativeEnergyResolution", energyResolution.at(0)},
    {"RelativeEnergyResolutionUncertUp", energyResolution.at(1)},
    {"RelativeEnergyResolutionUncertDown", energyResolution.at(2)}
  };
  ParticleWeightingBinLimits* energyBinValues = new ParticleWeightingBinLimits(binning.at(0).first, binning.at(0).second);
  ParticleWeightingBinLimits* thetaBinValues = new ParticleWeightingBinLimits(binning.at(1).first, binning.at(1).second);
  ParticleWeightingBinLimits* phiBinValues = new ParticleWeightingBinLimits(binning.at(2).first, binning.at(2).second);
 
  NDBin energyBinning = {{"Energy", energyBinValues}, {"Theta", thetaBinValues}, {"Phi", phiBinValues}};
  m_resolutionBinningTable.addEntry(energyResolutionInfo, energyBinning);
}

◆ addThetaPhiResolution()

void addThetaPhiResolution	(	std::vector< double >	thetaResolution,
		std::vector< double >	phiResolution,
		Binning	binning )

For neutral hadrons Add theta resolution [vector of three values] and phi resolution [vector of three values] in bins of (clusterUncEnergy, theta, phi) [vector of three pairs] thetaResolution structure: [value, uncertainty up, uncertainty down] phiResolution structure: [value, uncertainty up, uncertainty down] bin structure for bin limit pairs: [energy limits, theta limits, phi limits].

Parameters

thetaResolution	theta resolution with uncertainties
phiResolution	phi resolution with uncertainties
binning	binning of theta/phi resolution valid region in bins of [energy, theta phi]

Build energy resolution with uncertainties as info map

Extract energy, theta and phi bin values

Add resolution value to binning table

Definition at line 74 of file ECLPhotonEnergyResolution.cc.

{
  WeightInfo resolutionInfo = {
    {"ThetaResolution", thetaResolution.at(0)},
    {"ThetaResolutionUncertUp", thetaResolution.at(1)},
    {"ThetaResolutionUncertDown", thetaResolution.at(2)},
    {"PhiResolution", phiResolution.at(0)},
    {"PhiResolutionUncertUp", phiResolution.at(1)},
    {"PhiResolutionUncertDown", phiResolution.at(2)}
  };
  ParticleWeightingBinLimits* energyBinValues = new ParticleWeightingBinLimits(binning.at(0).first, binning.at(0).second);
  ParticleWeightingBinLimits* thetaBinValues = new ParticleWeightingBinLimits(binning.at(1).first, binning.at(1).second);
  ParticleWeightingBinLimits* phiBinValues = new ParticleWeightingBinLimits(binning.at(2).first, binning.at(2).second);
 
  NDBin energyBinning = {{"Energy", energyBinValues}, {"Theta", thetaBinValues}, {"Phi", phiBinValues}};
  m_resolutionBinningTable.addEntry(resolutionInfo, energyBinning);
}

◆ getFullEnergyResolution()

ParticleWeightingLookUpTable getFullEnergyResolution ( ) const

inline

Get ParticleWeightingLookUpTable with binned energy resolution.

Returns: ParticleWeightingLookUpTable of the energy resolution in energy, theta, and phi bins

Definition at line 48 of file ECLPhotonEnergyResolution.h.

48{ return m_resolutionBinningTable; }

◆ getRelativeEnergyResolution()

double getRelativeEnergyResolution	(	double	energy,
		double	theta,
		double	phi ) const

Get energy resolution for given energy, theta, phi.

Parameters

energy	photon energy value
theta	photon theta value
phi	photon phi value

Returns: Energy resolution

Construct energy, theta, phi map corresponding to binning structure

Extract energy resolution information

Definition at line 55 of file ECLPhotonEnergyResolution.cc.

{
  std::map<std::string, double> binning = {{"Energy", energy}, {"Theta", theta}, {"Phi", phi}};
  WeightInfo energyResolutionInfo = m_resolutionBinningTable.getInfo(binning);
 
  if (energyResolutionInfo.at("RelativeEnergyResolution") == -1.) {
    B2DEBUG(20, "Energy resolution was returned as -1 since provided energy was not within energy resolution binning.");
  }
 
  return energyResolutionInfo.at("RelativeEnergyResolution");
}

◆ getThetaPhiResolution()

double getThetaPhiResolution	(	double	energy,
		double	theta,
		double	phi,
		bool	thetaOrPhi ) const

Get theta or phi resolution for given energy, theta, phi.

Parameters

energy	- for photons, energy value, for hadrons clusterUnCorr energy value
theta	photon theta value
phi	photon phi value
thetaOrPhi	flag to pick theta or phi, 0 - return theta resolution, 1 - return phi resolution

Returns: Theta or Phiresolution

Construct energy, theta, phi map corresponding to binning structure

Extract resolution information

Definition at line 104 of file ECLPhotonEnergyResolution.cc.

{
  std::map<std::string, double> binning = {{"Energy", energy}, {"Theta", theta}, {"Phi", phi}};
  WeightInfo resolutionInfo = m_resolutionBinningTable.getInfo(binning);
 
  if (thetaOrPhi) {
    if (resolutionInfo.at("PhiResolution") == -1.) {
      B2DEBUG(20, "Phi resolution was returned as -1 since provided energy was not within theta resolution binning.");
    }
    return resolutionInfo.at("PhiResolution");
  } else {
    if (resolutionInfo.at("ThetaResolution") == -1.) {
      B2DEBUG(20, "Theta resolution was returned as -1 since provided energy was not within theta resolution binning.");
    }
    return resolutionInfo.at("ThetaResolution");
  }
}

Member Data Documentation

◆ m_resolutionBinningTable

ParticleWeightingLookUpTable m_resolutionBinningTable

private

Utilise ParticleWeighting modules as holders for binning of energy resolution.

Definition at line 95 of file ECLPhotonEnergyResolution.h.

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

analysis/dbobjects/include/ECLPhotonEnergyResolution.h
analysis/dbobjects/src/ECLPhotonEnergyResolution.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ECLPhotonEnergyResolution()

◆ ~ECLPhotonEnergyResolution()

Member Function Documentation

◆ addRelativeEnergyResolution()

◆ addThetaPhiResolution()

◆ getFullEnergyResolution()

◆ getRelativeEnergyResolution()

◆ getThetaPhiResolution()

Member Data Documentation

◆ m_resolutionBinningTable