#include <CATopJetProducer.h>

Inheritance diagram for cms::CATopJetProducer:

Public Member Functions
	CATopJetProducer (const edm::ParameterSet &ps)

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

virtual void	runAlgorithm (edm::Event &iEvent, const edm::EventSetup &iSetup)

virtual	~CATopJetProducer ()

Public Member Functions inherited from FastjetJetProducer
	FastjetJetProducer (const edm::ParameterSet &iConfig)

virtual	~FastjetJetProducer ()

Public Member Functions inherited from VirtualJetProducer
std::string	jetType () const

	VirtualJetProducer (const edm::ParameterSet &iConfig)

virtual	~VirtualJetProducer ()

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

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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 Types
enum	tagalgos { CA_TOPTAGGER, FJ_CMS_TOPTAG, FJ_HEP_TOPTAG, FJ_JHU_TOPTAG, FJ_NSUB_TAG }

Private Attributes
double	centralEtaCut_

std::auto_ptr < fastjet::CMSTopTagger >	fjCMSTopTagger_
	The algorithm to do the work. More...

std::auto_ptr < fastjet::HEPTopTagger >	fjHEPTopTagger_

std::auto_ptr < fastjet::JHTopTagger >	fjJHUTopTagger_

std::auto_ptr < fastjet::RestFrameNSubjettinessTagger >	fjNSUBTagger_

std::auto_ptr< CATopJetAlgorithm >	legacyCMSTopTagger_

double	ptMin_

int	tagAlgo_

bool	verbose_

Additional Inherited Members
Public Types inherited from FastjetJetProducer
typedef boost::shared_ptr < DynamicRfilt >	DynamicRfiltPtr

typedef fastjet::Transformer	transformer

typedef std::vector < transformer_ptr >	transformer_coll

typedef std::unique_ptr < transformer >	transformer_ptr

Public Types inherited from VirtualJetProducer
typedef boost::shared_ptr < fastjet::GhostedAreaSpec >	ActiveAreaSpecPtr

typedef boost::shared_ptr < fastjet::AreaDefinition >	AreaDefinitionPtr

typedef boost::shared_ptr < fastjet::ClusterSequence >	ClusterSequencePtr

typedef boost::shared_ptr < fastjet::JetDefinition >	JetDefPtr

typedef boost::shared_ptr < fastjet::JetDefinition::Plugin >	PluginPtr

typedef boost::shared_ptr < fastjet::RangeDefinition >	RangeDefPtr

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 Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from FastjetJetProducer
virtual void	produceTrackJets (edm::Event &iEvent, const edm::EventSetup &iSetup)

Protected Member Functions inherited from VirtualJetProducer
virtual void	copyConstituents (const std::vector< fastjet::PseudoJet > &fjConstituents, reco::Jet *jet)

virtual std::vector < reco::CandidatePtr >	getConstituents (const std::vector< fastjet::PseudoJet > &fjConstituents)

virtual void	inputTowers ()

virtual bool	isAnomalousTower (reco::CandidatePtr input)

bool	makeBasicJet (const JetType::Type &fTag)

bool	makeCaloJet (const JetType::Type &fTag)

bool	makeGenJet (const JetType::Type &fTag)

bool	makePFClusterJet (const JetType::Type &fTag)

bool	makePFJet (const JetType::Type &fTag)

virtual void	makeProduces (std::string s, std::string tag="")

bool	makeTrackJet (const JetType::Type &fTag)

void	offsetCorrectJets (std::vector< fastjet::PseudoJet > &orphanInput)

virtual void	output (edm::Event &iEvent, edm::EventSetup const &iSetup)

template<typename T >
void	writeCompoundJets (edm::Event &iEvent, edm::EventSetup const &iSetup)
	function template to write out the outputs More...

template<typename T >
void	writeJets (edm::Event &iEvent, edm::EventSetup const &iSetup)

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)

Protected Attributes inherited from VirtualJetProducer
bool	doAreaDiskApprox_

bool	doAreaFastjet_

bool	doFastJetNonUniform_

bool	doPUOffsetCorr_

bool	doPVCorrection_

bool	doRhoFastjet_

ActiveAreaSpecPtr	fjActiveArea_

