#include <HTTTopJetProducer.h>

Inheritance diagram for cms::HTTTopJetProducer:

Public Member Functions
virtual void	addHTTTopJetTagInfoCollection (edm::Event &iEvent, const edm::EventSetup &iSetup, edm::OrphanHandle< reco::BasicJetCollection > &oh)

	HTTTopJetProducer (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	~HTTTopJetProducer ()

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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
int	filtN_

double	filtR_

std::auto_ptr< fastjet::HEPTopTaggerV2 >	fjHEPTopTagger_

double	massRatioWidth_

double	maxCandMass_

double	maxFatjetAbsEta_

double	maxM13Cut_

double	maxR_

double	minCandMass_

double	minCandPt_

double	minFatjetPt_

double	minM13Cut_

double	minM23Cut_

double	minR_

double	minSubjetPt_

int	mode_

double	muCut_

bool	optimalR_

bool	qJets_

bool	rejectMinR_

double	subjetMass_

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from FastjetJetProducer
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

static void	fillDescriptionsFromFastJetProducer (edm::ParameterSetDescription &desc)

Static Public Member Functions inherited from VirtualJetProducer
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

static void	fillDescriptionsFromVirtualJetProducer (edm::ParameterSetDescription &desc)

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::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
int	activeAreaRepeats_

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_

bool	fromHTTTopJetProducer_

double	ghostArea_

double	ghostEtaMax_

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	rhoEtaMax_

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 77 of file HTTTopJetProducer.h.

Constructor & Destructor Documentation

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

Definition at line 11 of file HTTTopJetProducer.cc.

References edm::ParameterSet::exists(), filtN_, filtR_, fjHEPTopTagger_, VirtualJetProducer::fromHTTTopJetProducer_, edm::ParameterSet::getParameter(), massRatioWidth_, maxCandMass_, maxFatjetAbsEta_, maxM13Cut_, maxR_, minCandMass_, minCandPt_, minFatjetPt_, minM13Cut_, minM23Cut_, minR_, minSubjetPt_, mode_, muCut_, optimalR_, qJets_, rejectMinR_, subjetMass_, and verbose_.

                                                                :
        FastjetJetProducer( conf ),
        optimalR_(false),
        qJets_(false),
        minFatjetPt_(200.),
        minSubjetPt_(20.),
        minCandPt_(200.),
        maxFatjetAbsEta_(2.5),
        subjetMass_(30.),
        muCut_(0.8),
        filtR_(0.3),
        filtN_(5),       
        mode_(0),
        minCandMass_(150.),
        maxCandMass_(200.),
        massRatioWidth_(15),
        minM23Cut_(0.35),
        minM13Cut_(0.2),
        maxM13Cut_(1.3),
        maxR_(1.5),
        minR_(0.5),
        rejectMinR_(false),
        verbose_(false )
 {
   
   // Read in all the options from the configuration
   if ( conf.exists("optimalR") ) 
     optimalR_ = conf.getParameter<bool>("optimalR");
 
   if ( conf.exists("qJets") ) 
     qJets_ = conf.getParameter<bool>("qJets");
 
   if ( conf.exists("minFatjetPt") ) 
     minFatjetPt_ = conf.getParameter<double>("minFatjetPt");
   
   if ( conf.exists("minSubjetPt") ) 
     minSubjetPt_ = conf.getParameter<double>("minSubjetPt");
   
   if ( conf.exists("minCandPt") ) 
     minCandPt_ = conf.getParameter<double>("minCandPt");
   
   if ( conf.exists("maxFatjetAbsEta") )
     maxFatjetAbsEta_ = conf.getParameter<double>("maxFatjetAbsEta");
   
   if ( conf.exists("subjetMass") )
     subjetMass_ = conf.getParameter<double>("subjetMass");
   
   if ( conf.exists("muCut") )
     muCut_ = conf.getParameter<double>("muCut");
 
   if ( conf.exists("filtR") )
     filtR_ = conf.getParameter<double>("filtR");
 
   if ( conf.exists("filtN") )
     filtN_ = conf.getParameter<int>("filtN");
   
   if ( conf.exists("mode") )
     mode_ = conf.getParameter<int>("mode");
   
   if ( conf.exists("minCandMass") )
     minCandMass_ = conf.getParameter<double>("minCandMass");
   
   if ( conf.exists("maxCandMass") )
     maxCandMass_ = conf.getParameter<double>("maxCandMass");
   
   if ( conf.exists("massRatioWidth") )
     massRatioWidth_ = conf.getParameter<double>("massRatioWidth");
   
   if ( conf.exists("minM23Cut") )
     minM23Cut_ = conf.getParameter<double>("minM23Cut");
 
   if ( conf.exists("minM13Cut") )
     minM13Cut_ = conf.getParameter<double>("minM13Cut");
   
   if ( conf.exists("maxM13Cut") )
     maxM13Cut_ = conf.getParameter<double>("maxM13Cut");
 
   if ( conf.exists("maxR") )
     maxR_ = conf.getParameter<double>("maxR");
 
   if ( conf.exists("minR") )
     minR_ = conf.getParameter<double>("minR");
 
   if ( conf.exists("rejectMinR") )
     rejectMinR_ = conf.getParameter<bool>("rejectMinR");
 
   if ( conf.exists("verbose") )
     verbose_ = conf.getParameter<bool>("verbose");
   
   // Create the tagger-wrapper
   produces<HTTTopJetTagInfoCollection>();
 
   // Signal to the VirtualJetProducer that we have to add HTT information
   fromHTTTopJetProducer_ = 1;
   
   fjHEPTopTagger_ = std::auto_ptr<fastjet::HEPTopTaggerV2>(new fastjet::HEPTopTaggerV2(
                                                optimalR_,
                                                qJets_,
                                                minSubjetPt_, 
                                                minCandPt_,
                                                subjetMass_,         
                                                muCut_,          
                                                filtR_,
                                                filtN_,
                                                mode_,           
                                                minCandMass_,        
                                                maxCandMass_,        
                                                massRatioWidth_,         
                                                minM23Cut_,      
                                                minM13Cut_,      
                                                maxM13Cut_,
                                                rejectMinR_
                                                )); 
   
 }

virtual cms::HTTTopJetProducer::~HTTTopJetProducer ( )

inlinevirtual

Definition at line 83 of file HTTTopJetProducer.h.

References addHTTTopJetTagInfoCollection(), iEvent, produce(), and runAlgorithm().

83 {}

Member Function Documentation

void HTTTopJetProducer::addHTTTopJetTagInfoCollection	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		edm::OrphanHandle< reco::BasicJetCollection > &	oh
	)

