#include <CATopJetProducer.h>

Inheritance diagram for cms::CATopJetProducer:

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

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

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

	~CATopJetProducer () 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 Types
enum	tagalgos { CA_TOPTAGGER, FJ_CMS_TOPTAG, FJ_JHU_TOPTAG, FJ_NSUB_TAG }

Private Attributes
double	centralEtaCut_

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

std::unique_ptr< fastjet::JHTopTagger >	fjJHUTopTagger_

std::unique_ptr< fastjet::RestFrameNSubjettinessTagger >	fjNSUBTagger_

std::unique_ptr< CATopJetAlgorithm >	legacyCMSTopTagger_

double	ptMin_

int	tagAlgo_

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)

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

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

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	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 71 of file CATopJetProducer.h.

Member Enumeration Documentation

◆ tagalgos

enum cms::CATopJetProducer::tagalgos

private

Enumerator
CA_TOPTAGGER
FJ_CMS_TOPTAG
FJ_JHU_TOPTAG
FJ_NSUB_TAG

Definition at line 92 of file CATopJetProducer.h.

92 { CA_TOPTAGGER, FJ_CMS_TOPTAG, FJ_JHU_TOPTAG, FJ_NSUB_TAG };

Constructor & Destructor Documentation

◆ CATopJetProducer()

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

Definition at line 13 of file CATopJetProducer.cc.

     : 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::make_unique<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::make_unique<fastjet::CMSTopTagger>(conf.getParameter<double>("ptFrac"),
                                                               conf.getParameter<double>("rFrac"),
                                                               conf.getParameter<double>("adjacencyParam"));
   } else if (tagAlgo_ == FJ_JHU_TOPTAG) {
     fjJHUTopTagger_ = std::make_unique<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::make_unique<fastjet::RestFrameNSubjettinessTagger>(NsubJetDef,
                                                                             conf.getParameter<double>("tau2Cut"),
                                                                             conf.getParameter<double>("cosThetaSCut"),
                                                                             conf.getParameter<bool>("useExclusive"));
   }
 }

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

◆ ~CATopJetProducer()

cms::CATopJetProducer::~CATopJetProducer ( )

inlineoverride

Definition at line 75 of file CATopJetProducer.h.

75 {}

Member Function Documentation

◆ fillDescriptions()

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

static

Cambridge-Aachen top jet producer parameters

Definition at line 120 of file CATopJetProducer.cc.

                                                                                   {
   edm::ParameterSetDescription desc;
   desc.add<int>("tagAlgo", 0);             // choice of top tagging algorithm
   desc.add<double>("centralEtaCut", 2.5);  // eta for defining "central" jets
   desc.add<bool>("verbose", false);
   desc.add<string>("jetCollInstanceName", "caTopSubJets");  // subjet collection
   desc.add<int>("algorithm", 1);                            // 0 = KT, 1 = CA, 2 = anti-KT
   desc.add<int>("useAdjacency", 2);                         // veto adjacent subjets:
                                                             // 0,   no adjacency
                                                             //  1,  deltar adjacency
                                                             //  2,  modified adjacency
                                                             //  3,  calotower neirest neigbor based adjacency (untested)
   vector<double> sumEtBinsDefault = {0., 1600., 2600.};
   desc.add<vector<double>>(
       "sumEtBins",
       sumEtBinsDefault);  // sumEt bins over which cuts vary. vector={bin 0 lower bound, bin 1 lower bound, ...}
   vector<double> rBinsDefault(3, 0.8);
   desc.add<vector<double>>("rBins", rBinsDefault);  // Jet distance paramter R. R values depend on sumEt bins.
   vector<double> ptFracBinsDefault(3, 0.05);
   desc.add<vector<double>>("ptFracBins", ptFracBinsDefault);  // minimum fraction of central jet pt for subjets (deltap)
   vector<double> deltarBinsDefault(3, 0.019);
   desc.add<vector<double>>(
       "deltarBins", deltarBinsDefault);  // Applicable only if useAdjacency=1. deltar adjacency values for each sumEtBin
   vector<double> nCellBinsDefault(3, 1.9);
   desc.add<vector<double>>(
       "nCellBins",
       nCellBinsDefault);  // Applicable only if useAdjacency=3. number of cells apart for two subjets to be considered "independent"
   desc.add<bool>("useMaxTower", false);  // use max tower in adjacency criterion, otherwise use centroid - NOT USED
   desc.add<double>("sumEtEtaCut", 3.0);  // eta for event SumEt - NOT USED
   desc.add<double>("etFrac", 0.7);       // fraction of event sumEt / 2 for a jet to be considered "hard" - NOT USED
   desc.add<double>("ptFrac", 0.05);
   desc.add<double>("rFrac", 0.);
   desc.add<double>("adjacencyParam", 0.);
   desc.add<double>("deltaRCut", 0.19);
   desc.add<double>("cosThetaWMax", 0.7);
   desc.add<double>("tau2Cut", 0.);
   desc.add<double>("cosThetaSCut", 0.);
   desc.add<bool>("useExclusive", false);
   FastjetJetProducer::fillDescriptionsFromFastJetProducer(desc);
   VirtualJetProducer::fillDescriptionsFromVirtualJetProducer(desc);
   descriptions.add("CATopJetProducer", desc);
 }

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

◆ produce()

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

override

Definition at line 63 of file CATopJetProducer.cc.

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

References c, MillePedeFileConverter_cfg::e, and FastjetJetProducer::produce().

◆ runAlgorithm()

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

overridevirtual

Implements VirtualJetProducer.

Definition at line 65 of file CATopJetProducer.cc.

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

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

Member Data Documentation

◆ centralEtaCut_

double cms::CATopJetProducer::centralEtaCut_

private

Definition at line 90 of file CATopJetProducer.h.

Referenced by runAlgorithm().

◆ fjCMSTopTagger_

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

private

The algorithm to do the work.

Definition at line 84 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

◆ fjJHUTopTagger_

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

private

Definition at line 85 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

◆ fjNSUBTagger_

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

private

Definition at line 86 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

◆ legacyCMSTopTagger_

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

private

Definition at line 83 of file CATopJetProducer.h.

Referenced by CATopJetProducer().

◆ ptMin_

double cms::CATopJetProducer::ptMin_

private

Definition at line 89 of file CATopJetProducer.h.

Referenced by runAlgorithm().

◆ tagAlgo_

int cms::CATopJetProducer::tagAlgo_

private

Definition at line 88 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

◆ verbose_

bool cms::CATopJetProducer::verbose_

private

Definition at line 91 of file CATopJetProducer.h.

Referenced by runAlgorithm().

Public Member Functions

Static Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ tagalgos

Constructor & Destructor Documentation

◆ CATopJetProducer()

◆ ~CATopJetProducer()

Member Function Documentation

◆ fillDescriptions()

◆ produce()

◆ runAlgorithm()

Member Data Documentation

◆ centralEtaCut_

◆ fjCMSTopTagger_

◆ fjJHUTopTagger_

◆ fjNSUBTagger_

◆ legacyCMSTopTagger_

◆ ptMin_

◆ tagAlgo_

◆ verbose_