NjettinessAdder Class Reference

#include <NjettinessAdder.h>

Inheritance diagram for NjettinessAdder:

Public Types
enum	AxesDefinition_t {   KT_Axes = 0, CA_Axes, AntiKT_Axes, WTA_KT_Axes,   WTA_CA_Axes, Manual_Axes, OnePass_KT_Axes, OnePass_CA_Axes,   OnePass_AntiKT_Axes, OnePass_WTA_KT_Axes, OnePass_WTA_CA_Axes, OnePass_Manual_Axes,   MultiPass_Axes, N_AXES_DEFINITIONS }
enum	MeasureDefinition_t {   NormalizedMeasure = 0, UnnormalizedMeasure, OriginalGeometricMeasure, NormalizedCutoffMeasure,   UnnormalizedCutoffMeasure, GeometricCutoffMeasure, N_MEASURE_DEFINITIONS }
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

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

Private Attributes
double	akAxesR0_
unsigned	axesDefinition_
double	beta_
edm::EDGetTokenT< edm::ValueMap< float > >	input_weights_token_
unsigned	measureDefinition_
std::vector< unsigned >	Njets_
int	nPass_
double	R0_
double	Rcutoff_
std::unique_ptr< fastjet::contrib::Njettiness >	routine_
edm::InputTag	src_
edm::EDGetTokenT< edm::View< reco::Jet > >	src_token_
edm::ValueMap< float > const *	weightsHandle_

Detailed Description

Definition at line 14 of file NjettinessAdder.h.

Member Enumeration Documentation

AxesDefinition_t

enum NjettinessAdder::AxesDefinition_t

Enumerator
KT_Axes
CA_Axes
AntiKT_Axes
WTA_KT_Axes
WTA_CA_Axes
Manual_Axes
OnePass_KT_Axes
OnePass_CA_Axes
OnePass_AntiKT_Axes
OnePass_WTA_KT_Axes
OnePass_WTA_CA_Axes
OnePass_Manual_Axes
MultiPass_Axes
N_AXES_DEFINITIONS

Definition at line 25 of file NjettinessAdder.h.

                        {
    KT_Axes = 0,
    CA_Axes,
    AntiKT_Axes,  // (axAxesR0)
    WTA_KT_Axes,
    WTA_CA_Axes,
    Manual_Axes,
    OnePass_KT_Axes,
    OnePass_CA_Axes,
    OnePass_AntiKT_Axes,  // (axAxesR0)
    OnePass_WTA_KT_Axes,
    OnePass_WTA_CA_Axes,
    OnePass_Manual_Axes,
    MultiPass_Axes,
    N_AXES_DEFINITIONS
  };

MeasureDefinition_t

enum NjettinessAdder::MeasureDefinition_t

Enumerator
NormalizedMeasure
UnnormalizedMeasure
OriginalGeometricMeasure
NormalizedCutoffMeasure
UnnormalizedCutoffMeasure
GeometricCutoffMeasure
N_MEASURE_DEFINITIONS

Definition at line 16 of file NjettinessAdder.h.

                           {
    NormalizedMeasure = 0,      // (beta,R0)
    UnnormalizedMeasure,        // (beta)
    OriginalGeometricMeasure,   // (beta)
    NormalizedCutoffMeasure,    // (beta,R0,Rcutoff)
    UnnormalizedCutoffMeasure,  // (beta,Rcutoff)
    GeometricCutoffMeasure,     // (beta,Rcutoff)
    N_MEASURE_DEFINITIONS
  };

Constructor & Destructor Documentation

NjettinessAdder()

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

explicit

Definition at line 7 of file NjettinessAdder.cc.

References akAxesR0_, AntiKT_Axes, axesDefinition_, beta_, CA_Axes, edm::ParameterSet::getParameter(), input_weights_token_, KT_Axes, measureDefinition_, MultiPass_Axes, create_idmaps::n, Njets_, NormalizedCutoffMeasure, NormalizedMeasure, nPass_, OnePass_AntiKT_Axes, OnePass_CA_Axes, OnePass_KT_Axes, OnePass_WTA_CA_Axes, OnePass_WTA_KT_Axes, OriginalGeometricMeasure, R0_, Rcutoff_, routine_, runBasic_cfi::srcWeights, UnnormalizedCutoffMeasure, UnnormalizedMeasure, WTA_CA_Axes, and WTA_KT_Axes.

    : src_(iConfig.getParameter<edm::InputTag>("src")),
      src_token_(consumes<edm::View<reco::Jet>>(src_)),
      Njets_(iConfig.getParameter<std::vector<unsigned>>("Njets")),
      measureDefinition_(iConfig.getParameter<unsigned>("measureDefinition")),
      beta_(iConfig.getParameter<double>("beta")),
      R0_(iConfig.getParameter<double>("R0")),
      Rcutoff_(iConfig.getParameter<double>("Rcutoff")),
      axesDefinition_(iConfig.getParameter<unsigned>("axesDefinition")),
      nPass_(iConfig.getParameter<int>("nPass")),
      akAxesR0_(iConfig.getParameter<double>("akAxesR0")) {
  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 tauN_str;
    tauN_str << "tau" << *n;

    produces<edm::ValueMap<float>>(tauN_str.str());
  }

  // Get the measure definition
  fastjet::contrib::NormalizedMeasure normalizedMeasure(beta_, R0_);
  fastjet::contrib::UnnormalizedMeasure unnormalizedMeasure(beta_);
  fastjet::contrib::OriginalGeometricMeasure geometricMeasure(beta_);  // changed in 1.020
  fastjet::contrib::NormalizedCutoffMeasure normalizedCutoffMeasure(beta_, R0_, Rcutoff_);
  fastjet::contrib::UnnormalizedCutoffMeasure unnormalizedCutoffMeasure(beta_, Rcutoff_);
  //fastjet::contrib::GeometricCutoffMeasure     geometricCutoffMeasure   (beta_,Rcutoff_); // removed in 1.020

  fastjet::contrib::MeasureDefinition const* measureDef = nullptr;
  switch (measureDefinition_) {
    case UnnormalizedMeasure:
      measureDef = &unnormalizedMeasure;
      break;
    case OriginalGeometricMeasure:
      measureDef = &geometricMeasure;
      break;  // changed in 1.020
    case NormalizedCutoffMeasure:
      measureDef = &normalizedCutoffMeasure;
      break;
    case UnnormalizedCutoffMeasure:
      measureDef = &unnormalizedCutoffMeasure;
      break;
    //case GeometricCutoffMeasure : measureDef = &geometricCutoffMeasure; break; // removed in 1.020
    case NormalizedMeasure:
    default:
      measureDef = &normalizedMeasure;
      break;
  }

  // Get the axes definition
  fastjet::contrib::KT_Axes kt_axes;
  fastjet::contrib::CA_Axes ca_axes;
  fastjet::contrib::AntiKT_Axes antikt_axes(akAxesR0_);
  fastjet::contrib::WTA_KT_Axes wta_kt_axes;
  fastjet::contrib::WTA_CA_Axes wta_ca_axes;
  fastjet::contrib::OnePass_KT_Axes onepass_kt_axes;
  fastjet::contrib::OnePass_CA_Axes onepass_ca_axes;
  fastjet::contrib::OnePass_AntiKT_Axes onepass_antikt_axes(akAxesR0_);
  fastjet::contrib::OnePass_WTA_KT_Axes onepass_wta_kt_axes;
  fastjet::contrib::OnePass_WTA_CA_Axes onepass_wta_ca_axes;
  fastjet::contrib::MultiPass_Axes multipass_axes(nPass_);

  fastjet::contrib::AxesDefinition const* axesDef = nullptr;
  switch (axesDefinition_) {
    case KT_Axes:
    default:
      axesDef = &kt_axes;
      break;
    case CA_Axes:
      axesDef = &ca_axes;
      break;
    case AntiKT_Axes:
      axesDef = &antikt_axes;
      break;
    case WTA_KT_Axes:
      axesDef = &wta_kt_axes;
      break;
    case WTA_CA_Axes:
      axesDef = &wta_ca_axes;
      break;
    case OnePass_KT_Axes:
      axesDef = &onepass_kt_axes;
      break;
    case OnePass_CA_Axes:
      axesDef = &onepass_ca_axes;
      break;
    case OnePass_AntiKT_Axes:
      axesDef = &onepass_antikt_axes;
      break;
    case OnePass_WTA_KT_Axes:
      axesDef = &onepass_wta_kt_axes;
      break;
    case OnePass_WTA_CA_Axes:
      axesDef = &onepass_wta_ca_axes;
      break;
    case MultiPass_Axes:
      axesDef = &multipass_axes;
      break;
  };

  routine_ = std::make_unique<fastjet::contrib::Njettiness>(*axesDef, *measureDef);
}

