#include <SubEventGenJetProducer.h>

Inheritance diagram for cms::SubEventGenJetProducer:

Public Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

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

	SubEventGenJetProducer (const edm::ParameterSet &ps)

	~SubEventGenJetProducer () 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

Protected Member Functions
void	inputTowers () 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 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
std::vector< int >	hydroTag_

bool	ignoreHydro_

std::unique_ptr< std::vector< reco::GenJet > >	jets_

std::vector< int >	nSubParticles_

std::vector< std::vector< fastjet::PseudoJet > >	subInputs_

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_

Private Attributes
edm::EDGetTokenT< reco::CandidateView >	input_cand_token_

Additional Inherited Members
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

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

static void	fillDescriptionsFromVirtualJetProducer (edm::ParameterSetDescription &desc)

Detailed Description

Definition at line 19 of file SubEventGenJetProducer.h.

Constructor & Destructor Documentation

◆ SubEventGenJetProducer()

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

Definition at line 36 of file SubEventGenJetProducer.cc.

                                                                           : VirtualJetProducer(conf) {
   ignoreHydro_ = conf.getUntrackedParameter<bool>("ignoreHydro", true);
  
   // the subjet collections are set through the config file in the "jetCollInstanceName" field.
  
   input_cand_token_ = consumes<reco::CandidateView>(src_);
 }

References edm::ParameterSet::getUntrackedParameter(), ignoreHydro_, input_cand_token_, and VirtualJetProducer::src_.

◆ ~SubEventGenJetProducer()

cms::SubEventGenJetProducer::~SubEventGenJetProducer ( )

inlineoverride

Definition at line 22 of file SubEventGenJetProducer.h.

22 {}

Member Function Documentation

◆ inputTowers()

void SubEventGenJetProducer::inputTowers ( )

overrideprotectedvirtual

Reimplemented from VirtualJetProducer.

Definition at line 44 of file SubEventGenJetProducer.cc.

                                          {
   std::vector<edm::Ptr<reco::Candidate>>::const_iterator inBegin = inputs_.begin(), inEnd = inputs_.end(), i = inBegin;
   for (; i != inEnd; ++i) {
     reco::CandidatePtr input = inputs_[i - inBegin];
     if (edm::isNotFinite(input->pt()))
       continue;
     if (input->et() < inputEtMin_)
       continue;
     if (input->energy() < inputEMin_)
       continue;
     if (isAnomalousTower(input))
       continue;
  
     edm::Ptr<reco::Candidate> p = inputs_[i - inBegin];
     const GenParticle* pref = dynamic_cast<const GenParticle*>(p.get());
     int subevent = pref->collisionId();
     LogDebug("SubEventContainers") << "SubEvent is : " << subevent << endl;
     LogDebug("SubEventContainers") << "SubSize is : " << subInputs_.size() << endl;
  
     if (subevent >= (int)subInputs_.size()) {
       hydroTag_.resize(subevent + 1, -1);
       subInputs_.resize(subevent + 1);
       LogDebug("SubEventContainers") << "SubSize is : " << subInputs_.size() << endl;
       LogDebug("SubEventContainers") << "HydroTagSize is : " << hydroTag_.size() << endl;
     }
  
     LogDebug("SubEventContainers") << "HydroTag is : " << hydroTag_[subevent] << endl;
     if (hydroTag_[subevent] != 0)
       hydroTag_[subevent] = (int)checkHydro(pref);
  
     subInputs_[subevent].push_back(fastjet::PseudoJet(input->px(), input->py(), input->pz(), input->energy()));
  
     subInputs_[subevent].back().set_user_index(i - inBegin);
   }
 }

References reco::GenParticle::collisionId(), hydroTag_, mps_fire::i, input, VirtualJetProducer::inputEMin_, VirtualJetProducer::inputEtMin_, VirtualJetProducer::inputs_, createfilelist::int, VirtualJetProducer::isAnomalousTower(), edm::isNotFinite(), LogDebug, AlCaHLTBitMon_ParallelJobs::p, and subInputs_.

Referenced by produce().

◆ produce()

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

override

