#include <VirtualJetProducer.h>

Inheritance diagram for VirtualJetProducer:

Classes
struct	JetType

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

Public Member Functions
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	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

static void	fillDescriptionsFromVirtualJetProducer (edm::ParameterSetDescription &desc)

Protected Member Functions
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)

virtual void	runAlgorithm (edm::Event &iEvent, const edm::EventSetup &iSetup)=0

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
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
std::unique_ptr< AnomalousTower >	anomalousTowerDef_

edm::EDGetTokenT< std::vector< edm::FwdPtr< reco::PFCandidate > > >	input_candidatefwdptr_token_

edm::EDGetTokenT< reco::CandidateView >	input_candidateview_token_

edm::EDGetTokenT< std::vector< edm::FwdPtr< reco::GenParticle > > >	input_gencandidatefwdptr_token_

edm::EDGetTokenT< std::vector< edm::FwdPtr< pat::PackedCandidate > > >	input_packedcandidatefwdptr_token_

edm::EDGetTokenT< std::vector< edm::FwdPtr< pat::PackedGenParticle > > >	input_packedgencandidatefwdptr_token_

Detailed Description

Definition at line 35 of file VirtualJetProducer.h.

Member Typedef Documentation

◆ ActiveAreaSpecPtr

typedef std::shared_ptr<fastjet::GhostedAreaSpec> VirtualJetProducer::ActiveAreaSpecPtr

Definition at line 76 of file VirtualJetProducer.h.

◆ AreaDefinitionPtr

typedef std::shared_ptr<fastjet::AreaDefinition> VirtualJetProducer::AreaDefinitionPtr

Definition at line 77 of file VirtualJetProducer.h.

◆ ClusterSequencePtr

typedef std::shared_ptr<fastjet::ClusterSequence> VirtualJetProducer::ClusterSequencePtr

Definition at line 73 of file VirtualJetProducer.h.

◆ JetDefPtr

typedef std::shared_ptr<fastjet::JetDefinition> VirtualJetProducer::JetDefPtr

Definition at line 75 of file VirtualJetProducer.h.

◆ PluginPtr

typedef std::shared_ptr<fastjet::JetDefinition::Plugin> VirtualJetProducer::PluginPtr

Definition at line 74 of file VirtualJetProducer.h.

◆ SelectorPtr

typedef std::shared_ptr<fastjet::Selector> VirtualJetProducer::SelectorPtr

Definition at line 78 of file VirtualJetProducer.h.

Constructor & Destructor Documentation

◆ VirtualJetProducer()

VirtualJetProducer::VirtualJetProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 113 of file VirtualJetProducer.cc.

                                                                      {
   moduleLabel_ = iConfig.getParameter<string>("@module_label");
   src_ = iConfig.getParameter<edm::InputTag>("src");
   srcPVs_ = iConfig.getParameter<edm::InputTag>("srcPVs");
   jetType_ = iConfig.getParameter<string>("jetType");
   jetAlgorithm_ = iConfig.getParameter<string>("jetAlgorithm");
   rParam_ = iConfig.getParameter<double>("rParam");
   inputEtMin_ = iConfig.getParameter<double>("inputEtMin");
   inputEMin_ = iConfig.getParameter<double>("inputEMin");
   jetPtMin_ = iConfig.getParameter<double>("jetPtMin");
   doPVCorrection_ = iConfig.getParameter<bool>("doPVCorrection");
   doAreaFastjet_ = iConfig.getParameter<bool>("doAreaFastjet");
   doRhoFastjet_ = iConfig.getParameter<bool>("doRhoFastjet");
   jetCollInstanceName_ = iConfig.getParameter<string>("jetCollInstanceName");
   doPUOffsetCorr_ = iConfig.getParameter<bool>("doPUOffsetCorr");
   puSubtractorName_ = iConfig.getParameter<string>("subtractorName");
   useExplicitGhosts_ =
       iConfig.getParameter<bool>("useExplicitGhosts");  // use explicit ghosts in the fastjet clustering sequence?
   doAreaDiskApprox_ = iConfig.getParameter<bool>("doAreaDiskApprox");
   voronoiRfact_ = iConfig.getParameter<double>(
       "voronoiRfact");  // Voronoi-based area calculation allows for an empirical scale factor
   rhoEtaMax_ = iConfig.getParameter<double>(
       "Rho_EtaMax");  // do fasjet area / rho calcluation? => accept corresponding parameters
   ghostEtaMax_ = iConfig.getParameter<double>("Ghost_EtaMax");
   activeAreaRepeats_ = iConfig.getParameter<int>("Active_Area_Repeats");
   ghostArea_ = iConfig.getParameter<double>("GhostArea");
   restrictInputs_ = iConfig.getParameter<bool>("restrictInputs");  // restrict inputs to first "maxInputs" towers?
   maxInputs_ = iConfig.getParameter<unsigned int>("maxInputs");
   writeCompound_ = iConfig.getParameter<bool>(
       "writeCompound");  // Check to see if we are writing compound jets for substructure and jet grooming
   writeJetsWithConst_ = iConfig.getParameter<bool>("writeJetsWithConst");  //write subtracted jet constituents
   doFastJetNonUniform_ = iConfig.getParameter<bool>("doFastJetNonUniform");
   puCenters_ = iConfig.getParameter<vector<double>>("puCenters");
   puWidth_ = iConfig.getParameter<double>("puWidth");
   nExclude_ = iConfig.getParameter<unsigned int>("nExclude");
   useDeterministicSeed_ = iConfig.getParameter<bool>("useDeterministicSeed");
   minSeed_ = iConfig.getParameter<unsigned int>("minSeed");
   verbosity_ = iConfig.getParameter<int>("verbosity");
   applyWeight_ = iConfig.getParameter<bool>("applyWeight");
   if (applyWeight_) {
     edm::InputTag srcWeights = iConfig.getParameter<edm::InputTag>("srcWeights");
     if (srcWeights.label() == src_.label())
       LogWarning("VirtualJetProducer")
           << "Particle and weights collection have the same label. You may be applying the same weights twice.\n";
     if (!srcWeights.label().empty())
       input_weights_token_ = consumes<edm::ValueMap<float>>(srcWeights);
   }
  
   anomalousTowerDef_ = std::make_unique<AnomalousTower>(iConfig);
  
   input_vertex_token_ = consumes<reco::VertexCollection>(srcPVs_);
   input_candidateview_token_ = consumes<reco::CandidateView>(src_);
   input_candidatefwdptr_token_ =
       consumes<vector<edm::FwdPtr<reco::PFCandidate>>>(iConfig.getParameter<edm::InputTag>("src"));
   input_packedcandidatefwdptr_token_ =
       consumes<vector<edm::FwdPtr<pat::PackedCandidate>>>(iConfig.getParameter<edm::InputTag>("src"));
   input_gencandidatefwdptr_token_ =
       consumes<vector<edm::FwdPtr<reco::GenParticle>>>(iConfig.getParameter<edm::InputTag>("src"));
   input_packedgencandidatefwdptr_token_ =
       consumes<vector<edm::FwdPtr<pat::PackedGenParticle>>>(iConfig.getParameter<edm::InputTag>("src"));
  
   //
   // additional parameters to think about:
   // - overlap threshold (set to 0.75 for the time being)
   // - p parameter for generalized kT (set to -2 for the time being)
   // - fastjet PU subtraction parameters (not yet considered)
   //
   if (jetAlgorithm_ == "Kt")
     fjJetDefinition_ = std::make_shared<fastjet::JetDefinition>(fastjet::kt_algorithm, rParam_);
  
   else if (jetAlgorithm_ == "CambridgeAachen")
     fjJetDefinition_ = std::make_shared<fastjet::JetDefinition>(fastjet::cambridge_algorithm, rParam_);
  
   else if (jetAlgorithm_ == "AntiKt")
     fjJetDefinition_ = std::make_shared<fastjet::JetDefinition>(fastjet::antikt_algorithm, rParam_);
  
   else if (jetAlgorithm_ == "GeneralizedKt")
     fjJetDefinition_ = std::make_shared<fastjet::JetDefinition>(fastjet::genkt_algorithm, rParam_, -2);
  
   else if (jetAlgorithm_ == "SISCone") {
     fjPlugin_ = PluginPtr(new fastjet::SISConePlugin(rParam_, 0.75, 0, 0.0, false, fastjet::SISConePlugin::SM_pttilde));
     fjJetDefinition_ = std::make_shared<fastjet::JetDefinition>(&*fjPlugin_);
  
   } else if (jetAlgorithm_ == "IterativeCone") {
     fjPlugin_ = PluginPtr(new fastjet::CMSIterativeConePlugin(rParam_, 1.0));
     fjJetDefinition_ = std::make_shared<fastjet::JetDefinition>(&*fjPlugin_);
  
   } else if (jetAlgorithm_ == "CDFMidPoint") {
     fjPlugin_ = PluginPtr(new fastjet::CDFMidPointPlugin(rParam_, 0.75));
     fjJetDefinition_ = std::make_shared<fastjet::JetDefinition>(&*fjPlugin_);
  
   } else if (jetAlgorithm_ == "ATLASCone") {
     fjPlugin_ = PluginPtr(new fastjet::ATLASConePlugin(rParam_));
     fjJetDefinition_ = std::make_shared<fastjet::JetDefinition>(&*fjPlugin_);
  
   } else {
     throw cms::Exception("Invalid jetAlgorithm") << "Jet algorithm for VirtualJetProducer is invalid, Abort!\n";
   }
  
   jetTypeE = JetType::byName(jetType_);
  
   if (doPUOffsetCorr_) {
     if (puSubtractorName_.empty()) {
       LogWarning("VirtualJetProducer")
           << "Pile Up correction on; however, pile up type is not specified. Using default... \n";
       subtractor_ = std::make_shared<PileUpSubtractor>(iConfig, consumesCollector());
     } else
       subtractor_ = std::shared_ptr<PileUpSubtractor>(
           PileUpSubtractorFactory::get()->create(puSubtractorName_, iConfig, consumesCollector()));
   }
  
   // do approximate disk-based area calculation => warn if conflicting request
   if (doAreaDiskApprox_ && doAreaFastjet_)
     throw cms::Exception("Conflicting area calculations")
         << "Both the calculation of jet area via fastjet and via an analytical disk approximation have been requested. "
            "Please decide on one.\n";
  
   if (doAreaFastjet_ || doRhoFastjet_) {
     if (voronoiRfact_ <= 0) {
       fjActiveArea_ = std::make_shared<fastjet::GhostedAreaSpec>(ghostEtaMax_, activeAreaRepeats_, ghostArea_);
  
       if (!useExplicitGhosts_) {
         fjAreaDefinition_ = std::make_shared<fastjet::AreaDefinition>(fastjet::active_area, *fjActiveArea_);
       } else {
         fjAreaDefinition_ =
             std::make_shared<fastjet::AreaDefinition>(fastjet::active_area_explicit_ghosts, *fjActiveArea_);
       }
     }
     fjSelector_ = std::make_shared<fastjet::Selector>(fastjet::SelectorAbsRapMax(rhoEtaMax_));
   }
  
   if ((doFastJetNonUniform_) && (puCenters_.empty()))
     throw cms::Exception("doFastJetNonUniform")
         << "Parameter puCenters for doFastJetNonUniform is not defined." << std::endl;
  
   // make the "produces" statements
   makeProduces(moduleLabel_, jetCollInstanceName_);
   produces<vector<double>>("rhos");
   produces<vector<double>>("sigmas");
   produces<double>("rho");
   produces<double>("sigma");
 }

