ZDCTriggerEfficiency Class Reference

#include <ZDCTriggerEfficiency.h>

Collaboration diagram for ZDCTriggerEfficiency:

Public Member Functions
void	UpdatelumiBlock (unsigned int lumiBlock)
	ZDCTriggerEfficiency ()
void	SetEffParamsAndErrors (std::array< std::vector< TSpline3 * >, 2 > effParams, std::array< std::vector< TSpline3 * >, 2 > effParamErrors)
void	SetEffParamCorrCoeffs (std::array< std::vector< TSpline3 * >, 2 > effParamsCorrCoeffs)
double	GetEfficiency (int side, float ADCSum)
std::pair< double, double >	GetEfficiencyAndError (MsgStream &msg, int side, float ADCSum)

Private Attributes
bool	m_haveParams
bool	m_haveCorrCoeffs
unsigned int	m_minLB
unsigned int	m_maxLB
unsigned int	m_currentLB
std::array< std::vector< TSpline3 * >, 2 >	m_effParams
std::array< std::vector< TSpline3 * >, 2 >	m_effParamErrors
std::array< std::vector< TSpline3 * >, 2 >	m_effParamCorrCoeffs
std::array< std::vector< double >, 2 >	m_currentParams
std::array< std::vector< double >, 2 >	m_currentParamErrors
std::array< std::vector< double >, 2 >	m_currentCorrCoefff

Detailed Description

Definition at line 16 of file ZDCTriggerEfficiency.h.

Constructor & Destructor Documentation

ZDCTriggerEfficiency()

ZDCTriggerEfficiency::ZDCTriggerEfficiency ( )

inline

Definition at line 59 of file ZDCTriggerEfficiency.h.

                         : // in order, alpha-beta, alpha-theta, beta-theta
    m_haveParams(false), m_haveCorrCoeffs(false),
    m_minLB(0), m_maxLB(0), m_currentLB(0) {
 
  }

Member Function Documentation

GetEfficiency()

double ZDCTriggerEfficiency::GetEfficiency	(	int	side,
		float	ADCSum
	)

Definition at line 9 of file ZDCTriggerEfficiency.cxx.

                                                                {
  double alpha = m_currentParams[side][0];
  double  beta = m_currentParams[side][1];
  double  theta = m_currentParams[side][2];
  if(alpha<10e-5||beta<10e-5||theta<10e-5)
    {
      return -1;
    }
  double m = std::exp(-ADCSum/theta);
  double p = std::exp(-(ADCSum - alpha) / beta);
  double efficiency = (1 - m) / (1 + p);
  return efficiency;
}

GetEfficiencyAndError()

std::pair< double, double > ZDCTriggerEfficiency::GetEfficiencyAndError	(	MsgStream &	msg,
		int	side,
		float	ADCSum
	)

Definition at line 23 of file ZDCTriggerEfficiency.cxx.

                                                                                                         {
  double alpha = m_currentParams[side][0];
  double  beta = m_currentParams[side][1];
  double  theta = m_currentParams[side][2];
  if(alpha<10e-5||beta<10e-5||theta<10e-5)
    {
      return std::make_pair(-1,-1);
    }
  
  double alphaErr = m_currentParamErrors[side][0];
  double betaErr = m_currentParamErrors[side][1];
  double thetaErr = m_currentParamErrors[side][2];
  
  double corrAlphaBeta = m_currentCorrCoefff[side][0];
  double corrAlphaTheta = m_currentCorrCoefff[side][1];
  double corrBetaTheta = m_currentCorrCoefff[side][2];
  
  double m = std::exp(-ADCSum/theta);
  double p = std::exp(-(ADCSum - alpha) / beta);
  double efficiency = (1 - m) / (1 + p);
  double dda = (1 - m) * -p / beta / (1 + p) / (1 + p);
  double ddb = (1 - m) * -p / std::pow(1 + p, 2.0) * ADCSum / beta / beta;
  double ddt = -ADCSum * m / theta / theta / (1 + p);
  
  double efficiencyErr = std::sqrt(
                   std::pow(alphaErr * dda, 2) +
                   std::pow(betaErr * ddb, 2) +
                   std::pow(thetaErr * ddt, 2) +
                   2 * corrAlphaBeta*alphaErr*betaErr * dda * ddb +
                   2 * corrAlphaTheta*alphaErr*thetaErr * dda * ddt +
                   2 * corrBetaTheta*betaErr*thetaErr * ddb * ddt);
  
  if (efficiencyErr<0)
    {
      efficiencyErr = 0;
      msg << MSG::INFO << "Efficiency Uncertainty calculation abnormal." << endmsg;
    }
  
  return std::make_pair(efficiency, efficiencyErr);
}

