CDCDedxCosineEdge Class Reference

dE/dx special large cosine calibration to fix bending shoulder at large costh More...

#include <CDCDedxCosineEdge.h>

Inheritance diagram for CDCDedxCosineEdge:

Public Member Functions
	CDCDedxCosineEdge ()
	Default constructor.
	CDCDedxCosineEdge (const std::vector< std::vector< double > > &largecosth)
	Constructor.
	~CDCDedxCosineEdge ()
	Destructor.
int	getSize (const int side) const
	Get the number of bins of requested side.
double	getMean (int side, unsigned int ibin) const
	Return calibration constant for given side and bin #.
double	getMean (double costh) const
	Return calibration constant for cosine value.
double	getCosEdgePar (int side, unsigned int ibin) const
	Return specific large cosine constants on give side.
void	printCosEdgePars (int side)
	Print large cosine constants array on requested side.
void	setCosthEdgePar (int side, unsigned int ibin, double value)
	Set specific hadron parameter.

Private Member Functions
	ClassDef (CDCDedxCosineEdge, 1)
	ClassDef.

Private Attributes
std::vector< std::vector< double > >	m_largeCos
	large cosine calibration constants

Detailed Description

dE/dx special large cosine calibration to fix bending shoulder at large costh

Definition at line 30 of file CDCDedxCosineEdge.h.

Constructor & Destructor Documentation

CDCDedxCosineEdge() [1/2]

CDCDedxCosineEdge ( )

inline

Default constructor.

Definition at line 36 of file CDCDedxCosineEdge.h.

: m_largeCos() {};

CDCDedxCosineEdge() [2/2]

CDCDedxCosineEdge ( const std::vector< std::vector< double > > & largecosth )

inlineexplicit

Constructor.

Parameters

largecosth

calibration constants

Definition at line 42 of file CDCDedxCosineEdge.h.

: m_largeCos(largecosth) {};

~CDCDedxCosineEdge()

~CDCDedxCosineEdge ( )

inline

Destructor.

Definition at line 47 of file CDCDedxCosineEdge.h.

{};

Member Function Documentation

getCosEdgePar()

double getCosEdgePar ( int side, unsigned int ibin ) const

inline

Return specific large cosine constants on give side.

Parameters

side	<0 for backward and >0 for forward
ibin	is bin number for given costh

Definition at line 146 of file CDCDedxCosineEdge.h.

    {
 
      std::vector<double> temp;
      double coslow = 0.0, coshigh = 0.0;
      if (side < 0) {
        if (m_largeCos.size() > 0) temp = m_largeCos[0];
        coslow = -0.870; coshigh = -0.850;
      } else if (side > 0) {
        if (m_largeCos.size() > 1) temp = m_largeCos[1];
        coslow = 0.950; coshigh = 0.960;
      } else {
        B2ERROR("CDCDedxCosineEdge:choose > 0 for forward and <0 for backward side");
        return  -99.0;
      }
 
      if (ibin >= temp.size()) {
        B2ERROR("CDCDedxCosineEdge:Problem with bin index: choose 0 and " << temp.size() - 1); //
        return -99.0;
      }
 
      if (temp.size() == 0) return -99.0;
      double bw = std::abs(coshigh - coslow) / temp.size();
      double bc = coslow  + (0.5 + ibin) * bw; //bin centre
      std::cout << "Par # " << ibin << ", costh bin centre = " << bc << ", const =" << temp[ibin] << std::endl;
      return temp[ibin];
    };

getMean() [1/2]

double getMean ( double costh ) const

inline

Return calibration constant for cosine value.

Parameters

costh

is a signed cosine theta value so by definition it recog it's side

Returns

calibration constant

