Inheritance diagram for RecoTauJetRegionProducer:

Public Types
typedef edm::Association < reco::PFJetCollection >	PFJetMatchMap

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

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

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

	RecoTauJetRegionProducer (const edm::ParameterSet &pset)

	~RecoTauJetRegionProducer ()

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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Attributes
float	deltaR2_

edm::InputTag	inputJets_

edm::InputTag	pfSrc_

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 26 of file RecoTauJetRegionProducer.cc.

Member Typedef Documentation

typedef edm::Association<reco::PFJetCollection> RecoTauJetRegionProducer::PFJetMatchMap

Definition at line 28 of file RecoTauJetRegionProducer.cc.

Constructor & Destructor Documentation

RecoTauJetRegionProducer::RecoTauJetRegionProducer ( const edm::ParameterSet & pset )

explicit

Definition at line 38 of file RecoTauJetRegionProducer.cc.

References deltaR2_, edm::ParameterSet::getParameter(), inputJets_, and pfSrc_.

                                  {
   deltaR2_ = pset.getParameter<double>("deltaR"); deltaR2_*=deltaR2_;
   inputJets_ = pset.getParameter<edm::InputTag>("src");
   pfSrc_ = pset.getParameter<edm::InputTag>("pfSrc");
   produces<reco::PFJetCollection>("jets");
   produces<PFJetMatchMap>();
 }

RecoTauJetRegionProducer::~RecoTauJetRegionProducer ( )

inline

Definition at line 30 of file RecoTauJetRegionProducer.cc.

30 {}

Member Function Documentation

void RecoTauJetRegionProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

virtual

Implements edm::EDProducer.

Definition at line 47 of file RecoTauJetRegionProducer.cc.

References reco::CompositePtrCandidate::addDaughter(), reco::CompositePtrCandidate::clearDaughters(), reco::deltaR2(), deltaR2_, alignCSCRings::e, edm::helper::Filler< Map >::fill(), edm::Event::get(), edm::Event::getByLabel(), edm::RefVector< C, T, F >::id(), inputJets_, edm::helper::Filler< Map >::insert(), fwrapper::jets, edm::Ref< C, T, F >::key(), pfSrc_, edm::Event::put(), and edm::RefVector< C, T, F >::size().

                              {
 
   edm::Handle<reco::PFCandidateCollection> pfCandsHandle;
   evt.getByLabel(pfSrc_, pfCandsHandle);
 
   // Build Ptrs for all the PFCandidates
   typedef edm::Ptr<reco::PFCandidate> PFCandPtr;
   std::vector<PFCandPtr> pfCands;
   pfCands.reserve(pfCandsHandle->size());
   for (size_t icand = 0; icand < pfCandsHandle->size(); ++icand) {
     pfCands.push_back(PFCandPtr(pfCandsHandle, icand));
   }
 
   // Get the jets
   edm::Handle<reco::CandidateView> jetView;
   evt.getByLabel(inputJets_, jetView);
   // Convert to a vector of PFJetRefs
   reco::PFJetRefVector jets =
       reco::tau::castView<reco::PFJetRefVector>(jetView);
   size_t nJets = jets.size();
 
   // Get the original product, so we can match against it - otherwise the
   // indices don't match up.
   edm::ProductID originalId = jets.id();
   edm::Handle<reco::PFJetCollection> originalJets;
   size_t nOriginalJets = 0;
   // We have to make sure that we have some selected jets, otherwise we don't
   // actually have a valid product ID to the original jets.
   if (nJets) {
     try {
       evt.get(originalId, originalJets);
     } catch(const cms::Exception &e) {
       edm::LogError("MissingOriginalCollection")
         << "Can't get the original jets that made: " << inputJets_
         << " that have product ID: " << originalId
         << " from the event!!";
       throw e;
     }
     nOriginalJets = originalJets->size();
   }
 
   std::auto_ptr<reco::PFJetCollection> newJets(new reco::PFJetCollection);
 
   // Keep track of the indices of the current jet and the old (original) jet
   // -1 indicats no match.
   std::vector<int> matchInfo(nOriginalJets, -1);
   newJets->reserve(nJets);
   for (size_t ijet = 0; ijet < nJets; ++ijet) {
     // Get a ref to jet
     reco::PFJetRef jetRef = jets[ijet];
     // Make an initial copy.
     newJets->emplace_back(*jetRef);
     reco::PFJet & newJet = newJets->back();
     // Clear out all the constituents
     newJet.clearDaughters();
     // Loop over all the PFCands
     for ( auto cand :  pfCands )
       if ( reco::deltaR2(*jetRef,*cand)<deltaR2_ ) newJet.addDaughter(cand);
     // Match the index of the jet we just made to the index into the original
     // collection.
     matchInfo[jetRef.key()] = ijet;
   }
 
   // Put our new jets into the event
   edm::OrphanHandle<reco::PFJetCollection> newJetsInEvent =
     evt.put(newJets, "jets");
 
   // Create a matching between original jets -> extra collection
   std::auto_ptr<PFJetMatchMap> matching(new PFJetMatchMap(newJetsInEvent));
   if (nJets) {
     PFJetMatchMap::Filler filler(*matching);
     filler.insert(originalJets, matchInfo.begin(), matchInfo.end());
     filler.fill();
   }
   evt.put(matching);
 }