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<>
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

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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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 5 of file NjettinessAdder.cc.

    : 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::unique_ptr<fastjet::contrib::Njettiness>(new fastjet::contrib::Njettiness(*axesDef, *measureDef));
}

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

~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 144 of file NjettinessAdder.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 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);
}

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

Referenced by produce().

produce()

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

override

Definition at line 110 of file NjettinessAdder.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_);
 
  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());
  }
}

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

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().