References beamerCreator::create(), Exception, get, edm::ParameterSet::getParameter(), and HLT_FULL_cff::srcWeights.

◆ ~VirtualJetProducer()

VirtualJetProducer::~VirtualJetProducer ( )

override

Definition at line 257 of file VirtualJetProducer.cc.

257 {}

Member Function Documentation

◆ addHTTTopJetTagInfoCollection()

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

inlineprotectedvirtual

Reimplemented in cms::HTTTopJetProducer.

Definition at line 113 of file VirtualJetProducer.h.

115 {};

◆ copyConstituents()

void VirtualJetProducer::copyConstituents	(	const std::vector< fastjet::PseudoJet > &	fjConstituents,
		reco::Jet *	jet
	)

protectedvirtual

Definition at line 513 of file VirtualJetProducer.cc.

                                                                                                       {
   for (unsigned int i = 0; i < fjConstituents.size(); ++i) {
     int index = fjConstituents[i].user_index();
     if (index >= 0 && static_cast<unsigned int>(index) < inputs_.size())
       jet->addDaughter(inputs_[index]);
   }
 }

References mps_fire::i, and metsig::jet.

◆ fillDescriptions()

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

static

Definition at line 1004 of file VirtualJetProducer.cc.

                                                                                     {
   edm::ParameterSetDescription desc;
   fillDescriptionsFromVirtualJetProducer(desc);
   desc.add<string>("jetCollInstanceName", "");
  
   // addDefault must be used here instead of add unless
   // all the classes that inherit from this class redefine
   // the fillDescriptions function. Otherwise, the autogenerated
   // cfi filenames are the same and conflict.
   descriptions.addDefault(desc);
 }

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

◆ fillDescriptionsFromVirtualJetProducer()

void VirtualJetProducer::fillDescriptionsFromVirtualJetProducer ( edm::ParameterSetDescription & desc )

static

Definition at line 1016 of file VirtualJetProducer.cc.

                                                                                                 {
   desc.add<edm::InputTag>("src", edm::InputTag("particleFlow"));
   desc.add<edm::InputTag>("srcPVs", edm::InputTag(""));
   desc.add<string>("jetType", "PFJet");
   desc.add<string>("jetAlgorithm", "AntiKt");
   desc.add<double>("rParam", 0.4);
   desc.add<double>("inputEtMin", 0.0);
   desc.add<double>("inputEMin", 0.0);
   desc.add<double>("jetPtMin", 5.);
   desc.add<bool>("doPVCorrection", false);
   desc.add<bool>("doAreaFastjet", false);
   desc.add<bool>("doRhoFastjet", false);
   desc.add<bool>("doPUOffsetCorr", false);
   desc.add<double>("puPtMin", 10.);
   desc.add<double>("nSigmaPU", 1.0);
   desc.add<double>("radiusPU", 0.5);
   desc.add<string>("subtractorName", "");
   desc.add<bool>("useExplicitGhosts", false);
   desc.add<bool>("doAreaDiskApprox", false);
   desc.add<double>("voronoiRfact", -0.9);
   desc.add<double>("Rho_EtaMax", 4.4);
   desc.add<double>("Ghost_EtaMax", 5.);
   desc.add<int>("Active_Area_Repeats", 1);
   desc.add<double>("GhostArea", 0.01);
   desc.add<bool>("restrictInputs", false);
   desc.add<unsigned int>("maxInputs", 1);
   desc.add<bool>("writeCompound", false);
   desc.add<bool>("writeJetsWithConst", false);
   desc.add<bool>("doFastJetNonUniform", false);
   desc.add<bool>("useDeterministicSeed", false);
   desc.add<unsigned int>("minSeed", 14327);
   desc.add<int>("verbosity", 0);
   desc.add<double>("puWidth", 0.);
   desc.add<unsigned int>("nExclude", 0);
   desc.add<unsigned int>("maxBadEcalCells", 9999999);
   desc.add<unsigned int>("maxBadHcalCells", 9999999);
   desc.add<unsigned int>("maxProblematicEcalCells", 9999999);
   desc.add<unsigned int>("maxProblematicHcalCells", 9999999);
   desc.add<unsigned int>("maxRecoveredEcalCells", 9999999);
   desc.add<unsigned int>("maxRecoveredHcalCells", 9999999);
   vector<double> puCentersDefault;
   desc.add<vector<double>>("puCenters", puCentersDefault);
   desc.add<bool>("applyWeight", false);
   desc.add<edm::InputTag>("srcWeights", edm::InputTag(""));
   desc.add<double>("minimumTowersFraction", 0.);
 }

References submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

Referenced by cms::CSJetProducer::fillDescriptions(), FastjetJetProducer::fillDescriptions(), cms::CATopJetProducer::fillDescriptions(), and cms::HTTTopJetProducer::fillDescriptions().

◆ getConstituents()

vector< reco::CandidatePtr > VirtualJetProducer::getConstituents ( const std::vector< fastjet::PseudoJet > & fjConstituents )

protectedvirtual

Definition at line 522 of file VirtualJetProducer.cc.

                                                                                                            {
   vector<reco::CandidatePtr> result;
   result.reserve(fjConstituents.size() / 2);
   for (unsigned int i = 0; i < fjConstituents.size(); i++) {
     auto index = fjConstituents[i].user_index();
     if (index >= 0 && static_cast<unsigned int>(index) < inputs_.size()) {
       result.emplace_back(inputs_[index]);
     }
   }
   return result;
 }

References mps_fire::i, and mps_fire::result.

Referenced by FastjetJetProducer::produceTrackJets(), and cms::SubEventGenJetProducer::runAlgorithm().

◆ inputTowers()

void VirtualJetProducer::inputTowers ( )

protectedvirtual

Reimplemented in cms::CompoundJetProducer, SubjetFilterJetProducer, and cms::SubEventGenJetProducer.

Definition at line 442 of file VirtualJetProducer.cc.

                                      {
   auto inBegin = inputs_.begin(), inEnd = inputs_.end(), i = inBegin;
   for (; i != inEnd; ++i) {
     auto const& input = **i;
     // std::cout << "CaloTowerVI jets " << input->pt() << " " << input->et() << ' '<< input->energy() << ' ' << (isAnomalousTower(input) ? " bad" : " ok") << std::endl;
     if (edm::isNotFinite(input.pt()))
       continue;
     if (input.et() < inputEtMin_)
       continue;
     if (input.energy() < inputEMin_)
       continue;
     if (isAnomalousTower(*i))
       continue;
     // Change by SRR : this is no longer an error nor warning, this can happen with PU mitigation algos.
     // Also switch to something more numerically safe. (VI: 10^-42GeV????)
     if (input.pt() < 100 * std::numeric_limits<double>::epsilon()) {
       continue;
     }
     if (makeCaloJet(jetTypeE) && doPVCorrection_) {
       const CaloTower& tower = dynamic_cast<const CaloTower&>(input);
       auto const& ct = tower.p4(vertex_);  // very expensive as computed in eta/phi
       fjInputs_.emplace_back(ct.px(), ct.py(), ct.pz(), ct.energy());
       fjInputs_.back().set_user_index(i - inBegin);
       //std::cout << "tower:" << *tower << '\n';
     } else {
       /*
       if(makePFJet(jetTypeE)) {
     reco::PFCandidate& pfc = (reco::PFCandidate&)input;
     std::cout << "PF cand:" << pfc << '\n';
       }
       */
       if (!applyWeight_) {
         fjInputs_.emplace_back(input.px(), input.py(), input.pz(), input.energy());
         fjInputs_.back().set_user_index(i - inBegin);
       } else {
         if (input_weights_token_.isUninitialized())
           throw cms::Exception("InvalidInput")
               << "applyWeight set to True, but no weights given in VirtualJetProducer\n";
         float w = weights_[*i];
         if (w > 0) {
           fjInputs_.emplace_back(input.px() * w, input.py() * w, input.pz() * w, input.energy() * w);
           fjInputs_.back().set_user_index(i - inBegin);
         }
       }
     }
   }
  
   if (restrictInputs_ && fjInputs_.size() > maxInputs_) {
     reco::helper::GreaterByPtPseudoJet pTComparator;
     std::sort(fjInputs_.begin(), fjInputs_.end(), pTComparator);
     fjInputs_.resize(maxInputs_);
     edm::LogWarning("JetRecoTooManyEntries")
         << "Too many inputs in the event, limiting to first " << maxInputs_ << ". Output is suspect.";
   }
 }

References geometryDiff::epsilon, mps_fire::i, input, edm::isNotFinite(), jetUpdater_cfi::sort, hgcalTowerProducer_cfi::tower, and w.

Referenced by SubjetFilterJetProducer::inputTowers(), cms::CompoundJetProducer::inputTowers(), and FastjetJetProducer::produceTrackJets().

◆ isAnomalousTower()

bool VirtualJetProducer::isAnomalousTower ( reco::CandidatePtr input )

protectedvirtual

Definition at line 499 of file VirtualJetProducer.cc.

                                                                 {
   if (!makeCaloJet(jetTypeE))
     return false;
   else
     return (*anomalousTowerDef_)(*input);
 }

Referenced by cms::SubEventGenJetProducer::inputTowers().

◆ jetType()

std::string VirtualJetProducer::jetType ( ) const

inline

Definition at line 85 of file VirtualJetProducer.h.

85 { return jetType_; }

◆ makeBasicJet()

bool VirtualJetProducer::makeBasicJet ( const JetType::Type & fTag )

inlineprotected

Definition at line 61 of file VirtualJetProducer.h.

61 { return fTag == JetType::BasicJet; }

◆ makeCaloJet()

bool VirtualJetProducer::makeCaloJet ( const JetType::Type & fTag )

inlineprotected

Definition at line 56 of file VirtualJetProducer.h.

56 { return fTag == JetType::CaloJet; }

References btagElecInJet_cfi::CaloJet.

◆ makeGenJet()

bool VirtualJetProducer::makeGenJet ( const JetType::Type & fTag )

inlineprotected

Definition at line 58 of file VirtualJetProducer.h.

58 { return fTag == JetType::GenJet; }

References nanoDQM_cfi::GenJet.

◆ makePFClusterJet()

bool VirtualJetProducer::makePFClusterJet ( const JetType::Type & fTag )

inlineprotected

Definition at line 60 of file VirtualJetProducer.h.

60 { return fTag == JetType::PFClusterJet; }

◆ makePFJet()

bool VirtualJetProducer::makePFJet ( const JetType::Type & fTag )

inlineprotected

Definition at line 57 of file VirtualJetProducer.h.

57 { return fTag == JetType::PFJet; }

◆ makeProduces()

void VirtualJetProducer::makeProduces	(	std::string	s,
		std::string	tag = `""`
	)

protectedvirtual

Definition at line 83 of file VirtualJetProducer.cc.

                                                                     {
   if (writeCompound_) {
     produces<reco::BasicJetCollection>();
   }
  
   if (writeJetsWithConst_) {
     produces<reco::PFCandidateCollection>(tag).setBranchAlias(alias);
     produces<reco::PFJetCollection>();
   } else {
     if (makeCaloJet(jetTypeE)) {
       produces<reco::CaloJetCollection>(tag).setBranchAlias(alias);
     } else if (makePFJet(jetTypeE)) {
       produces<reco::PFJetCollection>(tag).setBranchAlias(alias);
     } else if (makeGenJet(jetTypeE)) {
       produces<reco::GenJetCollection>(tag).setBranchAlias(alias);
     } else if (makeTrackJet(jetTypeE)) {
       produces<reco::TrackJetCollection>(tag).setBranchAlias(alias);
     } else if (makePFClusterJet(jetTypeE)) {
       produces<reco::PFClusterJetCollection>(tag).setBranchAlias(alias);
     } else if (makeBasicJet(jetTypeE)) {
       produces<reco::BasicJetCollection>(tag).setBranchAlias(alias);
     }
   }
 }

References SiStripOfflineCRack_cfg::alias, and GlobalPosition_Frontier_DevDB_cff::tag.

Referenced by SubjetFilterJetProducer::SubjetFilterJetProducer().

◆ makeTrackJet()

bool VirtualJetProducer::makeTrackJet ( const JetType::Type & fTag )

inlineprotected

Definition at line 59 of file VirtualJetProducer.h.

59 { return fTag == JetType::TrackJet; }

Referenced by FastjetJetProducer::produce().

◆ offsetCorrectJets()

void VirtualJetProducer::offsetCorrectJets ( std::vector< fastjet::PseudoJet > & orphanInput )

protected

◆ output()

void VirtualJetProducer::output	(	edm::Event &	iEvent,
		edm::EventSetup const &	iSetup
	)

protectedvirtual

Reimplemented in cms::CompoundJetProducer, and SubjetFilterJetProducer.

Definition at line 536 of file VirtualJetProducer.cc.

                                                                              {
   // Write jets and constitutents. Will use fjJets_, inputs_
   // and fjClusterSeq_
  
   if (writeCompound_) {
     // Write jets and subjets
     switch (jetTypeE) {
       case JetType::CaloJet:
         writeCompoundJets<reco::CaloJet>(iEvent, iSetup);
         break;
       case JetType::PFJet:
         writeCompoundJets<reco::PFJet>(iEvent, iSetup);
         break;
       case JetType::GenJet:
         writeCompoundJets<reco::GenJet>(iEvent, iSetup);
         break;
       case JetType::BasicJet:
         writeCompoundJets<reco::BasicJet>(iEvent, iSetup);
         break;
       default:
         throw cms::Exception("InvalidInput") << "invalid jet type in CompoundJetProducer\n";
         break;
     };
   } else if (writeJetsWithConst_) {
     // Write jets and new constituents.
     writeJetsWithConstituents<reco::PFJet>(iEvent, iSetup);
   } else {
     switch (jetTypeE) {
       case JetType::CaloJet:
         writeJets<reco::CaloJet>(iEvent, iSetup);
         break;
       case JetType::PFJet:
         writeJets<reco::PFJet>(iEvent, iSetup);
         break;
       case JetType::GenJet:
         writeJets<reco::GenJet>(iEvent, iSetup);
         break;
       case JetType::TrackJet:
         writeJets<reco::TrackJet>(iEvent, iSetup);
         break;
       case JetType::PFClusterJet:
         writeJets<reco::PFClusterJet>(iEvent, iSetup);
         break;
       case JetType::BasicJet:
         writeJets<reco::BasicJet>(iEvent, iSetup);
         break;
       default:
         throw cms::Exception("InvalidInput") << "invalid jet type in VirtualJetProducer\n";
         break;
     };
   }
 }

References btagElecInJet_cfi::CaloJet, Exception, nanoDQM_cfi::GenJet, and iEvent.