SetEffParamCorrCoeffs()

void ZDCTriggerEfficiency::SetEffParamCorrCoeffs ( std::array< std::vector< TSpline3 * >, 2 >  effParamsCorrCoeffs )

inline

Definition at line 77 of file ZDCTriggerEfficiency.h.

                                                                                    {
    for (int side : {0, 1}) {
      m_effParamCorrCoeffs[side] = effParamsCorrCoeffs[side];
    }
 
    m_haveCorrCoeffs = true;
  }

SetEffParamsAndErrors()

void ZDCTriggerEfficiency::SetEffParamsAndErrors ( std::array< std::vector< TSpline3 * >, 2 >  effParams, std::array< std::vector< TSpline3 * >, 2 >  effParamErrors  )

inline

Definition at line 65 of file ZDCTriggerEfficiency.h.

                                                                               {
    for (int side : {0, 1}) {
      m_effParams[side] = effParams[side];
      m_effParamErrors[side] = effParamErrors[side];
    }
 
    m_minLB = m_effParams[0][0]->GetXmin();
    m_maxLB = m_effParams[0][0]->GetXmax();
    m_haveParams = true;
  }

UpdatelumiBlock()

void ZDCTriggerEfficiency::UpdatelumiBlock ( unsigned int  lumiBlock )

inline

Definition at line 35 of file ZDCTriggerEfficiency.h.

  {
    if (!m_haveParams) throw std::runtime_error("No params for trigger efficiency.");
 
    if (lumiBlock != m_currentLB) {
      int newLB = std::max(std::min(lumiBlock, m_maxLB), m_minLB);
      for (int side : {0, 1}) {
        m_currentParams[side].clear();
        m_currentParamErrors[side].clear();
        m_currentCorrCoefff[side].clear();
 
        for (size_t ipar = 0; ipar < m_effParams[side].size(); ipar++) {
          m_currentParams[side].push_back(m_effParams[side][ipar]->Eval(newLB));
          m_currentParamErrors[side].push_back(m_effParamErrors[side][ipar]->Eval(newLB));
 
          if (m_haveCorrCoeffs) {
            m_currentCorrCoefff[side].push_back(m_effParamCorrCoeffs[side][ipar]->Eval(newLB));
          }
        }
        m_currentLB = lumiBlock;
      }
    }
  }

Member Data Documentation

m_currentCorrCoefff

std::array<std::vector<double>, 2> ZDCTriggerEfficiency::m_currentCorrCoefff

private

Definition at line 32 of file ZDCTriggerEfficiency.h.

m_currentLB

unsigned int ZDCTriggerEfficiency::m_currentLB

private

Definition at line 24 of file ZDCTriggerEfficiency.h.

m_currentParamErrors

std::array<std::vector<double>, 2> ZDCTriggerEfficiency::m_currentParamErrors

private

Definition at line 31 of file ZDCTriggerEfficiency.h.

m_currentParams

std::array<std::vector<double>, 2> ZDCTriggerEfficiency::m_currentParams

private

Definition at line 30 of file ZDCTriggerEfficiency.h.

m_effParamCorrCoeffs

std::array<std::vector<TSpline3*>, 2> ZDCTriggerEfficiency::m_effParamCorrCoeffs

private

Definition at line 28 of file ZDCTriggerEfficiency.h.

m_effParamErrors

std::array<std::vector<TSpline3*>, 2> ZDCTriggerEfficiency::m_effParamErrors

private

Definition at line 27 of file ZDCTriggerEfficiency.h.

m_effParams

std::array<std::vector<TSpline3*>, 2> ZDCTriggerEfficiency::m_effParams

private

Definition at line 26 of file ZDCTriggerEfficiency.h.

m_haveCorrCoeffs

bool ZDCTriggerEfficiency::m_haveCorrCoeffs

private

Definition at line 19 of file ZDCTriggerEfficiency.h.

m_haveParams

bool ZDCTriggerEfficiency::m_haveParams

private

Definition at line 18 of file ZDCTriggerEfficiency.h.

m_maxLB

unsigned int ZDCTriggerEfficiency::m_maxLB

private

Definition at line 22 of file ZDCTriggerEfficiency.h.

m_minLB

unsigned int ZDCTriggerEfficiency::m_minLB

private

Definition at line 21 of file ZDCTriggerEfficiency.h.

---

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

• ZDCTriggerEfficiency.h
• ZDCTriggerEfficiency.cxx