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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	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.

    : 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");
}

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

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.

                                                                             {
  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();
}

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().

readHcalDigis_()

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

private

Definition at line 237 of file L1TPFCaloProducer.cc.

                                                                                    {
  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);
    }
  }
}

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().

readHcalHGCTowers_()

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

private

Definition at line 281 of file L1TPFCaloProducer.cc.

                                                                                        {
  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());
    }
  }
}

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

Referenced by produce().

readPhase2BarrelCaloTowers_()

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

private

Definition at line 261 of file L1TPFCaloProducer.cc.

                                                                                          {
  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());
    }
  }
}

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

Referenced by produce().

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().