Definition at line 80 of file SubEventGenJetProducer.cc.

                                                                                   {
   LogDebug("VirtualJetProducer") << "Entered produce\n";
  
   fjJets_.clear();
   subInputs_.clear();
   nSubParticles_.clear();
   hydroTag_.clear();
   inputs_.clear();
  
   // get inputs and convert them to the fastjet format (fastjet::PeudoJet)
   edm::Handle<reco::CandidateView> inputsHandle;
   iEvent.getByToken(input_cand_token_, inputsHandle);
   for (size_t i = 0; i < inputsHandle->size(); ++i) {
     inputs_.push_back(inputsHandle->ptrAt(i));
   }
   LogDebug("VirtualJetProducer") << "Got inputs\n";
  
   inputTowers();
   // Convert candidates to fastjet::PseudoJets.
   // Also correct to Primary Vertex. Will modify fjInputs_
   // and use inputs_
  
  
   jets_ = std::make_unique<std::vector<GenJet>>();
  
   LogDebug("VirtualJetProducer") << "Inputted towers\n";
  
   size_t nsub = subInputs_.size();
  
   for (size_t isub = 0; isub < nsub; ++isub) {
     if (ignoreHydro_ && hydroTag_[isub])
       continue;
     fjJets_.clear();
     fjInputs_.clear();
     fjInputs_ = subInputs_[isub];
     runAlgorithm(iEvent, iSetup);
   }
  
   //Finalize
   LogDebug("SubEventJetProducer") << "Wrote jets\n";
  
   iEvent.put(std::move(jets_));
   return;
 }

References VirtualJetProducer::fjInputs_, VirtualJetProducer::fjJets_, hydroTag_, mps_fire::i, iEvent, ignoreHydro_, input_cand_token_, VirtualJetProducer::inputs_, inputTowers(), jets_, LogDebug, eostools::move(), nsub, nSubParticles_, edm::View< T >::ptrAt(), runAlgorithm(), edm::View< T >::size(), and subInputs_.

◆ runAlgorithm()

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

overridevirtual

Implements VirtualJetProducer.

Definition at line 126 of file SubEventGenJetProducer.cc.

                                                                                        {
   // run algorithm
   fjJets_.clear();
  
   fjClusterSeq_ = std::make_shared<fastjet::ClusterSequence>(fjInputs_, *fjJetDefinition_);
   fjJets_ = fastjet::sorted_by_pt(fjClusterSeq_->inclusive_jets(jetPtMin_));
  
   using namespace reco;
  
   for (unsigned int ijet = 0; ijet < fjJets_.size(); ++ijet) {
     GenJet jet;
     const fastjet::PseudoJet& fjJet = fjJets_[ijet];
  
     std::vector<fastjet::PseudoJet> fjConstituents = sorted_by_pt(fjClusterSeq_->constituents(fjJet));
  
     std::vector<CandidatePtr> constituents = getConstituents(fjConstituents);
  
     double px = fjJet.px();
     double py = fjJet.py();
     double pz = fjJet.pz();
     double E = fjJet.E();
     double jetArea = 0.0;
     double pu = 0.;
  
     writeSpecific(jet, Particle::LorentzVector(px, py, pz, E), vertex_, constituents, iSetup);
  
     jet.setJetArea(jetArea);
     jet.setPileup(pu);
  
     jets_->push_back(jet);
   }
 }

References VirtualJetProducer::fjClusterSeq_, VirtualJetProducer::fjInputs_, VirtualJetProducer::fjJetDefinition_, VirtualJetProducer::fjJets_, VirtualJetProducer::getConstituents(), metsig::jet, VirtualJetProducer::jetPtMin_, jets_, muons2muons_cfi::pu, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, VirtualJetProducer::vertex_, and reco::writeSpecific().

Referenced by produce().

Member Data Documentation

◆ hydroTag_

std::vector<int> cms::SubEventGenJetProducer::hydroTag_

protected

Definition at line 29 of file SubEventGenJetProducer.h.

Referenced by inputTowers(), and produce().

◆ ignoreHydro_

bool cms::SubEventGenJetProducer::ignoreHydro_

protected

Definition at line 31 of file SubEventGenJetProducer.h.

Referenced by produce(), and SubEventGenJetProducer().

◆ input_cand_token_

edm::EDGetTokenT<reco::CandidateView> cms::SubEventGenJetProducer::input_cand_token_

private

Definition at line 39 of file SubEventGenJetProducer.h.

Referenced by produce(), and SubEventGenJetProducer().

◆ jets_

std::unique_ptr<std::vector<reco::GenJet> > cms::SubEventGenJetProducer::jets_

protected

Definition at line 28 of file SubEventGenJetProducer.h.

Referenced by produce(), and runAlgorithm().

◆ nSubParticles_

std::vector<int> cms::SubEventGenJetProducer::nSubParticles_

protected

Definition at line 30 of file SubEventGenJetProducer.h.

Referenced by produce().

◆ subInputs_

std::vector<std::vector<fastjet::PseudoJet> > cms::SubEventGenJetProducer::subInputs_

protected

Definition at line 27 of file SubEventGenJetProducer.h.

Referenced by inputTowers(), and produce().

Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SubEventGenJetProducer()

◆ ~SubEventGenJetProducer()

Member Function Documentation

◆ inputTowers()

◆ produce()

◆ runAlgorithm()

Member Data Documentation

◆ hydroTag_

◆ ignoreHydro_

◆ input_cand_token_

◆ jets_

◆ nSubParticles_

◆ subInputs_