◆ produce()

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

override

Definition at line 264 of file VirtualJetProducer.cc.

                                                                               {
   // If requested, set the fastjet random seed to a deterministic function
   // of the run/lumi/event.
   // NOTE!!! The fastjet random number sequence is a global singleton.
   // Thus, we have to create an object and get access to the global singleton
   // in order to change it.
   if (useDeterministicSeed_) {
     fastjet::GhostedAreaSpec gas;
     std::vector<int> seeds(2);
     unsigned int runNum_uint = static_cast<unsigned int>(iEvent.id().run());
     unsigned int evNum_uint = static_cast<unsigned int>(iEvent.id().event());
     seeds[0] = std::max(runNum_uint, minSeed_ + 3) + 3 * evNum_uint;
     seeds[1] = std::max(runNum_uint, minSeed_ + 5) + 5 * evNum_uint;
     gas.set_random_status(seeds);
   }
  
   LogDebug("VirtualJetProducer") << "Entered produce\n";
   //determine signal vertex2
   vertex_ = reco::Jet::Point(0, 0, 0);
   if ((makeCaloJet(jetTypeE) || makePFJet(jetTypeE)) && doPVCorrection_) {
     LogDebug("VirtualJetProducer") << "Adding PV info\n";
     edm::Handle<reco::VertexCollection> pvCollection;
     iEvent.getByToken(input_vertex_token_, pvCollection);
     if (!pvCollection->empty())
       vertex_ = pvCollection->begin()->position();
   }
  
   // Get Weights Collection
   if ((applyWeight_) && (!input_weights_token_.isUninitialized()))
     weights_ = iEvent.get(input_weights_token_);
  
   // For Pileup subtraction using offset correction:
   // set up geometry map
   if (doPUOffsetCorr_) {
     subtractor_->setupGeometryMap(iEvent, iSetup);
   }
  
   // clear data
   LogDebug("VirtualJetProducer") << "Clear data\n";
   fjInputs_.clear();
   fjJets_.clear();
   inputs_.clear();
  
   // get inputs and convert them to the fastjet format (fastjet::PeudoJet)
   edm::Handle<reco::CandidateView> inputsHandle;
  
   edm::Handle<std::vector<edm::FwdPtr<reco::PFCandidate>>> pfinputsHandleAsFwdPtr;
   edm::Handle<std::vector<edm::FwdPtr<pat::PackedCandidate>>> packedinputsHandleAsFwdPtr;
   edm::Handle<std::vector<edm::FwdPtr<reco::GenParticle>>> geninputsHandleAsFwdPtr;
   edm::Handle<std::vector<edm::FwdPtr<pat::PackedGenParticle>>> packedgeninputsHandleAsFwdPtr;
  
   bool isView = iEvent.getByToken(input_candidateview_token_, inputsHandle);
   if (isView) {
     if (inputsHandle->empty()) {
       output(iEvent, iSetup);
       return;
     }
     for (size_t i = 0; i < inputsHandle->size(); ++i) {
       inputs_.push_back(inputsHandle->ptrAt(i));
     }
   } else {
     bool isPF = iEvent.getByToken(input_candidatefwdptr_token_, pfinputsHandleAsFwdPtr);
     bool isPFFwdPtr = iEvent.getByToken(input_packedcandidatefwdptr_token_, packedinputsHandleAsFwdPtr);
     bool isGen = iEvent.getByToken(input_gencandidatefwdptr_token_, geninputsHandleAsFwdPtr);
     bool isGenFwdPtr = iEvent.getByToken(input_packedgencandidatefwdptr_token_, packedgeninputsHandleAsFwdPtr);
  
     if (isPF) {
       if (pfinputsHandleAsFwdPtr->empty()) {
         output(iEvent, iSetup);
         return;
       }
       for (size_t i = 0; i < pfinputsHandleAsFwdPtr->size(); ++i) {
         if ((*pfinputsHandleAsFwdPtr)[i].ptr().isAvailable()) {
           inputs_.push_back((*pfinputsHandleAsFwdPtr)[i].ptr());
         } else if ((*pfinputsHandleAsFwdPtr)[i].backPtr().isAvailable()) {
           inputs_.push_back((*pfinputsHandleAsFwdPtr)[i].backPtr());
         }
       }
     } else if (isPFFwdPtr) {
       if (packedinputsHandleAsFwdPtr->empty()) {
         output(iEvent, iSetup);
         return;
       }
       for (size_t i = 0; i < packedinputsHandleAsFwdPtr->size(); ++i) {
         if ((*packedinputsHandleAsFwdPtr)[i].ptr().isAvailable()) {
           inputs_.push_back((*packedinputsHandleAsFwdPtr)[i].ptr());
         } else if ((*packedinputsHandleAsFwdPtr)[i].backPtr().isAvailable()) {
           inputs_.push_back((*packedinputsHandleAsFwdPtr)[i].backPtr());
         }
       }
     } else if (isGen) {
       if (geninputsHandleAsFwdPtr->empty()) {
         output(iEvent, iSetup);
         return;
       }
       for (size_t i = 0; i < geninputsHandleAsFwdPtr->size(); ++i) {
         if ((*geninputsHandleAsFwdPtr)[i].ptr().isAvailable()) {
           inputs_.push_back((*geninputsHandleAsFwdPtr)[i].ptr());
         } else if ((*geninputsHandleAsFwdPtr)[i].backPtr().isAvailable()) {
           inputs_.push_back((*geninputsHandleAsFwdPtr)[i].backPtr());
         }
       }
     } else if (isGenFwdPtr) {
       if (geninputsHandleAsFwdPtr->empty()) {
         output(iEvent, iSetup);
         return;
       }
       for (size_t i = 0; i < packedgeninputsHandleAsFwdPtr->size(); ++i) {
         if ((*packedgeninputsHandleAsFwdPtr)[i].ptr().isAvailable()) {
           inputs_.push_back((*packedgeninputsHandleAsFwdPtr)[i].ptr());
         } else if ((*packedgeninputsHandleAsFwdPtr)[i].backPtr().isAvailable()) {
           inputs_.push_back((*packedgeninputsHandleAsFwdPtr)[i].backPtr());
         }
       }
     } else {
       throw cms::Exception("Invalid Jet Inputs")
           << "Did not specify appropriate inputs for VirtualJetProducer, Abort!\n";
     }
   }
   LogDebug("VirtualJetProducer") << "Got inputs\n";
  
   // Convert candidates to fastjet::PseudoJets.
   // Also correct to Primary Vertex. Will modify fjInputs_
   // and use inputs_
   fjInputs_.reserve(inputs_.size());
   inputTowers();
   LogDebug("VirtualJetProducer") << "Inputted towers\n";
  
   // For Pileup subtraction using offset correction:
   // Subtract pedestal.
   if (doPUOffsetCorr_) {
     subtractor_->setDefinition(fjJetDefinition_);
     subtractor_->reset(inputs_, fjInputs_, fjJets_);
     subtractor_->calculatePedestal(fjInputs_);
     subtractor_->subtractPedestal(fjInputs_);
     LogDebug("VirtualJetProducer") << "Subtracted pedestal\n";
   }
   // Run algorithm. Will modify fjJets_ and allocate fjClusterSeq_.
   // This will use fjInputs_
   runAlgorithm(iEvent, iSetup);
  
   // if ( doPUOffsetCorr_ ) {
   //    subtractor_->setAlgorithm(fjClusterSeq_);
   // }
  
   LogDebug("VirtualJetProducer") << "Ran algorithm\n";
   // For Pileup subtraction using offset correction:
   // Now we find jets and need to recalculate their energy,
   // mark towers participated in jet,
   // remove occupied towers from the list and recalculate mean and sigma
   // put the initial towers collection to the jet,
   // and subtract from initial towers in jet recalculated mean and sigma of towers
   if (doPUOffsetCorr_) {
     LogDebug("VirtualJetProducer") << "Do PUOffsetCorr\n";
     vector<fastjet::PseudoJet> orphanInput;
     subtractor_->calculateOrphanInput(orphanInput);
     subtractor_->calculatePedestal(orphanInput);
     subtractor_->offsetCorrectJets();
   }
   // Write the output jets.
   // This will (by default) call the member function template
   // "writeJets", but can be overridden.
   // this will use inputs_
   output(iEvent, iSetup);
   LogDebug("VirtualJetProducer") << "Wrote jets\n";
  
   // Clear the work vectors so that memory is free for other modules.
   // Use the trick of swapping with an empty vector so that the memory
   // is actually given back rather than silently kept.
   decltype(fjInputs_)().swap(fjInputs_);
   decltype(fjJets_)().swap(fjJets_);
   decltype(inputs_)().swap(inputs_);
  
   return;
 }

