Skip to content
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.

Sign up for GitHub

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

Run3-205D Last step for modernize the HCAL AlCaReco producers #35509

Merged
merged 2 commits into from
Oct 2, 2021
Merged

Run3-205D Last step for modernize the HCAL AlCaReco producers #35509

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
31c743f
Modernize the remaining HCAL AlCaReco producers
Oct 2, 2021
5ec10cc
Code check
Oct 2, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
75 changes: 43 additions & 32 deletions Calibration/HcalAlCaRecoProducers/plugins/AlCaGammaJetProducer.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,10 +5,11 @@
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EgammaCandidates/interface/Photon.h"
Expand Down Expand Up @@ -39,7 +40,9 @@
class AlCaGammaJetProducer : public edm::global::EDProducer<> {
public:
explicit AlCaGammaJetProducer(const edm::ParameterSet&);
~AlCaGammaJetProducer() override = default;
void produce(edm::StreamID, edm::Event&, const edm::EventSetup&) const override;
static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

private:
bool select(const reco::PhotonCollection&, const reco::PFJetCollection&) const;
Expand Down Expand Up @@ -156,107 +159,95 @@ bool AlCaGammaJetProducer::select(const reco::PhotonCollection& ph, const reco::
}
return false;
}

// ------------ method called to produce the data ------------
void AlCaGammaJetProducer::produce(edm::StreamID, edm::Event& iEvent, const edm::EventSetup&) const {
// Access the collections from iEvent
edm::Handle<reco::PhotonCollection> phoHandle;
iEvent.getByToken(tok_Photon_, phoHandle);
auto const& phoHandle = iEvent.getHandle(tok_Photon_);
if (!phoHandle.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get the product " << labelPhoton_;
return;
}
const reco::PhotonCollection& photon = *(phoHandle.product());

edm::Handle<reco::PFJetCollection> pfjet;
iEvent.getByToken(tok_PFJet_, pfjet);
auto const& pfjet = iEvent.getHandle(tok_PFJet_);
if (!pfjet.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelPFJet_;
return;
}
const reco::PFJetCollection& pfjets = *(pfjet.product());

edm::Handle<reco::PFCandidateCollection> pfc;
iEvent.getByToken(tok_PFCand_, pfc);
auto const& pfc = iEvent.getHandle(tok_PFCand_);
if (!pfc.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelPFCandidate_;
return;
}
const reco::PFCandidateCollection& pfcand = *(pfc.product());

edm::Handle<reco::VertexCollection> vt;
iEvent.getByToken(tok_Vertex_, vt);
auto const& vt = iEvent.getHandle(tok_Vertex_);
if (!vt.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelVertex_;
return;
}
const reco::VertexCollection& vtx = *(vt.product());

edm::Handle<reco::PFMETCollection> pfmt;
iEvent.getByToken(tok_PFMET_, pfmt);
auto const& pfmt = iEvent.getHandle(tok_PFMET_);
if (!pfmt.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelPFMET_;
return;
}
const reco::PFMETCollection& pfmet = *(pfmt.product());

edm::Handle<edm::SortedCollection<HBHERecHit, edm::StrictWeakOrdering<HBHERecHit>>> hbhe;
iEvent.getByToken(tok_HBHE_, hbhe);
auto const& hbhe = iEvent.getHandle(tok_HBHE_);
if (!hbhe.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelHBHE_;
return;
}
const edm::SortedCollection<HBHERecHit, edm::StrictWeakOrdering<HBHERecHit>>& Hithbhe = *(hbhe.product());

edm::Handle<edm::SortedCollection<HORecHit, edm::StrictWeakOrdering<HORecHit>>> ho;
iEvent.getByToken(tok_HO_, ho);
auto const& ho = iEvent.getHandle(tok_HO_);
if (!ho.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelHO_;
return;
}
const edm::SortedCollection<HORecHit, edm::StrictWeakOrdering<HORecHit>>& Hitho = *(ho.product());

edm::Handle<edm::SortedCollection<HFRecHit, edm::StrictWeakOrdering<HFRecHit>>> hf;
iEvent.getByToken(tok_HF_, hf);
auto const& hf = iEvent.getHandle(tok_HF_);
if (!hf.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelHF_;
return;
}
const edm::SortedCollection<HFRecHit, edm::StrictWeakOrdering<HFRecHit>>& Hithf = *(hf.product());

edm::Handle<edm::TriggerResults> trig;
iEvent.getByToken(tok_TrigRes_, trig);
auto const& trig = iEvent.getHandle(tok_TrigRes_);
if (!trig.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelTrigger_;
return;
}
const edm::TriggerResults& trigres = *(trig.product());

edm::Handle<double> rh;
iEvent.getByToken(tok_Rho_, rh);
auto const& rh = iEvent.getHandle(tok_Rho_);
if (!rh.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelRho_;
return;
}
const double rho_val = *(rh.product());

edm::Handle<reco::GsfElectronCollection> gsf;
iEvent.getByToken(tok_GsfElec_, gsf);
auto const& gsf = iEvent.getHandle(tok_GsfElec_);
if (!gsf.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelGsfEle_;
return;
}
const reco::GsfElectronCollection& gsfele = *(gsf.product());

edm::Handle<reco::ConversionCollection> con;
iEvent.getByToken(tok_Conv_, con);
auto const& con = iEvent.getHandle(tok_Conv_);
if (!con.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelConv_;
return;
}
const reco::ConversionCollection& conv = *(con.product());

edm::Handle<reco::BeamSpot> bs;
iEvent.getByToken(tok_BS_, bs);
auto const& bs = iEvent.getHandle(tok_BS_);
if (!bs.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelBeamSpot_;
return;
Expand Down Expand Up @@ -287,10 +278,8 @@ void AlCaGammaJetProducer::produce(edm::StreamID, edm::Event& iEvent, const edm:
iEvent.emplace(put_rho_, rho_val);
iEvent.emplace(put_conv_, conv);

edm::Handle<edm::ValueMap<Bool_t>> loosePhotonQual;
iEvent.getByToken(tok_loosePhoton_, loosePhotonQual);
edm::Handle<edm::ValueMap<Bool_t>> tightPhotonQual;
iEvent.getByToken(tok_tightPhoton_, tightPhotonQual);
auto const& loosePhotonQual = iEvent.getHandle(tok_loosePhoton_);
auto const& tightPhotonQual = iEvent.getHandle(tok_tightPhoton_);
if (loosePhotonQual.isValid() && tightPhotonQual.isValid()) {
miniLoosePhoton.reserve(photon.size());
miniTightPhoton.reserve(photon.size());
Expand Down Expand Up @@ -329,4 +318,26 @@ void AlCaGammaJetProducer::produce(edm::StreamID, edm::Event& iEvent, const edm:
return;
}

void AlCaGammaJetProducer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
edm::ParameterSetDescription desc;
desc.add<edm::InputTag>("PhoInput", edm::InputTag("gedPhotons"));
desc.add<edm::InputTag>("PFjetInput", edm::InputTag("ak4PFJetsCHS"));
desc.add<edm::InputTag>("HBHEInput", edm::InputTag("hbhereco"));
desc.add<edm::InputTag>("HFInput", edm::InputTag("hfreco"));
desc.add<edm::InputTag>("HOInput", edm::InputTag("horeco"));
desc.add<edm::InputTag>("METInput", edm::InputTag("pfMet"));
desc.add<edm::InputTag>("TriggerResults", edm::InputTag("TriggerResults::HLT"));
desc.add<edm::InputTag>("gsfeleInput", edm::InputTag("gedGsfElectrons"));
desc.add<edm::InputTag>("particleFlowInput", edm::InputTag("particleFlow"));
desc.add<edm::InputTag>("VertexInput", edm::InputTag("offlinePrimaryVertices"));
desc.add<edm::InputTag>("ConversionsInput", edm::InputTag("allConversions"));
desc.add<edm::InputTag>("rhoInput", edm::InputTag("fixedGridRhoFastjetAll"));
desc.add<edm::InputTag>("BeamSpotInput", edm::InputTag("offlineBeamSpot"));
desc.add<edm::InputTag>("PhoLoose", edm::InputTag("PhotonIDProdGED", "PhotonCutBasedIDLoose"));
desc.add<edm::InputTag>("PhoTight", edm::InputTag("PhotonIDProdGED", "PhotonCutBasedIDTight"));
desc.add<double>("MinPtJet", 10.0);
desc.add<double>("MinPtPhoton", 10.0);
descriptions.add("alcaGammaJetProducer", desc);
}

DEFINE_FWK_MODULE(AlCaGammaJetProducer);
13 changes: 5 additions & 8 deletions Calibration/HcalAlCaRecoProducers/plugins/AlCaGammaJetSelector.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -27,6 +27,7 @@
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"

#include "DataFormats/EgammaCandidates/interface/Photon.h"
#include "DataFormats/EgammaCandidates/interface/PhotonFwd.h"
Expand All @@ -47,7 +48,7 @@ namespace AlCaGammaJet {
class AlCaGammaJetSelector : public edm::stream::EDFilter<edm::GlobalCache<AlCaGammaJet::Counters> > {
public:
explicit AlCaGammaJetSelector(const edm::ParameterSet&, const AlCaGammaJet::Counters* count);
~AlCaGammaJetSelector() override;
~AlCaGammaJetSelector() override = default;

static std::unique_ptr<AlCaGammaJet::Counters> initializeGlobalCache(edm::ParameterSet const&) {
return std::make_unique<AlCaGammaJet::Counters>();
Expand Down Expand Up @@ -101,8 +102,6 @@ AlCaGammaJetSelector::AlCaGammaJetSelector(const edm::ParameterSet& iConfig, con
tok_PFJet_ = consumes<reco::PFJetCollection>(labelPFJet_);
}

AlCaGammaJetSelector::~AlCaGammaJetSelector() {}

//
// member functions
//
Expand All @@ -116,24 +115,22 @@ void AlCaGammaJetSelector::fillDescriptions(edm::ConfigurationDescriptions& desc
desc.add<edm::InputTag>("PFjetInput", edm::InputTag("ak4PFJetsCHS"));
desc.add<double>("MinPtJet", 10.0);
desc.add<double>("MinPtPhoton", 10.0);
descriptions.addDefault(desc);
descriptions.add("alcaGammaJetSelector", desc);
}

// ------------ method called on each new Event ------------
bool AlCaGammaJetSelector::filter(edm::Event& iEvent, const edm::EventSetup& iSetup) {
nProcessed_++;

// Access the collections from iEvent
edm::Handle<reco::PhotonCollection> phoHandle;
iEvent.getByToken(tok_Photon_, phoHandle);
auto const& phoHandle = iEvent.getHandle(tok_Photon_);
if (!phoHandle.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get the product " << labelPhoton_;
return false; // do not filter
}
const reco::PhotonCollection photons = *(phoHandle.product());

edm::Handle<reco::PFJetCollection> pfjetHandle;
iEvent.getByToken(tok_PFJet_, pfjetHandle);
auto const& pfjetHandle = iEvent.getHandle(tok_PFJet_);
if (!pfjetHandle.isValid()) {
edm::LogWarning("AlCaGammaJet") << "AlCaGammaJetProducer: Error! can't get product " << labelPFJet_;
return false; // do not filter
Expand Down
Loading