~NjettinessAdder()

NjettinessAdder::~NjettinessAdder ( )

inlineoverride

Definition at line 44 of file NjettinessAdder.h.

{}

Member Function Documentation

getTau()

float NjettinessAdder::getTau	(	unsigned	num,
		const edm::Ptr< reco::Jet > &	object
	)		const

Definition at line 146 of file NjettinessAdder.cc.

References Options::const, Calorimetry_cff::dp, input_weights_token_, edm::EDGetTokenT< T >::isUninitialized(), isotrackApplyRegressor::k, MainPageGenerator::l, EgammaValidation_cff::num, default_text_conditions_cfi::object, routine_, and w().

Referenced by produce().

                                                                               {
  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 N-jettiness computation is missing!";
          }
        }
      }  // end if basic jet
    }  // end if daughter pointer is nonnull and available
    else
      edm::LogWarning("MissingJetConstituent") << "Jet constituent required for N-jettiness computation is missing!";
  }

  return routine_->getTau(num, fjParticles);
}

produce()

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

override

Definition at line 112 of file NjettinessAdder.cc.

References getTau(), iEvent, input_weights_token_, edm::EDGetTokenT< T >::isUninitialized(), PDWG_EXODelayedJetMET_cff::jets, eostools::move(), create_idmaps::n, Njets_, src_token_, submitPVValidationJobs::t, and weightsHandle_.

                                                                           {
  // 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_);

  for (std::vector<unsigned>::const_iterator n = Njets_.begin(); n != Njets_.end(); ++n) {
    std::ostringstream tauN_str;
    tauN_str << "tau" << *n;

    // prepare room for output
    std::vector<float> tauN;
    tauN.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 = getTau(*n, jetPtr);

      tauN.push_back(t);
    }

    auto outT = std::make_unique<edm::ValueMap<float>>();
    edm::ValueMap<float>::Filler fillerT(*outT);
    fillerT.insert(jets, tauN.begin(), tauN.end());
    fillerT.fill();

    iEvent.put(std::move(outT), tauN_str.str());
  }
}

Member Data Documentation

akAxesR0_

double NjettinessAdder::akAxesR0_

private

Definition at line 63 of file NjettinessAdder.h.

Referenced by NjettinessAdder().

axesDefinition_

unsigned NjettinessAdder::axesDefinition_

private

Definition at line 61 of file NjettinessAdder.h.

Referenced by NjettinessAdder().

beta_

double NjettinessAdder::beta_

private

Definition at line 56 of file NjettinessAdder.h.

Referenced by NjettinessAdder().

input_weights_token_

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

private

Definition at line 65 of file NjettinessAdder.h.

Referenced by getTau(), NjettinessAdder(), and produce().

measureDefinition_

unsigned NjettinessAdder::measureDefinition_

private

Definition at line 55 of file NjettinessAdder.h.

Referenced by NjettinessAdder().

Njets_

std::vector<unsigned> NjettinessAdder::Njets_

private

Definition at line 52 of file NjettinessAdder.h.

Referenced by NjettinessAdder(), and produce().

nPass_

int NjettinessAdder::nPass_

private

Definition at line 62 of file NjettinessAdder.h.

Referenced by NjettinessAdder().

R0_

double NjettinessAdder::R0_

private

Definition at line 57 of file NjettinessAdder.h.

Referenced by NjettinessAdder().

Rcutoff_

double NjettinessAdder::Rcutoff_

private

Definition at line 58 of file NjettinessAdder.h.

Referenced by NjettinessAdder().

routine_

std::unique_ptr<fastjet::contrib::Njettiness> NjettinessAdder::routine_

private

Definition at line 68 of file NjettinessAdder.h.

Referenced by getTau(), and NjettinessAdder().

src_

edm::InputTag NjettinessAdder::src_

private

Definition at line 50 of file NjettinessAdder.h.

src_token_

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

private

Definition at line 51 of file NjettinessAdder.h.

Referenced by produce().

weightsHandle_

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

private

Definition at line 66 of file NjettinessAdder.h.

Referenced by produce().