References edm::View< T >::empty(), Exception, mps_fire::i, iEvent, LogDebug, SiStripPI::max, convertSQLitetoXML_cfg::output, edm::View< T >::ptrAt(), anotherprimaryvertexanalyzer_cfi::pvCollection, InitialStep_cff::seeds, edm::View< T >::size(), and std::swap().

Referenced by cms::SubJetProducer::produce(), CMSInsideOutJetProducer::produce(), cms::CSJetProducer::produce(), SubjetFilterJetProducer::produce(), and FastjetJetProducer::produce().

◆ runAlgorithm()

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

protectedpure virtual

Implemented in cms::HTTTopJetProducer, cms::CATopJetProducer, FastjetJetProducer, SubjetFilterJetProducer, cms::CSJetProducer, CMSInsideOutJetProducer, cms::SubJetProducer, and cms::SubEventGenJetProducer.

◆ writeCompoundJets()

template<class T >

void VirtualJetProducer::writeCompoundJets	(	edm::Event &	iEvent,
		edm::EventSetup const &	iSetup
	)

protected

function template to write out the outputs

Definition at line 767 of file VirtualJetProducer.cc.

                                                                                         {
   if (verbosity_ >= 1) {
     std::cout << "<VirtualJetProducer::writeCompoundJets (moduleLabel = " << moduleLabel_ << ")>:" << std::endl;
   }
  
   // get a list of output jets
   auto jetCollection = std::make_unique<reco::BasicJetCollection>();
   // get a list of output subjets
   auto subjetCollection = std::make_unique<std::vector<T>>();
  
   // This will store the handle for the subjets after we write them
   edm::OrphanHandle<std::vector<T>> subjetHandleAfterPut;
   // this is the mapping of subjet to hard jet
   std::vector<std::vector<int>> indices;
   // this is the list of hardjet 4-momenta
   std::vector<math::XYZTLorentzVector> p4_hardJets;
   // this is the hardjet areas
   std::vector<double> area_hardJets;
  
   // Loop over the hard jets
   std::vector<fastjet::PseudoJet>::const_iterator it = fjJets_.begin(), iEnd = fjJets_.end(), iBegin = fjJets_.begin();
   indices.resize(fjJets_.size());
   for (; it != iEnd; ++it) {
     fastjet::PseudoJet const& localJet = *it;
     unsigned int jetIndex = it - iBegin;
     // Get the 4-vector for the hard jet
     p4_hardJets.push_back(math::XYZTLorentzVector(localJet.px(), localJet.py(), localJet.pz(), localJet.e()));
     double localJetArea = 0.0;
     if (doAreaFastjet_ && localJet.has_area()) {
       localJetArea = localJet.area();
     }
     area_hardJets.push_back(localJetArea);
  
     // create the subjet list
     std::vector<fastjet::PseudoJet> constituents;
     if (it->has_pieces()) {
       constituents = it->pieces();
     } else if (it->has_constituents()) {
       constituents = it->constituents();
     }
  
     std::vector<fastjet::PseudoJet>::const_iterator itSubJetBegin = constituents.begin(), itSubJet = itSubJetBegin,
                                                     itSubJetEnd = constituents.end();
     for (; itSubJet != itSubJetEnd; ++itSubJet) {
       fastjet::PseudoJet const& subjet = *itSubJet;
       if (verbosity_ >= 1) {
         std::cout << "subjet #" << (itSubJet - itSubJetBegin) << ": Pt = " << subjet.pt() << ", eta = " << subjet.eta()
                   << ", phi = " << subjet.phi() << ", mass = " << subjet.m()
                   << " (#constituents = " << subjet.constituents().size() << ")" << std::endl;
         std::vector<fastjet::PseudoJet> subjet_constituents = subjet.constituents();
         int idx_constituent = 0;
         for (std::vector<fastjet::PseudoJet>::const_iterator constituent = subjet_constituents.begin();
              constituent != subjet_constituents.end();
              ++constituent) {
           if (constituent->pt() < 1.e-3)
             continue;  // CV: skip ghosts
           std::cout << "  constituent #" << idx_constituent << ": Pt = " << constituent->pt()
                     << ", eta = " << constituent->eta() << ", phi = " << constituent->phi() << ","
                     << " mass = " << constituent->m() << std::endl;
           ++idx_constituent;
         }
       }
  
       if (verbosity_ >= 1) {
         std::cout << "subjet #" << (itSubJet - itSubJetBegin) << ": Pt = " << subjet.pt() << ", eta = " << subjet.eta()
                   << ", phi = " << subjet.phi() << ", mass = " << subjet.m()
                   << " (#constituents = " << subjet.constituents().size() << ")" << std::endl;
         std::vector<fastjet::PseudoJet> subjet_constituents = subjet.constituents();
         int idx_constituent = 0;
         for (std::vector<fastjet::PseudoJet>::const_iterator constituent = subjet_constituents.begin();
              constituent != subjet_constituents.end();
              ++constituent) {
           if (constituent->pt() < 1.e-3)
             continue;  // CV: skip ghosts
           std::cout << " constituent #" << idx_constituent << ": Pt = " << constituent->pt()
                     << ", eta = " << constituent->eta() << ", phi = " << constituent->phi() << ","
                     << " mass = " << constituent->m() << std::endl;
           ++idx_constituent;
         }
       }
  
       math::XYZTLorentzVector p4Subjet(subjet.px(), subjet.py(), subjet.pz(), subjet.e());
       reco::Particle::Point point(0, 0, 0);
  
       // This will hold ptr's to the subjets
       std::vector<reco::CandidatePtr> subjetConstituents;
  
       // Get the transient subjet constituents from fastjet
       std::vector<fastjet::PseudoJet> subjetFastjetConstituents = subjet.constituents();
       std::vector<reco::CandidatePtr> constituents = getConstituents(subjetFastjetConstituents);
  
       indices[jetIndex].push_back(subjetCollection->size());
  
       // Add the concrete subjet type to the subjet list to write to event record
       subjetCollection->push_back(*std::make_unique<T>());
       auto& jet = subjetCollection->back();
       if ((applyWeight_) && (makePFJet(jetTypeE)))
         reco::writeSpecific(dynamic_cast<reco::PFJet&>(jet), p4Subjet, point, constituents, iSetup, &weights_);
       else
         reco::writeSpecific(jet, p4Subjet, point, constituents, iSetup);
       jet.setIsWeighted(applyWeight_);
       double subjetArea = 0.0;
       if (doAreaFastjet_ && itSubJet->has_area()) {
         subjetArea = itSubJet->area();
       }
       jet.setJetArea(subjetArea);
     }
   }
  
   // put subjets into event record
   subjetHandleAfterPut = iEvent.put(std::move(subjetCollection), jetCollInstanceName_);
  
   // Now create the hard jets with ptr's to the subjets as constituents
   std::vector<math::XYZTLorentzVector>::const_iterator ip4 = p4_hardJets.begin(), ip4Begin = p4_hardJets.begin(),
                                                        ip4End = p4_hardJets.end();
  
   for (; ip4 != ip4End; ++ip4) {
     int p4_index = ip4 - ip4Begin;
     std::vector<int>& ind = indices[p4_index];
     std::vector<reco::CandidatePtr> i_hardJetConstituents;
     // Add the subjets to the hard jet
     for (std::vector<int>::const_iterator isub = ind.begin(); isub != ind.end(); ++isub) {
       reco::CandidatePtr candPtr(subjetHandleAfterPut, *isub, false);
       i_hardJetConstituents.push_back(candPtr);
     }
     reco::Particle::Point point(0, 0, 0);
     reco::BasicJet toput(*ip4, point, i_hardJetConstituents);
     toput.setJetArea(area_hardJets[ip4 - ip4Begin]);
     jetCollection->push_back(toput);
   }
  
   // put hard jets into event record
   // Store the Orphan handle for adding HTT information
   edm::OrphanHandle<reco::BasicJetCollection> oh = iEvent.put(std::move(jetCollection));
  
   if (fromHTTTopJetProducer_) {
     addHTTTopJetTagInfoCollection(iEvent, iSetup, oh);
   }
 }

References gather_cfg::cout, iEvent, bTagCombinedSVVariables_cff::indices, metsig::jet, jetfilter_cfi::jetCollection, eostools::move(), point, reco::Jet::setJetArea(), and reco::writeSpecific().

◆ writeJets()

template<typename T >

void VirtualJetProducer::writeJets	(	edm::Event &	iEvent,
		edm::EventSetup const &	iSetup
	)

protected

