#include <SubEventGenJetProducer.h>

Inheritance diagram for cms::SubEventGenJetProducer:

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

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

	SubEventGenJetProducer (const edm::ParameterSet &ps)

virtual	~SubEventGenJetProducer ()

Public Member Functions inherited from VirtualJetProducer
std::string	jetType () const

	VirtualJetProducer (const edm::ParameterSet &iConfig)

virtual	~VirtualJetProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
virtual void	inputTowers ()

Protected Member Functions inherited from VirtualJetProducer
virtual void	copyConstituents (const std::vector< fastjet::PseudoJet > &fjConstituents, reco::Jet *jet)

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

virtual 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)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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

RangeDefPtr	fjRangeDef_

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

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::RangeDefinition >	RangeDefPtr

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

virtual cms::SubEventGenJetProducer::~SubEventGenJetProducer ( )

inlinevirtual

Definition at line 24 of file SubEventGenJetProducer.h.

24 {}

Member Function Documentation

void SubEventGenJetProducer::inputTowers ( )

protectedvirtual

Reimplemented from VirtualJetProducer.

Definition at line 47 of file SubEventGenJetProducer.cc.

References reco::GenParticle::collisionId(), edm::Ptr< T >::get(), hydroTag_, i, input, VirtualJetProducer::inputEMin_, VirtualJetProducer::inputEtMin_, VirtualJetProducer::inputs_, VirtualJetProducer::isAnomalousTower(), edm::isNotFinite(), LogDebug, 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
	)

virtual

Reimplemented from VirtualJetProducer.

Definition at line 83 of file SubEventGenJetProducer.cc.

References VirtualJetProducer::fjInputs_, VirtualJetProducer::fjJets_, edm::Event::getByToken(), hydroTag_, i, ignoreHydro_, input_cand_token_, VirtualJetProducer::inputs_, inputTowers(), fwrapper::jets, LogDebug, nSubParticles_, edm::Event::put(), runAlgorithm(), subInputs_, and subJets_.

                                                                                 {
    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_
 
 
    std::auto_ptr<std::vector<GenJet> > jets(new 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(jets);  
    return;
 }

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

virtual

Implements VirtualJetProducer.

Definition at line 129 of file SubEventGenJetProducer.cc.

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

Referenced by produce().

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