#include <L1EGammaEECalibrator.h>

Public Member Functions
float	calibrationFactor (const float &pt, const float &eta) const

	L1EGammaEECalibrator (const edm::ParameterSet &)

Private Member Functions
int	bin (const std::set< float > &container, float value) const

int	etaBin (float eta) const

int	ptBin (float pt) const

Private Attributes
std::vector< float >	calib_factors

std::set< float >	eta_bins

std::set< float >	pt_bins

Detailed Description

Definition at line 9 of file L1EGammaEECalibrator.h.

Constructor & Destructor Documentation

◆ L1EGammaEECalibrator()

L1EGammaEECalibrator::L1EGammaEECalibrator ( const edm::ParameterSet & pset )

explicit

Definition at line 18 of file L1EGammaEECalibrator.cc.

References l1tp2::as_vector(), calib_factors, filterCSVwithJSON::copy, eta_bins, etaBin(), edm::FileInPath::fullPath(), muonDTDigis_cfi::pset, pt_bins, and ptBin().

                                                                       {
   //read the JSON file and populate the eta - pt bins and the value container
   boost::property_tree::ptree calibration_map;
   read_json(pset.getParameter<edm::FileInPath>("calibrationFile").fullPath(), calibration_map);
 
   auto eta_l = l1tp2::as_vector(calibration_map, "eta_l");
   std::copy(eta_l.begin(), eta_l.end(), std::inserter(eta_bins, eta_bins.end()));
   auto eta_h = l1tp2::as_vector(calibration_map, "eta_h");
   eta_bins.insert(eta_h.back());
 
   auto pt_l = l1tp2::as_vector(calibration_map, "pt_l");
   std::copy(pt_l.begin(), pt_l.end(), std::inserter(pt_bins, pt_bins.end()));
   auto pt_h = l1tp2::as_vector(calibration_map, "pt_h");
   pt_bins.insert(pt_h.back());
 
   auto calib_data = l1tp2::as_vector(calibration_map, "calib");
   auto n_bins_eta = eta_bins.size();
   auto n_bins_pt = pt_bins.size();
   calib_factors.reserve(n_bins_eta * n_bins_pt);
   for (auto calib_f = calib_data.begin(); calib_f != calib_data.end(); ++calib_f) {
     auto index = calib_f - calib_data.begin();
     int eta_bin = etaBin(eta_l[index]);
     int pt_bin = ptBin(pt_l[index]);
     calib_factors[(eta_bin * n_bins_pt) + pt_bin] = *calib_f;
   }
 }

Member Function Documentation

◆ bin()

int L1EGammaEECalibrator::bin	(	const std::set< float > &	container,
		float	value
	)		const

private

Definition at line 45 of file L1EGammaEECalibrator.cc.

References HLT_2022v12_cff::distance.

Referenced by etaBin(), and ptBin().

                                                                                {
   auto bin_l = container.upper_bound(value);
   if (bin_l == container.end()) {
     // value not mapped to any bin
     return -1;
   }
   return std::distance(container.begin(), bin_l) - 1;
 }

◆ calibrationFactor()

float L1EGammaEECalibrator::calibrationFactor	(	const float &	pt,
		const float &	eta
	)		const

Definition at line 54 of file L1EGammaEECalibrator.cc.

References calib_factors, PVValHelper::eta, etaBin(), DiDispStaMuonMonitor_cfi::pt, pt_bins, and ptBin().

Referenced by L1EGammaEEProducer::produce().

                                                                                      {
   int bin_eta = etaBin(eta);
   int bin_pt = ptBin(pt);
   if (bin_eta == -1 || bin_pt == -1)
     return 1.;
   auto n_bins_pt = pt_bins.size();
   return calib_factors[(bin_eta * n_bins_pt) + bin_pt];
 }