AreaDefinitionPtr	fjAreaDefinition_

ClusterSequencePtr	fjClusterSeq_

std::vector< fastjet::PseudoJet >	fjInputs_

JetDefPtr	fjJetDefinition_

std::vector< fastjet::PseudoJet >	fjJets_

PluginPtr	fjPlugin_

RangeDefPtr	fjRangeDef_

edm::EDGetTokenT < reco::VertexCollection >	input_vertex_token_

double	inputEMin_

double	inputEtMin_

std::vector< edm::Ptr < reco::Candidate > >	inputs_

std::string	jetAlgorithm_

std::string	jetCollInstanceName_

double	jetPtMin_

std::string	jetType_

JetType::Type	jetTypeE

unsigned int	maxInputs_

unsigned int	minSeed_

std::string	moduleLabel_

unsigned int	nExclude_

std::vector< double >	puCenters_

std::string	puSubtractorName_

double	puWidth_

bool	restrictInputs_

double	rParam_

edm::InputTag	src_

edm::InputTag	srcPVs_

boost::shared_ptr < PileUpSubtractor >	subtractor_

bool	useDeterministicSeed_

bool	useExplicitGhosts_

int	verbosity_

reco::Particle::Point	vertex_

double	voronoiRfact_

bool	writeCompound_

Detailed Description

Definition at line 76 of file CATopJetProducer.h.

Member Enumeration Documentation

enum cms::CATopJetProducer::tagalgos

private

Enumerator
CA_TOPTAGGER
FJ_CMS_TOPTAG
FJ_HEP_TOPTAG
FJ_JHU_TOPTAG
FJ_NSUB_TAG

Definition at line 101 of file CATopJetProducer.h.

                   {
     CA_TOPTAGGER,
     FJ_CMS_TOPTAG,
     FJ_HEP_TOPTAG,
     FJ_JHU_TOPTAG,
     FJ_NSUB_TAG
     };

Constructor & Destructor Documentation

CATopJetProducer::CATopJetProducer ( const edm::ParameterSet & ps )

Definition at line 10 of file CATopJetProducer.cc.

References CA_TOPTAGGER, FJ_CMS_TOPTAG, FJ_HEP_TOPTAG, FJ_JHU_TOPTAG, FJ_NSUB_TAG, fjCMSTopTagger_, fjHEPTopTagger_, fjJHUTopTagger_, fjNSUBTagger_, edm::ParameterSet::getParameter(), legacyCMSTopTagger_, fwrapper::plugin, VirtualJetProducer::src_, and tagAlgo_.

                                                              :
        FastjetJetProducer( conf ),
        tagAlgo_(conf.getParameter<int>("tagAlgo")),
        ptMin_(conf.getParameter<double>("jetPtMin")),
        centralEtaCut_(conf.getParameter<double>("centralEtaCut")),
        verbose_(conf.getParameter<bool>("verbose"))
 {
 
     if (tagAlgo_ == CA_TOPTAGGER ) {
         
         legacyCMSTopTagger_ = std::auto_ptr<CATopJetAlgorithm>(
             new CATopJetAlgorithm(src_,
             conf.getParameter<bool>  ("verbose"),              
                 conf.getParameter<int>  ("algorithm"),                  // 0 = KT, 1 = CA, 2 = anti-KT
                 conf.getParameter<int>   ("useAdjacency"),                  // choose adjacency requirement:
                                                                 //  0 = no adjacency
                                                                 //  1 = deltar adjacency 
                                                                 //  2 = modified adjacency
                                                                 //  3 = calotower neirest neigbor based adjacency (untested)
                conf.getParameter<double>("centralEtaCut"),              // eta for defining "central" jets
                conf.getParameter<double>("jetPtMin"),                   // min jet pt
                conf.getParameter<std::vector<double> >("sumEtBins"),    // sumEt bins over which cuts may vary. vector={bin 0 lower bound, bin 1 lower bound, ...} 
                conf.getParameter<std::vector<double> >("rBins"),        // Jet distance paramter R. R values depend on sumEt bins.
                conf.getParameter<std::vector<double> >("ptFracBins"),   // fraction of hard jet pt that subjet must have (deltap)
                conf.getParameter<std::vector<double> >("deltarBins"),   // Applicable only if useAdjacency=1. deltar adjacency values for each sumEtBin
                conf.getParameter<std::vector<double> >("nCellBins"),    // Applicable only if useAdjacency=3. number of cells to consider two subjets adjacent
                conf.getParameter<double>("inputEtMin"),                 // seed threshold - NOT USED
                conf.getParameter<bool>  ("useMaxTower"),                // use max tower as adjacency criterion, otherwise use centroid - NOT USED
                conf.getParameter<double>("sumEtEtaCut"),                // eta for event SumEt - NOT USED
             conf.getParameter<double>("etFrac")                     // fraction of event sumEt / 2 for a jet to be considered "hard" -NOT USED
             ));
     }
     else if (tagAlgo_ == FJ_CMS_TOPTAG ) {
         fjCMSTopTagger_ = std::auto_ptr<fastjet::CMSTopTagger>(
             new fastjet::CMSTopTagger(conf.getParameter<double> ("ptFrac"),
                           conf.getParameter<double> ("rFrac"),
                           conf.getParameter<double> ("adjacencyParam"))
         );
     }
     else if (tagAlgo_ == FJ_HEP_TOPTAG ) {
         fjHEPTopTagger_ = std::auto_ptr<fastjet::HEPTopTagger>(
             new fastjet::HEPTopTagger(conf.getParameter<double>("muCut"),
                           conf.getParameter<double>("maxSubjetMass"),
                           conf.getParameter<bool>("useSubjetMass")
                           )
         );
     }
     else if (tagAlgo_ == FJ_JHU_TOPTAG ) {
         fjJHUTopTagger_ = std::auto_ptr<fastjet::JHTopTagger>(
             new fastjet::JHTopTagger(conf.getParameter<double>("ptFrac"),
                          conf.getParameter<double>("deltaRCut"),
                          conf.getParameter<double>("cosThetaWMax")
                          )
         );
     }
     else if (tagAlgo_ == FJ_NSUB_TAG ) {
         
         fastjet::JetDefinition::Plugin *plugin = new fastjet::SISConePlugin(0.6, 0.75);
         fastjet::JetDefinition NsubJetDef(plugin);
         fjNSUBTagger_ = std::auto_ptr<fastjet::RestFrameNSubjettinessTagger>(
             new fastjet::RestFrameNSubjettinessTagger(NsubJetDef,
                                   conf.getParameter<double>("tau2Cut"),
                                   conf.getParameter<double>("cosThetaSCut"),
                                   conf.getParameter<bool>("useExclusive")
                                   )
         );
     }
                 
         
 
 
 
 }

virtual cms::CATopJetProducer::~CATopJetProducer ( )

inlinevirtual

Definition at line 82 of file CATopJetProducer.h.

82 {}

Member Function Documentation

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

virtual

Reimplemented from FastjetJetProducer.

Definition at line 84 of file CATopJetProducer.cc.

References FastjetJetProducer::produce().

 {
   FastjetJetProducer::produce(e, c);
 }

void CATopJetProducer::runAlgorithm	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Reimplemented from FastjetJetProducer.

Definition at line 89 of file CATopJetProducer.cc.

