#include <HLTrigger/JetMET/plugins/HLTScoutingPFProducer.cc>

Inheritance diagram for HLTScoutingPFProducer:

Public Member Functions
	HLTScoutingPFProducer (const edm::ParameterSet &)

	~HLTScoutingPFProducer () override

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::StreamID sid, edm::Event &iEvent, edm::EventSetup const &setup) const final

Private Attributes
const bool	doCandidates

const bool	doJetTags

const bool	doMet

const edm::EDGetTokenT< reco::PFMETCollection >	metCollection_

const edm::EDGetTokenT< reco::PFCandidateCollection >	pfCandidateCollection_

const double	pfCandidatePtCut

const edm::EDGetTokenT< reco::PFJetCollection >	pfJetCollection_

const double	pfJetEtaCut

const double	pfJetPtCut

const edm::EDGetTokenT< reco::JetTagCollection >	pfJetTagCollection_

const edm::EDGetTokenT< double >	rho_

const edm::EDGetTokenT< reco::VertexCollection >	vertexCollection_

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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)

Detailed Description

Description: Producer for ScoutingPFJets from reco::PFJet objects, ScoutingVertexs from reco::Vertexs and ScoutingParticles from reco::PFCandidates

Definition at line 42 of file HLTScoutingPFProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 71 of file HLTScoutingPFProducer.cc.

References ~HLTScoutingPFProducer().

                                                                           :
     pfJetCollection_(consumes<reco::PFJetCollection>(iConfig.getParameter<edm::InputTag>("pfJetCollection"))),
     pfJetTagCollection_(consumes<reco::JetTagCollection>(iConfig.getParameter<edm::InputTag>("pfJetTagCollection"))),
     pfCandidateCollection_(consumes<reco::PFCandidateCollection>(iConfig.getParameter<edm::InputTag>("pfCandidateCollection"))),
     vertexCollection_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("vertexCollection"))),
     metCollection_(consumes<reco::PFMETCollection>(iConfig.getParameter<edm::InputTag>("metCollection"))),
     rho_(consumes<double>(iConfig.getParameter<edm::InputTag>("rho"))),
     pfJetPtCut(iConfig.getParameter<double>("pfJetPtCut")),
     pfJetEtaCut(iConfig.getParameter<double>("pfJetEtaCut")),
     pfCandidatePtCut(iConfig.getParameter<double>("pfCandidatePtCut")),
     doJetTags(iConfig.getParameter<bool>("doJetTags")),
     doCandidates(iConfig.getParameter<bool>("doCandidates")),
     doMet(iConfig.getParameter<bool>("doMet"))
 {
     //register products
     produces<ScoutingPFJetCollection>();
     produces<ScoutingParticleCollection>();
     produces<double>("rho");
     produces<double>("pfMetPt");
     produces<double>("pfMetPhi");
 }

HLTScoutingPFProducer::~HLTScoutingPFProducer ( )

overridedefault

Referenced by HLTScoutingPFProducer().

Member Function Documentation

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

static

Definition at line 224 of file HLTScoutingPFProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and DEFINE_FWK_MODULE.

                                                                                        {
     edm::ParameterSetDescription desc;
     desc.add<edm::InputTag>("pfJetCollection",edm::InputTag("hltAK4PFJets"));
     desc.add<edm::InputTag>("pfJetTagCollection",edm::InputTag("hltCombinedSecondaryVertexBJetTagsPF"));
     desc.add<edm::InputTag>("pfCandidateCollection", edm::InputTag("hltParticleFlow"));
     desc.add<edm::InputTag>("vertexCollection", edm::InputTag("hltPixelVertices"));
     desc.add<edm::InputTag>("metCollection", edm::InputTag("hltPFMETProducer"));
     desc.add<edm::InputTag>("rho", edm::InputTag("hltFixedGridRhoFastjetAll"));
     desc.add<double>("pfJetPtCut", 20.0);
     desc.add<double>("pfJetEtaCut", 3.0);
     desc.add<double>("pfCandidatePtCut", 0.6);
     desc.add<bool>("doJetTags", true);
     desc.add<bool>("doCandidates", true);
     desc.add<bool>("doMet", true);
     descriptions.add("hltScoutingPFProducer", desc);
 }