Definition at line 602 of file VirtualJetProducer.cc.

                                                                                 {
   // std::cout << "writeJets " << typeid(T).name()
   //          << (doRhoFastjet_ ? " doRhoFastjet " : "")
   //          << (doAreaFastjet_ ? " doAreaFastjet " : "")
   //          << (doAreaDiskApprox_ ? " doAreaDiskApprox " : "")
   //          << std::endl;
  
   if (doRhoFastjet_) {
     // declare jet collection without the two jets,
     // for unbiased background estimation.
     std::vector<fastjet::PseudoJet> fjexcluded_jets;
     fjexcluded_jets = fjJets_;
  
     if (fjexcluded_jets.size() > 2)
       fjexcluded_jets.resize(nExclude_);
  
     if (doFastJetNonUniform_) {
       auto rhos = std::make_unique<std::vector<double>>();
       auto sigmas = std::make_unique<std::vector<double>>();
       int nEta = puCenters_.size();
       rhos->reserve(nEta);
       sigmas->reserve(nEta);
       fastjet::ClusterSequenceAreaBase const* clusterSequenceWithArea =
           dynamic_cast<fastjet::ClusterSequenceAreaBase const*>(fjClusterSeq_.get());
  
       if (clusterSequenceWithArea == nullptr) {
         if (!fjJets_.empty()) {
           throw cms::Exception("LogicError") << "fjClusterSeq is not initialized while inputs are present\n ";
         }
       } else {
         for (int ie = 0; ie < nEta; ++ie) {
           double eta = puCenters_[ie];
           double etamin = eta - puWidth_;
           double etamax = eta + puWidth_;
           fastjet::RangeDefinition range_rho(etamin, etamax);
           fastjet::BackgroundEstimator bkgestim(*clusterSequenceWithArea, range_rho);
           bkgestim.set_excluded_jets(fjexcluded_jets);
           rhos->push_back(bkgestim.rho());
           sigmas->push_back(bkgestim.sigma());
         }
       }
       iEvent.put(std::move(rhos), "rhos");
       iEvent.put(std::move(sigmas), "sigmas");
     } else {
       auto rho = std::make_unique<double>(0.0);
       auto sigma = std::make_unique<double>(0.0);
       double mean_area = 0;
  
       fastjet::ClusterSequenceAreaBase const* clusterSequenceWithArea =
           dynamic_cast<fastjet::ClusterSequenceAreaBase const*>(fjClusterSeq_.get());
       if (clusterSequenceWithArea == nullptr) {
         if (!fjJets_.empty()) {
           throw cms::Exception("LogicError") << "fjClusterSeq is not initialized while inputs are present\n ";
         }
       } else {
         clusterSequenceWithArea->get_median_rho_and_sigma(*fjSelector_, false, *rho, *sigma, mean_area);
         if ((*rho < 0) || (edm::isNotFinite(*rho))) {
           edm::LogError("BadRho") << "rho value is " << *rho << " area:" << mean_area
                                   << " and n_empty_jets: " << clusterSequenceWithArea->n_empty_jets(*fjSelector_)
                                   << " with range " << fjSelector_->description() << ". Setting rho to rezo.";
           *rho = 0;
         }
       }
       iEvent.put(std::move(rho), "rho");
       iEvent.put(std::move(sigma), "sigma");
     }
   }  // doRhoFastjet_
  
   // produce output jet collection
  
   using namespace reco;
  
   // allocate fjJets_.size() Ts in vector
   auto jets = std::make_unique<std::vector<T>>(fjJets_.size());
  
   if (!fjJets_.empty()) {
     // Distance between jet centers and overlap area -- for disk-based area calculation
     using RIJ = std::pair<double, double>;
     std::vector<RIJ> rijStorage(fjJets_.size() * (fjJets_.size() / 2));
     RIJ* rij[fjJets_.size()];
     unsigned int k = 0;
     for (unsigned int ijet = 0; ijet < fjJets_.size(); ++ijet) {
       rij[ijet] = &rijStorage[k];
       k += ijet;
     }
  
     float etaJ[fjJets_.size()], phiJ[fjJets_.size()];
  
     auto orParam_ = 1. / rParam_;
     // fill jets
     for (unsigned int ijet = 0; ijet < fjJets_.size(); ++ijet) {
       auto& jet = (*jets)[ijet];
  
       // get the fastjet jet
       const fastjet::PseudoJet& fjJet = fjJets_[ijet];
       // get the constituents from fastjet
       std::vector<fastjet::PseudoJet> const& fjConstituents = fastjet::sorted_by_pt(fjJet.constituents());
       // convert them to CandidatePtr vector
       std::vector<CandidatePtr> const& constituents = getConstituents(fjConstituents);
  
       // write the specifics to the jet (simultaneously sets 4-vector, vertex).
       // These are overridden functions that will call the appropriate
       // specific allocator.
       if ((applyWeight_) && (makePFJet(jetTypeE)))
         writeSpecific(dynamic_cast<reco::PFJet&>(jet),
                       Particle::LorentzVector(fjJet.px(), fjJet.py(), fjJet.pz(), fjJet.E()),
                       vertex_,
                       constituents,
                       iSetup,
                       &weights_);
       else
         writeSpecific(
             jet, Particle::LorentzVector(fjJet.px(), fjJet.py(), fjJet.pz(), fjJet.E()), vertex_, constituents, iSetup);
       phiJ[ijet] = jet.phi();
       etaJ[ijet] = jet.eta();
       jet.setIsWeighted(applyWeight_);
     }
  
     // calcuate the jet area
     for (unsigned int ijet = 0; ijet < fjJets_.size(); ++ijet) {
       // calcuate the jet area
       double jetArea = 0.0;
       // get the fastjet jet
       const auto& fjJet = fjJets_[ijet];
       if (doAreaFastjet_ && fjJet.has_area()) {
         jetArea = fjJet.area();
       } else if (doAreaDiskApprox_) {
         // Here it is assumed that fjJets_ is in decreasing order of pT,
         // which should happen in FastjetJetProducer::runAlgorithm()
         jetArea = M_PI;
         RIJ* distance = rij[ijet];
         for (unsigned jJet = 0; jJet < ijet; ++jJet) {
           distance[jJet].first = std::sqrt(reco::deltaR2(etaJ[ijet], phiJ[ijet], etaJ[jJet], phiJ[jJet])) * orParam_;
           distance[jJet].second = reco::helper::VirtualJetProducerHelper::intersection(distance[jJet].first);
           jetArea -= distance[jJet].second;
           for (unsigned kJet = 0; kJet < jJet; ++kJet) {
             jetArea += reco::helper::VirtualJetProducerHelper::intersection(distance[jJet].first,
                                                                             distance[kJet].first,
                                                                             rij[jJet][kJet].first,
                                                                             distance[jJet].second,
                                                                             distance[kJet].second,
                                                                             rij[jJet][kJet].second);
           }  // end loop over harder jets
         }    // end loop over harder jets
         jetArea *= (rParam_ * rParam_);
       }
       auto& jet = (*jets)[ijet];
       jet.setJetArea(jetArea);
  
       if (doPUOffsetCorr_) {
         jet.setPileup(subtractor_->getPileUpEnergy(ijet));
       } else {
         jet.setPileup(0.0);
       }
  
       // std::cout << "area " << ijet << " " << jetArea << " " << Area<T>::get(jet) << std::endl;
       // std::cout << "JetVI " << ijet << ' '<< jet.pt() << " " << jet.et() << ' '<< jet.energy() << ' '<< jet.mass() << std::endl;
     }
   }
   // put the jets in the collection
   iEvent.put(std::move(jets), jetCollInstanceName_);
 }

References reco::deltaR2(), HLT_FULL_cff::distance, PVValHelper::eta, muonTiming_cfi::etamax, muonTiming_cfi::etamin, Exception, first, iEvent, reco::helper::VirtualJetProducerHelper::intersection(), edm::isNotFinite(), metsig::jet, singleTopDQM_cfi::jets, dqmdumpme::k, M_PI, eostools::move(), HLT_FULL_cff::nEta, fastjet::BackgroundEstimator::rho(), edm::second(), fastjet::BackgroundEstimator::set_excluded_jets(), fastjet::BackgroundEstimator::sigma(), mathSSE::sqrt(), and reco::writeSpecific().

◆ writeJetsWithConstituents()

template<class T >

void VirtualJetProducer::writeJetsWithConstituents	(	edm::Event &	iEvent,
		edm::EventSetup const &	iSetup
	)

protected

function template to write out the outputs

