dE/dx 1D cell correction calibration constants More...

#include <CDCDedx1DCell.h>

Inheritance diagram for CDCDedx1DCell:

Public Member Functions
	CDCDedx1DCell ()
	Default constructor.

	CDCDedx1DCell (short version, const std::vector< std::vector< double > > &groupGains, const std::vector< unsigned int > &layerToGroup)
	Constructor with grouped constants.

	~CDCDedx1DCell ()
	Destructor.

CDCDedx1DCell &	operator*= (CDCDedx1DCell const &rhs)
	Combine payloads.

short	getVersion () const
	Get the version for the 1D cleanup.

unsigned int	getNGroups () const
	Number of groups.

unsigned int	getNLayers () const
	Number of layers.

unsigned int	getGroup (unsigned int layer) const
	Get group index for layer.

unsigned int	getNBins (unsigned int layer) const
	Get the number of bins for the entrance angle correction.

const std::vector< std::vector< double > > &	getOneDCell () const
	Get the grouped calibration constants.

const std::vector< unsigned int > &	getLayerMap () const
	Get the Layer-to-group map.

double	getMean (unsigned int layer, unsigned int bin) const
	Return dE/dx mean value for the given bin.

void	setMean (unsigned int layer, unsigned int bin, double value)
	Reset dE/dx mean value for the given bin.

double	getMean (unsigned int layer, double enta) const
	Return dE/dx mean value for given entrance angle.

Private Member Functions
bool	isValidGroupedPayload () const
	Payload validation.

	ClassDef (CDCDedx1DCell, 6)
	ClassDef.

Private Attributes
short	m_version
	dE/dx cleanup correction versus entrance angle may be different for different layers, so store as a vector of vectors keep a version number to identify which layers are valid

std::vector< std::vector< double > >	m_onedgains
	dE/dx means in entrance angle bins

std::vector< unsigned int >	m_layerToGroup
	layer to group map

Detailed Description

dE/dx 1D cell correction calibration constants

Definition at line 26 of file CDCDedx1DCell.h.

Constructor & Destructor Documentation

◆ CDCDedx1DCell() [1/2]

CDCDedx1DCell ( )

inline

Default constructor.

Definition at line 33 of file CDCDedx1DCell.h.

33: m_version(0), m_onedgains(), m_layerToGroup() {};

◆ CDCDedx1DCell() [2/2]

CDCDedx1DCell	(	short	version,
		const std::vector< std::vector< double > > &	groupGains,
		const std::vector< unsigned int > &	layerToGroup )

Constructor with grouped constants.

Definition at line 14 of file CDCDedx1DCell.cc.

  : m_version(version),
    m_onedgains(groupGains),
    m_layerToGroup(layerToGroup)
{
  // 1. groups exist
  if (groupGains.empty()) {
    B2ERROR("CDCDedx1DCell: groupGains is empty.");
    return;
  }
 
  // 2. mapping exists
  if (layerToGroup.empty()) {
    B2ERROR("CDCDedx1DCell: layerToGroup is empty.");
    return;
  }
 
  // 3. check groups not empty
  for (size_t g = 0; g < groupGains.size(); ++g) {
    if (groupGains[g].empty()) {
      B2ERROR("CDCDedx1DCell: group " << g << " is empty.");
      return;
    }
  }
 
  // 4. check Layer -> group mapping
  for (size_t layer = 0; layer < layerToGroup.size(); ++layer) {
    if (layerToGroup[layer] >= groupGains.size()) {
      B2ERROR("CDCDedx1DCell: layer " << layer
              << " refers to invalid group "
              << layerToGroup[layer]);
      return;
    }
  }
 
  m_onedgains = groupGains;
  m_layerToGroup = layerToGroup;
}

◆ ~CDCDedx1DCell()

~CDCDedx1DCell ( )

inline

Destructor.

Definition at line 43 of file CDCDedx1DCell.h.

43{};

Member Function Documentation

◆ getGroup()

unsigned int getGroup ( unsigned int layer ) const

inline

Get group index for layer.

Definition at line 69 of file CDCDedx1DCell.h.

    {
      if (layer >= getNLayers()) {
        B2ERROR("Layer out of range");
        return std::numeric_limits<unsigned int>::max();
      }
 
      if (!m_layerToGroup.empty()) {
        return m_layerToGroup[layer];
      }
 
      if (m_onedgains.size() == 2) {
        return (layer < 8) ? 0 : 1;
      }
 
      B2ERROR("Invalid 1DCell payload state");
      return std::numeric_limits<unsigned int>::max();
    };

◆ getLayerMap()

const std::vector< unsigned int > & getLayerMap ( ) const

inline

Get the Layer-to-group map.

Returns: layer map

Definition at line 113 of file CDCDedx1DCell.h.

    {
      return m_layerToGroup;
    }

◆ getMean() [1/2]

double getMean	(	unsigned int	layer,
		double	enta ) const

Return dE/dx mean value for given entrance angle.

Parameters

layer	continuous layer number
enta	entrance angle (-pi/2 to pi/2)

Definition at line 113 of file CDCDedx1DCell.cc.

{
  const unsigned int group = getGroup(layer);
 
  if (group >= m_onedgains.size()) return 1.0;
  if (m_onedgains[group].empty()) return 1.0;
 
  const double piby2 = M_PI / 2.0;
 
  if (enta < -piby2) enta += piby2;
  if (enta >  piby2) enta -= piby2;
 
  unsigned int nbins = m_onedgains[group].size();
 
  double binsize = 2.0 * piby2 / nbins;
 
  int bin = std::floor((enta + piby2) / binsize);
 
  if (bin < 0 || static_cast<unsigned int>(bin) >= nbins) {
    B2WARNING("Problem with CDC dE/dx 1D binning!");
    return 1.0;
  }
 
  return m_onedgains[group][bin];
}

◆ getMean() [2/2]

double getMean	(	unsigned int	layer,
		unsigned int	bin ) const

Return dE/dx mean value for the given bin.

Parameters

layer	is layer number between 0-55
bin	is enta bin number

Definition at line 91 of file CDCDedx1DCell.cc.

{
  const unsigned int group = getGroup(layer);
  if (group >= m_onedgains.size()) return 1.0;
 
  if (bin >= m_onedgains[group].size()) return 1.0;
 
  return m_onedgains[group][bin];
}

◆ getNBins()

unsigned int getNBins ( unsigned int layer ) const

inline

Get the number of bins for the entrance angle correction.

Definition at line 90 of file CDCDedx1DCell.h.

    {
      const unsigned int group = getGroup(layer);
 
      if (group >= m_onedgains.size()) return 0;
 
      return m_onedgains[group].size();
    };

◆ getNGroups()

unsigned int getNGroups ( ) const

inline

Number of groups.

Definition at line 56 of file CDCDedx1DCell.h.

56{return m_onedgains.size();};

◆ getNLayers()

unsigned int getNLayers ( ) const

inline

Number of layers.

Definition at line 59 of file CDCDedx1DCell.h.

    {
      if (!m_layerToGroup.empty()) return m_layerToGroup.size();
 
      if (m_onedgains.size() == 2) return 56;
 
      return 0;
    };

◆ getOneDCell()

const std::vector< std::vector< double > > & getOneDCell ( ) const

inline

Get the grouped calibration constants.

Returns: grouped calibration constants

Definition at line 103 of file CDCDedx1DCell.h.

    {
      return m_onedgains;
    }

◆ getVersion()

short getVersion ( ) const

inline

Get the version for the 1D cleanup.

Definition at line 53 of file CDCDedx1DCell.h.

53{ return m_version; };

◆ isValidGroupedPayload()

bool isValidGroupedPayload ( ) const

private

Payload validation.

Definition at line 55 of file CDCDedx1DCell.cc.

{
  if (m_onedgains.empty()) {
    B2ERROR("CDCDedx1DCell: no gain groups");
    return false;
  }
 
  for (unsigned int g = 0; g < m_onedgains.size(); ++g) {
    if (m_onedgains[g].empty()) {
      B2ERROR("CDCDedx1DCell: group " << g << " is empty");
      return false;
    }
  }
 
  // new-style payload
  if (!m_layerToGroup.empty()) {
 
    for (unsigned int layer = 0; layer < m_layerToGroup.size(); ++layer) {
      if (m_layerToGroup[layer] >= m_onedgains.size()) {
        B2ERROR("1DCell: invalid group index in layer map");
        return false;
      }
    }
 
    return true;
  }
 
  // old-style payload
  if (m_onedgains.size() == 2) {
    return true;
  }
 
  B2ERROR("CDCDedx1DCell: invalid legacy payload shape");
  return false;
}

◆ operator*=()

CDCDedx1DCell & operator*= ( CDCDedx1DCell const & rhs )

Combine payloads.

Definition at line 139 of file CDCDedx1DCell.cc.

{
  if (m_version != rhs.getVersion()) {
    B2WARNING("1DCell versions do not match");
    return *this;
  }
 
  if (!isValidGroupedPayload() || !rhs.isValidGroupedPayload()) {
    B2WARNING("Invalid payload");
    return *this;
  }
 
  if (m_layerToGroup != rhs.m_layerToGroup) {
    B2WARNING("Layer mapping mismatch");
    return *this;
  }
 
  if (m_onedgains.size() != rhs.m_onedgains.size()) {
    B2WARNING("Group count mismatch");
    return *this;
  }
 
  for (unsigned int group = 0; group < m_onedgains.size(); ++group) {
 
    if (m_onedgains[group].size() != rhs.m_onedgains[group].size()) {
      B2WARNING("Bin count mismatch");
      return *this;
    }
 
    for (unsigned int bin = 0; bin < m_onedgains[group].size(); ++bin) {
      m_onedgains[group][bin] *= rhs.m_onedgains[group][bin];
    }
  }
 
  return *this;
}

◆ setMean()

void setMean	(	unsigned int	layer,
		unsigned int	bin,
		double	value )

Reset dE/dx mean value for the given bin.

Parameters

layer	is layer number between 0-55
bin	is enta bin number
value	is constant for requested entabin and layer

Definition at line 101 of file CDCDedx1DCell.cc.

{
  const unsigned int group = getGroup(layer);
 
  if (group >= m_onedgains.size()) return;
  if (bin >= m_onedgains[group].size()) return;
 
  m_onedgains[group][bin] = value;
}

Member Data Documentation

◆ m_layerToGroup

std::vector<unsigned int> m_layerToGroup

private

layer to group map

Definition at line 150 of file CDCDedx1DCell.h.

◆ m_onedgains

std::vector<std::vector<double> > m_onedgains

private

dE/dx means in entrance angle bins

Definition at line 149 of file CDCDedx1DCell.h.

◆ m_version

short m_version

private

dE/dx cleanup correction versus entrance angle may be different for different layers, so store as a vector of vectors keep a version number to identify which layers are valid

version number for 1D cleanup correction

Definition at line 148 of file CDCDedx1DCell.h.

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

cdc/dbobjects/include/CDCDedx1DCell.h
cdc/dbobjects/src/CDCDedx1DCell.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CDCDedx1DCell() [1/2]

◆ CDCDedx1DCell() [2/2]

◆ ~CDCDedx1DCell()

Member Function Documentation

◆ getGroup()

◆ getLayerMap()

◆ getMean() [1/2]

◆ getMean() [2/2]

◆ getNBins()

◆ getNGroups()

◆ getNLayers()

◆ getOneDCell()

◆ getVersion()

◆ isValidGroupedPayload()

◆ operator*=()

◆ setMean()

Member Data Documentation

◆ m_layerToGroup

◆ m_onedgains

◆ m_version