void HLTScoutingPFProducer::produce	(	edm::StreamID	sid,
		edm::Event &	iEvent,
		edm::EventSetup const &	setup
	)		const

finalprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 96 of file HLTScoutingPFProducer.cc.

References reco::deltaR2(), doCandidates, doJetTags, doMet, edm::Event::getByToken(), metsig::jet, MT2Analyzer::metCollection, metCollection_, eostools::move(), pfCandidateCollection_, pfCandidatePtCut, pfJetCollection_, pfJetEtaCut, pfJetPtCut, pfJetTagCollection_, funct::pow(), edm::Event::put(), rho, rho_, GlobalPosition_Frontier_DevDB_cff::tag, particleFlowSuperClusterECAL_cfi::vertexCollection, vertexCollection_, and badGlobalMuonTaggersAOD_cff::vtx.

 {
     using namespace edm;
 
     //get vertices
      Handle<reco::VertexCollection> vertexCollection;
      std::unique_ptr<ScoutingVertexCollection> outVertices(new ScoutingVertexCollection());
      if(iEvent.getByToken(vertexCollection_, vertexCollection)){
         for(auto &vtx : *vertexCollection){
             outVertices->emplace_back(
                           vtx.x(), vtx.y(), vtx.z(), vtx.zError(), vtx.xError(), vtx.yError(), vtx.tracksSize(), vtx.chi2(), vtx.ndof(), vtx.isValid()
                     );
         }
       }
 
     //get rho
     Handle<double> rho;
     std::unique_ptr<double> outRho(new double(-999));
     if(iEvent.getByToken(rho_, rho)){
         outRho.reset(new double(*rho));
     }
 
     //get MET
     Handle<reco::PFMETCollection> metCollection;
     std::unique_ptr<double> outMetPt(new double(-999));
     std::unique_ptr<double> outMetPhi(new double(-999));
     if(doMet && iEvent.getByToken(metCollection_, metCollection)){
         outMetPt.reset(new double(metCollection->front().pt()));
         outMetPhi.reset(new double(metCollection->front().phi()));
     }
 
     //get PF candidates
     Handle<reco::PFCandidateCollection> pfCandidateCollection;
     std::unique_ptr<ScoutingParticleCollection> outPFCandidates(new ScoutingParticleCollection());
     if(doCandidates && iEvent.getByToken(pfCandidateCollection_, pfCandidateCollection)){
         for(auto &cand : *pfCandidateCollection){
             if(cand.pt() > pfCandidatePtCut){
                 int vertex_index = -1;
                 int index_counter = 0;
                 double dr2 = 0.0001;
                 for (auto &vtx: *outVertices) {
                     double tmp_dr2 = pow(vtx.x() - cand.vx(), 2) + pow(vtx.y() - cand.vy(), 2)
                         + pow(vtx.z() - cand.vz(), 2);
                     if (tmp_dr2 < dr2) {
                         dr2 = tmp_dr2;
                         vertex_index = index_counter;
                     }
                     if (dr2 == 0.0)
                         break;
                     ++index_counter;
                 }
                 outPFCandidates->emplace_back(
                         cand.pt(), cand.eta(), cand.phi(), cand.mass(), cand.pdgId(), vertex_index
                         );
             }
         }
     }
 
     //get PF jets
     Handle<reco::PFJetCollection> pfJetCollection;
     std::unique_ptr<ScoutingPFJetCollection> outPFJets(new ScoutingPFJetCollection());
     if(iEvent.getByToken(pfJetCollection_, pfJetCollection)){
         //get PF jet tags
         Handle<reco::JetTagCollection> pfJetTagCollection;
         bool haveJetTags = false;
         if(doJetTags && iEvent.getByToken(pfJetTagCollection_, pfJetTagCollection)){
             haveJetTags = true;
         }
 
         for(auto &jet : *pfJetCollection){
             if(jet.pt() < pfJetPtCut || fabs(jet.eta()) > pfJetEtaCut) continue;
             //find the jet tag corresponding to the jet
             float tagValue = -20;
             float minDR2 = 0.01;
             if(haveJetTags){
                 for(auto &tag : *pfJetTagCollection){
                     float dR2 = reco::deltaR2(jet, *(tag.first));
                     if(dR2 < minDR2){
                         minDR2 = dR2;
                         tagValue = tag.second;
                     }
                 }
             }
             //get the PF constituents of the jet
             std::vector<int> candIndices;
             if(doCandidates){
                 for(auto &cand : jet.getPFConstituents()){
                     if(cand->pt() > pfCandidatePtCut){
                         //search for the candidate in the collection
                         float minDR2 = 0.0001;
                         int matchIndex = -1;
                         int outIndex = 0;
                         for(auto &outCand : *outPFCandidates){
                             float dR2 = pow(cand->eta() - outCand.eta(), 2) + pow(cand->phi() - outCand.phi(), 2);
                             if(dR2 < minDR2){
                                 minDR2 = dR2;
                                 matchIndex = outIndex;
                             }
                             if(minDR2 == 0){
                                 break;
                             }
                             outIndex++;
                         }
                         candIndices.push_back(matchIndex);
                     }
                 }
             }
             outPFJets->emplace_back(
                     jet.pt(), jet.eta(), jet.phi(), jet.mass(), jet.jetArea(),
                     jet.chargedHadronEnergy(), jet.neutralHadronEnergy(), jet.photonEnergy(),
                     jet.electronEnergy(), jet.muonEnergy(), jet.HFHadronEnergy(), jet.HFEMEnergy(),
                     jet.chargedHadronMultiplicity(), jet.neutralHadronMultiplicity(), jet.photonMultiplicity(),
                     jet.electronMultiplicity(), jet.muonMultiplicity(), 
                     jet.HFHadronMultiplicity(), jet.HFEMMultiplicity(), 
                     jet.hoEnergy(), tagValue, 0.0, candIndices
                     );
         }
     }
 
     //put output
     iEvent.put(std::move(outPFCandidates));
     iEvent.put(std::move(outPFJets));
     iEvent.put(std::move(outRho), "rho");
     iEvent.put(std::move(outMetPt), "pfMetPt");
     iEvent.put(std::move(outMetPhi), "pfMetPhi");
 }

Member Data Documentation

const bool HLTScoutingPFProducer::doCandidates

private

Definition at line 64 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const bool HLTScoutingPFProducer::doJetTags

private

Definition at line 63 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const bool HLTScoutingPFProducer::doMet

private

Definition at line 65 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<reco::PFMETCollection> HLTScoutingPFProducer::metCollection_

private

Definition at line 56 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<reco::PFCandidateCollection> HLTScoutingPFProducer::pfCandidateCollection_

private

Definition at line 54 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const double HLTScoutingPFProducer::pfCandidatePtCut

private

Definition at line 61 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<reco::PFJetCollection> HLTScoutingPFProducer::pfJetCollection_

private

Definition at line 52 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const double HLTScoutingPFProducer::pfJetEtaCut

private

Definition at line 60 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const double HLTScoutingPFProducer::pfJetPtCut

private

Definition at line 59 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<reco::JetTagCollection> HLTScoutingPFProducer::pfJetTagCollection_

private

Definition at line 53 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<double> HLTScoutingPFProducer::rho_

private

Definition at line 57 of file HLTScoutingPFProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<reco::VertexCollection> HLTScoutingPFProducer::vertexCollection_

private

Definition at line 55 of file HLTScoutingPFProducer.cc.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation