DipzLikelihood Class Reference

#include <DipzLikelihoodCmp.h>

Collaboration diagram for DipzLikelihood:

Public Member Functions
	DipzLikelihood (const std::string &decName_z, const std::string &decName_negLogSigma2)
	~DipzLikelihood ()=default
double	operator() (const HypoJetVector &) const

Private Member Functions
double	getDipzMLPLDecValue (const pHypoJet &ip, const std::string &decName) const
double	calcNum (double acmlt, const pHypoJet &ip) const
double	calcDenom (double acmlt, const pHypoJet &ip) const
double	calcLogTerm (double acmlt, const pHypoJet &ip, double zhat) const
double	checkedRatio (double num, double den) const

Private Attributes
const std::string	m_decName_z
const std::string	m_decName_negLogSigma2

Detailed Description

Definition at line 40 of file DipzLikelihoodCmp.h.

Constructor & Destructor Documentation

DipzLikelihood()

DipzLikelihood::DipzLikelihood	(	const std::string &	decName_z,
		const std::string &	decName_negLogSigma2
	)

Definition at line 11 of file DipzLikelihoodCmp.cxx.

                                                         :
  m_decName_z(decName_z),
  m_decName_negLogSigma2(decName_negLogSigma2){
}

~DipzLikelihood()

DipzLikelihood::~DipzLikelihood ( )

default

Member Function Documentation

calcDenom()

double DipzLikelihood::calcDenom ( double  acmlt, const pHypoJet &  ip  ) const

private

Definition at line 54 of file DipzLikelihoodCmp.cxx.

                                                                      { 
  double sigma_squared =
    std::exp(-1 * getDipzMLPLDecValue(ip, m_decName_negLogSigma2));
 
  double oneoversigmasq = checkedRatio(1, sigma_squared);
  
  return acml + oneoversigmasq;
}

calcLogTerm()

double DipzLikelihood::calcLogTerm ( double  acmlt, const pHypoJet &  ip, double  zhat  ) const

private

Definition at line 63 of file DipzLikelihoodCmp.cxx.

                                      {
 
  double dipz_mu = getDipzMLPLDecValue(ip, m_decName_z);
 
  double dipz_negLogSigmaSq =
    getDipzMLPLDecValue(ip, m_decName_negLogSigma2);
 
  double sigma_squared = std::exp(-1 * dipz_negLogSigmaSq); 
  
  double logterm =
    -0.5 * std::log(2.0 * M_PI)
    + 0.5 * dipz_negLogSigmaSq
    - checkedRatio(std::pow(zhat - dipz_mu, 2), (2.0 * sigma_squared) );
  
  return acml + logterm;
 
}

calcNum()

double DipzLikelihood::calcNum ( double  acmlt, const pHypoJet &  ip  ) const

private

Definition at line 42 of file DipzLikelihoodCmp.cxx.

                                                                    {
  
  double sigma_squared =
    std::exp(-1 * getDipzMLPLDecValue(ip, m_decName_negLogSigma2));
  
  double muoversigmasq =
    checkedRatio( getDipzMLPLDecValue(ip, m_decName_z), sigma_squared);
  
  return acml + muoversigmasq;
}

checkedRatio()

double DipzLikelihood::checkedRatio ( double  num, double  den  ) const

private

Definition at line 17 of file DipzLikelihoodCmp.cxx.

                                                                {
  if (den == 0.) {
    // dividing x/0  is picked up by FPEAditor. C++ simply returns
    // +inf or -inf if x != 0, or a nan otherwise.
    throw std::runtime_error("DipzLikelihood::checkedRatio dividing by 0");
  }
  
  return num/den;
}

getDipzMLPLDecValue()

double DipzLikelihood::getDipzMLPLDecValue ( const pHypoJet &  ip, const std::string &  decName  ) const

private

Definition at line 27 of file DipzLikelihoodCmp.cxx.

{
  
  float momentValue;
  if (!(ip->getAttribute(decName, momentValue))) {
    throw std::runtime_error("Impossible to retrieve decorator \'" +
                 decName + "\' for jet hypo");
  }
  
  // momentValue is retrieved as a float, but will be heavily used in
  // further calculations. Convert to a double here
  return momentValue;
}

operator()()

double DipzLikelihood::operator()

(

const HypoJetVector &

ips

)

const

Definition at line 84 of file DipzLikelihoodCmp.cxx.

                                                                {
  
  auto zhat_num = std::accumulate(ips.begin(),
                  ips.end(),
                  0.0,
                  [this](double sum, const pHypoJet& jp) {
                return this->calcNum(sum, jp);});
  
  auto zhat_den = std::accumulate(ips.begin(),
                  ips.end(),
                  0.0,
                  [this](double sum, const pHypoJet& jp) {
                    return this->calcDenom(sum, jp);});
  
  auto zhat = checkedRatio(zhat_num, zhat_den);
  
  auto logproduct =
    std::accumulate(ips.begin(), 
            ips.end(), 
            0.0, 
            [&zhat,this](double sum, const pHypoJet& jp) {
              return this->calcLogTerm(sum, jp, zhat);});
  
  return logproduct;
}

Member Data Documentation

m_decName_negLogSigma2

const std::string DipzLikelihood::m_decName_negLogSigma2

private

Definition at line 51 of file DipzLikelihoodCmp.h.

m_decName_z

const std::string DipzLikelihood::m_decName_z

private

Definition at line 50 of file DipzLikelihoodCmp.h.

---

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

• DipzLikelihoodCmp.h
• DipzLikelihoodCmp.cxx