References CA_TOPTAGGER, centralEtaCut_, gather_cfg::cout, VirtualJetProducer::doAreaFastjet_, VirtualJetProducer::doRhoFastjet_, FJ_CMS_TOPTAG, FJ_HEP_TOPTAG, FJ_JHU_TOPTAG, FJ_NSUB_TAG, VirtualJetProducer::fjAreaDefinition_, VirtualJetProducer::fjClusterSeq_, fjCMSTopTagger_, fjHEPTopTagger_, VirtualJetProducer::fjInputs_, VirtualJetProducer::fjJetDefinition_, VirtualJetProducer::fjJets_, fjJHUTopTagger_, fjNSUBTagger_, i, perp(), ptMin_, tagAlgo_, verbose_, and VirtualJetProducer::voronoiRfact_.

 {
   if ( !doAreaFastjet_ && !doRhoFastjet_) {
     fjClusterSeq_ = ClusterSequencePtr( new fastjet::ClusterSequence( fjInputs_, *fjJetDefinition_ ) );
   } else if (voronoiRfact_ <= 0) {
     fjClusterSeq_ = ClusterSequencePtr( new fastjet::ClusterSequenceArea( fjInputs_, *fjJetDefinition_ , *fjAreaDefinition_ ) );
   } else {
     fjClusterSeq_ = ClusterSequencePtr( new fastjet::ClusterSequenceVoronoiArea( fjInputs_, *fjJetDefinition_ , fastjet::VoronoiAreaSpec(voronoiRfact_) ) );
   }
 
   if (tagAlgo_ == CA_TOPTAGGER){
     (*legacyCMSTopTagger_).run( fjInputs_, fjJets_, fjClusterSeq_ );
     
   }
   else {
     
     //Run the jet clustering
     vector<fastjet::PseudoJet> inclusiveJets = fjClusterSeq_->inclusive_jets(ptMin_);
 
     if ( verbose_ ) cout << "Getting central jets" << endl;
     // Find the transient central jets
     vector<fastjet::PseudoJet> centralJets;
     for (unsigned int i = 0; i < inclusiveJets.size(); i++) {
         
         if (inclusiveJets[i].perp() > ptMin_ && fabs(inclusiveJets[i].rapidity()) < centralEtaCut_) {
             centralJets.push_back(inclusiveJets[i]);
         }
     }
 
     fastjet::CMSTopTagger & CMSTagger = *fjCMSTopTagger_;
     fastjet::HEPTopTagger & HEPTagger = *fjHEPTopTagger_;
     fastjet::JHTopTagger & JHUTagger = *fjJHUTopTagger_;
     fastjet::RestFrameNSubjettinessTagger & NSUBTagger = *fjNSUBTagger_;
 
 
     vector<fastjet::PseudoJet>::iterator jetIt = centralJets.begin(), centralJetsEnd = centralJets.end();
     if ( verbose_ )cout<<"Loop over jets"<<endl;
     for ( ; jetIt != centralJetsEnd; ++jetIt ) {
         
         if (verbose_) cout << "CMS FJ jet pt: " << (*jetIt).perp() << endl;
 
         fastjet::PseudoJet taggedJet;
         if (tagAlgo_ == FJ_CMS_TOPTAG) taggedJet = CMSTagger.result(*jetIt);
         else if (tagAlgo_ == FJ_HEP_TOPTAG) taggedJet = HEPTagger.result(*jetIt);
         else if (tagAlgo_ == FJ_JHU_TOPTAG) taggedJet = JHUTagger.result(*jetIt);
         else if (tagAlgo_ == FJ_NSUB_TAG) taggedJet = NSUBTagger.result(*jetIt);
         else cout << "NOT A VALID TAGGING ALGORITHM CHOICE!" << endl;
 
         if (taggedJet != 0) fjJets_.push_back(taggedJet);
     }
   }
 } 

Member Data Documentation

double cms::CATopJetProducer::centralEtaCut_

private

Definition at line 99 of file CATopJetProducer.h.

Referenced by runAlgorithm().

std::auto_ptr<fastjet::CMSTopTagger> cms::CATopJetProducer::fjCMSTopTagger_

private

The algorithm to do the work.

Definition at line 90 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

std::auto_ptr<fastjet::HEPTopTagger> cms::CATopJetProducer::fjHEPTopTagger_

private

Definition at line 91 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

std::auto_ptr<fastjet::JHTopTagger> cms::CATopJetProducer::fjJHUTopTagger_

private

Definition at line 92 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

std::auto_ptr<fastjet::RestFrameNSubjettinessTagger> cms::CATopJetProducer::fjNSUBTagger_

private

Definition at line 93 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

std::auto_ptr<CATopJetAlgorithm> cms::CATopJetProducer::legacyCMSTopTagger_

private

Definition at line 89 of file CATopJetProducer.h.

Referenced by CATopJetProducer().

double cms::CATopJetProducer::ptMin_

private

Definition at line 98 of file CATopJetProducer.h.

Referenced by runAlgorithm().

int cms::CATopJetProducer::tagAlgo_

private

Definition at line 97 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

bool cms::CATopJetProducer::verbose_

private

Definition at line 100 of file CATopJetProducer.h.

Referenced by runAlgorithm().

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation