JME::JetResolutionVsEta Class Reference

Inheritance diagram for JME::JetResolutionVsEta:

Public Member Functions
bool	fill () override
	JetResolutionVsEta ()
Public Member Functions inherited from cond::payloadInspector::Histogram1< AxisType, PayloadType, IOV_M >
bool	fill () override
void	fillWithBinAndValue (size_t bin, AxisType weight=1)
void	fillWithValue (AxisType value, AxisType weight=1)
virtual AxisType	getFromPayload (PayloadType &payload)
	Histogram1 (const std::string &title, const std::string &xLabel, size_t nbins, float min, float max, const std::string &yLabel="entries")
void	init () override
Public Member Functions inherited from cond::payloadInspector::Plot2D< PayloadType, AxisType, AxisType, IOV_M, 1 >
std::shared_ptr< PayloadType >	fetchPayload (const cond::Hash &payloadHash)
	Plot2D (const std::string &type, const std::string &title, const std::string xLabel, const std::string &yLabel)
std::string	serializeData () override
	~Plot2D () override=default
Public Member Functions inherited from cond::payloadInspector::PlotImpl< IOV_M, NTAGS >
	PlotImpl (const std::string &type, const std::string &title)
std::string	processData () override
	~PlotImpl () override=default
Public Member Functions inherited from cond::payloadInspector::PlotBase
void	addInputParam (const std::string &paramName)
cond::persistency::Session	dbSession ()
template<typename PayloadType >
std::shared_ptr< PayloadType >	fetchPayload (const cond::Hash &payloadHash)
template<int index>
TagReference	getTag ()
cond::Tag_t	getTagInfo (const std::string &tag)
const std::map< std::string, std::string > &	inputParamValues () const
bool	isSingleIov () const
bool	isTwoTags () const
unsigned int	ntags () const
std::string	payloadType () const
	PlotBase ()
std::string	title () const
std::string	type () const
virtual	~PlotBase ()=default

Additional Inherited Members
Public Types inherited from cond::payloadInspector::Histogram1< AxisType, PayloadType, IOV_M >
typedef Plot2D< PayloadType, AxisType, AxisType, IOV_M, 1 >	Base
Public Types inherited from cond::payloadInspector::Plot2D< PayloadType, AxisType, AxisType, IOV_M, 1 >
typedef PlotImpl< IOV_M, NTAGS >	Base
Protected Attributes inherited from cond::payloadInspector::Plot2D< PayloadType, AxisType, AxisType, IOV_M, 1 >
std::vector< std::tuple< AxisType, AxisType > >	m_plotData
Protected Attributes inherited from cond::payloadInspector::PlotBase
std::set< std::string >	m_inputParams
std::map< std::string, std::string >	m_inputParamValues
PlotAnnotations	m_plotAnnotations
std::vector< std::pair< cond::Time_t, cond::Time_t > >	m_tagBoundaries
std::vector< std::vector< std::tuple< cond::Time_t, cond::Hash > > >	m_tagIovs
std::vector< std::string >	m_tagNames

Detailed Description

Definition at line 48 of file JetResolution_PayloadInspector.cc.

Constructor & Destructor Documentation

JetResolutionVsEta()

JME::JetResolutionVsEta::JetResolutionVsEta ( )

inline

Definition at line 50 of file JetResolution_PayloadInspector.cc.

References cond::payloadInspector::PlotBase::addInputParam().

        : cond::payloadInspector::Histogram1D<JetResolutionObject, SINGLE_IOV>(
              "Jet Resolution", "#eta", NBIN_ETA, MIN_ETA, MAX_ETA, "Resolution") {
      cond::payloadInspector::PlotBase::addInputParam("Jet_Pt");
      cond::payloadInspector::PlotBase::addInputParam("Jet_Rho");
    }

Member Function Documentation

fill()

bool JME::JetResolutionVsEta::fill ( )

inlineoverridevirtual

Implements cond::payloadInspector::PlotImpl< IOV_M, NTAGS >.

Definition at line 57 of file JetResolution_PayloadInspector.cc.

References f, mps_fire::i, heavyIonCSV_trainingSettings::idx, cond::payloadInspector::PlotBase::inputParamValues(), MAX_ETA, MIN_ETA, NBIN_ETA, jetsAK4_Puppi_cff::payload, AlCaHarvesting_cff::record, makeGlobalPositionRcd_cfg::tag, and trigObjTnPSource_cfi::var.

                         {
      double par_Pt = 100.;
      double par_Eta = 1.;
      double par_Rho = 20.;

      // Default values will be used if no input parameters
      auto paramValues = cond::payloadInspector::PlotBase::inputParamValues();
      auto ip = paramValues.find("Jet_Pt");
      if (ip != paramValues.end()) {
        par_Pt = std::stod(ip->second);
        edm::LogPrint("JER_PI") << "Jet Pt: " << par_Pt;
      }
      ip = paramValues.find("Jet_Rho");
      if (ip != paramValues.end()) {
        par_Rho = std::stod(ip->second);
        edm::LogPrint("JER_PI") << "Jet Rho: " << par_Rho;
      }

      auto tag = PlotBase::getTag<0>();
      for (auto const& iov : tag.iovs) {
        std::shared_ptr<JetResolutionObject> payload = Base::fetchPayload(std::get<1>(iov));
        if (payload.get()) {
          if (!payload->getRecords().empty() &&  // No formula for SF
              payload->getDefinition().getFormulaString().compare("") == 0)
            return false;

          for (const auto& record : payload->getRecords()) {
            // Check Pt & Rho
            if (!record.getVariablesRange().empty() && payload->getDefinition().getVariableName(0) == "JetPt" &&
                record.getVariablesRange()[0].is_inside(par_Pt)) {
              if (record.getBinsRange().size() > 1 && payload->getDefinition().getBinName(1) == "Rho" &&
                  record.getBinsRange()[1].is_inside(par_Rho)) {
                if (!record.getBinsRange().empty() && payload->getDefinition().getBinName(0) == "JetEta") {
                  reco::FormulaEvaluator f(payload->getDefinition().getFormulaString());

                  for (size_t idx = 0; idx <= NBIN_ETA; idx++) {
                    par_Eta = (idx + 0.5) * (MAX_ETA - MIN_ETA) / NBIN_ETA + MIN_ETA;
                    if (record.getBinsRange()[0].is_inside(par_Eta)) {
                      std::vector<double> var = {par_Pt};
                      std::vector<double> param;
                      for (size_t i = 0; i < record.getParametersValues().size(); i++) {
                        double par = record.getParametersValues()[i];
                        param.push_back(par);
                      }
                      float res = f.evaluate(var, param);
                      fillWithBinAndValue(idx + 1, res);
                    }
                  }
                }
              }
            }
          }  // records
          return true;
        }
      }
      return false;
    }