#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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<>
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

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)

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	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_

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_

bool	writeCompound_

bool	writeJetsWithConst_

Detailed Description

Definition at line 73 of file HTTTopJetProducer.h.

Constructor & Destructor Documentation

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

Definition at line 11 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::unique_ptr<fastjet::HEPTopTaggerV2>(new fastjet::HEPTopTaggerV2(optimalR_,
                                                                                          qJets_,
                                                                                          minSubjetPt_,
                                                                                          minCandPt_,
                                                                                          subjetMass_,
                                                                                          muCut_,
                                                                                          filtR_,
                                                                                          filtN_,
                                                                                          mode_,
                                                                                          minCandMass_,
                                                                                          maxCandMass_,
                                                                                          massRatioWidth_,
                                                                                          minM23Cut_,
                                                                                          minM13Cut_,
                                                                                          maxM13Cut_,
                                                                                          rejectMinR_));
 }

cms::HTTTopJetProducer::~HTTTopJetProducer ( )

inlineoverride

Definition at line 77 of file HTTTopJetProducer.h.

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

77 {}

Member Function Documentation

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

overridevirtual

Reimplemented from VirtualJetProducer.

Definition at line 134 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::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 186 of file HTTTopJetProducer.cc.

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

Referenced by ~HTTTopJetProducer().

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

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

override

Definition at line 81 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
	)

overridevirtual

Reimplemented from FastjetJetProducer.

Definition at line 91 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);
     }
   }
 }