#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, GeometricMeasure, NormalizedCutoffMeasure, UnnormalizedCutoffMeasure, GeometricCutoffMeasure, N_MEASURE_DEFINITIONS }

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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)

virtual	~NjettinessAdder ()

Public Member Functions inherited from edm::one::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
double	akAxesR0_

unsigned	axesDefinition_

double	beta_

unsigned	measureDefinition_

std::vector< unsigned >	Njets_

int	nPass_

double	R0_

double	Rcutoff_

std::auto_ptr < fastjet::contrib::Njettiness >	routine_

edm::InputTag	src_

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

Additional Inherited Members
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 15 of file NjettinessAdder.h.

Member Enumeration Documentation

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 27 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
     };

enum NjettinessAdder::MeasureDefinition_t

Enumerator
NormalizedMeasure
UnnormalizedMeasure
GeometricMeasure
NormalizedCutoffMeasure
UnnormalizedCutoffMeasure
GeometricCutoffMeasure
N_MEASURE_DEFINITIONS

Definition at line 18 of file NjettinessAdder.h.

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

Constructor & Destructor Documentation

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

explicit

Definition at line 6 of file NjettinessAdder.cc.

References akAxesR0_, AntiKT_Axes, axesDefinition_, beta_, CA_Axes, GeometricCutoffMeasure, GeometricMeasure, KT_Axes, measureDefinition_, MultiPass_Axes, gen::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, R0_, Rcutoff_, routine_, 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"))
 {
   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().c_str());
     }
 
 
   // Get the measure definition
   fastjet::contrib::NormalizedMeasure          normalizedMeasure        (beta_,R0_);
   fastjet::contrib::UnnormalizedMeasure        unnormalizedMeasure      (beta_);
   fastjet::contrib::GeometricMeasure           geometricMeasure         (beta_);
   fastjet::contrib::NormalizedCutoffMeasure    normalizedCutoffMeasure  (beta_,R0_,Rcutoff_);
   fastjet::contrib::UnnormalizedCutoffMeasure  unnormalizedCutoffMeasure(beta_,Rcutoff_);
   fastjet::contrib::GeometricCutoffMeasure     geometricCutoffMeasure   (beta_,Rcutoff_);
 
   fastjet::contrib::MeasureDefinition const * measureDef = 0;
   switch ( measureDefinition_ ) {
   case UnnormalizedMeasure : measureDef = &unnormalizedMeasure; break;
   case GeometricMeasure    : measureDef = &geometricMeasure; break;
   case NormalizedCutoffMeasure : measureDef = &normalizedCutoffMeasure; break;
   case UnnormalizedCutoffMeasure : measureDef = &unnormalizedCutoffMeasure; break;
   case GeometricCutoffMeasure : measureDef = &geometricCutoffMeasure; break;
   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 = 0;
   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::auto_ptr<fastjet::contrib::Njettiness> ( new fastjet::contrib::Njettiness( *axesDef, *measureDef ) );
       
 }

virtual NjettinessAdder::~NjettinessAdder ( )

inlinevirtual

Definition at line 48 of file NjettinessAdder.h.

48 {}

Member Function Documentation

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

Definition at line 109 of file NjettinessAdder.cc.

References edm::Ptr< T >::isAvailable(), edm::Ptr< T >::isNonnull(), relval_steps::k, and routine_.

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() )
     FJparticles.push_back( fastjet::PseudoJet( dp->px(), dp->py(), dp->pz(), dp->energy() ) );
       else
     edm::LogWarning("MissingJetConstituent") << "Jet constituent required for N-subjettiness computation is missing!";
     }
 
   return routine_->getTau(num, FJparticles); 
 }

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

virtual

Implements edm::one::EDProducerBase.

Definition at line 77 of file NjettinessAdder.cc.

References edm::helper::Filler< Map >::fill(), edm::Event::getByToken(), getTau(), edm::helper::Filler< Map >::insert(), fwrapper::jets, gen::n, Njets_, edm::Event::put(), src_token_, and tree::t.

                                                                              {
   // read input collection
   edm::Handle<edm::View<reco::Jet> > jets;
   iEvent.getByToken(src_token_, jets);
   
   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);
       }
 
       std::auto_ptr<edm::ValueMap<float> > outT(new edm::ValueMap<float>());
       edm::ValueMap<float>::Filler fillerT(*outT);
       fillerT.insert(jets, tauN.begin(), tauN.end());
       fillerT.fill();
 
       iEvent.put(outT,tauN_str.str().c_str());
     }
 }