Definition at line 909 of file VirtualJetProducer.cc.

                                                                                                 {
   if (verbosity_ >= 1) {
     std::cout << "<VirtualJetProducer::writeJetsWithConstituents (moduleLabel = " << moduleLabel_ << ")>:" << std::endl;
   }
  
   // get a list of output jets  MV: make this compatible with template
   auto jetCollection = std::make_unique<reco::PFJetCollection>();
  
   // this is the mapping of jet to constituents
   std::vector<std::vector<int>> indices;
   // this is the list of jet 4-momenta
   std::vector<math::XYZTLorentzVector> p4_Jets;
   // this is the jet areas
   std::vector<double> area_Jets;
  
   // get a list of output constituents
   auto constituentCollection = std::make_unique<reco::PFCandidateCollection>();
  
   // This will store the handle for the constituents after we write them
   edm::OrphanHandle<reco::PFCandidateCollection> constituentHandleAfterPut;
  
   // Loop over the jets and extract constituents
   std::vector<fastjet::PseudoJet> constituentsSub;
   std::vector<fastjet::PseudoJet>::const_iterator it = fjJets_.begin(), iEnd = fjJets_.end(), iBegin = fjJets_.begin();
   indices.resize(fjJets_.size());
  
   for (; it != iEnd; ++it) {
     fastjet::PseudoJet const& localJet = *it;
     unsigned int jetIndex = it - iBegin;
     // Get the 4-vector for the hard jet
     p4_Jets.push_back(math::XYZTLorentzVector(localJet.px(), localJet.py(), localJet.pz(), localJet.e()));
     double localJetArea = 0.0;
     if (doAreaFastjet_ && localJet.has_area()) {
       localJetArea = localJet.area();
     }
     area_Jets.push_back(localJetArea);
  
     // create the constituent list
     std::vector<fastjet::PseudoJet> constituents, ghosts;
     if (it->has_pieces())
       constituents = it->pieces();
     else if (it->has_constituents())
       fastjet::SelectorIsPureGhost().sift(it->constituents(), ghosts, constituents);  //filter out ghosts
  
     //loop over constituents of jet (can be subjets or normal constituents)
     indices[jetIndex].reserve(constituents.size());
     constituentsSub.reserve(constituentsSub.size() + constituents.size());
     for (fastjet::PseudoJet const& constit : constituents) {
       indices[jetIndex].push_back(constituentsSub.size());
       constituentsSub.push_back(constit);
     }
   }
  
   //Loop over constituents and store in the event
   static const reco::PFCandidate dummySinceTranslateIsNotStatic;
   for (fastjet::PseudoJet const& constit : constituentsSub) {
     auto orig = inputs_[constit.user_index()];
     auto id = dummySinceTranslateIsNotStatic.translatePdgIdToType(orig->pdgId());
     reco::PFCandidate pCand(reco::PFCandidate(orig->charge(), orig->p4(), id));
     math::XYZTLorentzVector pVec;
     pVec.SetPxPyPzE(constit.px(), constit.py(), constit.pz(), constit.e());
     pCand.setP4(pVec);
     pCand.setSourceCandidatePtr(orig->sourceCandidatePtr(0));
     constituentCollection->push_back(pCand);
   }
   // put constituents into event record
   constituentHandleAfterPut = iEvent.put(std::move(constituentCollection), jetCollInstanceName_);
  
   // Now create the jets with ptr's to the constituents
   std::vector<math::XYZTLorentzVector>::const_iterator ip4 = p4_Jets.begin(), ip4Begin = p4_Jets.begin(),
                                                        ip4End = p4_Jets.end();
  
   for (; ip4 != ip4End; ++ip4) {
     int p4_index = ip4 - ip4Begin;
     std::vector<int>& ind = indices[p4_index];
     std::vector<reco::CandidatePtr> i_jetConstituents;
     // Add the constituents to the jet
     for (std::vector<int>::const_iterator iconst = ind.begin(); iconst != ind.end(); ++iconst) {
       reco::CandidatePtr candPtr(constituentHandleAfterPut, *iconst, false);
       i_jetConstituents.push_back(candPtr);
     }
     if (!i_jetConstituents.empty()) {  //only keep jets which have constituents after subtraction
       reco::Particle::Point point(0, 0, 0);
       reco::PFJet jet;
       reco::writeSpecific(jet, *ip4, point, i_jetConstituents, iSetup);
       jet.setJetArea(area_Jets[ip4 - ip4Begin]);
       jetCollection->emplace_back(jet);
     }
   }
  
   // put jets into event record
   iEvent.put(std::move(jetCollection));
 }

References gather_cfg::cout, triggerObjects_cff::id, iEvent, bTagCombinedSVVariables_cff::indices, metsig::jet, jetfilter_cfi::jetCollection, eostools::move(), point, reco::LeafCandidate::setP4(), reco::PFCandidate::setSourceCandidatePtr(), reco::PFCandidate::translatePdgIdToType(), and reco::writeSpecific().

Member Data Documentation

◆ activeAreaRepeats_

int VirtualJetProducer::activeAreaRepeats_

protected

Definition at line 175 of file VirtualJetProducer.h.

◆ anomalousTowerDef_

std::unique_ptr<AnomalousTower> VirtualJetProducer::anomalousTowerDef_

private

Definition at line 213 of file VirtualJetProducer.h.

◆ applyWeight_

bool VirtualJetProducer::applyWeight_

protected

Definition at line 209 of file VirtualJetProducer.h.

◆ doAreaDiskApprox_

bool VirtualJetProducer::doAreaDiskApprox_

protected

Definition at line 165 of file VirtualJetProducer.h.

◆ doAreaFastjet_

bool VirtualJetProducer::doAreaFastjet_

protected

Definition at line 163 of file VirtualJetProducer.h.

Referenced by cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

◆ doFastJetNonUniform_

bool VirtualJetProducer::doFastJetNonUniform_

protected

Definition at line 168 of file VirtualJetProducer.h.

◆ doPUOffsetCorr_

bool VirtualJetProducer::doPUOffsetCorr_

protected

Definition at line 179 of file VirtualJetProducer.h.

◆ doPVCorrection_

bool VirtualJetProducer::doPVCorrection_

protected

Definition at line 156 of file VirtualJetProducer.h.

◆ doRhoFastjet_

bool VirtualJetProducer::doRhoFastjet_

protected

Definition at line 167 of file VirtualJetProducer.h.

Referenced by cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

◆ fjActiveArea_

ActiveAreaSpecPtr VirtualJetProducer::fjActiveArea_

protected

Definition at line 188 of file VirtualJetProducer.h.

◆ fjAreaDefinition_

AreaDefinitionPtr VirtualJetProducer::fjAreaDefinition_

protected

Definition at line 189 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::FastjetJetProducer(), cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

◆ fjClusterSeq_

ClusterSequencePtr VirtualJetProducer::fjClusterSeq_

protected

Definition at line 185 of file VirtualJetProducer.h.

Referenced by cms::CSJetProducer::produce(), FastjetJetProducer::produce(), FastjetJetProducer::produceTrackJets(), cms::SubEventGenJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

◆ fjInputs_

std::vector<fastjet::PseudoJet> VirtualJetProducer::fjInputs_

protected

Definition at line 191 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::produce(), FastjetJetProducer::produceTrackJets(), cms::SubJetProducer::runAlgorithm(), cms::SubEventGenJetProducer::runAlgorithm(), CMSInsideOutJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), SubjetFilterJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

◆ fjJetDefinition_

JetDefPtr VirtualJetProducer::fjJetDefinition_

protected

Definition at line 186 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

◆ fjJets_

std::vector<fastjet::PseudoJet> VirtualJetProducer::fjJets_

protected

Definition at line 192 of file VirtualJetProducer.h.

Referenced by cms::HTTTopJetProducer::addHTTTopJetTagInfoCollection(), cms::SubEventGenJetProducer::produce(), FastjetJetProducer::produceTrackJets(), cms::SubEventGenJetProducer::runAlgorithm(), CMSInsideOutJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

◆ fjPlugin_

PluginPtr VirtualJetProducer::fjPlugin_

protected

Definition at line 187 of file VirtualJetProducer.h.

◆ fjSelector_

SelectorPtr VirtualJetProducer::fjSelector_

protected

Definition at line 190 of file VirtualJetProducer.h.

◆ fromHTTTopJetProducer_

bool VirtualJetProducer::fromHTTTopJetProducer_ = false

protected

Definition at line 208 of file VirtualJetProducer.h.

Referenced by cms::HTTTopJetProducer::HTTTopJetProducer().

◆ ghostArea_

double VirtualJetProducer::ghostArea_

protected

Definition at line 176 of file VirtualJetProducer.h.

◆ ghostEtaMax_

double VirtualJetProducer::ghostEtaMax_

protected

Definition at line 174 of file VirtualJetProducer.h.

◆ input_candidatefwdptr_token_

edm::EDGetTokenT<std::vector<edm::FwdPtr<reco::PFCandidate> > > VirtualJetProducer::input_candidatefwdptr_token_

private

Definition at line 217 of file VirtualJetProducer.h.

◆ input_candidateview_token_

edm::EDGetTokenT<reco::CandidateView> VirtualJetProducer::input_candidateview_token_

private

Definition at line 216 of file VirtualJetProducer.h.

◆ input_gencandidatefwdptr_token_