virtual

Reimplemented from VirtualJetProducer.

Definition at line 186 of file HTTTopJetProducer.cc.

Referenced by ~HTTTopJetProducer().

                                                                             {
 
 
   // Set up output list
   auto tagInfos = std::make_unique<HTTTopJetTagInfoCollection>();
 
   // Loop over jets
   for (size_t ij=0; ij != fjJets_.size(); ij++){
 
     HTTTopJetProperties properties;
     HTTTopJetTagInfo tagInfo;
 
     // Black magic:
     // In the standard CA treatment the RefToBase is made from the handle directly
     // Since we only have a OrphanHandle (the JetCollection is created by this process) 
     // we have to take the detour via the Ref
     edm::Ref<reco::BasicJetCollection> ref(oh, ij);  
     edm::RefToBase<reco::Jet> rtb(ref);  
     
     fastjet::HEPTopTaggerV2Structure *s = (fastjet::HEPTopTaggerV2Structure*) fjJets_[ij].structure_non_const_ptr();
       
     properties.fjMass           = s->fj_mass();
     properties.fjPt             = s->fj_pt();
     properties.fjEta            = s->fj_eta();
     properties.fjPhi            = s->fj_phi();
     
     properties.topMass           = s->top_mass();
     properties.unfilteredMass    = s->unfiltered_mass();
     properties.prunedMass    = s->pruned_mass();
     properties.fRec      = s->fRec();
     properties.massRatioPassed   = s->mass_ratio_passed();
     
     properties.ropt           = s->ropt();
     properties.roptCalc           = s->roptCalc();
     properties.ptForRoptCalc      = s->ptForRoptCalc();
     
     properties.tau1Unfiltered     = s->Tau1Unfiltered();
     properties.tau2Unfiltered     = s->Tau2Unfiltered();
     properties.tau3Unfiltered     = s->Tau3Unfiltered();
     properties.tau1Filtered   = s->Tau1Filtered();
     properties.tau2Filtered   = s->Tau2Filtered();
     properties.tau3Filtered   = s->Tau3Filtered();
     properties.qWeight        = s->qweight();
     properties.qEpsilon       = s->qepsilon(); 
     properties.qSigmaM            = s->qsigmaM();
 
     tagInfo.insert(rtb, properties );
     tagInfos->push_back( tagInfo );   
   }
 
   iEvent.put(std::move(tagInfos));
   
 };

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

virtual

Reimplemented from FastjetJetProducer.

Definition at line 131 of file HTTTopJetProducer.cc.

References fjHEPTopTagger_, edm::RandomNumberGenerator::getEngine(), FastjetJetProducer::produce(), qJets_, and edm::Event::streamID().

Referenced by ~HTTTopJetProducer().

 {
 
   if (qJets_){
     edm::Service<edm::RandomNumberGenerator> rng;
     CLHEP::HepRandomEngine* engine = &rng->getEngine(e.streamID());
     fjHEPTopTagger_->set_rng(engine);
   }
 
   FastjetJetProducer::produce(e, c);
 }

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

virtual

Reimplemented from FastjetJetProducer.

Definition at line 143 of file HTTTopJetProducer.cc.

References gather_cfg::cout, VirtualJetProducer::doAreaFastjet_, VirtualJetProducer::doRhoFastjet_, VirtualJetProducer::fjAreaDefinition_, VirtualJetProducer::fjClusterSeq_, fjHEPTopTagger_, VirtualJetProducer::fjInputs_, VirtualJetProducer::fjJetDefinition_, VirtualJetProducer::fjJets_, mps_fire::i, maxFatjetAbsEta_, minFatjetPt_, perp(), verbose_, and VirtualJetProducer::voronoiRfact_.

Referenced by ~HTTTopJetProducer().

 {
 
   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_) ) );
   }
 
   //Run the jet clustering
   vector<fastjet::PseudoJet> inclusiveJets = fjClusterSeq_->inclusive_jets(minFatjetPt_);
 
   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() > minFatjetPt_ && fabs(inclusiveJets[i].rapidity()) < maxFatjetAbsEta_) {
       centralJets.push_back(inclusiveJets[i]);
     }
   }
 
   fastjet::HEPTopTaggerV2 & HEPTagger = *fjHEPTopTagger_;
 
   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;
 
     taggedJet = HEPTagger.result(*jetIt);
 
     if (taggedJet != 0){
       fjJets_.push_back(taggedJet);           
     }
   }
   
 }