dE/dx wire gain 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 > > &onedgains)
	Constructor.

	~CDCDedx1DCell ()
	Destructor.

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

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

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

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

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
	ClassDef (CDCDedx1DCell, 5)
	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

Detailed Description

dE/dx wire gain 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() {};

◆ CDCDedx1DCell() [2/2]

CDCDedx1DCell	(	short	version,
		const std::vector< std::vector< double > > &	onedgains )

inline

Constructor.

Definition at line 38 of file CDCDedx1DCell.h.

38: m_version(version), m_onedgains(onedgains) {};

◆ ~CDCDedx1DCell()

~CDCDedx1DCell ( )

inline

Destructor.

Definition at line 43 of file CDCDedx1DCell.h.

43{};

Member Function Documentation

◆ getMean() [1/2]

double getMean	(	unsigned int	layer,
		double	enta ) const

inline

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 132 of file CDCDedx1DCell.h.

    {
      if (layer > 56) {
        B2ERROR("No such layer!");
        return 0;
      }
 
      unsigned int mylayer = 0;
      if (layer >= 8 && m_onedgains.size() == 2) mylayer = 1;
      else if (m_onedgains.size() == 56) mylayer = layer;
 
      if (mylayer >= m_onedgains.size()) return 1.0;
 
      double piby2 =  M_PI / 2.0;
      // assume rotational symmetry
      if (enta < -piby2) enta += piby2;
      if (enta > piby2) enta -= piby2;
 
      double binsize = (m_onedgains[mylayer].size() != 0) ? 2.0 * piby2 / m_onedgains[mylayer].size() : 0.0;
      int bin = (binsize != 0.0) ? std::floor((enta + piby2) / binsize) : -1;
      if (bin < 0 || (unsigned)bin >= m_onedgains[mylayer].size()) {
        B2WARNING("Problem with CDC dE/dx 1D binning!");
        return 1.0;
      }
 
      return m_onedgains[mylayer][bin];
    };

◆ getMean() [2/2]

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

inline

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 103 of file CDCDedx1DCell.h.

    {
      unsigned int mylayer = 0;
      if (layer >= 8 && m_onedgains.size() == 2) mylayer = 1;
      else if (m_onedgains.size() == 56) mylayer = layer;
 
      if (mylayer < m_onedgains.size() and bin < m_onedgains[mylayer].size()) return m_onedgains[mylayer][bin];
 
      return 1.0;
    };

◆ getNBins()

unsigned int getNBins ( unsigned int layer ) const

inline

Get the number of bins for the entrance angle correction.

Definition at line 76 of file CDCDedx1DCell.h.

    {
      if (m_onedgains.size() == 0) {
        B2ERROR("ERROR!");
        return 0;
      }
 
      if (m_onedgains.size() == 2) {
        if (layer == 0) return m_onedgains[0].size();
        else if (layer == 1) return m_onedgains[1].size();
        else {
          B2ERROR("ERROR! const vector not found");
          return 0;
        }
      } else if (m_onedgains.size() == 56) {
        if (layer < m_onedgains.size()) return m_onedgains[layer].size();
      } else {
        B2ERROR("ERROR! Wrong # of constants vector array: getNBins()");
        return 0;
      }
      return 0;
    };

◆ getSize()

unsigned int getSize ( ) const

inline

Get the number of bins for the entrance angle correction.

Definition at line 72 of file CDCDedx1DCell.h.

72{ return m_onedgains.size(); };

◆ getVersion()

short getVersion ( ) const

inline

Get the version for the 1D cleanup.

Definition at line 68 of file CDCDedx1DCell.h.

68{ return m_version; };

◆ operator*=()

CDCDedx1DCell & operator*= ( CDCDedx1DCell const & rhs )

inline

Combine payloads.

Definition at line 48 of file CDCDedx1DCell.h.

    {
      if (m_version != rhs.getVersion()) {
        B2WARNING("1D cell gain parameters do not match, cannot merge!");
        return *this;
      }
      for (unsigned int layer = 0; layer < getSize(); ++layer) {
        double tempLayer = -1; //out of bound layer
        if (getSize() == 2)tempLayer = layer * 20; // 20 to make <8 and >=8 difference
        else if (getSize() == 56)tempLayer = layer;
        else B2ERROR("ERROR! Wrong # of constants vector array");
        for (unsigned int bin = 0; bin < getNBins(layer); ++bin) {
          m_onedgains[layer][bin] *= rhs.getMean(tempLayer, bin);
        }
      }
      return *this;
    }

◆ setMean()

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

inline

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 119 of file CDCDedx1DCell.h.

    {
      unsigned int mylayer = 0;
      if (layer >= 8 && m_onedgains.size() == 2) mylayer = 1;
      else if (m_onedgains.size() == 56) mylayer = layer;
 
      if (mylayer < m_onedgains.size() and bin < m_onedgains[mylayer].size()) m_onedgains[mylayer][bin] = value;
    };

Member Data Documentation

◆ m_onedgains

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

private

dE/dx means in entrance angle bins

Definition at line 165 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 164 of file CDCDedx1DCell.h.

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

cdc/dbobjects/include/CDCDedx1DCell.h

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CDCDedx1DCell() [1/2]

◆ CDCDedx1DCell() [2/2]

◆ ~CDCDedx1DCell()

Member Function Documentation

◆ getMean() [1/2]

◆ getMean() [2/2]

◆ getNBins()

◆ getSize()

◆ getVersion()

◆ operator*=()

◆ setMean()

Member Data Documentation

◆ m_onedgains

◆ m_version