Skip to content

Commit

Permalink
Merge pull request #37565 from arsahasransu/ForTSGv2Menu
Browse files Browse the repository at this point in the history
Change tracks to optional requirement in HLTDisplacedEgammaFilter
  • Loading branch information
cmsbuild committed Apr 14, 2022
2 parents 1290e52 + d8ad893 commit a0b62d8
Show file tree
Hide file tree
Showing 2 changed files with 22 additions and 13 deletions.
34 changes: 21 additions & 13 deletions HLTrigger/Egamma/plugins/HLTDisplacedEgammaFilter.cc
Expand Up @@ -36,11 +36,14 @@ HLTDisplacedEgammaFilter::HLTDisplacedEgammaFilter(const edm::ParameterSet& iCon
sMaj_max = iConfig.getParameter<double>("sMaj_max");
seedTimeMin = iConfig.getParameter<double>("seedTimeMin");
seedTimeMax = iConfig.getParameter<double>("seedTimeMax");
useTrackVeto = iConfig.getParameter<bool>("useTrackVeto");

inputToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(inputTag_);
rechitsEBToken_ = consumes<EcalRecHitCollection>(rechitsEB);
rechitsEEToken_ = consumes<EcalRecHitCollection>(rechitsEE);
inputTrkToken_ = consumes<reco::TrackCollection>(inputTrk);
if (useTrackVeto) {
inputTrkToken_ = consumes<reco::TrackCollection>(inputTrk);
}
}

HLTDisplacedEgammaFilter::~HLTDisplacedEgammaFilter() = default;
Expand All @@ -61,6 +64,7 @@ void HLTDisplacedEgammaFilter::fillDescriptions(edm::ConfigurationDescriptions&
desc.add<double>("sMaj_max", 999.0);
desc.add<double>("seedTimeMin", -25.0);
desc.add<double>("seedTimeMax", 25.0);
desc.add<bool>("useTrackVeto", true);
desc.add<int>("maxTrackCut", 0);
desc.add<double>("trackPtCut", 3.0);
desc.add<double>("trackdRCut", 0.5);
Expand All @@ -86,7 +90,9 @@ bool HLTDisplacedEgammaFilter::hltFilter(edm::Event& iEvent,

// get hold of collection of objects
edm::Handle<reco::TrackCollection> tracks;
iEvent.getByToken(inputTrkToken_, tracks);
if (useTrackVeto) {
iEvent.getByToken(inputTrkToken_, tracks);
}

// get the EcalRecHit
edm::Handle<EcalRecHitCollection> rechitsEB_;
Expand Down Expand Up @@ -128,19 +134,21 @@ bool HLTDisplacedEgammaFilter::hltFilter(edm::Event& iEvent,

//Track Veto

int nTrk = 0;
for (auto const& it : *tracks) {
if (it.pt() < trkPtCut)
continue;
LorentzVector trkP4(it.px(), it.py(), it.pz(), it.p());
double dR = ROOT::Math::VectorUtil::DeltaR(trkP4, ref->p4());
if (dR < trkdRCut)
nTrk++;
if (useTrackVeto) {
int nTrk = 0;
for (auto const& it : *tracks) {
if (it.pt() < trkPtCut)
continue;
LorentzVector trkP4(it.px(), it.py(), it.pz(), it.p());
double dR = ROOT::Math::VectorUtil::DeltaR(trkP4, ref->p4());
if (dR < trkdRCut)
nTrk++;
if (nTrk > maxTrkCut)
break;
}
if (nTrk > maxTrkCut)
break;
continue;
}
if (nTrk > maxTrkCut)
continue;

n++;
// std::cout << "Passed eta: " << ref->eta() << std::endl;
Expand Down
1 change: 1 addition & 0 deletions HLTrigger/Egamma/plugins/HLTDisplacedEgammaFilter.h
Expand Up @@ -46,6 +46,7 @@ class HLTDisplacedEgammaFilter : public HLTFilter {
double sMaj_max;
double seedTimeMin;
double seedTimeMax;
bool useTrackVeto;

edm::InputTag inputTrk;
edm::EDGetTokenT<reco::TrackCollection> inputTrkToken_;
Expand Down

0 comments on commit a0b62d8

Please sign in to comment.