Class to hold the prediction of resolution depending dE/dx, nhit, and cos(theta) More...

#include <CDCDedxSigmaPred.h>

Collaboration diagram for CDCDedxSigmaPred:

[legend]

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	getSigma (double dedx, double nhit, double cos, double timereso)
	Return the predicted resolution depending on dE/dx, nhit, and cos(theta)

double	nhitPrediction (double nhit)
	Return sigma from the nhit parameterization.

double	ionzPrediction (double dedx)
	Return sigma from the ionization parameterization.

double	cosPrediction (double cos)
	Return sigma from the cos parameterization.

double	getDedxPars (int i)
	get the dedx parameters

void	setDedxPars (int i, double val)
	set the dedx parameters

double	getNHitPars (int i)
	get the nhit parameters

void	setNHitPars (int i, double val)
	set the nhit parameters

double	getCosPars (int i)
	get the cos(theta) parameters

void	setCosPars (int i, double val)
	set the cos(theta) parameters

Private Attributes
double	m_dedxpars [2]
	parameters for sigma vs.

double	m_cospars [10]
	parameters for sigma vs.

double	m_nhitpars [5]
	parameters for sigma vs.

const DBObjPtr< CDCDedxSigmaPars >	m_DBSigmaPars
	db object for dE/dx resolution parameters

Detailed Description

Class to hold the prediction of resolution depending dE/dx, nhit, and cos(theta)

Definition at line 36 of file CDCDedxSigmaPred.h.

Member Function Documentation

◆ cosPrediction()

double cosPrediction ( double cos )

Return sigma from the cos parameterization.

Definition at line 141 of file CDCDedxSigmaPred.cc.

{
  // Create an instance of the CDCDedxWidgetSigma class
  CDCDedxWidgetSigma gs;
 
  // Define the parameter array and initialize it
  double cospar[11];
  cospar[0] = 3;
  for (int i = 0; i < 10; ++i)  cospar[i + 1] = m_cospars[i];
 
  // Initialize x array to pass to sigmaCurve
  double x[1];
  x[0] = cos;
 
  double corCos;
  corCos = gs.sigmaCurve(x, std::vector<double>(cospar, cospar + 11));
 
  return corCos;
}

◆ getCosPars()

double getCosPars ( int i )

inline

get the cos(theta) parameters

Definition at line 98 of file CDCDedxSigmaPred.h.

98{ return m_cospars[i]; };

◆ getDedxPars()

double getDedxPars ( int i )

inline

get the dedx parameters

Definition at line 78 of file CDCDedxSigmaPred.h.

78{ return m_dedxpars[i]; };

◆ getNHitPars()

double getNHitPars ( int i )

inline

get the nhit parameters

Definition at line 88 of file CDCDedxSigmaPred.h.

88{ return m_nhitpars[i]; };

◆ getSigma()

double getSigma	(	double	dedx,
		double	nhit,
		double	cos,
		double	timereso )

Return the predicted resolution depending on dE/dx, nhit, and cos(theta)

Definition at line 81 of file CDCDedxSigmaPred.cc.

{
  double correction  = cosPrediction(cos) * nhitPrediction(nhit) * ionzPrediction(dedx) * timereso;
  return correction;
}

◆ ionzPrediction()

double ionzPrediction ( double dedx )

Return sigma from the ionization parameterization.

Definition at line 122 of file CDCDedxSigmaPred.cc.

{
  // Define the parameter array and initialize it
  double dedxpar[3];
  dedxpar[0] = 1;
  for (int i = 0; i < 2; ++i) dedxpar[i + 1] = m_dedxpars[i];
 
  // Create an instance of the CDCDedxWidgetSigma class
  CDCDedxWidgetSigma gs;
 
// Initialize x array to pass to sigmaCurve
  double x[1];
  x[0] = dedx;
 
  double corDedx = gs.sigmaCurve(x, std::vector<double>(dedxpar, dedxpar + 3));
 
  return corDedx;
}

◆ nhitPrediction()

double nhitPrediction ( double nhit )

Return sigma from the nhit parameterization.

Definition at line 88 of file CDCDedxSigmaPred.cc.

{
 
// Define minimum and maximum nhit values
  int nhit_min = 8, nhit_max = 37;
 
  // Define the parameter array and initialize it
  double nhitpar[6];
  nhitpar[0] = 2;
  for (int i = 0; i < 5; ++i) nhitpar[i + 1] = m_nhitpars[i];
 
  // Create an instance of the CDCDedxWidgetSigma class
  CDCDedxWidgetSigma gs;
 
  // Initialize x array to pass to sigmaCurve
  double x[1];
 
  x[0] = nhit;
 
  double corNHit;
 
  if (nhit < nhit_min) {
    x[0] = nhit_min;
    corNHit = gs.sigmaCurve(x, std::vector<double>(nhitpar, nhitpar + 6)) * std::sqrt(nhit_min / nhit);
  } else if (nhit > nhit_max) {
    x[0] = nhit_max;
    corNHit = gs.sigmaCurve(x, std::vector<double>(nhitpar, nhitpar + 6)) * std::sqrt(nhit_max / nhit);
  } else {
    corNHit = gs.sigmaCurve(x, std::vector<double>(nhitpar, nhitpar + 6));
  }
 
  return corNHit;
}

◆ printParameters()

void printParameters ( std::string infile )

write the parameters in file

Definition at line 60 of file CDCDedxSigmaPred.cc.

{
 
  B2INFO("\n\tCDCDedxSigmaPred: Printing parameters to file --> " << outfile.c_str());
 
  // write out the parameters to file
  std::ofstream fout(outfile.c_str());
 
  for (int i = 1; i < 3; ++i) fout << i << "\t" << m_dedxpars[i - 1] << std::endl;
 
  fout << std::endl;
  for (int i = 3; i < 8; ++i) fout << i << "\t" << m_nhitpars[i - 3] << std::endl;
 
  fout << std::endl;
  for (int i = 8; i < 18; ++i) fout << i << "\t" << m_cospars[i - 8] << std::endl;
 
  fout << std::endl;
 
  fout.close();
}

◆ setCosPars()

void setCosPars	(	int	i,
		double	val )

inline

set the cos(theta) parameters

Definition at line 103 of file CDCDedxSigmaPred.h.

103{ m_cospars[i] = val; };

◆ setDedxPars()

void setDedxPars	(	int	i,
		double	val )

inline

set the dedx parameters

Definition at line 83 of file CDCDedxSigmaPred.h.

83{ m_dedxpars[i] = val; };

◆ setNHitPars()

void setNHitPars	(	int	i,
		double	val )

inline

set the nhit parameters

Definition at line 93 of file CDCDedxSigmaPred.h.

93{ m_nhitpars[i] = val; };

◆ setParameters() [1/2]

void setParameters ( )

set the parameters

Definition at line 13 of file CDCDedxSigmaPred.cc.

{
 
  // make sure the resolution parameters are reasonable
  if (!m_DBSigmaPars || m_DBSigmaPars->getSize() == 0)
    B2FATAL("No dE/dx sigma parameters!");
 
  std::vector<double> sigmapar;
  sigmapar = m_DBSigmaPars->getSigmaPars();
  for (int i = 0; i < 2; ++i)  m_dedxpars[i] = sigmapar[i];
  for (int i = 0; i < 5; ++i) m_nhitpars[i] = sigmapar[i + 2];
  for (int i = 0; i < 10; ++i) m_cospars[i] = sigmapar[i + 7];
}

◆ setParameters() [2/2]

void setParameters ( std::string infile )

set the parameters from file

Definition at line 27 of file CDCDedxSigmaPred.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 --> dedx parameters");
  for (int i = 0; i < 2; ++i) {
    fin >> par >> m_dedxpars[i];
    B2INFO("\t\t (" << i << ")" << m_dedxpars[i]);
  }
 
  B2INFO("\t --> nhit parameters");
  for (int i = 0; i <= 4; ++i) {
    fin >> par >> m_nhitpars[i];
    B2INFO("\t\t (" << i << ")" << m_nhitpars[i]);
  }
 
  B2INFO("\t --> cos parameters");
  for (int i = 0; i < 10; ++i) {   //4
    fin >> par >> m_cospars[i];
    B2INFO("\t\t (" << i << ")" << m_cospars[i]);
  }
 
  fin.close();
}

Member Data Documentation

◆ m_cospars

double m_cospars[10]

private

parameters for sigma vs.

cos(theta) curve

Definition at line 108 of file CDCDedxSigmaPred.h.

◆ m_DBSigmaPars

const DBObjPtr<CDCDedxSigmaPars> m_DBSigmaPars

private

db object for dE/dx resolution parameters

Definition at line 111 of file CDCDedxSigmaPred.h.

◆ m_dedxpars

double m_dedxpars[2]

private

parameters for sigma vs.

dE/dx curve

Definition at line 107 of file CDCDedxSigmaPred.h.

◆ m_nhitpars

double m_nhitpars[5]

private

parameters for sigma vs.

nhit curve

Definition at line 109 of file CDCDedxSigmaPred.h.

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

cdc/utilities/include/CDCDedxSigmaPred.h
cdc/utilities/src/CDCDedxSigmaPred.cc

Public Member Functions

Private Attributes

Detailed Description

Member Function Documentation

◆ cosPrediction()

◆ getCosPars()

◆ getDedxPars()

◆ getNHitPars()

◆ getSigma()

◆ ionzPrediction()

◆ nhitPrediction()

◆ printParameters()

◆ setCosPars()

◆ setDedxPars()

◆ setNHitPars()

◆ setParameters() [1/2]

◆ setParameters() [2/2]

Member Data Documentation

◆ m_cospars

◆ m_DBSigmaPars

◆ m_dedxpars

◆ m_nhitpars