Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge pull request #3321 from jstupak/jetToolboxPart2
Reco -- Code for grooming value maps and pileup jet ID
- Loading branch information
Showing
9 changed files
with
1,029 additions
and
3 deletions.
There are no files selected for viewing
122 changes: 122 additions & 0 deletions
122
CommonTools/RecoAlgos/plugins/JetDeltaRValueMapProducer.cc
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,122 @@ | ||
| /* \class PFJetSelector | ||
| * | ||
| * Selects jets with a configurable string-based cut, | ||
| * and also writes out the constituents of the jet | ||
| * into a separate collection. | ||
| * | ||
| * \author: Sal Rappoccio | ||
| * | ||
| * | ||
| * for more details about the cut syntax, see the documentation | ||
| * page below: | ||
| * | ||
| * https://twiki.cern.ch/twiki/bin/view/CMS/SWGuidePhysicsCutParser | ||
| * | ||
| * | ||
| */ | ||
|
|
||
|
|
||
| #include "FWCore/Framework/interface/EDProducer.h" | ||
|
|
||
| #include "DataFormats/JetReco/interface/Jet.h" | ||
| #include "DataFormats/JetReco/interface/PFJet.h" | ||
| #include "DataFormats/JetReco/interface/CaloJet.h" | ||
|
|
||
| #include "FWCore/Framework/interface/MakerMacros.h" | ||
| #include "CommonTools/Utils/interface/StringObjectFunction.h" | ||
| #include "FWCore/Framework/interface/Event.h" | ||
| #include "DataFormats/Common/interface/ValueMap.h" | ||
| #include "DataFormats/Math/interface/deltaR.h" | ||
|
|
||
| template < class T > | ||
| class JetDeltaRValueMapProducer : public edm::EDProducer { | ||
|
|
||
| public: | ||
|
|
||
| typedef std::vector<T> JetsInput; | ||
| typedef edm::ValueMap<float> JetValueMap; | ||
|
|
||
| JetDeltaRValueMapProducer ( edm::ParameterSet const & params ) : | ||
| srcToken_( consumes< typename edm::View<T> >( params.getParameter<edm::InputTag>("src") ) ), | ||
| matchedToken_( consumes< typename edm::View<T> >( params.getParameter<edm::InputTag>( "matched" ) ) ), | ||
| distMin_( params.getParameter<double>( "distMin" ) ), | ||
| value_( params.getParameter<std::string>("value") ), | ||
| evaluation_( value_ ) | ||
| { | ||
| produces< JetValueMap >(); | ||
| } | ||
|
|
||
| virtual ~JetDeltaRValueMapProducer() {} | ||
|
|
||
| virtual void beginJob() override {} | ||
| virtual void endJob() override {} | ||
|
|
||
| virtual void produce(edm::Event& iEvent, const edm::EventSetup& iSetup) override { | ||
|
|
||
|
|
||
|
|
||
| std::auto_ptr< JetValueMap > jetValueMap ( new JetValueMap() ); | ||
| edm::ValueMap<float>::Filler filler(*jetValueMap); | ||
|
|
||
| edm::Handle< typename edm::View<T> > h_jets1; | ||
| iEvent.getByToken( srcToken_, h_jets1 ); | ||
| edm::Handle< typename edm::View<T> > h_jets2; | ||
| iEvent.getByToken( matchedToken_, h_jets2 ); | ||
|
|
||
| std::vector<float> values; values.reserve( h_jets1->size() ); | ||
|
|
||
| // Now set the Ptrs with the orphan handles. | ||
| std::vector<float> v_jets2_eta, v_jets2_phi; | ||
| float jet1_eta, jet1_phi; | ||
| for ( typename edm::View<T>::const_iterator ibegin = h_jets1->begin(), | ||
| iend = h_jets1->end(), ijet = ibegin; | ||
| ijet != iend; ++ijet ) { | ||
| float minDR2=9999; | ||
| float value=-9999; | ||
|
|
||
| jet1_eta=ijet->eta(); | ||
| jet1_phi=ijet->phi(); | ||
|
|
||
| for ( typename edm::View<T>::const_iterator jbegin = h_jets2->begin(), | ||
| jend = h_jets2->end(), jjet = jbegin; | ||
| jjet != jend; ++jjet ) { | ||
|
|
||
| if(ijet==ibegin){ | ||
| v_jets2_eta.push_back(jjet->eta()); | ||
| v_jets2_phi.push_back(jjet->phi()); | ||
| } | ||
|
|
||
| int index=jjet - jbegin; | ||
| float dR2=reco::deltaR2(jet1_eta,jet1_phi,v_jets2_eta.at(index),v_jets2_phi.at(index)); | ||
| if ( dR2 < distMin_*distMin_ && dR2 < minDR2) { | ||
| // Check the selection | ||
| value = evaluation_(*jjet); | ||
| minDR2 = dR2; | ||
| } | ||
| }// end loop over matched jets | ||
|
|
||
| // Fill to the vector | ||
| values.push_back( value ); | ||
|
|
||
| }// end loop over src jets | ||
|
|
||
| filler.insert(h_jets1, values.begin(), values.end()); | ||
| filler.fill(); | ||
|
|
||
| // put in Event | ||
| iEvent.put(jetValueMap); | ||
|
|
||
| } | ||
|
|
||
| protected: | ||
| edm::EDGetTokenT< typename edm::View<T> > srcToken_; | ||
| edm::EDGetTokenT< typename edm::View<T> > matchedToken_; | ||
| double distMin_; | ||
| std::string value_; | ||
| StringObjectFunction<T> evaluation_; | ||
|
|
||
| }; | ||
|
|
||
| typedef JetDeltaRValueMapProducer<reco::Jet> RecoJetDeltaRValueMapProducer; | ||
|
|
||
| DEFINE_FWK_MODULE( RecoJetDeltaRValueMapProducer ); |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Oops, something went wrong.