Inheritance diagram for L1TPFCaloProducer:

Classes
struct	SimpleHGCTC

Public Member Functions
	L1TPFCaloProducer (const edm::ParameterSet &)

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

void	readHcalDigis_ (edm::Event &event, const edm::EventSetup &)

void	readHcalHGCTowers_ (edm::Event &event, const edm::EventSetup &)

void	readPhase2BarrelCaloTowers_ (edm::Event &event, const edm::EventSetup &)

Private Attributes
std::unique_ptr< l1tpf_calo::SimpleCaloLinkerBase >	caloLinker_

bool	debug_

edm::ESGetToken< CaloTPGTranscoder, CaloTPGRecord >	decoderTag_

std::vector< edm::EDGetTokenT< reco::CandidateView > >	ecalCands_

l1tpf_calo::SingleCaloClusterer	ecalClusterer_

bool	ecalOnly_

l1tpf::corrector	emCorrector_

l1tpf::corrector	hadCorrector_

std::vector< edm::EDGetTokenT< reco::CandidateView > >	hcalCands_

l1tpf_calo::SingleCaloClusterer	hcalClusterer_

std::vector< edm::EDGetTokenT< HcalTrigPrimDigiCollection > >	hcalDigis_

bool	hcalDigisBarrel_

bool	hcalDigisHF_

std::vector< edm::EDGetTokenT< l1t::HGCalTowerBxCollection > >	hcalHGCTowers_

bool	hcalHGCTowersHadOnly_

l1tpf::corrector	hcCorrector_

std::vector< edm::EDGetTokenT< l1tp2::CaloTowerCollection > >	phase2barrelTowers_

l1tpf::ParametricResolution	resol_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 30 of file L1TPFCaloProducer.cc.

Constructor & Destructor Documentation

◆ L1TPFCaloProducer()

L1TPFCaloProducer::L1TPFCaloProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 69 of file L1TPFCaloProducer.cc.

References ecalCands_, ecalOnly_, edm::ParameterSet::getParameter(), hcalCands_, hcalDigis_, hcalDigisBarrel_, hcalDigisHF_, hcalHGCTowers_, hcalHGCTowersHadOnly_, phase2barrelTowers_, and makeGlobalPositionRcd_cfg::tag.

     : ecalOnly_(iConfig.existsAs<bool>("ecalOnly") ? iConfig.getParameter<bool>("ecalOnly") : false),
       debug_(iConfig.getUntrackedParameter<int>("debug", 0)),
       decoderTag_(esConsumes<CaloTPGTranscoder, CaloTPGRecord>(edm::ESInputTag("", ""))),
       emCorrector_(iConfig.getParameter<std::string>("emCorrector"), -1, debug_),
       hcCorrector_(iConfig.getParameter<std::string>("hcCorrector"), -1, debug_),
       hadCorrector_(iConfig.getParameter<std::string>("hadCorrector"),
                     iConfig.getParameter<double>("hadCorrectorEmfMax"),
                     debug_),
       ecalClusterer_(iConfig.getParameter<edm::ParameterSet>("ecalClusterer")),
       hcalClusterer_(iConfig.getParameter<edm::ParameterSet>("hcalClusterer")),
       caloLinker_(l1tpf_calo::makeCaloLinker(
           iConfig.getParameter<edm::ParameterSet>("linker"), ecalClusterer_, hcalClusterer_)),
       resol_(iConfig.getParameter<edm::ParameterSet>("resol")) {
   produces<l1t::PFClusterCollection>("ecalCells");
 
   produces<l1t::PFClusterCollection>("emCalibrated");
   produces<l1t::PFClusterCollection>("emUncalibrated");
 
   for (auto &tag : iConfig.getParameter<std::vector<edm::InputTag>>("ecalCandidates")) {
     ecalCands_.push_back(consumes<reco::CandidateView>(tag));
   }
 
   if (ecalOnly_)
     return;
 
   produces<l1t::PFClusterCollection>("hcalCells");
 
   produces<l1t::PFClusterCollection>("hcalUnclustered");
   produces<l1t::PFClusterCollection>("hcalUncalibrated");
   produces<l1t::PFClusterCollection>("hcalCalibrated");
 
   produces<l1t::PFClusterCollection>("uncalibrated");
   produces<l1t::PFClusterCollection>("calibrated");
 
   for (auto &tag : iConfig.getParameter<std::vector<edm::InputTag>>("hcalCandidates")) {
     hcalCands_.push_back(consumes<reco::CandidateView>(tag));
   }
 
   for (auto &tag : iConfig.getParameter<std::vector<edm::InputTag>>("hcalDigis")) {
     hcalDigis_.push_back(consumes<HcalTrigPrimDigiCollection>(tag));
   }
   if (!hcalDigis_.empty()) {
     hcalDigisBarrel_ = iConfig.getParameter<bool>("hcalDigisBarrel");
     hcalDigisHF_ = iConfig.getParameter<bool>("hcalDigisHF");
   }
 
   for (auto &tag : iConfig.getParameter<std::vector<edm::InputTag>>("phase2barrelCaloTowers")) {
     phase2barrelTowers_.push_back(consumes<l1tp2::CaloTowerCollection>(tag));
   }
 
   for (auto &tag : iConfig.getParameter<std::vector<edm::InputTag>>("hcalHGCTowers")) {
     hcalHGCTowers_.push_back(consumes<l1t::HGCalTowerBxCollection>(tag));
   }
   if (!hcalHGCTowers_.empty())
     hcalHGCTowersHadOnly_ = iConfig.getParameter<bool>("hcalHGCTowersHadOnly");
 }

