#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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 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 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::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)

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

boost::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 77 of file CATopJetProducer.h.

Member Enumeration Documentation

enum cms::CATopJetProducer::tagalgos

private

Enumerator
CA_TOPTAGGER
FJ_CMS_TOPTAG
FJ_JHU_TOPTAG
FJ_NSUB_TAG

Definition at line 102 of file CATopJetProducer.h.

                   {
     CA_TOPTAGGER,
     FJ_CMS_TOPTAG,
     FJ_JHU_TOPTAG,
     FJ_NSUB_TAG
     };

Constructor & Destructor Documentation

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

Definition at line 11 of file CATopJetProducer.cc.

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

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

cms::CATopJetProducer::~CATopJetProducer ( )

inlineoverride

Definition at line 83 of file CATopJetProducer.h.

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

83 {}

Member Function Documentation

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

static

Cambridge-Aachen top jet producer parameters

Definition at line 133 of file CATopJetProducer.cc.

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

Referenced by ~CATopJetProducer().

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

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

override

Definition at line 77 of file CATopJetProducer.cc.

References FastjetJetProducer::produce().

Referenced by ~CATopJetProducer().

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

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

overridevirtual

Implements VirtualJetProducer.

Definition at line 82 of file CATopJetProducer.cc.

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

Referenced by ~CATopJetProducer().

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

Member Data Documentation

double cms::CATopJetProducer::centralEtaCut_

private

Definition at line 100 of file CATopJetProducer.h.

Referenced by runAlgorithm().

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

private

The algorithm to do the work.

Definition at line 92 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

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

private

Definition at line 93 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

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

private

Definition at line 94 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

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

private

Definition at line 91 of file CATopJetProducer.h.

Referenced by CATopJetProducer().

double cms::CATopJetProducer::ptMin_

private

Definition at line 99 of file CATopJetProducer.h.

Referenced by runAlgorithm().

int cms::CATopJetProducer::tagAlgo_

private

Definition at line 98 of file CATopJetProducer.h.

Referenced by CATopJetProducer(), and runAlgorithm().

bool cms::CATopJetProducer::verbose_

private

Definition at line 101 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

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation