ECFAdder Class Reference

#include <ECFAdder.h>

Inheritance diagram for ECFAdder:

Public Member Functions
	ECFAdder (const edm::ParameterSet &iConfig)
float	getECF (unsigned index, const edm::Ptr< reco::Jet > &object) const
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
double	alpha_
double	beta_
std::vector< std::string >	cuts_
std::string	ecftype_
edm::EDGetTokenT< edm::ValueMap< float > >	input_weights_token_
std::vector< unsigned >	Njets_
std::vector< std::shared_ptr< fastjet::FunctionOfPseudoJet< double > > >	routine_
std::vector< StringCutObjectSelector< reco::Jet > >	selectors_
edm::InputTag	src_
edm::EDGetTokenT< edm::View< reco::Jet > >	src_token_
std::vector< std::string >	variables_
edm::ValueMap< float > const *	weightsHandle_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 15 of file ECFAdder.h.

Constructor & Destructor Documentation

ECFAdder()

ECFAdder::ECFAdder ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 9 of file ECFAdder.cc.

    : src_(iConfig.getParameter<edm::InputTag>("src")),
      src_token_(consumes<edm::View<reco::Jet>>(src_)),
      Njets_(iConfig.getParameter<std::vector<unsigned>>("Njets")),
      cuts_(iConfig.getParameter<std::vector<std::string>>("cuts")),
      ecftype_(iConfig.getParameter<std::string>("ecftype")),
      alpha_(iConfig.getParameter<double>("alpha")),
      beta_(iConfig.getParameter<double>("beta")) {
  if (cuts_.size() != Njets_.size()) {
    throw cms::Exception("ConfigurationError") << "cuts and Njets must be the same size in ECFAdder" << std::endl;
  }
 
  edm::InputTag srcWeights = iConfig.getParameter<edm::InputTag>("srcWeights");
  if (!srcWeights.label().empty())
    input_weights_token_ = consumes<edm::ValueMap<float>>(srcWeights);
 
  for (std::vector<unsigned>::const_iterator n = Njets_.begin(); n != Njets_.end(); ++n) {
    std::ostringstream ecfN_str;
    std::shared_ptr<fastjet::FunctionOfPseudoJet<double>> pfunc;
 
    if (ecftype_ == "ECF" || ecftype_.empty()) {
      ecfN_str << "ecf" << *n;
      pfunc.reset(new fastjet::contrib::EnergyCorrelator(*n, beta_, fastjet::contrib::EnergyCorrelator::pt_R));
    } else if (ecftype_ == "C") {
      ecfN_str << "ecfC" << *n;
      pfunc.reset(new fastjet::contrib::EnergyCorrelatorCseries(*n, beta_, fastjet::contrib::EnergyCorrelator::pt_R));
    } else if (ecftype_ == "D") {
      ecfN_str << "ecfD" << *n;
      pfunc.reset(
          new fastjet::contrib::EnergyCorrelatorGeneralizedD2(alpha_, beta_, fastjet::contrib::EnergyCorrelator::pt_R));
    } else if (ecftype_ == "N") {
      ecfN_str << "ecfN" << *n;
      pfunc.reset(new fastjet::contrib::EnergyCorrelatorNseries(*n, beta_, fastjet::contrib::EnergyCorrelator::pt_R));
    } else if (ecftype_ == "M") {
      ecfN_str << "ecfM" << *n;
      pfunc.reset(new fastjet::contrib::EnergyCorrelatorMseries(*n, beta_, fastjet::contrib::EnergyCorrelator::pt_R));
    } else if (ecftype_ == "U") {
      ecfN_str << "ecfU" << *n;
      pfunc.reset(new fastjet::contrib::EnergyCorrelatorUseries(*n, beta_, fastjet::contrib::EnergyCorrelator::pt_R));
    }
    variables_.push_back(ecfN_str.str());
    produces<edm::ValueMap<float>>(ecfN_str.str());
    routine_.push_back(pfunc);
 
    selectors_.push_back(StringCutObjectSelector<reco::Jet>(cuts_[n - Njets_.begin()]));
  }
}