Member Function Documentation

◆ produce()

void L1TPFCaloProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

-------———ECAL INFO----------———

----------——HCAL INFO-------------——

------------— CLUSTERING ---------------—

Definition at line 128 of file L1TPFCaloProducer.cc.

References funct::abs(), l1tpf_calo::SingleCaloClusterer::add(), HltBtagPostValidation_cff::c, caloLinker_, l1tpf_calo::SingleCaloClusterer::clear(), bsc_activity_cfg::clusters, l1tpf_calo::SingleCaloClusterer::correct(), l1tpf::corrector::correctedPt(), debug_, ecalCands_, ecalClusterer_, ecalOnly_, emCorrector_, l1tpf_calo::SingleCaloClusterer::fetch(), l1tpf_calo::SingleCaloClusterer::fetchCells(), hadCorrector_, hcalCands_, hcalClusterer_, hcalDigis_, hcalHGCTowers_, hcCorrector_, iEvent, eostools::move(), phase2barrelTowers_, readHcalDigis_(), readHcalHGCTowers_(), readPhase2BarrelCaloTowers_(), resol_, l1tpf_calo::SingleCaloClusterer::run(), unpackBuffers-CaloStage2::token, and l1tpf::corrector::valid().

                                                                              {
   edm::Handle<reco::CandidateView> ecals;
   for (const auto &token : ecalCands_) {
     iEvent.getByToken(token, ecals);
     for (const reco::Candidate &it : *ecals) {
       if (debug_)
         edm::LogWarning("L1TPFCaloProducer")
             << "adding ECal input pt " << it.pt() << ", eta " << it.eta() << ", phi " << it.phi() << "\n";
       ecalClusterer_.add(it);
     }
   }
 
   if (!ecalOnly_) {
     edm::Handle<reco::CandidateView> hcals;
     for (const auto &token : hcalCands_) {
       iEvent.getByToken(token, hcals);
       for (const reco::Candidate &it : *hcals) {
         if (debug_)
           edm::LogWarning("L1TPFCaloProducer")
               << "adding HCal cand input pt " << it.pt() << ", eta " << it.eta() << ", phi " << it.phi() << "\n";
         hcalClusterer_.add(it);
       }
     }
     if (!hcalDigis_.empty()) {
       readHcalDigis_(iEvent, iSetup);
     }
     if (!phase2barrelTowers_.empty()) {
       readPhase2BarrelCaloTowers_(iEvent, iSetup);
     }
     if (!hcalHGCTowers_.empty()) {
       readHcalHGCTowers_(iEvent, iSetup);
     }
   }
 
   ecalClusterer_.run();
 
   auto ecalCellsH = iEvent.put(ecalClusterer_.fetchCells(), "ecalCells");
 
   iEvent.put(ecalClusterer_.fetch(ecalCellsH), "emUncalibrated");
 
   if (emCorrector_.valid()) {
     ecalClusterer_.correct(
         [&](const l1tpf_calo::Cluster &c) -> float { return emCorrector_.correctedPt(0., c.et, std::abs(c.eta)); });
   }
 
   std::unique_ptr<l1t::PFClusterCollection> corrEcal = ecalClusterer_.fetch(ecalCellsH);
 
   if (debug_) {
     for (const l1t::PFCluster &it : *corrEcal) {
       edm::LogWarning("L1TPFCaloProducer")
           << "corrected ECal cluster pt " << it.pt() << ", eta " << it.eta() << ", phi " << it.phi() << "\n";
     }
   }
 
   auto ecalClustH = iEvent.put(std::move(corrEcal), "emCalibrated");
 
   if (ecalOnly_) {
     ecalClusterer_.clear();
     return;
   }
 
   hcalClusterer_.run();
 
   auto hcalCellsH = iEvent.put(hcalClusterer_.fetchCells(), "hcalCells");
 
   // this we put separately for debugging
   iEvent.put(hcalClusterer_.fetchCells(/*unclustered=*/true), "hcalUnclustered");
 
   iEvent.put(hcalClusterer_.fetch(hcalCellsH), "hcalUncalibrated");
 
   if (hcCorrector_.valid()) {
     hcalClusterer_.correct(
         [&](const l1tpf_calo::Cluster &c) -> float { return hcCorrector_.correctedPt(c.et, 0., std::abs(c.eta)); });
   }
 
   auto hcalClustH = iEvent.put(hcalClusterer_.fetch(hcalCellsH), "hcalCalibrated");
 
   // Calorimeter linking
   caloLinker_->run();
 
   iEvent.put(caloLinker_->fetch(ecalClustH, hcalClustH), "uncalibrated");
 
   if (hadCorrector_.valid()) {
     caloLinker_->correct([&](const l1tpf_calo::CombinedCluster &c) -> float {
       if (debug_)
         edm::LogWarning("L1TPFCaloProducer") << "raw linked cluster pt " << c.et << ", eta " << c.eta << ", phi "
                                              << c.phi << ", emPt " << c.ecal_et << "\n";
       return hadCorrector_.correctedPt(c.et, c.ecal_et, std::abs(c.eta));
     });
   }
 
   std::unique_ptr<l1t::PFClusterCollection> clusters = caloLinker_->fetch(ecalClustH, hcalClustH);
   for (l1t::PFCluster &c : *clusters) {
     c.setPtError(resol_(c.pt(), std::abs(c.eta())));
     if (debug_)
       edm::LogWarning("L1TPFCaloProducer") << "calibrated linked cluster pt " << c.pt() << ", eta " << c.eta()
                                            << ", phi " << c.phi() << ", emPt " << c.emEt() << "\n";
   }
   iEvent.put(std::move(clusters), "calibrated");
 
   ecalClusterer_.clear();
   hcalClusterer_.clear();
   caloLinker_->clear();
 }

