#include <HTTTopJetProducer.h>

Inheritance diagram for cms::HTTTopJetProducer:

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

	HTTTopJetProducer (const edm::ParameterSet &ps)

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

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

	~HTTTopJetProducer () override

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

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

	~FastjetJetProducer () override

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

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

	VirtualJetProducer (const edm::ParameterSet &iConfig)

	~VirtualJetProducer () override

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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

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)

Private Attributes
int	filtN_

double	filtR_

std::unique_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 std::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 std::shared_ptr< fastjet::GhostedAreaSpec >	ActiveAreaSpecPtr

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

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

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

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

typedef std::shared_ptr< fastjet::Selector >	SelectorPtr

Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

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)

CaloGeometry const &	getGeometry (edm::EventSetup const &) const

HcalTopology const &	getTopology (edm::EventSetup const &) const

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)

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

Protected Attributes inherited from VirtualJetProducer
int	activeAreaRepeats_

bool	applyWeight_

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_

SelectorPtr	fjSelector_

bool	fromHTTTopJetProducer_ = false

double	ghostArea_

double	ghostEtaMax_

edm::EDGetTokenT< reco::VertexCollection >	input_vertex_token_

edm::EDGetTokenT< edm::ValueMap< float > >	input_weights_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_

std::shared_ptr< PileUpSubtractor >	subtractor_

bool	useDeterministicSeed_

bool	useExplicitGhosts_

int	verbosity_

reco::Particle::Point	vertex_

double	voronoiRfact_

edm::ValueMap< float >	weights_

bool	writeCompound_

bool	writeJetsWithConst_

Detailed Description

Definition at line 73 of file HTTTopJetProducer.h.

Constructor & Destructor Documentation

◆ HTTTopJetProducer()

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

Definition at line 13 of file HTTTopJetProducer.cc.

References 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
   optimalR_ = conf.getParameter<bool>("optimalR");
   qJets_ = conf.getParameter<bool>("qJets");
   minFatjetPt_ = conf.getParameter<double>("minFatjetPt");
   minSubjetPt_ = conf.getParameter<double>("minSubjetPt");
   minCandPt_ = conf.getParameter<double>("minCandPt");
   maxFatjetAbsEta_ = conf.getParameter<double>("maxFatjetAbsEta");
   subjetMass_ = conf.getParameter<double>("subjetMass");
   muCut_ = conf.getParameter<double>("muCut");
   filtR_ = conf.getParameter<double>("filtR");
   filtN_ = conf.getParameter<int>("filtN");
   mode_ = conf.getParameter<int>("mode");
   minCandMass_ = conf.getParameter<double>("minCandMass");
   maxCandMass_ = conf.getParameter<double>("maxCandMass");
   massRatioWidth_ = conf.getParameter<double>("massRatioWidth");
   minM23Cut_ = conf.getParameter<double>("minM23Cut");
   minM13Cut_ = conf.getParameter<double>("minM13Cut");
   maxM13Cut_ = conf.getParameter<double>("maxM13Cut");
   maxR_ = conf.getParameter<double>("maxR");
   minR_ = conf.getParameter<double>("minR");
   rejectMinR_ = conf.getParameter<bool>("rejectMinR");
   verbose_ = conf.getParameter<bool>("verbose");
 
   // Create the tagger-wrapper
   produces<HTTTopJetTagInfoCollection>();
 
   // Signal to the VirtualJetProducer that we have to add HTT information
   fromHTTTopJetProducer_ = true;
 
   fjHEPTopTagger_ = std::make_unique<fastjet::HEPTopTaggerV2>(optimalR_,
                                                               qJets_,
                                                               minSubjetPt_,
                                                               minCandPt_,
                                                               subjetMass_,
                                                               muCut_,
                                                               filtR_,
                                                               filtN_,
                                                               mode_,
                                                               minCandMass_,
                                                               maxCandMass_,
                                                               massRatioWidth_,
                                                               minM23Cut_,
                                                               minM13Cut_,
                                                               maxM13Cut_,
                                                               rejectMinR_);
 }

◆ ~HTTTopJetProducer()

cms::HTTTopJetProducer::~HTTTopJetProducer ( )

inlineoverride

Definition at line 77 of file HTTTopJetProducer.h.

77 {}

Member Function Documentation

◆ addHTTTopJetTagInfoCollection()

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

overridevirtual

Reimplemented from VirtualJetProducer.

Definition at line 136 of file HTTTopJetProducer.cc.

                                                                                                    {
   // 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));
 };

◆ fillDescriptions()

void HTTTopJetProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 188 of file HTTTopJetProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, FastjetJetProducer::fillDescriptionsFromFastJetProducer(), VirtualJetProducer::fillDescriptionsFromVirtualJetProducer(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                    {
   edm::ParameterSetDescription desc;
 
   desc.add<bool>("optimalR", true);
   desc.add<bool>("qJets", false);
   desc.add<double>("minFatjetPt", 200.);
   desc.add<double>("minSubjetPt", 0.);
   desc.add<double>("minCandPt", 0.);
   desc.add<double>("maxFatjetAbsEta", 99.);
   desc.add<double>("subjetMass", 30.);
   desc.add<double>("filtR", 0.3);
   desc.add<int>("filtN", 5);
   desc.add<int>("mode", 4);
   desc.add<double>("minCandMass", 0.);
   desc.add<double>("maxCandMass", 999999.);
   desc.add<double>("massRatioWidth", 999999.);
   desc.add<double>("minM23Cut", 0.);
   desc.add<double>("minM13Cut", 0.);
   desc.add<double>("maxM13Cut", 999999.);
   desc.add<double>("maxR", 1.5);
   desc.add<double>("minR", 0.5);
   desc.add<bool>("rejectMinR", false);
   desc.add<bool>("verbose", false);
 
   desc.add<int>("algorithm", 1);  // where is it needed?
 
   FastjetJetProducer::fillDescriptionsFromFastJetProducer(desc);
   desc.add<std::string>("jetCollInstanceName", "SubJets");
   VirtualJetProducer::fillDescriptionsFromVirtualJetProducer(desc);
 
   descriptions.add("HTTTopJetProducer", desc);
 }

◆ produce()

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

override

Definition at line 83 of file HTTTopJetProducer.cc.

References c, MillePedeFileConverter_cfg::e, fjHEPTopTagger_, edm::RandomNumberGenerator::getEngine(), FastjetJetProducer::produce(), and qJets_.

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

◆ runAlgorithm()

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

overridevirtual

Reimplemented from FastjetJetProducer.

Definition at line 93 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_.

                                                                                   {
   if (!doAreaFastjet_ && !doRhoFastjet_) {
     fjClusterSeq_ = std::make_shared<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);
     }
   }
 }