#include <ParametricResolution.h>

Public Member Functions
float	operator() (const float pt, const float abseta) const

	ParametricResolution ()

	ParametricResolution (const edm::ParameterSet &cpset)

Static Public Member Functions
static std::vector< float >	getVFloat (const edm::ParameterSet &cpset, const std::string &name)

Protected Types
enum	Kind { Kind::Calo, Kind::Track }

Protected Attributes
std::vector< float >	etas_

Kind	kind_

std::vector< float >	offsets_

std::vector< float >	ptMaxs_

std::vector< float >	ptMins_

std::vector< float >	scales_

Detailed Description

Definition at line 10 of file ParametricResolution.h.

Member Enumeration Documentation

enum l1tpf::ParametricResolution::Kind

strongprotected

Enumerator
Calo
Track

Definition at line 21 of file ParametricResolution.h.

21 { Calo, Track };

dumpRecoGeometry_cfg.Calo

tuple Calo

Definition: dumpRecoGeometry_cfg.py:191

Constructor & Destructor Documentation

l1tpf::ParametricResolution::ParametricResolution ( )

inline

Definition at line 14 of file ParametricResolution.h.

14 {}

l1tpf::ParametricResolution::ParametricResolution ( const edm::ParameterSet & cpset )

Definition at line 8 of file ParametricResolution.cc.

References Calo, etas_, Exception, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), getVFloat(), kind_, ptMaxs_, ptMin, ptMins_, AlCaHLTBitMon_QueryRunRegistry::string, and Track.

     : etas_(getVFloat(cpset, "etaBins")), offsets_(getVFloat(cpset, "offset")), scales_(getVFloat(cpset, "scale")) {
   if (cpset.existsAs<std::vector<double>>("ptMin")) {
     ptMins_ = getVFloat(cpset, "ptMin");
   } else {
     float ptMin = cpset.existsAs<double>("ptMin") ? cpset.getParameter<double>("ptMin") : 0;
     ptMins_ = std::vector<float>(etas_.size(), ptMin);
   }
   if (cpset.existsAs<std::vector<double>>("ptMax")) {
     ptMaxs_ = getVFloat(cpset, "ptMax");
   } else {
     ptMaxs_ = std::vector<float>(etas_.size(), 1e6);
   }
 
   std::string skind = cpset.getParameter<std::string>("kind");
   if (skind == "track")
     kind_ = Kind::Track;
   else if (skind == "calo")
     kind_ = Kind::Calo;
   else
     throw cms::Exception("Configuration", "Bad kind of resolution: " + skind);
 }

Member Function Documentation

std::vector< float > l1tpf::ParametricResolution::getVFloat	(	const edm::ParameterSet &	cpset,
		const std::string &	name
	)

static

Definition at line 3 of file ParametricResolution.cc.

References edm::ParameterSet::getParameter(), and mergeVDriftHistosByStation::name.

Referenced by ParametricResolution().

                                                                                                          {
   std::vector<double> vd = cpset.getParameter<std::vector<double>>(name);
   return std::vector<float>(vd.begin(), vd.end());
 }

float l1tpf::ParametricResolution::operator()	(	const float	pt,
		const float	abseta
	)		const

Definition at line 31 of file ParametricResolution.cc.

References dumpRecoGeometry_cfg::Calo, mps_fire::i, dqmiodumpmetadata::n, and DiDispStaMuonMonitor_cfi::pt.

                                                                                     {
   for (unsigned int i = 0, n = etas_.size(); i < n; ++i) {
     if (pt > ptMaxs_[i])
       continue;
     if (abseta < etas_[i]) {
       switch (kind_) {
         case Kind::Track:
           return pt * std::min<float>(1.f, std::hypot(pt * scales_[i] * 0.001, offsets_[i]));
         case Kind::Calo:
           return std::min<float>(pt, pt * scales_[i] + offsets_[i]);
           if (pt < ptMins_[i])
             return pt * std::min<float>(1, scales_[i] + offsets_[i] / ptMins_[i]);
           return std::min<float>(pt, pt * scales_[i] + offsets_[i]);
       }
     }
   }
   return std::min<float>(pt, 0.3 * pt + 7);  // saturate to 100% at 10 GeV, and to 30% at high pt
 }