◆ readHcalDigis_()

void L1TPFCaloProducer::readHcalDigis_	(	edm::Event &	event,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 236 of file L1TPFCaloProducer.cc.

References funct::abs(), l1tpf_calo::SingleCaloClusterer::add(), debug_, decoderTag_, EgHLTOffHistBins_cfi::et, edm::EventSetup::getData(), hcalClusterer_, hcalDigis_, hcalDigisBarrel_, hcalDigisHF_, hcalRecHitTable_cff::ieta, iEvent, hcalRecHitTable_cff::iphi, unpackBuffers-CaloStage2::token, l1t::CaloTools::towerEta(), and l1t::CaloTools::towerPhi().

Referenced by produce().

                                                                                     {
   const auto &decoder = iSetup.getData(decoderTag_);
   edm::Handle<HcalTrigPrimDigiCollection> hcalTPs;
   for (const auto &token : hcalDigis_) {
     iEvent.getByToken(token, hcalTPs);
     for (const auto &itr : *hcalTPs) {
       HcalTrigTowerDetId id = itr.id();
       double et = decoder.hcaletValue(itr.id(), itr.t0());
       if (et <= 0)
         continue;
       float towerEta = l1t::CaloTools::towerEta(id.ieta());
       float towerPhi = l1t::CaloTools::towerPhi(id.ieta(), id.iphi());
       if (!hcalDigisBarrel_ && std::abs(towerEta) < 2)  // |eta| < 2 => barrel (there's no HE in Phase2)
         continue;
       if (!hcalDigisHF_ && std::abs(towerEta) > 2)  // |eta| > 2 => HF
         continue;
       if (debug_)
         edm::LogWarning("L1TPFCaloProducer")
             << "adding HCal digi input pt " << et << ", eta " << towerEta << ", phi " << towerPhi << "\n";
       hcalClusterer_.add(et, towerEta, towerPhi);
     }
   }
 }

◆ readHcalHGCTowers_()

void L1TPFCaloProducer::readHcalHGCTowers_	(	edm::Event &	event,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 280 of file L1TPFCaloProducer.cc.

References l1tpf_calo::SingleCaloClusterer::add(), BXVector< T >::begin(), debug_, ecalClusterer_, BXVector< T >::end(), hcalClusterer_, hcalHGCTowers_, hcalHGCTowersHadOnly_, iEvent, and unpackBuffers-CaloStage2::token.

Referenced by produce().

                                                                                         {
   edm::Handle<l1t::HGCalTowerBxCollection> hgcTowers;
 
   for (const auto &token : hcalHGCTowers_) {
     iEvent.getByToken(token, hgcTowers);
     for (auto it = hgcTowers->begin(0), ed = hgcTowers->end(0); it != ed; ++it) {
       if (debug_)
         edm::LogWarning("L1TPFCaloProducer")
             << "adding HGC Tower hadEt " << it->etHad() << ", emEt " << it->etEm() << ", pt " << it->pt() << ", eta "
             << it->eta() << ", phi " << it->phi() << "\n";
       hcalClusterer_.add(it->etHad(), it->eta(), it->phi());
       if (!hcalHGCTowersHadOnly_)
         ecalClusterer_.add(it->etEm(), it->eta(), it->phi());
     }
   }
 }

◆ readPhase2BarrelCaloTowers_()

void L1TPFCaloProducer::readPhase2BarrelCaloTowers_	(	edm::Event &	event,
		const edm::EventSetup &
	)

private

Definition at line 260 of file L1TPFCaloProducer.cc.

References l1tpf_calo::SingleCaloClusterer::add(), debug_, ecalClusterer_, hcalClusterer_, createfilelist::int, phase2barrelTowers_, submitPVValidationJobs::t, unpackBuffers-CaloStage2::token, and HLT_2023v12_cff::towers.

Referenced by produce().

                                                                                           {
   edm::Handle<l1tp2::CaloTowerCollection> towers;
   for (const auto &token : phase2barrelTowers_) {
     event.getByToken(token, towers);
     for (const auto &t : *towers) {
       // sanity check from https://github.com/cms-l1t-offline/cmssw/blob/l1t-phase2-v3.0.2/L1Trigger/L1CaloTrigger/plugins/L1TowerCalibrator.cc#L259-L263
       if ((int)t.towerIEta() == -1016 && (int)t.towerIPhi() == -962)
         continue;
       if (debug_ && (t.hcalTowerEt() > 0 || t.ecalTowerEt() > 0)) {
         edm::LogWarning("L1TPFCaloProducer")
             << "adding phase2 L1 CaloTower eta " << t.towerEta() << "   phi " << t.towerPhi() << "   ieta "
             << t.towerIEta() << "   iphi " << t.towerIPhi() << "   ecal " << t.ecalTowerEt() << "    hcal "
             << t.hcalTowerEt() << "\n";
       }
       hcalClusterer_.add(t.hcalTowerEt(), t.towerEta(), t.towerPhi());
       ecalClusterer_.add(t.ecalTowerEt(), t.towerEta(), t.towerPhi());
     }
   }
 }

Member Data Documentation

◆ caloLinker_

std::unique_ptr<l1tpf_calo::SimpleCaloLinkerBase> L1TPFCaloProducer::caloLinker_

private

Definition at line 51 of file L1TPFCaloProducer.cc.

Referenced by produce().

◆ debug_

bool L1TPFCaloProducer::debug_

private

Definition at line 35 of file L1TPFCaloProducer.cc.

Referenced by produce(), readHcalDigis_(), readHcalHGCTowers_(), and readPhase2BarrelCaloTowers_().

◆ decoderTag_

edm::ESGetToken<CaloTPGTranscoder, CaloTPGRecord> L1TPFCaloProducer::decoderTag_

private

Definition at line 40 of file L1TPFCaloProducer.cc.

Referenced by readHcalDigis_().

◆ ecalCands_

std::vector<edm::EDGetTokenT<reco::CandidateView> > L1TPFCaloProducer::ecalCands_

private

Definition at line 36 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and produce().

◆ ecalClusterer_

l1tpf_calo::SingleCaloClusterer L1TPFCaloProducer::ecalClusterer_

private

Definition at line 50 of file L1TPFCaloProducer.cc.

Referenced by produce(), readHcalHGCTowers_(), and readPhase2BarrelCaloTowers_().

◆ ecalOnly_

bool L1TPFCaloProducer::ecalOnly_

private

Definition at line 35 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and produce().

◆ emCorrector_

l1tpf::corrector L1TPFCaloProducer::emCorrector_

private

Definition at line 46 of file L1TPFCaloProducer.cc.

Referenced by produce().

◆ hadCorrector_

l1tpf::corrector L1TPFCaloProducer::hadCorrector_

private

Definition at line 48 of file L1TPFCaloProducer.cc.

Referenced by produce().

◆ hcalCands_

std::vector<edm::EDGetTokenT<reco::CandidateView> > L1TPFCaloProducer::hcalCands_

private

Definition at line 37 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and produce().

◆ hcalClusterer_

l1tpf_calo::SingleCaloClusterer L1TPFCaloProducer::hcalClusterer_

private

Definition at line 50 of file L1TPFCaloProducer.cc.

Referenced by produce(), readHcalDigis_(), readHcalHGCTowers_(), and readPhase2BarrelCaloTowers_().

◆ hcalDigis_

std::vector<edm::EDGetTokenT<HcalTrigPrimDigiCollection> > L1TPFCaloProducer::hcalDigis_

private

Definition at line 39 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), produce(), and readHcalDigis_().

◆ hcalDigisBarrel_

bool L1TPFCaloProducer::hcalDigisBarrel_

private

Definition at line 41 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and readHcalDigis_().

◆ hcalDigisHF_

bool L1TPFCaloProducer::hcalDigisHF_

private

Definition at line 41 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and readHcalDigis_().

◆ hcalHGCTowers_

std::vector<edm::EDGetTokenT<l1t::HGCalTowerBxCollection> > L1TPFCaloProducer::hcalHGCTowers_

private

Definition at line 43 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), produce(), and readHcalHGCTowers_().

