#include <NjettinessAdder.h>

Inheritance diagram for NjettinessAdder:

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::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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
	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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
const double	cone_

const std::vector< unsigned >	Njets_

const edm::InputTag	src_

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

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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 13 of file NjettinessAdder.h.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 15 of file NjettinessAdder.h.

References n, and Njets_.

                                                              :
       src_(iConfig.getParameter<edm::InputTag>("src")),
       src_token_(consumes<edm::View<reco::Jet>>(src_)),
       cone_(iConfig.getParameter<double>("cone")),
       Njets_(iConfig.getParameter<std::vector<unsigned> >("Njets"))
     {
       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());
       }
     }

virtual NjettinessAdder::~NjettinessAdder ( )

inlinevirtual

Definition at line 30 of file NjettinessAdder.h.

30 {}

Member Function Documentation

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

Definition at line 38 of file NjettinessAdder.cc.

References cone_, edm::Ptr< T >::isAvailable(), edm::Ptr< T >::isNonnull(), gen::k, and nJettinessAdder_cfi::Njettiness.

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!";
   }
 
   fastjet::contrib::NsubParameters paraNsub = fastjet::contrib::NsubParameters(1.0, cone_);
   fastjet::contrib::Njettiness routine(fastjet::contrib::Njettiness::onepass_kt_axes, paraNsub);
   return routine.getTau(num, FJparticles); 
 }

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

virtual

Implements edm::EDProducer.

Definition at line 6 of file NjettinessAdder.cc.

References edm::helper::Filler< Map >::fill(), edm::Event::getByToken(), getTau(), edm::helper::Filler< Map >::insert(), fwrapper::jets, n, Njets_, edm::Event::put(), src_token_, and edmStreamStallGrapher::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());
   }
 }

Member Data Documentation

const double NjettinessAdder::cone_

private

Definition at line 38 of file NjettinessAdder.h.

Referenced by getTau().

const std::vector<unsigned> NjettinessAdder::Njets_

private

Definition at line 39 of file NjettinessAdder.h.

Referenced by NjettinessAdder(), and produce().

const edm::InputTag NjettinessAdder::src_

private

Definition at line 36 of file NjettinessAdder.h.

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

private

Definition at line 37 of file NjettinessAdder.h.

Referenced by produce().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation