From 3db457f666cfb4ba55d0b08a29be0e4f6f96ae7a Mon Sep 17 00:00:00 2001 From: swmukher Date: Thu, 19 Nov 2020 10:20:09 +0100 Subject: [PATCH 1/2] produce variable nLayerIT --- .../plugins/EgammaHLTGsfTrackVarProducer.cc | 11 +++++++++++ 1 file changed, 11 insertions(+) diff --git a/RecoEgamma/EgammaHLTProducers/plugins/EgammaHLTGsfTrackVarProducer.cc b/RecoEgamma/EgammaHLTProducers/plugins/EgammaHLTGsfTrackVarProducer.cc index a9f2645073e24..b74818d0088cf 100644 --- a/RecoEgamma/EgammaHLTProducers/plugins/EgammaHLTGsfTrackVarProducer.cc +++ b/RecoEgamma/EgammaHLTProducers/plugins/EgammaHLTGsfTrackVarProducer.cc @@ -65,6 +65,7 @@ class EgammaHLTGsfTrackVarProducer : public edm::global::EDProducer<> { const edm::EDPutTokenT oneOverESeedMinusOneOverPMapPutToken_; const edm::EDPutTokenT missingHitsMapPutToken_; const edm::EDPutTokenT validHitsMapPutToken_; + const edm::EDPutTokenT nLayerITMapPutToken_; const edm::EDPutTokenT chi2MapPutToken_; }; @@ -88,6 +89,7 @@ EgammaHLTGsfTrackVarProducer::EgammaHLTGsfTrackVarProducer(const edm::ParameterS missingHitsMapPutToken_{ produces("MissingHits").setBranchAlias("missinghits")}, validHitsMapPutToken_{produces("ValidHits").setBranchAlias("validhits")}, + nLayerITMapPutToken_{produces("NLayerIT").setBranchAlias("nlayerit")}, chi2MapPutToken_{produces("Chi2").setBranchAlias("chi2")} {} void EgammaHLTGsfTrackVarProducer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) { @@ -120,6 +122,7 @@ void EgammaHLTGsfTrackVarProducer::produce(edm::StreamID, edm::Event& iEvent, co reco::RecoEcalCandidateIsolationMap oneOverESeedMinusOneOverPMap(recoEcalCandHandle); reco::RecoEcalCandidateIsolationMap missingHitsMap(recoEcalCandHandle); reco::RecoEcalCandidateIsolationMap validHitsMap(recoEcalCandHandle); + reco::RecoEcalCandidateIsolationMap nLayerITMap(recoEcalCandHandle); reco::RecoEcalCandidateIsolationMap chi2Map(recoEcalCandHandle); for (unsigned int iRecoEcalCand = 0; iRecoEcalCand < recoEcalCandHandle->size(); ++iRecoEcalCand) { @@ -144,6 +147,7 @@ void EgammaHLTGsfTrackVarProducer::produce(edm::StreamID, edm::Event& iEvent, co } } + int nLayerITValue = -1; int validHitsValue = 0; float chi2Value = 9999999.; float missingHitsValue = 9999999; @@ -161,6 +165,7 @@ void EgammaHLTGsfTrackVarProducer::produce(edm::StreamID, edm::Event& iEvent, co : useDefaultValuesForEndcap_ && nrTracks >= 1; if (rmCutsDueToNrTracks || useDefaultValues) { + nLayerITValue = 100; dEtaInValue = 0; dEtaSeedInValue = 0; dPhiInValue = 0; @@ -205,6 +210,10 @@ void EgammaHLTGsfTrackVarProducer::produce(edm::StreamID, edm::Event& iEvent, co validHitsValue = gsfTracks[trkNr]->numberOfValidHits(); } + if (gsfTracks[trkNr]->hitPattern().pixelLayersWithMeasurement() > nLayerITValue) { + nLayerITValue = gsfTracks[trkNr]->hitPattern().pixelLayersWithMeasurement(); + } + if (gsfTracks[trkNr]->normalizedChi2() < chi2Value) { chi2Value = gsfTracks[trkNr]->normalizedChi2(); } @@ -232,6 +241,7 @@ void EgammaHLTGsfTrackVarProducer::produce(edm::StreamID, edm::Event& iEvent, co oneOverESeedMinusOneOverPMap.insert(recoEcalCandRef, oneOverESeedMinusOneOverPValue); missingHitsMap.insert(recoEcalCandRef, missingHitsValue); validHitsMap.insert(recoEcalCandRef, validHitsValue); + nLayerITMap.insert(recoEcalCandRef, nLayerITValue); chi2Map.insert(recoEcalCandRef, chi2Value); } @@ -242,6 +252,7 @@ void EgammaHLTGsfTrackVarProducer::produce(edm::StreamID, edm::Event& iEvent, co iEvent.emplace(oneOverESeedMinusOneOverPMapPutToken_, oneOverESeedMinusOneOverPMap); iEvent.emplace(missingHitsMapPutToken_, missingHitsMap); iEvent.emplace(validHitsMapPutToken_, validHitsMap); + iEvent.emplace(nLayerITMapPutToken_, nLayerITMap); iEvent.emplace(chi2MapPutToken_, chi2Map); } From 526910e259b65475f5a4cc45b18bc76365a0a17f Mon Sep 17 00:00:00 2001 From: swmukher Date: Thu, 19 Nov 2020 10:59:00 +0100 Subject: [PATCH 2/2] minor adjustment for code-checks test --- .../EgammaHLTProducers/plugins/EgammaHLTGsfTrackVarProducer.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/RecoEgamma/EgammaHLTProducers/plugins/EgammaHLTGsfTrackVarProducer.cc b/RecoEgamma/EgammaHLTProducers/plugins/EgammaHLTGsfTrackVarProducer.cc index b74818d0088cf..ce66e70a0daab 100644 --- a/RecoEgamma/EgammaHLTProducers/plugins/EgammaHLTGsfTrackVarProducer.cc +++ b/RecoEgamma/EgammaHLTProducers/plugins/EgammaHLTGsfTrackVarProducer.cc @@ -210,7 +210,7 @@ void EgammaHLTGsfTrackVarProducer::produce(edm::StreamID, edm::Event& iEvent, co validHitsValue = gsfTracks[trkNr]->numberOfValidHits(); } - if (gsfTracks[trkNr]->hitPattern().pixelLayersWithMeasurement() > nLayerITValue) { + if (gsfTracks[trkNr]->hitPattern().pixelLayersWithMeasurement() > nLayerITValue) { nLayerITValue = gsfTracks[trkNr]->hitPattern().pixelLayersWithMeasurement(); }