CP::SSVWeightsAlg::nFMethodPileupBJetBasedClass Class Reference

Collaboration diagram for CP::SSVWeightsAlg::nFMethodPileupBJetBasedClass:

Public Member Functions
	nFMethodPileupBJetBasedClass (const nlohmann::json &jsonConfig)
double	getPFake (const double muactual, const int b_jet_count, const int N_fake, const double SF_fake_low, const double SF_fake_high) const

Private Attributes
std::map< std::string, std::map< std::string, double > >	m_nFPileupBJetMap {}
int	m_upperboundNbjets = -999
double	m_lowMuHighMuThreshold = -999

Detailed Description

Definition at line 95 of file SSVWeightsAlg.h.

Constructor & Destructor Documentation

nFMethodPileupBJetBasedClass()

CP::SSVWeightsAlg::nFMethodPileupBJetBasedClass::nFMethodPileupBJetBasedClass

(

const nlohmann::json &

jsonConfig

)

Definition at line 640 of file SSVWeightsAlg.cxx.

    : m_nFPileupBJetMap(jsonConfig["nF_pileup_bjet_based"])
  {
    // Extract information from JSON file for nFMethodPileupBJetBased
    std::map<std::string, double>::iterator lastItem = std::prev(m_nFPileupBJetMap.at("high_muactual").end());
    std::string lastItemKey = lastItem->first;
    m_upperboundNbjets = std::stoi(lastItemKey);
    m_lowMuHighMuThreshold = jsonConfig["CalibrationInformation"]["lowMuHighMuThreshold"];
  }

Member Function Documentation

getPFake()

double CP::SSVWeightsAlg::nFMethodPileupBJetBasedClass::getPFake	(	const double	muactual,
		const int	b_jet_count,
		const int	N_fake,
		const double	SF_fake_low,
		const double	SF_fake_high ) const

Definition at line 651 of file SSVWeightsAlg.cxx.

                                    {
 
    double P_fake = 1;  
    // 2D map muactual and Nbjets
    std::string mu_key = (muactual >= m_lowMuHighMuThreshold) ? "high_muactual" : "low_muactual";
 
    // Get the value
    double n_F_value = 0;
 
    if (b_jet_count < m_upperboundNbjets){
      // Build the bjets key string
      std::string bjets_key = std::to_string(b_jet_count) + "_bjets";
      n_F_value = m_nFPileupBJetMap.at(mu_key).at(bjets_key);
    }
    else{
      // Build the bjets key string
      std::string bjets_key = std::to_string(m_upperboundNbjets) + "p_bjets";
      n_F_value = m_nFPileupBJetMap.at(mu_key).at(bjets_key);
    }
    auto denom = poisson_pmf(N_fake, n_F_value);
    if (denom == 0. )[[unlikely]]{
      throw std::runtime_error("nFMethodPileupBJetBasedClass::getPFake: divide-by-zero");
    }
    if (muactual >= m_lowMuHighMuThreshold){
      P_fake = (poisson_pmf(N_fake, SF_fake_high*n_F_value))/denom;
    }
    else {
      P_fake = (poisson_pmf(N_fake, SF_fake_low*n_F_value))/denom;
    }
 
    return P_fake;
  }

Member Data Documentation

m_lowMuHighMuThreshold

double CP::SSVWeightsAlg::nFMethodPileupBJetBasedClass::m_lowMuHighMuThreshold = -999

private

Definition at line 107 of file SSVWeightsAlg.h.

m_nFPileupBJetMap

std::map<std::string, std::map<std::string, double> > CP::SSVWeightsAlg::nFMethodPileupBJetBasedClass::m_nFPileupBJetMap {}

private

Definition at line 105 of file SSVWeightsAlg.h.

{};

m_upperboundNbjets

int CP::SSVWeightsAlg::nFMethodPileupBJetBasedClass::m_upperboundNbjets = -999

private

Definition at line 106 of file SSVWeightsAlg.h.

---

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

• SSVWeightsAlg.h
• SSVWeightsAlg.cxx