◆ hcalHGCTowersHadOnly_

bool L1TPFCaloProducer::hcalHGCTowersHadOnly_

private

Definition at line 44 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and readHcalHGCTowers_().

◆ hcCorrector_

l1tpf::corrector L1TPFCaloProducer::hcCorrector_

private

Definition at line 47 of file L1TPFCaloProducer.cc.

Referenced by produce().

◆ phase2barrelTowers_

std::vector<edm::EDGetTokenT<l1tp2::CaloTowerCollection> > L1TPFCaloProducer::phase2barrelTowers_

private

Definition at line 42 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), produce(), and readPhase2BarrelCaloTowers_().

◆ resol_

l1tpf::ParametricResolution L1TPFCaloProducer::resol_

private

Definition at line 53 of file L1TPFCaloProducer.cc.

Referenced by produce().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ L1TPFCaloProducer()

Member Function Documentation

◆ produce()

◆ readHcalDigis_()

◆ readHcalHGCTowers_()

◆ readPhase2BarrelCaloTowers_()

Member Data Documentation

◆ caloLinker_

◆ debug_

◆ decoderTag_

◆ ecalCands_

◆ ecalClusterer_

◆ ecalOnly_

◆ emCorrector_

◆ hadCorrector_

◆ hcalCands_

◆ hcalClusterer_

◆ hcalDigis_

◆ hcalDigisBarrel_

◆ hcalDigisHF_

◆ hcalHGCTowers_

◆ hcalHGCTowersHadOnly_

◆ hcCorrector_

◆ phase2barrelTowers_

◆ resol_