#include <PixeldEdxData.h>

Collaboration diagram for PixeldEdxData:

Public Member Functions
void	setPar (const int i, const double param)
void	setPosNeg (const bool posneg)
void	setFunctionType (const std::string &fun)
void	setBetheBlochType (const std::string &bb)
void	setMinimumdEdxForMass (const double mindedxMass)
double	dEdxPdf (const double dedx, const double signedP, const double mass, const std::array< double, 9 > &par, const int offset) const
double	fdEdxZero (const double x, const std::array< double, 9 > &par) const
double	getPar (const int i, const int j) const
std::array< double, 3 >	getP (const double dedxArg, const double signedP, const int nGoodPixels) const
std::array< double, 9 >	getFirstNPar (const double p, const int nGoodPixels, const int np) const
double	getMass (const double dedx, const double signedP, const int nGoodPixels) const
double	getdEdx (const double p, const double mass, const int nGoodPixels) const
double	getdEdx (const double p, const double mass, const std::array< double, 9 > &par) const
std::vector< float >	getLikelihoods (const double dedx2, const double p2, const int nGoodPixels) const

Static Public Member Functions
static double	crystalBall (const double x, const double x0, const double sig, const double alp, const double n)
static double	asymGaus (const double x, const double x0, const double sig, const double asym)
static double	moyal (const double x, const double Ep, const double R)
static double	dEdx_5p_BG_aleph (const double xbg, const std::array< double, 9 > &pp)
static double	dEdx_5p_aleph (const double p, double mass, const std::array< double, 9 > &pp)
static double	dEdx_5p_BG (const double xbg, const std::array< double, 9 > &pp)
static double	dEdx_5p (const double p, const double mass, const std::array< double, 9 > &pp)
static double	dEdx_BG (const double xbg, const std::array< double, 9 > &pp)
static double	dEdx_def (const double p, const double mass, const std::array< double, 9 > &pp)
static double	dEdx_3p (const double p, const double mass, const std::array< double, 9 > &pp)

Private Attributes
std::unordered_map< uint32_t, std::vector< double > >	m_par
bool	m_posneg {}
std::string	m_fun
std::string	m_bb
double	m_mindedxformass {}

Static Private Attributes
static constexpr double	m_piMass = 0.13957
static constexpr double	m_kMass = 0.49368
static constexpr double	m_pMass = 0.93827

Detailed Description

Definition at line 23 of file PixeldEdxData.h.

Member Function Documentation

◆ asymGaus()

double PixeldEdxData::asymGaus	(	const double	x,
		const double	x0,
		const double	sig,
		const double	asym )

static

Definition at line 164 of file PixeldEdxData.cxx.

                                                                                                {
  double sigp = sig*asym;
  double sigm = sig/asym;
  double fun;
  if (x>x0) { fun = TMath::Gaus(x,x0,sigp,false); } 
  else      { fun = TMath::Gaus(x,x0,sigm,false); }
  return fun/(std::sqrt(2*M_PI)*((sigp*0.5) + (sigm*0.5)));
}

◆ crystalBall()

double PixeldEdxData::crystalBall	(	const double	x,
		const double	x0,
		const double	sig,
		const double	alp,
		const double	n )

static

Definition at line 152 of file PixeldEdxData.cxx.

                                                                                                                 {
  double A = std::pow(n/std::abs(alp),n)*std::exp(-std::pow(std::abs(alp),2)*0.5);
  double B = n/std::abs(alp)-std::abs(alp);
  double pull = (x-x0)/sig;
  double N = sig*(n/alp*(1/(n-1))*std::exp(-alp*alp*0.5)+std::sqrt(M_PI_2)*(1+std::erf(alp*M_SQRT1_2)));
  if (pull<alp){
    return std::exp(-pull*pull*0.5)/N;
  }
  return A*std::pow(B+pull,-n)/N;
}

◆ dEdx_3p()

double PixeldEdxData::dEdx_3p	(	const double	p,
		const double	mass,
		const std::array< double, 9 > &	pp )

static

Definition at line 216 of file PixeldEdxData.cxx.

                                                                                            {
  double bg = p/mass;
  return dEdx_BG(bg,pp);
}

◆ dEdx_5p()

double PixeldEdxData::dEdx_5p	(	const double	p,
		const double	mass,
		const std::array< double, 9 > &	pp )

static

Definition at line 197 of file PixeldEdxData.cxx.

                                                                                            {
  double bg = p/mass;
  return dEdx_5p_BG(bg,pp);
}

◆ dEdx_5p_aleph()

double PixeldEdxData::dEdx_5p_aleph	(	const double	p,
		double	mass,
		const std::array< double, 9 > &	pp )

static

Definition at line 185 of file PixeldEdxData.cxx.

                                                                                                  {
  double bg = p/mass;
  return dEdx_5p_BG_aleph(bg,pp);
}

◆ dEdx_5p_BG()

double PixeldEdxData::dEdx_5p_BG	(	const double	xbg,
		const std::array< double, 9 > &	pp )

static

Definition at line 190 of file PixeldEdxData.cxx.

                                                                               {
  double gamma   = std::sqrt(xbg*xbg+1);
  double beta    = std::sqrt(1.-1./(gamma*gamma));
  double fdedx   = pp[0]/std::pow(beta,pp[2])*std::log(1+std::pow(std::abs(pp[1])*beta*gamma,pp[4]))-pp[3];
  return fdedx;
}

◆ dEdx_5p_BG_aleph()

double PixeldEdxData::dEdx_5p_BG_aleph	(	const double	xbg,
		const std::array< double, 9 > &	pp )

static

Definition at line 179 of file PixeldEdxData.cxx.

                                                                                     {
  double gamma = std::sqrt(xbg*xbg+1.);
  double beta  = std::sqrt(1.-1./(gamma*gamma));
  return pp[0]/std::pow(beta,pp[3])*(pp[1]-std::pow(beta,pp[3])-std::log(pp[2]+1./std::pow(beta*gamma,pp[4])));
}

◆ dEdx_BG()

double PixeldEdxData::dEdx_BG	(	const double	xbg,
		const std::array< double, 9 > &	pp )

static

Definition at line 202 of file PixeldEdxData.cxx.

                                                                            {
  double gamma     = std::sqrt(xbg*xbg+1);
  double beta      = std::sqrt(1-1/(gamma*gamma));
  double beta2     = beta*beta;
  double gamma2    = gamma*gamma;
  double fdedx = (pp[0]/std::pow(beta2,pp[2]))*(std::log(std::abs(pp[1])*beta2*gamma2))-pp[0];
  return fdedx;
}

◆ dEdx_def()

double PixeldEdxData::dEdx_def	(	const double	p,
		const double	mass,
		const std::array< double, 9 > &	pp )

static

Definition at line 211 of file PixeldEdxData.cxx.

                                                                                             {
  double bg = p/mass;
  return dEdx_BG(bg,pp);
}

◆ dEdxPdf()

double PixeldEdxData::dEdxPdf	(	const double	dedx,
		const double	signedP,
		const double	mass,
		const std::array< double, 9 > &	par,
		const int	offset ) const

Definition at line 30 of file PixeldEdxData.cxx.

                                                                                                                                             {
  const auto p = std::abs(signedP);
 
  double x0 = getdEdx(p,mass,par);
 
  double pLG[4]={1,1,1,1};
 
  if (m_fun=="AG") {
    x0     = std::log(x0)   - 9.4;
    pLG[1] = x0;
    pLG[0] = par[offset]+par[offset+2];
    pLG[3] = par[offset+1]+par[offset+3];
    return asymGaus(dedx, pLG[1], pLG[0], pLG[3]);
  } 
  else if (m_fun=="CB") {
    x0     = std::log(x0)   - 9.4;
    pLG[0] = x0;
    pLG[1] = par[offset];
    pLG[2] = par[offset+1];
    pLG[3] = par[offset+2];
    return crystalBall(dedx, pLG[0], pLG[1], pLG[2], pLG[3]);
  } 
  return 0;
}

◆ fdEdxZero()

double PixeldEdxData::fdEdxZero	(	const double	x,
		const std::array< double, 9 > &	par ) const

Definition at line 55 of file PixeldEdxData.cxx.

                                                                                    {
  return getdEdx(par[6],x,par)-par[5];
}

◆ getdEdx() [1/2]

double PixeldEdxData::getdEdx	(	const double	p,
		const double	mass,
		const int	nGoodPixels ) const

Definition at line 133 of file PixeldEdxData.cxx.

                                                                                                  {
  if (nGoodPixels<=0) { return -2; }
  const auto & par = getFirstNPar(signedP,nGoodPixels,5);
  const auto p = std::abs(signedP);
  if      (m_bb=="3p")       { return dEdx_3p(p,mass,par); } 
  else if (m_bb=="5p")       { return dEdx_5p(p,mass,par); } 
  else if (m_bb=="5p_aleph") { return dEdx_5p_aleph(p,mass,par); }
  return 0;
}

◆ getdEdx() [2/2]

double PixeldEdxData::getdEdx	(	const double	p,
		const double	mass,
		const std::array< double, 9 > &	par ) const

Definition at line 143 of file PixeldEdxData.cxx.

                                                                                                           {
  const auto p = std::abs(signedP);
  if      (m_bb=="3p")       { return dEdx_3p(p,mass,par); } 
  else if (m_bb=="5p")       { return dEdx_5p(p,mass,par); } 
  else if (m_bb=="5p_aleph") { return dEdx_5p_aleph(p,mass,par); }
  return 0;
}

◆ getFirstNPar()

std::array< double, 9 > PixeldEdxData::getFirstNPar	(	const double	p,
		const int	nGoodPixels,
		const int	np ) const

Definition at line 93 of file PixeldEdxData.cxx.

                                                                                       {
  std::array<double,9> par{};
  auto setNGoodPixels{nGoodPixels};
  if (setNGoodPixels>5) { setNGoodPixels=5; }
  if (p<0 && m_posneg) { setNGoodPixels += 5; }
  for (int i=0;i<npar;i++){
    par[i] = getPar(setNGoodPixels-1,i);
  }
  return par;
}

◆ getLikelihoods()

std::vector< float > PixeldEdxData::getLikelihoods	(	const double	dedx2,
		const double	p2,
		const int	nGoodPixels ) const

Definition at line 222 of file PixeldEdxData.cxx.

                                                                                              {
  double dedx=dedx2; 
  double p=p2*0.001; 
  std::vector<float> vecResult(3,-1); 
  if (nGoodPixels>0){
    const auto & array { getP(dedx,p,nGoodPixels)};
    std::copy(array.begin(), array.end(),vecResult.begin());
  }
  return vecResult;
}

◆ getMass()

double PixeldEdxData::getMass	(	const double	dedx,
		const double	signedP,
		const int	nGoodPixels ) const

Definition at line 104 of file PixeldEdxData.cxx.

                                                                                                  {
 
  if (dedx<m_mindedxformass) { return m_piMass; }
 
  if (nGoodPixels<=0) { return -2; }
  auto  par = getFirstNPar(signedP,nGoodPixels,5);
 
  const auto p = std::abs(signedP);
 
  par[5] = dedx;
  par[6] = p;
 
  double xmin=0;
  double xmax=14000;
 
  double xmed = (xmin+xmax)*0.5;
  double prec=1e-6;
 
  if (dedx<15000){ return m_piMass; }
 
  do {
    if (fdEdxZero(xmed,par)==0) { return xmed; }
    if (fdEdxZero(xmed,par)*fdEdxZero(xmax,par)>0) { xmax=xmed; }
    else { xmin=xmed; }
    xmed=(xmin+xmax)*0.5;
  } while (std::abs(xmax-xmin)>prec);
  return xmed;
}

◆ getP()

std::array< double, 3 > PixeldEdxData::getP	(	const double	dedxArg,
		const double	signedP,
		const int	nGoodPixels ) const

Definition at line 71 of file PixeldEdxData.cxx.

                                                                                           {
  std::array<double,3> vhypo{-2,-2,-2};
  if (nGoodPixels<=0) { return vhypo; }
  const auto & par = getFirstNPar(signedP,nGoodPixels,9);
  const auto p = std::fabs(signedP);
 
  const auto dedx = std::log(dedxArg)-9.4;
 
  double pk  = dEdxPdf(dedx,p,m_kMass ,par,5);
  double pp  = dEdxPdf(dedx,p,m_pMass ,par,5);
  double ppi = dEdxPdf(dedx,p,m_piMass,par,5);
 
  if (ppi+pk+pp>=0) {
    vhypo = {ppi, pk,pp};
  } 
  else { 
    vhypo ={-1,-1,-1};
  }
  return vhypo;
}

◆ getPar()

double PixeldEdxData::getPar	(	const int	i,
		const int	j ) const

Definition at line 59 of file PixeldEdxData.cxx.

                                                           {
  auto itr = m_par.find(i);
  if (itr!=m_par.end()) {
    const std::vector<double> &param = itr->second;
    if (j<(int)param.size()) {
      return param[j];
    }
  }
  return 0;
}

◆ moyal()

double PixeldEdxData::moyal	(	const double	x,
		const double	Ep,
		const double	R )

static

Definition at line 174 of file PixeldEdxData.cxx.

                                                                         {
  double lambda = (x-Ep)/R;
  return std::sqrt(std::exp(-lambda-std::exp(-lambda))*M_1_PI*0.5);
}

◆ setBetheBlochType()

void PixeldEdxData::setBetheBlochType ( const std::string & bb )

Definition at line 22 of file PixeldEdxData.cxx.

                                                         {
  m_bb = bb;
}

◆ setFunctionType()

void PixeldEdxData::setFunctionType ( const std::string & fun )

Definition at line 18 of file PixeldEdxData.cxx.

                                                        {
  m_fun = fun;
}

◆ setMinimumdEdxForMass()

void PixeldEdxData::setMinimumdEdxForMass ( const double mindedxMass )

Definition at line 26 of file PixeldEdxData.cxx.

                                                                  {
  m_mindedxformass = mindedxMass;
}

◆ setPar()

void PixeldEdxData::setPar	(	const int	i,
		const double	param )

Definition at line 10 of file PixeldEdxData.cxx.

                                                          {
  m_par[i].push_back(param);
}

◆ setPosNeg()

void PixeldEdxData::setPosNeg ( const bool posneg )

Definition at line 14 of file PixeldEdxData.cxx.

                                               {
  m_posneg = posneg;
}

Member Data Documentation

◆ m_bb

std::string PixeldEdxData::m_bb

private

Definition at line 66 of file PixeldEdxData.h.

◆ m_fun

std::string PixeldEdxData::m_fun

private

Definition at line 65 of file PixeldEdxData.h.

◆ m_kMass

double PixeldEdxData::m_kMass = 0.49368

staticconstexprprivate

Definition at line 60 of file PixeldEdxData.h.

◆ m_mindedxformass

double PixeldEdxData::m_mindedxformass {}

private

Definition at line 67 of file PixeldEdxData.h.

67{};

◆ m_par

std::unordered_map<uint32_t,std::vector<double> > PixeldEdxData::m_par

private

Definition at line 63 of file PixeldEdxData.h.

◆ m_piMass

double PixeldEdxData::m_piMass = 0.13957

staticconstexprprivate

Definition at line 59 of file PixeldEdxData.h.

◆ m_pMass

double PixeldEdxData::m_pMass = 0.93827

staticconstexprprivate

Definition at line 61 of file PixeldEdxData.h.

◆ m_posneg

bool PixeldEdxData::m_posneg {}

private

Definition at line 64 of file PixeldEdxData.h.

64{};

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

Public Member Functions

Static Public Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Function Documentation

◆ asymGaus()

◆ crystalBall()

◆ dEdx_3p()

◆ dEdx_5p()

◆ dEdx_5p_aleph()

◆ dEdx_5p_BG()

◆ dEdx_5p_BG_aleph()

◆ dEdx_BG()

◆ dEdx_def()

◆ dEdxPdf()

◆ fdEdxZero()

◆ getdEdx() [1/2]

◆ getdEdx() [2/2]

◆ getFirstNPar()

◆ getLikelihoods()

◆ getMass()

◆ getP()

◆ getPar()

◆ moyal()

◆ setBetheBlochType()

◆ setFunctionType()

◆ setMinimumdEdxForMass()

◆ setPar()

◆ setPosNeg()

Member Data Documentation

◆ m_bb

◆ m_fun

◆ m_kMass

◆ m_mindedxformass

◆ m_par

◆ m_piMass

◆ m_pMass

◆ m_posneg