L1TPFCaloProducer Class Reference

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 31 of file L1TPFCaloProducer.cc.

Constructor & Destructor Documentation

L1TPFCaloProducer()

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

explicit

Definition at line 70 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 129 of file L1TPFCaloProducer.cc.

References funct::abs(), l1tpf_calo::SingleCaloClusterer::add(), 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 237 of file L1TPFCaloProducer.cc.

References funct::abs(), l1tpf_calo::SingleCaloClusterer::add(), debug_, decoderTag_, EgHLTOffHistBins_cfi::et, edm::EventSetup::getData(), hcalClusterer_, hcalDigis_, hcalDigisBarrel_, hcalDigisHF_, LEDCalibrationChannels::ieta, iEvent, LEDCalibrationChannels::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 281 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 261 of file L1TPFCaloProducer.cc.

References l1tpf_calo::SingleCaloClusterer::add(), debug_, ecalClusterer_, hcalClusterer_, createfilelist::int, phase2barrelTowers_, submitPVValidationJobs::t, unpackBuffers-CaloStage2::token, and HLT_2022v12_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 52 of file L1TPFCaloProducer.cc.

Referenced by produce().

debug_

bool L1TPFCaloProducer::debug_

private

Definition at line 36 of file L1TPFCaloProducer.cc.

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

decoderTag_

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

private

Definition at line 41 of file L1TPFCaloProducer.cc.

Referenced by readHcalDigis_().

ecalCands_

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

private

Definition at line 37 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and produce().

ecalClusterer_

l1tpf_calo::SingleCaloClusterer L1TPFCaloProducer::ecalClusterer_

private

Definition at line 51 of file L1TPFCaloProducer.cc.

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

ecalOnly_

bool L1TPFCaloProducer::ecalOnly_

private

Definition at line 36 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and produce().

emCorrector_

l1tpf::corrector L1TPFCaloProducer::emCorrector_

private

Definition at line 47 of file L1TPFCaloProducer.cc.

Referenced by produce().

hadCorrector_

l1tpf::corrector L1TPFCaloProducer::hadCorrector_

private

Definition at line 49 of file L1TPFCaloProducer.cc.

Referenced by produce().

hcalCands_

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

private

Definition at line 38 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and produce().

hcalClusterer_

l1tpf_calo::SingleCaloClusterer L1TPFCaloProducer::hcalClusterer_

private

Definition at line 51 of file L1TPFCaloProducer.cc.

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

hcalDigis_

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

private

Definition at line 40 of file L1TPFCaloProducer.cc.

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

hcalDigisBarrel_

bool L1TPFCaloProducer::hcalDigisBarrel_

private

Definition at line 42 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and readHcalDigis_().

hcalDigisHF_

bool L1TPFCaloProducer::hcalDigisHF_

private

Definition at line 42 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and readHcalDigis_().

hcalHGCTowers_

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

private

Definition at line 44 of file L1TPFCaloProducer.cc.

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

hcalHGCTowersHadOnly_

bool L1TPFCaloProducer::hcalHGCTowersHadOnly_

private

Definition at line 45 of file L1TPFCaloProducer.cc.

Referenced by L1TPFCaloProducer(), and readHcalHGCTowers_().

hcCorrector_

l1tpf::corrector L1TPFCaloProducer::hcCorrector_

private

Definition at line 48 of file L1TPFCaloProducer.cc.

Referenced by produce().

phase2barrelTowers_

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

private

Definition at line 43 of file L1TPFCaloProducer.cc.

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

resol_

l1tpf::ParametricResolution L1TPFCaloProducer::resol_

private

Definition at line 54 of file L1TPFCaloProducer.cc.

Referenced by produce().