#include <IPProducer.h>

Public Member Functions
	FromJetAndCands (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iC, const std::string &jets="jets")

std::vector< reco::JetTagInfo >	makeBaseVector (const edm::Event &iEvent)

const std::vector < reco::CandidatePtr > &	tracks (const reco::JetTagInfo &it)

Public Attributes
bool	explicitJTA

std::vector< std::vector < reco::CandidatePtr > >	m_map

double	maxDeltaR

edm::EDGetTokenT< edm::View < reco::Candidate > >	token_cands

edm::EDGetTokenT< edm::View < reco::Jet > >	token_jets

Detailed Description

Definition at line 82 of file IPProducer.h.

Constructor & Destructor Documentation

IPProducerHelpers::FromJetAndCands::FromJetAndCands	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC,
		const std::string &	jets = `"jets"`
	)

inline

Definition at line 84 of file IPProducer.h.

         : token_jets(iC.consumes<edm::View<reco::Jet> >(iConfig.getParameter<edm::InputTag>(jets))),
           token_cands(iC.consumes<edm::View<reco::Candidate> >(iConfig.getParameter<edm::InputTag>("candidates"))),
           maxDeltaR(iConfig.getParameter<double>("maxDeltaR")),
           explicitJTA(iConfig.existsAs<bool>("explicitJTA") ? iConfig.getParameter<bool>("explicitJTA") : false) {}

Member Function Documentation

std::vector<reco::JetTagInfo> IPProducerHelpers::FromJetAndCands::makeBaseVector ( const edm::Event & iEvent )

inline

Definition at line 91 of file IPProducer.h.

References HLT_FULL_cff::cands, reco::deltaR2(), explicitJTA, edm::Event::getByToken(), mps_fire::i, dqmiolumiharvest::j, fwrapper::jets, m_map, maxDeltaR, token_cands, and token_jets.

                                                                      {
       edm::Handle<edm::View<reco::Jet> > jets;
       iEvent.getByToken(token_jets, jets);
       std::vector<reco::JetTagInfo> bases;
 
       edm::Handle<edm::View<reco::Candidate> > cands;
       iEvent.getByToken(token_cands, cands);
       m_map.clear();
       m_map.resize(jets->size());
       double maxDeltaR2 = maxDeltaR * maxDeltaR;
       size_t i = 0;
       for (edm::View<reco::Jet>::const_iterator it = jets->begin(); it != jets->end(); it++, i++) {
         edm::RefToBase<reco::Jet> jRef(jets, i);
         bases.push_back(reco::JetTagInfo(jRef));
         if (explicitJTA) {
           for (size_t j = 0; j < it->numberOfDaughters(); ++j) {
             if (it->daughterPtr(j)->bestTrack() != nullptr && it->daughterPtr(j)->charge() != 0) {
               m_map[i].push_back(it->daughterPtr(j));
             }
           }
         } else {
           for (size_t j = 0; j < cands->size(); ++j) {
             if ((*cands)[j].bestTrack() != nullptr && (*cands)[j].charge() != 0 && (*cands)[j].pt() > 0 &&
                 Geom::deltaR2((*cands)[j], (*jets)[i]) < maxDeltaR2) {
               m_map[i].push_back(cands->ptrAt(j));
             }
           }
         }
       }
       return bases;
     }