References alpha_, beta_, cuts_, ecftype_, Exception, edm::ParameterSet::getParameter(), input_weights_token_, dqmiodumpmetadata::n, Njets_, routine_, selectors_, HLT_FULL_cff::srcWeights, and variables_.

Member Function Documentation

fillDescriptions()

void ECFAdder::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 137 of file ECFAdder.cc.

                                                                          {
  edm::ParameterSetDescription iDesc;
  iDesc.setComment("Energy Correlation Functions adder");
 
  iDesc.add<edm::InputTag>("src", edm::InputTag("no default"))->setComment("input collection");
  iDesc.add<std::vector<unsigned>>("Njets", {1, 2, 3})->setComment("Number of jets to emulate");
  iDesc.add<std::vector<std::string>>("cuts", {"", "", ""})
      ->setComment("Jet selections for each N value. Size must match Njets.");
  iDesc.add<double>("alpha", 1.0)->setComment("alpha factor, only valid for N2");
  iDesc.add<double>("beta", 1.0)->setComment("angularity factor");
  iDesc.add<std::string>("ecftype", "")->setComment("ECF type: ECF or empty; C; D; N; M; U;");
  iDesc.add<edm::InputTag>("srcWeights", edm::InputTag("puppi"));
  descriptions.add("ECFAdder", iDesc);
}

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HLT_FULL_cff::InputTag, edm::ParameterSetDescription::setComment(), edm::ParameterDescriptionNode::setComment(), and AlCaHLTBitMon_QueryRunRegistry::string.

getECF()

float ECFAdder::getECF	(	unsigned	index,
		const edm::Ptr< reco::Jet > &	object
	)		const

Definition at line 92 of file ECFAdder.cc.

                                                                          {
  std::vector<fastjet::PseudoJet> fjParticles;
  for (unsigned k = 0; k < object->numberOfDaughters(); ++k) {
    const reco::CandidatePtr& dp = object->daughterPtr(k);
    if (dp.isNonnull() && dp.isAvailable()) {
      // Here, the daughters are the "end" node, so this is a PFJet
      if (dp->numberOfDaughters() == 0) {
        if (object->isWeighted()) {
          if (input_weights_token_.isUninitialized())
            throw cms::Exception("MissingConstituentWeight")
                << "ECFAdder: No weights (e.g. PUPPI) given for weighted jet collection" << std::endl;
          float w = (*weightsHandle_)[dp];
          fjParticles.push_back(fastjet::PseudoJet(dp->px() * w, dp->py() * w, dp->pz() * w, dp->energy() * w));
        } else
          fjParticles.push_back(fastjet::PseudoJet(dp->px(), dp->py(), dp->pz(), dp->energy()));
      } else {  // Otherwise, this is a BasicJet, so you need to descend further.
        auto subjet = dynamic_cast<reco::Jet const*>(dp.get());
        for (unsigned l = 0; l < subjet->numberOfDaughters(); ++l) {
          if (subjet != nullptr) {
            const reco::CandidatePtr& ddp = subjet->daughterPtr(l);
            if (subjet->isWeighted()) {
              if (input_weights_token_.isUninitialized())
                throw cms::Exception("MissingConstituentWeight")
                    << "ECFAdder: No weights (e.g. PUPPI) given for weighted jet collection" << std::endl;
              float w = (*weightsHandle_)[ddp];
              fjParticles.push_back(fastjet::PseudoJet(ddp->px() * w, ddp->py() * w, ddp->pz() * w, ddp->energy() * w));
            } else
              fjParticles.push_back(fastjet::PseudoJet(ddp->px(), ddp->py(), ddp->pz(), ddp->energy()));
          } else {
            edm::LogWarning("MissingJetConstituent") << "BasicJet constituent required for ECF computation is missing!";
          }
        }
      }  // end if basic jet
    }    // end if daughter pointer is nonnull and available
    else
      edm::LogWarning("MissingJetConstituent") << "Jet constituent required for ECF computation is missing!";
  }
  if (fjParticles.size() > Njets_[index]) {
    return routine_[index]->result(join(fjParticles));
  } else {
    return -1.0;
  }
}

References Calorimetry_cff::dp, input_weights_token_, edm::EDGetTokenT< T >::isUninitialized(), join(), dqmdumpme::k, cmsLHEtoEOSManager::l, Njets_, resolutioncreator_cfi::object, routine_, and w.

Referenced by produce().

produce()

void ECFAdder::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 57 of file ECFAdder.cc.

                                                                    {
  // read input collection
  edm::Handle<edm::View<reco::Jet>> jets;
  iEvent.getByToken(src_token_, jets);
 
  // Get Weights Collection
  if (!input_weights_token_.isUninitialized())
    weightsHandle_ = &iEvent.get(input_weights_token_);
 
  unsigned i = 0;
  for (std::vector<unsigned>::const_iterator n = Njets_.begin(); n != Njets_.end(); ++n) {
    // prepare room for output
    std::vector<float> ecfN;
    ecfN.reserve(jets->size());
 
    for (typename edm::View<reco::Jet>::const_iterator jetIt = jets->begin(); jetIt != jets->end(); ++jetIt) {
      edm::Ptr<reco::Jet> jetPtr = jets->ptrAt(jetIt - jets->begin());
 
      float t = -1.0;
      if (selectors_[n - Njets_.begin()](*jetIt))
        t = getECF(i, jetPtr);
 
      ecfN.push_back(t);
    }
 
    auto outT = std::make_unique<edm::ValueMap<float>>();
    edm::ValueMap<float>::Filler fillerT(*outT);
    fillerT.insert(jets, ecfN.begin(), ecfN.end());
    fillerT.fill();
 
    iEvent.put(std::move(outT), variables_[i]);
    ++i;
  }
}

References getECF(), mps_fire::i, iEvent, input_weights_token_, edm::EDGetTokenT< T >::isUninitialized(), singleTopDQM_cfi::jets, eostools::move(), dqmiodumpmetadata::n, Njets_, selectors_, src_token_, submitPVValidationJobs::t, variables_, and weightsHandle_.

Member Data Documentation

alpha_

double ECFAdder::alpha_

private

Definition at line 31 of file ECFAdder.h.

Referenced by ECFAdder().

beta_

double ECFAdder::beta_

private

Definition at line 32 of file ECFAdder.h.

Referenced by ECFAdder().

cuts_

std::vector<std::string> ECFAdder::cuts_

private

Definition at line 28 of file ECFAdder.h.

Referenced by ECFAdder().

ecftype_

std::string ECFAdder::ecftype_

private

Definition at line 29 of file ECFAdder.h.

Referenced by ECFAdder().

input_weights_token_

edm::EDGetTokenT<edm::ValueMap<float> > ECFAdder::input_weights_token_

private

Definition at line 36 of file ECFAdder.h.

Referenced by ECFAdder(), getECF(), and produce().

Njets_

std::vector<unsigned> ECFAdder::Njets_

private

Definition at line 27 of file ECFAdder.h.

Referenced by ECFAdder(), getECF(), and produce().

routine_

std::vector<std::shared_ptr<fastjet::FunctionOfPseudoJet<double> > > ECFAdder::routine_

private

Definition at line 34 of file ECFAdder.h.

Referenced by ECFAdder(), and getECF().

selectors_

std::vector<StringCutObjectSelector<reco::Jet> > ECFAdder::selectors_

private

Definition at line 35 of file ECFAdder.h.

Referenced by ECFAdder(), and produce().

src_

edm::InputTag ECFAdder::src_

private

Definition at line 25 of file ECFAdder.h.

src_token_

edm::EDGetTokenT<edm::View<reco::Jet> > ECFAdder::src_token_

private

Definition at line 26 of file ECFAdder.h.

Referenced by produce().

variables_

std::vector<std::string> ECFAdder::variables_

private

Definition at line 30 of file ECFAdder.h.

Referenced by ECFAdder(), and produce().

weightsHandle_

edm::ValueMap<float> const* ECFAdder::weightsHandle_

private

Definition at line 37 of file ECFAdder.h.

Referenced by produce().