#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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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::vector< int >	nSubParticles_

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

std::vector< reco::GenJet > *	subJets_

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_

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

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

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 ( const edm::ParameterSet & ps )

Definition at line 34 of file SubEventGenJetProducer.cc.

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

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

cms::SubEventGenJetProducer::~SubEventGenJetProducer ( )

inlineoverride

Definition at line 24 of file SubEventGenJetProducer.h.

References produce(), and runAlgorithm().

24 {}

Member Function Documentation

void SubEventGenJetProducer::inputTowers ( )

overrideprotectedvirtual

Reimplemented from VirtualJetProducer.

Definition at line 47 of file SubEventGenJetProducer.cc.

References reco::GenParticle::collisionId(), reco::Candidate::energy(), reco::Candidate::et(), edm::Ptr< T >::get(), hydroTag_, mps_fire::i, input, VirtualJetProducer::inputEMin_, VirtualJetProducer::inputEtMin_, VirtualJetProducer::inputs_, createfilelist::int, VirtualJetProducer::isAnomalousTower(), edm::isNotFinite(), LogDebug, reco::Candidate::pt(), reco::Candidate::px(), reco::Candidate::py(), reco::Candidate::pz(), and subInputs_.

Referenced by produce().

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

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

override

Definition at line 83 of file SubEventGenJetProducer.cc.

References VirtualJetProducer::fjInputs_, VirtualJetProducer::fjJets_, edm::Event::getByToken(), hydroTag_, mps_fire::i, ignoreHydro_, input_cand_token_, VirtualJetProducer::inputs_, inputTowers(), fwrapper::jets, LogDebug, eostools::move(), nSubParticles_, edm::View< T >::ptrAt(), edm::Event::put(), runAlgorithm(), edm::View< T >::size(), subInputs_, and subJets_.

Referenced by ~SubEventGenJetProducer().

                                                                                 {
    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_
 
 
    auto jets = std::make_unique<std::vector<GenJet>>();
    subJets_ = jets.get();
 
    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;
 }

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

overridevirtual

Implements VirtualJetProducer.

Definition at line 129 of file SubEventGenJetProducer.cc.

References DEFINE_FWK_MODULE, VirtualJetProducer::fjClusterSeq_, VirtualJetProducer::fjInputs_, VirtualJetProducer::fjJetDefinition_, VirtualJetProducer::fjJets_, VirtualJetProducer::getConstituents(), metsig::jet, VirtualJetProducer::jetPtMin_, muons2muons_cfi::pu, reco::Jet::setJetArea(), reco::Jet::setPileup(), subJets_, VirtualJetProducer::vertex_, and reco::writeSpecific().

Referenced by produce(), and ~SubEventGenJetProducer().

 {
    // run algorithm
    fjJets_.clear();
 
    fjClusterSeq_ = ClusterSequencePtr( new 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);
     
     subJets_->push_back(jet);
    }   
 }