CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
PFTauL1TJetsMatching.cc
Go to the documentation of this file.
2 #include "Math/GenVector/VectorUtil.h"
7 
8 //
9 // class declaration
10 //
12  : tauSrc_(consumes<reco::PFTauCollection>(iConfig.getParameter<edm::InputTag>("TauSrc"))),
13  L1JetSrc_(consumes<trigger::TriggerFilterObjectWithRefs>(iConfig.getParameter<edm::InputTag>("L1JetSrc"))),
14  matchingR2_(iConfig.getParameter<double>("MatchingdR") * iConfig.getParameter<double>("MatchingdR")),
15  minTauPt_(iConfig.getParameter<double>("MinTauPt")),
16  minL1TPt_(iConfig.getParameter<double>("MinL1TPt")) {
17  produces<reco::PFTauCollection>();
18 }
20 
22  std::unique_ptr<reco::PFTauCollection> L1TmatchedPFTau(new reco::PFTauCollection);
23 
25  iEvent.getByToken(tauSrc_, taus);
26 
28  iEvent.getByToken(L1JetSrc_, L1Jets);
29 
30  l1t::JetVectorRef jetCandRefVec;
31  L1Jets->getObjects(trigger::TriggerL1Jet, jetCandRefVec);
32 
33  /* Loop over taus that must pass a certain minTauPt_ cut */
34  /* then loop over L1T jets that must pass minL1TPt_ */
35  /* and check whether they match, if yes -> include the taus in */
36  /* the new L1T matched PFTau collection */
37 
38  for (unsigned int iTau = 0; iTau < taus->size(); iTau++) {
39  bool isMatched = false;
40  if ((*taus)[iTau].pt() > minTauPt_) {
41  for (unsigned int iJet = 0; iJet < jetCandRefVec.size(); iJet++) {
42  if (jetCandRefVec[iJet]->pt() > minL1TPt_) {
43  if (reco::deltaR2((*taus)[iTau].p4(), jetCandRefVec[iJet]->p4()) < matchingR2_) {
44  isMatched = true;
45  break;
46  }
47  }
48  }
49  }
50  if (isMatched)
51  L1TmatchedPFTau->push_back((*taus)[iTau]);
52  }
53  iEvent.put(std::move(L1TmatchedPFTau));
54 }
55 
58  desc.add<edm::InputTag>("L1JetSrc", edm::InputTag("hltL1VBFDiJetOR"))
59  ->setComment("Input filter objects passing L1 seed");
60  desc.add<edm::InputTag>("TauSrc", edm::InputTag("hltSelectedPFTausTrackFindingLooseChargedIsolationAgainstMuon"))
61  ->setComment("Input collection of PFTaus");
62  desc.add<double>("MatchingdR", 0.5)->setComment("Maximum dR for matching between PFTaus and L1 filter jets");
63  desc.add<double>("MinTauPt", 20.0)->setComment("PFTaus above this pt will be considered");
64  desc.add<double>("MinL1TPt", 115.0)->setComment("L1T Objects above this pt will be considered");
65  descriptions.setComment(
66  "This module produces a collection of PFTaus matched to the leading jet passing the L1 seed filter.");
67  descriptions.add("PFTauL1TJetsMatching", desc);
68 }
void setComment(std::string const &value)
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:133
std::vector< PFTau > PFTauCollection
collection of PFTau objects
Definition: PFTauFwd.h:9
std::vector< l1t::PFTau > PFTauCollection
Definition: PFTau.h:61
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > L1JetSrc_
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
void produce(edm::StreamID, edm::Event &, const edm::EventSetup &) const override
int iEvent
Definition: GenABIO.cc:224
def move
Definition: eostools.py:511
bool isMatched(TrackingRecHit const &hit)
ParameterDescriptionBase * add(U const &iLabel, T const &value)
void setComment(std::string const &value)
constexpr auto deltaR2(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:16
void add(std::string const &label, ParameterSetDescription const &psetDescription)
PFTauL1TJetsMatching(const edm::ParameterSet &)
const edm::EDGetTokenT< reco::PFTauCollection > tauSrc_
std::vector< JetRef > JetVectorRef
Definition: Jet.h:14