Class to hold the prediction of mean as a function of beta-gamma (bg) More...

#include <CDCDedxMeanPred.h>

Public Member Functions
std::vector< double >	getMeanVector () const
	Return the mean vector from payload.

double	getMean (double bg)
	Return the predicted mean value as a function of beta-gamma (bg)

double	meanCurve (double x, double par, int version) const
	beta-gamma (bg) curve function

Private Attributes
std::vector< double >	m_meanpars
	dE/dx mean parameters

const DBObjPtr< CDCDedxMeanPars >	m_DBMeanPars
	db object for dE/dx mean parameters

Detailed Description

Class to hold the prediction of mean as a function of beta-gamma (bg)

Definition at line 26 of file CDCDedxMeanPred.h.

Member Function Documentation

◆ getMean()

double getMean ( double bg )

Return the predicted mean value as a function of beta-gamma (bg)

Definition at line 32 of file CDCDedxMeanPred.cc.

{
  m_meanpars = getMeanVector();
 
  // define the section of the mean to use
  double A = 0, B = 0, C = 0;
  if (bg < 4.5)
    A = 1;
  else if (bg < 10)
    B = 1;
  else
    C = 1;
 
  double x[1]; x[0] = bg;
  double parsA[9];
  double parsB[5];
  double parsC[5];
 
  parsA[0] = 1; parsB[0] = 2; parsC[0] = 3;
  for (int i = 0; i < 15; ++i) {
    if (i < 7) parsA[i + 1] = m_meanpars[i];
    else if (i < 11) parsB[i % 7 + 1] = m_meanpars[i];
    else parsC[i % 11 + 1] = m_meanpars[i];
  }
 
  // calculate dE/dx from the Bethe-Bloch mean
  double partA = meanCurve(x, parsA, 0);
  double partB = meanCurve(x, parsB, 0);
  double partC = meanCurve(x, parsC, 0);
 
  return (A * partA + B * partB + C * partC);
}

◆ getMeanVector()

std::vector< double > getMeanVector ( ) const

inline

Return the mean vector from payload.

Definition at line 33 of file CDCDedxMeanPred.h.

    {
      if (!m_DBMeanPars || m_DBMeanPars->getSize() == 0) {
        B2WARNING("No dE/dx mean parameters!");
        std::vector<double> meanpar;
        for (int i = 0; i < 15; ++i)
          meanpar.push_back(1.0);
        return meanpar;
      } else
        return m_DBMeanPars->getMeanPars();
    }

◆ meanCurve()

double meanCurve	(	double *	x,
		double *	par,
		int	version
	)		const

beta-gamma (bg) curve function

Definition at line 13 of file CDCDedxMeanPred.cc.

{
  // calculate the predicted mean value as a function of beta-gamma (bg)
  // this is done with a different function depending on the value of bg
  double f = 0;
 
  if (version == 0) {
    if (par[0] == 1)
      f = par[1] * std::pow(std::sqrt(x[0] * x[0] + 1), par[3]) / std::pow(x[0], par[3]) *
          (par[2] - par[5] * std::log(1 / x[0])) - par[4] + std::exp(par[6] + par[7] * x[0]);
    else if (par[0] == 2)
      f = par[1] * std::pow(x[0], 3) + par[2] * x[0] * x[0] + par[3] * x[0] + par[4];
    else if (par[0] == 3)
      f = -1.0 * par[1] * std::log(par[4] + std::pow(1 / x[0], par[2])) + par[3];
  }
 
  return f;
}

Member Data Documentation

◆ m_DBMeanPars

const DBObjPtr<CDCDedxMeanPars> m_DBMeanPars

private

db object for dE/dx mean parameters

Definition at line 60 of file CDCDedxMeanPred.h.

◆ m_meanpars

std::vector<double> m_meanpars

private

dE/dx mean parameters

Definition at line 58 of file CDCDedxMeanPred.h.

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

reconstruction/utility/include/CDCDedxMeanPred.h
reconstruction/utility/src/CDCDedxMeanPred.cc

Public Member Functions

Private Attributes