edm::EDGetTokenT<std::vector<edm::FwdPtr<reco::GenParticle> > > VirtualJetProducer::input_gencandidatefwdptr_token_

private

Definition at line 219 of file VirtualJetProducer.h.

◆ input_packedcandidatefwdptr_token_

edm::EDGetTokenT<std::vector<edm::FwdPtr<pat::PackedCandidate> > > VirtualJetProducer::input_packedcandidatefwdptr_token_

private

Definition at line 218 of file VirtualJetProducer.h.

◆ input_packedgencandidatefwdptr_token_

edm::EDGetTokenT<std::vector<edm::FwdPtr<pat::PackedGenParticle> > > VirtualJetProducer::input_packedgencandidatefwdptr_token_

private

Definition at line 220 of file VirtualJetProducer.h.

◆ input_vertex_token_

edm::EDGetTokenT<reco::VertexCollection> VirtualJetProducer::input_vertex_token_

protected

Definition at line 223 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::produceTrackJets().

◆ input_weights_token_

edm::EDGetTokenT<edm::ValueMap<float> > VirtualJetProducer::input_weights_token_

protected

Definition at line 224 of file VirtualJetProducer.h.

◆ inputEMin_

double VirtualJetProducer::inputEMin_

protected

Definition at line 154 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::inputTowers().

◆ inputEtMin_

double VirtualJetProducer::inputEtMin_

protected

Definition at line 153 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::inputTowers().

◆ inputs_

std::vector<edm::Ptr<reco::Candidate> > VirtualJetProducer::inputs_

protected

Definition at line 183 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::inputTowers(), cms::SubEventGenJetProducer::produce(), FastjetJetProducer::produceTrackJets(), SubjetFilterJetProducer::writeCompoundJets(), and cms::CompoundJetProducer::writeCompoundJets().

◆ jetAlgorithm_

std::string VirtualJetProducer::jetAlgorithm_

protected

Definition at line 151 of file VirtualJetProducer.h.

◆ jetCollInstanceName_

std::string VirtualJetProducer::jetCollInstanceName_

protected

Definition at line 199 of file VirtualJetProducer.h.

Referenced by cms::CompoundJetProducer::writeCompoundJets().

◆ jetPtMin_

double VirtualJetProducer::jetPtMin_

protected

Definition at line 155 of file VirtualJetProducer.h.

Referenced by cms::SubEventGenJetProducer::runAlgorithm(), cms::CSJetProducer::runAlgorithm(), and FastjetJetProducer::runAlgorithm().

◆ jetType_

std::string VirtualJetProducer::jetType_

protected

Definition at line 150 of file VirtualJetProducer.h.

◆ jetTypeE

JetType::Type VirtualJetProducer::jetTypeE

protected

Definition at line 54 of file VirtualJetProducer.h.

Referenced by SubjetFilterJetProducer::output(), cms::CompoundJetProducer::output(), and FastjetJetProducer::produce().

◆ maxInputs_

unsigned int VirtualJetProducer::maxInputs_

protected

Definition at line 160 of file VirtualJetProducer.h.

◆ minSeed_

unsigned int VirtualJetProducer::minSeed_

protected

Definition at line 205 of file VirtualJetProducer.h.

◆ moduleLabel_

std::string VirtualJetProducer::moduleLabel_

protected

Definition at line 147 of file VirtualJetProducer.h.

Referenced by ExternalGeneratorFilter.ExternalGeneratorFilter::insertInto(), Modules.SwitchProducer::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), Modules.ESProducer::nameInProcessDesc_(), and SubjetFilterJetProducer::SubjetFilterJetProducer().

◆ nExclude_

unsigned int VirtualJetProducer::nExclude_

protected

Definition at line 197 of file VirtualJetProducer.h.

◆ puCenters_

std::vector<double> VirtualJetProducer::puCenters_

protected

Definition at line 195 of file VirtualJetProducer.h.

◆ puSubtractorName_

std::string VirtualJetProducer::puSubtractorName_

protected

Definition at line 180 of file VirtualJetProducer.h.

◆ puWidth_

double VirtualJetProducer::puWidth_

protected

Definition at line 196 of file VirtualJetProducer.h.

◆ restrictInputs_

bool VirtualJetProducer::restrictInputs_

protected

Definition at line 159 of file VirtualJetProducer.h.

◆ rhoEtaMax_

double VirtualJetProducer::rhoEtaMax_

protected

Definition at line 173 of file VirtualJetProducer.h.

◆ rParam_

double VirtualJetProducer::rParam_

protected

Definition at line 152 of file VirtualJetProducer.h.

◆ src_

edm::InputTag VirtualJetProducer::src_

protected

Definition at line 148 of file VirtualJetProducer.h.

Referenced by cms::CATopJetProducer::CATopJetProducer(), and cms::SubEventGenJetProducer::SubEventGenJetProducer().

◆ srcPVs_

edm::InputTag VirtualJetProducer::srcPVs_

protected

Definition at line 149 of file VirtualJetProducer.h.

◆ subtractor_

std::shared_ptr<PileUpSubtractor> VirtualJetProducer::subtractor_

protected

Definition at line 202 of file VirtualJetProducer.h.

◆ useDeterministicSeed_

bool VirtualJetProducer::useDeterministicSeed_

protected

Definition at line 204 of file VirtualJetProducer.h.

◆ useExplicitGhosts_

bool VirtualJetProducer::useExplicitGhosts_

protected

Definition at line 164 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::FastjetJetProducer().

◆ verbosity_

int VirtualJetProducer::verbosity_

protected

Definition at line 207 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::runAlgorithm().

◆ vertex_

reco::Particle::Point VirtualJetProducer::vertex_

protected

Definition at line 184 of file VirtualJetProducer.h.

Referenced by FastjetJetProducer::produceTrackJets(), and cms::SubEventGenJetProducer::runAlgorithm().

◆ voronoiRfact_

double VirtualJetProducer::voronoiRfact_

protected

Definition at line 170 of file VirtualJetProducer.h.

Referenced by cms::CSJetProducer::runAlgorithm(), FastjetJetProducer::runAlgorithm(), cms::CATopJetProducer::runAlgorithm(), and cms::HTTTopJetProducer::runAlgorithm().

◆ weights_

edm::ValueMap<float> VirtualJetProducer::weights_

protected

Definition at line 210 of file VirtualJetProducer.h.

◆ writeCompound_

bool VirtualJetProducer::writeCompound_

protected

Definition at line 200 of file VirtualJetProducer.h.

◆ writeJetsWithConst_

bool VirtualJetProducer::writeJetsWithConst_

protected

Definition at line 201 of file VirtualJetProducer.h.

Classes

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Member Typedef Documentation

◆ ActiveAreaSpecPtr

◆ AreaDefinitionPtr

◆ ClusterSequencePtr

◆ JetDefPtr

◆ PluginPtr

◆ SelectorPtr

Constructor & Destructor Documentation

◆ VirtualJetProducer()

◆ ~VirtualJetProducer()

Member Function Documentation

◆ addHTTTopJetTagInfoCollection()

◆ copyConstituents()

◆ fillDescriptions()

◆ fillDescriptionsFromVirtualJetProducer()

◆ getConstituents()

◆ inputTowers()

◆ isAnomalousTower()

◆ jetType()

◆ makeBasicJet()

◆ makeCaloJet()

◆ makeGenJet()

◆ makePFClusterJet()

◆ makePFJet()

◆ makeProduces()

◆ makeTrackJet()

◆ offsetCorrectJets()

◆ output()

◆ produce()

◆ runAlgorithm()

◆ writeCompoundJets()

◆ writeJets()

◆ writeJetsWithConstituents()

Member Data Documentation

◆ activeAreaRepeats_

◆ anomalousTowerDef_

◆ applyWeight_

◆ doAreaDiskApprox_

◆ doAreaFastjet_

◆ doFastJetNonUniform_

◆ doPUOffsetCorr_

◆ doPVCorrection_

◆ doRhoFastjet_

◆ fjActiveArea_

◆ fjAreaDefinition_

◆ fjClusterSeq_

◆ fjInputs_

◆ fjJetDefinition_

◆ fjJets_

◆ fjPlugin_

◆ fjSelector_

◆ fromHTTTopJetProducer_

◆ ghostArea_

◆ ghostEtaMax_

◆ input_candidatefwdptr_token_

◆ input_candidateview_token_

◆ input_gencandidatefwdptr_token_

◆ input_packedcandidatefwdptr_token_

◆ input_packedgencandidatefwdptr_token_

◆ input_vertex_token_

◆ input_weights_token_

◆ inputEMin_

◆ inputEtMin_

◆ inputs_

◆ jetAlgorithm_

◆ jetCollInstanceName_

◆ jetPtMin_

◆ jetType_

◆ jetTypeE

◆ maxInputs_

◆ minSeed_

◆ moduleLabel_