Definition at line 91 of file CDCDedxCosineEdge.h.

    {
      if (m_largeCos.size() < 2) {
        B2ERROR("CDCDedxCosineEdge: vector-of-vectors too short, returning 1");
        return 1.0;
      }
 
      const std::vector<double>& temp = costh < 0 ?  m_largeCos[0] : m_largeCos[1];
      if (temp.size() < 2) {
        B2ERROR("CDCDedxCosineEdge: vector too short, returning 1");
        return 1.0;
      }
 
      double coslow = 0.0, coshigh = 0.0;
      if (costh < 0) {
        coslow = -0.870; coshigh = -0.850; //this is hardcoded and fixed
      } else {
        coslow = 0.950; coshigh = 0.960; //this is hardcoded and fixed
      }
 
      //don't do anything for other cosine range
      if (costh < coslow or costh > coshigh) {
        return 1.0;
      }
 
      if (costh <= -0.866 or costh >= 0.9575) return 1.0;
 
      // gains are stored at the center of the bins
      // find the bin center immediately preceding this value of costh
      int nbins = int(temp.size());
      double binsize = (coshigh - coslow) / nbins;
      int bin = std::floor((costh - 0.5 * binsize - coslow) / binsize);
      if (bin < 0) bin = 0;
      else if (bin > nbins - 2) bin = nbins - 2;
 
      // extrapolate backward for lowest half-bin and center positive half-bin
      // extrapolate forward for highest half-bin and center negative half-bin
      // MS: not clear why this is needed
      if (temp[bin + 1] - temp[bin] < -0.6) bin++;
      else if (temp[bin + 1] - temp[bin] > 0.6) bin--;
 
      if (bin < 0 or bin + 1 >= nbins) {
        B2WARNING("CDCDedxCosineEdge:no constants for costh: " << costh << " as it is not in range");
        return 1.0;
      }
 
      double frac = ((costh - 0.5 * binsize - coslow) / binsize) - bin;
      return ((temp[bin + 1] - temp[bin]) * frac + temp[bin]);
    }

getMean() [2/2]

double getMean ( int side, unsigned int ibin ) const

inline

Return calibration constant for given side and bin #.

Parameters

side	<0 for backward and >0 for forward
ibin	is bin number for given costh

Returns

calibration constant

Definition at line 72 of file CDCDedxCosineEdge.h.

    {
      if (m_largeCos.size() < 2) {
        B2ERROR("CDCDedxCosineEdge: vector-of-vectors too short, returning 1");
        return 1.0;
      }
 
      const std::vector<double>& temp = side < 0 ?  m_largeCos[0] : m_largeCos[1];
      if (ibin < temp.size()) return temp[ibin];
 
      B2WARNING("CDCDedxCosineEdge: invalid bin number, returning 1");
      return 1.0;
    };

getSize()

int getSize ( const int side ) const

inline

Get the number of bins of requested side.

Parameters

side	<0 for backward and >0 for forward

Returns

number of bins

Definition at line 54 of file CDCDedxCosineEdge.h.

    {
      if (side < 0) {
        if (m_largeCos.empty()) return 0;
        return m_largeCos[0].size();
      } else {
        if (m_largeCos.size() < 2) return 0;
        return m_largeCos[1].size();
      }
    };

printCosEdgePars()

void printCosEdgePars ( int side )

inline

Print large cosine constants array on requested side.

Parameters

side	<0 for backward and >0 for forward constants

Definition at line 179 of file CDCDedxCosineEdge.h.

    {
      std::vector<double> temp;
      double coslow = 0.0, coshigh = 0.0;
      if (side < 0) {
        if (m_largeCos.size() > 0) temp = m_largeCos[0];
        coslow = -0.870; coshigh = -0.850;
      } else if (side > 0) {
        if (m_largeCos.size() > 1) temp = m_largeCos[1];
        coslow = 0.950; coshigh = 0.960;
      } else {
        B2ERROR("CDCDedxCosineEdge:choose > 0 for forward and <0 for backward side");
        return;
      }
 
      if (temp.size() == 0) return;
      double bw = std::abs(coshigh - coslow) / temp.size();
      B2INFO("Printing parameters (0=backward and 1=forward): " << side << ", nPars = " << temp.size());
      for (unsigned int ibin = 0; ibin < temp.size(); ibin++) {
        double bc = coslow  + (0.5 + ibin) * bw; //bin centre
        std::cout << "Par # " << ibin << ", costh bin centre = " << bc << ", const = " << temp[ibin] << std::endl;
      }
    };

setCosthEdgePar()

void setCosthEdgePar ( int side, unsigned int ibin, double value )

inline

Set specific hadron parameter.

Parameters

side	0 for backward and 1 for forward
ibin	number starts from 0
value	of parameter to set

Definition at line 210 of file CDCDedxCosineEdge.h.

    {
      if (m_largeCos.size() < 2) {
        B2ERROR("CDCDedxCosineEdge: vector-of-vectors too short, value not set");
        return;
      }
 
      int iside = side < 0 ? 0 : 1;
      if (ibin >= m_largeCos[iside].size()) {
        B2ERROR("CDCDedxCosineEdge: invalid bin number, value not set");
        return;
      }
 
      m_largeCos[iside][ibin] = value;
      B2DEBUG(20, "Par # " << ibin << ", const = " << m_largeCos[iside][ibin]);
    };

Member Data Documentation

m_largeCos

std::vector<std::vector<double> > m_largeCos

private

large cosine calibration constants

Definition at line 229 of file CDCDedxCosineEdge.h.

---

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

• cdc/dbobjects/include/CDCDedxCosineEdge.h