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	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () 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 7 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::make_unique<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_, 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_FULL_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 146 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 112 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_, submitPVValidationJobs::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().