#include <NjettinessAdder.h>

Inheritance diagram for NjettinessAdder:

Public Member Functions
float	getTau (int num, edm::Ptr< reco::PFJet > 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
double	cone_

edm::InputTag	src_

edm::EDGetTokenT< edm::View < reco::PFJet > >	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 12 of file NjettinessAdder.h.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 14 of file NjettinessAdder.h.

                                                            :
     src_(iConfig.getParameter<edm::InputTag>("src")),
     src_token_(consumes<edm::View<reco::PFJet>>(src_)),
     cone_(iConfig.getParameter<double>("cone"))
       {
     produces<edm::ValueMap<float> >("tau1");
     produces<edm::ValueMap<float> >("tau2");
     produces<edm::ValueMap<float> >("tau3");
       }

virtual NjettinessAdder::~NjettinessAdder ( )

inlinevirtual

Definition at line 24 of file NjettinessAdder.h.

24 {}

Member Function Documentation

float NjettinessAdder::getTau	(	int	num,
		edm::Ptr< reco::PFJet >	object
	)		const

Definition at line 47 of file NjettinessAdder.cc.

References cone_, reco::LeafCandidate::energy(), edm::Ptr< T >::get(), gen::k, nJettinessAdder_cfi::Njettiness, reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::LeafCandidate::pz().

Referenced by produce().

 {
   std::vector<const reco::PFCandidate*> all_particles;
   for (unsigned k =0; k < object->getPFConstituents().size(); k++)
     all_particles.push_back( object->getPFConstituent(k).get() );
 
   std::vector<fastjet::PseudoJet> FJparticles;
   for (unsigned particle = 0; particle < all_particles.size(); particle++){
     const reco::PFCandidate *thisParticle = all_particles.at(particle);
     FJparticles.push_back( fastjet::PseudoJet( thisParticle->px(), thisParticle->py(), thisParticle->pz(), thisParticle->energy() ) );  
   }
   fastjet::contrib::NsubParameters paraNsub = fastjet::contrib::NsubParameters(1.0, cone_); //assume R=0.7 jet clusering used
   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, edm::Event::put(), src_token_, and reco::t2.

                                                                              {
   // read input collection
   edm::Handle<edm::View<reco::PFJet> > jets;
   iEvent.getByToken(src_token_, jets);
   
   // prepare room for output
   std::vector<float> tau1;         tau1.reserve(jets->size());
   std::vector<float> tau2;         tau2.reserve(jets->size());
   std::vector<float> tau3;         tau3.reserve(jets->size());
 
   for ( typename edm::View<reco::PFJet>::const_iterator jetIt = jets->begin() ; jetIt != jets->end() ; ++jetIt ) {
     reco::PFJet newCand(*jetIt);
     edm::Ptr<reco::PFJet> jetPtr = jets->ptrAt(jetIt - jets->begin());
 
     float t1=getTau(1, jetPtr );
     float t2=getTau(2, jetPtr );
     float t3=getTau(3, jetPtr );
 
     tau1.push_back(t1);
     tau2.push_back(t2);
     tau3.push_back(t3);
   }
 
   std::auto_ptr<edm::ValueMap<float> > outT1(new edm::ValueMap<float>());
   std::auto_ptr<edm::ValueMap<float> > outT2(new edm::ValueMap<float>());
   std::auto_ptr<edm::ValueMap<float> > outT3(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerT1(*outT1);
   edm::ValueMap<float>::Filler fillerT2(*outT2);
   edm::ValueMap<float>::Filler fillerT3(*outT3);
   fillerT1.insert(jets, tau1.begin(), tau1.end());
   fillerT2.insert(jets, tau2.begin(), tau2.end());
   fillerT3.insert(jets, tau3.begin(), tau3.end());
   fillerT1.fill();
   fillerT2.fill();
   fillerT3.fill();
   
   iEvent.put(outT1,"tau1");
   iEvent.put(outT2,"tau2");
   iEvent.put(outT3,"tau3");
 }

Member Data Documentation

double NjettinessAdder::cone_

private

Definition at line 32 of file NjettinessAdder.h.

Referenced by getTau().

edm::InputTag NjettinessAdder::src_

private

Definition at line 30 of file NjettinessAdder.h.

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

private

Definition at line 31 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