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

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 79 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::auto_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 85 of file HTTTopJetProducer.h.

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

85 {}

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 145 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 201 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 90 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 102 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);           
     }
   }
   
 }