-
Notifications
You must be signed in to change notification settings - Fork 4.3k
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
adding a PFJets-Tau overlap removal module for HLT #18997
Changes from 6 commits
81e2091
73364f2
610e68a
a60bc8e
8a15e70
6b772d8
9a4e110
33af49b
74b124e
File filter
Filter by extension
Conversations
Jump to
Diff view
Diff view
There are no files selected for viewing
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,35 @@ | ||
#ifndef RecoTauTag_HLTProducers_PFJetsTauOverlapRemoval_H | ||
#define RecoTauTag_HLTProducers_PFJetsTauOverlapRemoval_H | ||
|
||
// user include files | ||
#include "FWCore/Framework/interface/Frameworkfwd.h" | ||
#include "FWCore/Framework/interface/global/EDProducer.h" | ||
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h" | ||
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h" | ||
#include "FWCore/Framework/interface/Event.h" | ||
#include "FWCore/ParameterSet/interface/ParameterSet.h" | ||
#include "FWCore/Utilities/interface/InputTag.h" | ||
#include "DataFormats/Common/interface/Handle.h" | ||
#include "DataFormats/L1Trigger/interface/Tau.h" | ||
#include "DataFormats/JetReco/interface/CaloJetCollection.h" | ||
#include "DataFormats/TauReco/interface/PFTauFwd.h" | ||
#include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h" | ||
#include "DataFormats/HLTReco/interface/TriggerObject.h" | ||
#include "DataFormats/HLTReco/interface/TriggerEvent.h" | ||
|
||
#include <map> | ||
#include <vector> | ||
class PFJetsTauOverlapRemoval: public edm::global::EDProducer<> { | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. could this not be a stream module? There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. I'm not sure, I haven't found anything that does exactly this.. |
||
public: | ||
explicit PFJetsTauOverlapRemoval(const edm::ParameterSet&); | ||
~PFJetsTauOverlapRemoval(); | ||
virtual void produce(edm::StreamID, edm::Event&, const edm::EventSetup&) const override; | ||
static void fillDescriptions(edm::ConfigurationDescriptions& descriptions); | ||
|
||
private: | ||
|
||
const edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> tauSrc_; | ||
const edm::EDGetTokenT<reco::PFJetCollection> pfJetSrc_; | ||
const double mindR_; | ||
}; | ||
#endif |
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,60 @@ | ||
#include "RecoTauTag/HLTProducers/interface/PFJetsTauOverlapRemoval.h" | ||
#include "Math/GenVector/VectorUtil.h" | ||
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h" | ||
#include "FWCore/Utilities/interface/EDMException.h" | ||
#include "DataFormats/TauReco/interface/PFTau.h" | ||
|
||
// | ||
// class declaration | ||
// | ||
PFJetsTauOverlapRemoval::PFJetsTauOverlapRemoval(const edm::ParameterSet& iConfig): | ||
tauSrc_ ( consumes<trigger::TriggerFilterObjectWithRefs>(iConfig.getParameter<edm::InputTag>("TauSrc" ) ) ), | ||
pfJetSrc_ ( consumes<reco::PFJetCollection>(iConfig.getParameter<edm::InputTag>("PFJetSrc") ) ), | ||
mindR_ ( iConfig.getParameter<double>("Min_dR") ) | ||
{ | ||
produces<reco::PFJetCollection>(); | ||
} | ||
PFJetsTauOverlapRemoval::~PFJetsTauOverlapRemoval(){ } | ||
|
||
void PFJetsTauOverlapRemoval::produce(edm::StreamID iSId, edm::Event& iEvent, const edm::EventSetup& iES) const | ||
{ | ||
|
||
std::unique_ptr<reco::PFJetCollection> cleanedPFJets(new reco::PFJetCollection); | ||
|
||
double matchingR2 = mindR_*mindR_; | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. better to do this in the constructor |
||
|
||
edm::Handle<trigger::TriggerFilterObjectWithRefs> tauJets; | ||
iEvent.getByToken(tauSrc_, tauJets); | ||
|
||
edm::Handle<reco::PFJetCollection> PFJets; | ||
iEvent.getByToken(pfJetSrc_,PFJets); | ||
|
||
trigger::VRpftau taus; | ||
tauJets->getObjects(trigger::TriggerTau,taus); | ||
|
||
if(PFJets->size() > 1){ | ||
for(unsigned int iJet = 0; iJet < 2; iJet++){ | ||
bool isMatched = false; | ||
const reco::PFJet & myPFJet = (*PFJets)[iJet]; | ||
for(unsigned int iTau = 0; iTau < taus.size(); iTau++){ | ||
if(ROOT::Math::VectorUtil::DeltaR2((taus[iTau]->p4()).Vect(), myPFJet.p4().Vect()) < matchingR2){ | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. this is still not reworked as suggested.... |
||
isMatched = true; | ||
break; | ||
} | ||
} | ||
if(isMatched == false) cleanedPFJets->push_back(myPFJet); | ||
} | ||
if(PFJets->size() > 2) cleanedPFJets->push_back((*PFJets)[2]); | ||
} | ||
iEvent.put(std::move(cleanedPFJets)); | ||
} | ||
|
||
void PFJetsTauOverlapRemoval::fillDescriptions(edm::ConfigurationDescriptions& descriptions) | ||
{ | ||
edm::ParameterSetDescription desc; | ||
desc.add<edm::InputTag>("PFJetSrc", edm::InputTag("hltAK4PFJetsCorrected"))->setComment("Input collection of PFJets" ); | ||
desc.add<edm::InputTag>("TauSrc", edm::InputTag("hltPFTau20TrackLooseIso"))->setComment("Input collection of PFTaus that have passed ID and isolation requirements"); | ||
desc.add<double> ("Min_dR",0.5)->setComment("Minimum dR outside of which PFJets will be saved"); | ||
descriptions.setComment("This module produces a collection of PFJets that are cross-cleaned with respect to PFTaus passing a HLT filter."); | ||
descriptions.add ("PFJetsTauOverlapRemoval",desc); | ||
} |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
are the map or vector includes needed here? It doesn't seem so