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

#include <CDCDedxMeanPred.h>

Public Member Functions
void	setParameters (std::string infile)
	set the parameters from file

void	setParameters ()
	set the parameters

void	printParameters (std::string infile)
	write the parameters in file

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

double	getCurvePars (int i)
	get the curve parameters

void	setCurvePars (int i, double val)
	set the curve parameters

Private Attributes
double	m_meanpars [15]
	parameters for beta-gamma curve

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 31 of file CDCDedxMeanPred.h.

Member Function Documentation

◆ getCurvePars()

double getCurvePars ( int i )

inline

get the curve parameters

Definition at line 58 of file CDCDedxMeanPred.h.

58{ return m_meanpars[i]; };

Belle2::CDCDedxMeanPred::m_meanpars

double m_meanpars[15]

parameters for beta-gamma curve

Definition: CDCDedxMeanPred.h:67

◆ getMean()

double getMean ( double bg )

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

Definition at line 61 of file CDCDedxMeanPred.cc.

{
 
  // 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[8];
  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
  CDCDedxWidgetCurve gc;
  double partA = gc.meanCurve(x, std::vector<double>(parsA, parsA + 8));
  double partB = gc.meanCurve(x, std::vector<double>(parsB, parsB + 5));
  double partC = gc.meanCurve(x, std::vector<double>(parsC, parsC + 5));
 
  return (A * partA + B * partB + C * partC);
}

◆ printParameters()

void printParameters ( std::string infile )

write the parameters in file

Definition at line 49 of file CDCDedxMeanPred.cc.

{
 
  B2INFO("\tCDCDedxMeanPred: Printing parameters to file --> " << outfile.c_str());
 
  // write out the parameters to file
  std::ofstream fout(outfile.c_str());
 
  for (int i = 1; i < 16; ++i) fout << i << "\t" << m_meanpars[i - 1] << std::endl;
  fout.close();
}

◆ setCurvePars()

void setCurvePars	(	int	i,
		double	val
	)

inline

set the curve parameters

Definition at line 63 of file CDCDedxMeanPred.h.

63{ m_meanpars[i] = val; };

◆ setParameters() [1/2]

void setParameters ( )

set the parameters

Definition at line 13 of file CDCDedxMeanPred.cc.

{
 
  // make sure the resolution parameters are reasonable
  if (!m_DBMeanPars || m_DBMeanPars->getSize() == 0)
    B2FATAL("No dE/dx mean parameters!");
 
  std::vector<double> meanpar;
  meanpar = m_DBMeanPars->getMeanPars();
 
  for (int i = 0; i < 15; ++i)
    m_meanpars[i] = meanpar[i];
}

◆ setParameters() [2/2]

void setParameters ( std::string infile )

set the parameters from file

Definition at line 28 of file CDCDedxMeanPred.cc.

{
 
  B2INFO("\n\tWidgetParameterization: Using parameters from file --> " << infile);
 
  std::ifstream fin;
  fin.open(infile.c_str());
 
  if (!fin.good()) B2FATAL("\tWARNING: CANNOT FIND" << infile);
 
  int par;
 
  B2INFO("\t --> Curve parameters");
  for (int i = 0; i < 15; ++i) {
    fin >> par >> m_meanpars[i];
    B2INFO("\t\t (" << i << ")" << m_meanpars[i]);
  }
 
  fin.close();
}

Member Data Documentation

◆ m_DBMeanPars

const DBObjPtr<CDCDedxMeanPars> m_DBMeanPars

private

db object for dE/dx mean parameters

Definition at line 69 of file CDCDedxMeanPred.h.

◆ m_meanpars