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

Inheritance diagram for HLTScoutingPFProducer:

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

	~HLTScoutingPFProducer ()

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
	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
std::vector< ConsumesInfo >	consumesInfo () const

	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

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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
virtual void	produce (edm::StreamID sid, edm::Event &iEvent, edm::EventSetup const &setup) const overridefinal

Private Attributes
const bool	doCandidates

const bool	doJetTags

const bool	doMet

const edm::EDGetTokenT < reco::METCollection >	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

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)

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.

                                                                           :
     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::METCollection>(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<ScoutingVertexCollection>();
     produces<double>("rho");
     produces<double>("pfMetPt");
     produces<double>("pfMetPhi");
 }

HLTScoutingPFProducer::~HLTScoutingPFProducer ( )

Definition at line 94 of file HLTScoutingPFProducer.cc.

95 { }

Member Function Documentation

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

static

Definition at line 245 of file HLTScoutingPFProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLT_25ns14e33_v1_cff::InputTag.

                                                                                        {
     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

finaloverrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 98 of file HLTScoutingPFProducer.cc.

References reco::deltaR2(), doCandidates, doJetTags, doMet, edm::Event::getByToken(), metsig::jet, HLT_FULL_cff::metCollection, metCollection_, HLT_FULL_cff::pfCandidateCollection, pfCandidateCollection_, pfCandidatePtCut, HLT_FULL_cff::pfJetCollection, pfJetCollection_, pfJetEtaCut, pfJetPtCut, HLT_FULL_cff::pfJetTagCollection, pfJetTagCollection_, funct::pow(), edm::Event::put(), rho, rho_, GlobalPosition_Frontier_DevDB_cff::tag, GoodVertex_cfg::vertexCollection, and vertexCollection_.

 {
     using namespace edm;
 
     //get PF jets
     Handle<reco::PFJetCollection> pfJetCollection;
     if(!iEvent.getByToken(pfJetCollection_, pfJetCollection)){
         edm::LogError ("HLTScoutingPFProducer") << "invalid collection: pfJetCollection" << "\n";
         return;
     }
     //get PF jet tags
     Handle<reco::JetTagCollection> pfJetTagCollection;
     if(doJetTags && !iEvent.getByToken(pfJetTagCollection_, pfJetTagCollection)){
         edm::LogError ("HLTScoutingPFProducer") << "invalid collection: pfJetTagCollection" << "\n";
         return;
     }
     //get PF candidates
     Handle<reco::PFCandidateCollection> pfCandidateCollection;
     if(doCandidates && !iEvent.getByToken(pfCandidateCollection_, pfCandidateCollection)){
         edm::LogError ("HLTScoutingPFProducer") << "invalid collection: pfCandidateCollection" << "\n";
         return;
     }
     //get vertices
     Handle<reco::VertexCollection> vertexCollection;
     if(!iEvent.getByToken(vertexCollection_, vertexCollection)){
         edm::LogError ("HLTScoutingPFProducer") << "invalid collection: vertexCollection" << "\n";
         return;
     }
     //get rho
     Handle<double> rho;
     if(!iEvent.getByToken(rho_, rho)){
         edm::LogError ("HLTScoutingPFProducer") << "invalid collection: rho" << "\n";
         return;
     }
     std::auto_ptr<double> outRho(new double(*rho));
     //get MET
     Handle<reco::METCollection> metCollection;
     if(doMet && !iEvent.getByToken(metCollection_, metCollection)){
         edm::LogError ("HLTScoutingPFProducer") << "invalid collection: metCollection" << "\n";
         return;
     }
 
     //produce vertices
     std::auto_ptr<ScoutingVertexCollection> outVertices(new ScoutingVertexCollection());
     for(auto &vtx : *vertexCollection){
         outVertices->emplace_back(
                 vtx.x(), vtx.y(), vtx.z(), vtx.zError()
                 );
     }
 
     //produce PF candidates
     std::auto_ptr<ScoutingParticleCollection> outPFCandidates(new ScoutingParticleCollection());
     if(doCandidates){
         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
                         );
             }
         }
     }
 
     //produce PF jets
     std::auto_ptr<ScoutingPFJetCollection> outPFJets(new ScoutingPFJetCollection());
     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(doJetTags){
             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
                     );
     }
 
     double metPt = -999;
     double metPhi = -999;
     if(doMet){
         metPt = metCollection->front().pt();
         metPhi = metCollection->front().phi();
     }
     std::auto_ptr<double> outMetPt(new double(metPt));
     std::auto_ptr<double> outMetPhi(new double(metPhi));
 
 
     //put output
     iEvent.put(outVertices);
     iEvent.put(outPFCandidates);
     iEvent.put(outPFJets);
     iEvent.put(outRho, "rho");
     iEvent.put(outMetPt, "pfMetPt");
     iEvent.put(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::METCollection> 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