#include <L1TowerCalibrator.cc>

Inheritance diagram for L1TowerCalibrator:

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

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

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

Private Attributes
std::map< std::string, std::map< std::string, TF1 > >	all_nvtx_to_PU_sub_funcs

const double	barrelSF

const bool	debug

const edm::ESGetToken< CaloTPGTranscoder, CaloTPGRecord >	decoderTag_

std::map< std::string, TF1 >	ecal_nvtx_to_PU_sub_funcs

const double	EcalTpEtMin

std::map< std::string, TF1 >	hcal_nvtx_to_PU_sub_funcs

const edm::EDGetTokenT< HcalTrigPrimDigiCollection >	hcalToken_

edm::Handle< HcalTrigPrimDigiCollection >	hcalTowerHandle

const double	HcalTpEtMin

std::map< std::string, TF1 >	hf_nvtx_to_PU_sub_funcs

const double	hfSF

const double	HFTpEtMin

std::map< std::string, TF1 >	hgcalEM_nvtx_to_PU_sub_funcs

const double	HGCalEmTpEtMin

std::map< std::string, TF1 >	hgcalHad_nvtx_to_PU_sub_funcs

const double	HGCalHadTpEtMin

const double	hgcalSF

l1t::HGCalTowerBxCollection	hgcalTowers

edm::Handle< l1t::HGCalTowerBxCollection >	hgcalTowersHandle

const edm::EDGetTokenT< l1t::HGCalTowerBxCollection >	hgcalTowersToken_

edm::Handle< l1tp2::CaloTowerCollection >	l1CaloTowerHandle

const edm::EDGetTokenT< l1tp2::CaloTowerCollection >	l1TowerToken_

std::map< std::string, TF1 >	nHits_to_nvtx_funcs

std::vector< edm::ParameterSet >	nHits_to_nvtx_params

std::vector< edm::ParameterSet >	nvtx_to_PU_sub_params

const double	puThreshold

const double	puThresholdEcalMax

const double	puThresholdEcalMin

const double	puThresholdHcalMax

const double	puThresholdHcalMin

const double	puThresholdHFMax

const double	puThresholdHFMin

const double	puThresholdHGCalEMMax

const double	puThresholdHGCalEMMin

const double	puThresholdHGCalHadMax

const double	puThresholdHGCalHadMin

const double	puThresholdL1eg

const bool	skipCalibrations

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: Take the calibrated unclustered ECAL energy and total HCAL energy associated with the L1CaloTower collection output from L1EGammaCrystalsEmulatorProducer: l1CaloTowerCollection, "L1CaloTowerCollection"

as well as HGCal Tower level inputs: BXVector<l1t::HGCalTower> "hgcalTriggerPrimitiveDigiProducer" "tower" "HLT"
and HCAL HF inputs from: edm::SortedCollection<HcalTriggerPrimitiveDigi,edm::StrictWeakOrdering<HcalTriggerPrimitiveDigi> > "simHcalTriggerPrimitiveDigis" "" "HLT"

Implement PU-based calibrations which scale down the ET in the towers based on mapping nTowers with ECAL(HCAL) ET <= defined PU threshold. This value has been shown to be similar between TTbar, QCD, and minBias samples. This allows a prediction of nvtx. Which can be mapped to the total minBias energy in an eta slice of the detector. Subtract the total estimated minBias energy per eta slice divided by nTowers in that eta slice from each trigger tower in that eta slice.

This is all ECAL / HCAL specific or EM vs. Hadronic for HGCal.

Implementation: [Notes on implementation]

Definition at line 63 of file L1TowerCalibrator.cc.

Constructor & Destructor Documentation

◆ L1TowerCalibrator()

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

explicit

Definition at line 118 of file L1TowerCalibrator.cc.

References all_nvtx_to_PU_sub_funcs, debug, ecal_nvtx_to_PU_sub_funcs, hcal_nvtx_to_PU_sub_funcs, hf_nvtx_to_PU_sub_funcs, hgcalEM_nvtx_to_PU_sub_funcs, hgcalHad_nvtx_to_PU_sub_funcs, mps_fire::i, l1tTowerCalibrationProducer_cfi::iEta, nHits_to_nvtx_funcs, nHits_to_nvtx_params, nvtx_to_PU_sub_params, LaserDQM_cfg::p1, SiStripOfflineCRack_cfg::p2, muonDTDigis_cfi::pset, AlCaHLTBitMon_QueryRunRegistry::string, and parallelization::uint.

     : HcalTpEtMin(iConfig.getParameter<double>("HcalTpEtMin")),
       EcalTpEtMin(iConfig.getParameter<double>("EcalTpEtMin")),
       HGCalHadTpEtMin(iConfig.getParameter<double>("HGCalHadTpEtMin")),
       HGCalEmTpEtMin(iConfig.getParameter<double>("HGCalEmTpEtMin")),
       HFTpEtMin(iConfig.getParameter<double>("HFTpEtMin")),
       puThreshold(iConfig.getParameter<double>("puThreshold")),
       puThresholdL1eg(iConfig.getParameter<double>("puThresholdL1eg")),
       puThresholdHcalMin(iConfig.getParameter<double>("puThresholdHcalMin")),
       puThresholdHcalMax(iConfig.getParameter<double>("puThresholdHcalMax")),
       puThresholdEcalMin(iConfig.getParameter<double>("puThresholdEcalMin")),
       puThresholdEcalMax(iConfig.getParameter<double>("puThresholdEcalMax")),
       puThresholdHGCalEMMin(iConfig.getParameter<double>("puThresholdHGCalEMMin")),
       puThresholdHGCalEMMax(iConfig.getParameter<double>("puThresholdHGCalEMMax")),
       puThresholdHGCalHadMin(iConfig.getParameter<double>("puThresholdHGCalHadMin")),
       puThresholdHGCalHadMax(iConfig.getParameter<double>("puThresholdHGCalHadMax")),
       puThresholdHFMin(iConfig.getParameter<double>("puThresholdHFMin")),
       puThresholdHFMax(iConfig.getParameter<double>("puThresholdHFMax")),
       barrelSF(iConfig.getParameter<double>("barrelSF")),
       hgcalSF(iConfig.getParameter<double>("hgcalSF")),
       hfSF(iConfig.getParameter<double>("hfSF")),
       debug(iConfig.getParameter<bool>("debug")),
       skipCalibrations(iConfig.getParameter<bool>("skipCalibrations")),
       l1TowerToken_(consumes<l1tp2::CaloTowerCollection>(iConfig.getParameter<edm::InputTag>("l1CaloTowers"))),
       hgcalTowersToken_(
           consumes<l1t::HGCalTowerBxCollection>(iConfig.getParameter<edm::InputTag>("L1HgcalTowersInputTag"))),
       hcalToken_(consumes<HcalTrigPrimDigiCollection>(iConfig.getParameter<edm::InputTag>("hcalDigis"))),
       decoderTag_(esConsumes<CaloTPGTranscoder, CaloTPGRecord>(edm::ESInputTag("", ""))),
       nHits_to_nvtx_params(iConfig.getParameter<std::vector<edm::ParameterSet> >("nHits_to_nvtx_params")),
       nvtx_to_PU_sub_params(iConfig.getParameter<std::vector<edm::ParameterSet> >("nvtx_to_PU_sub_params")) {
   // Initialize the nHits --> nvtx functions
   for (uint i = 0; i < nHits_to_nvtx_params.size(); i++) {
     edm::ParameterSet* pset = &nHits_to_nvtx_params.at(i);
     std::string calo = pset->getParameter<std::string>("fit");
     nHits_to_nvtx_funcs[calo.c_str()] = TF1(calo.c_str(), "[0] + [1] * x");
     nHits_to_nvtx_funcs[calo.c_str()].SetParameter(0, pset->getParameter<std::vector<double> >("params").at(0));
     nHits_to_nvtx_funcs[calo.c_str()].SetParameter(1, pset->getParameter<std::vector<double> >("params").at(1));
 
     if (debug) {
       printf(
           "nHits_to_nvtx_params[%i]\n \
                 fit: %s \n \
                 p1: %f \n \
                 p2 %f \n",
           i,
           calo.c_str(),
           nHits_to_nvtx_funcs[calo.c_str()].GetParameter(0),
           nHits_to_nvtx_funcs[calo.c_str()].GetParameter(1));
     }
   }
 
   // Initialize the nvtx --> PU subtraction functions
   all_nvtx_to_PU_sub_funcs["ecal"] = ecal_nvtx_to_PU_sub_funcs;
   all_nvtx_to_PU_sub_funcs["hcal"] = hcal_nvtx_to_PU_sub_funcs;
   all_nvtx_to_PU_sub_funcs["hgcalEM"] = hgcalEM_nvtx_to_PU_sub_funcs;
   all_nvtx_to_PU_sub_funcs["hgcalHad"] = hgcalHad_nvtx_to_PU_sub_funcs;
   all_nvtx_to_PU_sub_funcs["hf"] = hf_nvtx_to_PU_sub_funcs;
 
   for (uint i = 0; i < nvtx_to_PU_sub_params.size(); i++) {
     edm::ParameterSet* pset = &nvtx_to_PU_sub_params.at(i);
     std::string calo = pset->getParameter<std::string>("calo");
     std::string iEta = pset->getParameter<std::string>("iEta");
     double p1 = pset->getParameter<std::vector<double> >("params").at(0);
     double p2 = pset->getParameter<std::vector<double> >("params").at(1);
 
     all_nvtx_to_PU_sub_funcs[calo.c_str()][iEta.c_str()] = TF1(calo.c_str(), "[0] + [1] * x");
     all_nvtx_to_PU_sub_funcs[calo.c_str()][iEta.c_str()].SetParameter(0, p1);
     all_nvtx_to_PU_sub_funcs[calo.c_str()][iEta.c_str()].SetParameter(1, p2);
 
     if (debug) {
       printf(
           "nvtx_to_PU_sub_params[%i]\n \
                 sub detector: %s \n \
                 iEta: %s \n \
                 p1: %f \n \
                 p2 %f \n",
           i,
           calo.c_str(),
           iEta.c_str(),
           p1,
           p2);
     }
   }
 
   // Our two outputs, calibrated towers and estimated nvtx for fun
   produces<l1tp2::CaloTowerCollection>("L1CaloTowerCalibratedCollection");
   produces<double>("EstimatedNvtx");
 }

Member Function Documentation

◆ produce()

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

overrideprivatevirtual

Implements edm::one::EDProducerBase.

Definition at line 207 of file L1TowerCalibrator.cc.

                                                                              {
   // Estimated number of vertices used for calibration estimattion
   std::unique_ptr<double> EstimatedNvtx(new double);
   // Calibrated output collection
   std::unique_ptr<l1tp2::CaloTowerCollection> L1CaloTowerCalibratedCollection(new l1tp2::CaloTowerCollection);
 
   // Load the ECAL+HCAL tower sums coming from L1EGammaCrystalsEmulatorProducer.cc
   iEvent.getByToken(l1TowerToken_, l1CaloTowerHandle);
 
   // HGCal info
   iEvent.getByToken(hgcalTowersToken_, hgcalTowersHandle);
   hgcalTowers = (*hgcalTowersHandle.product());
 
   // HF Tower info
   iEvent.getByToken(hcalToken_, hcalTowerHandle);
 
   // Barrel ECAL (unclustered) and HCAL
   for (auto& hit : *l1CaloTowerHandle.product()) {
     l1tp2::CaloTower l1Hit;
     l1Hit.setEcalTowerEt(hit.ecalTowerEt());
     l1Hit.setHcalTowerEt(hit.hcalTowerEt());
     l1Hit.setL1egTowerEt(hit.l1egTowerEt());
     // Add min ET thresholds for tower ET
     if (l1Hit.ecalTowerEt() < EcalTpEtMin)
       l1Hit.setEcalTowerEt(0.0);
     if (l1Hit.hcalTowerEt() < HcalTpEtMin)
       l1Hit.setHcalTowerEt(0.0);
     l1Hit.setTowerIEta(hit.towerIEta());
     l1Hit.setTowerIPhi(hit.towerIPhi());
     l1Hit.setTowerEta(hit.towerEta());
     l1Hit.setTowerPhi(hit.towerPhi());
     l1Hit.setIsBarrel(hit.isBarrel());
     l1Hit.setNL1eg(hit.nL1eg());
     l1Hit.setL1egTrkSS(hit.l1egTrkSS());
     l1Hit.setL1egTrkIso(hit.l1egTrkIso());
     l1Hit.setL1egStandaloneSS(hit.l1egStandaloneSS());
     l1Hit.setL1egStandaloneIso(hit.l1egStandaloneIso());
 
     // FIXME There is an error in the L1EGammaCrystalsEmulatorProducer.cc which is
     // returning towers with minimal ECAL energy, and no HCAL energy with these
     // iEta/iPhi coordinates and eta = -88.653152 and phi = -99.000000.
     // Skip these for the time being until the upstream code has been debugged
     if ((int)l1Hit.towerIEta() == -1016 && (int)l1Hit.towerIPhi() == -962)
       continue;
 
     (*L1CaloTowerCalibratedCollection).push_back(l1Hit);
     if (debug)
       printf("Barrel tower iEta %i iPhi %i eta %f phi %f ecal_et %f hcal_et_sum %f\n",
              (int)l1Hit.towerIEta(),
              (int)l1Hit.towerIPhi(),
              l1Hit.towerEta(),
              l1Hit.towerPhi(),
              l1Hit.ecalTowerEt(),
              l1Hit.hcalTowerEt());
   }
 
   // Loop over HGCalTowers and create L1CaloTowers for them and add to collection
   // This iterator is taken from the PF P2 group
   // https://github.com/p2l1pfp/cmssw/blob/170808db68038d53794bc65fdc962f8fc337a24d/L1Trigger/Phase2L1ParticleFlow/plugins/L1TPFCaloProducer.cc#L278-L289
   for (auto it = hgcalTowers.begin(0), ed = hgcalTowers.end(0); it != ed; ++it) {
     // skip lowest ET towers
     if (it->etEm() < HGCalEmTpEtMin && it->etHad() < HGCalHadTpEtMin)
       continue;
 
     l1tp2::CaloTower l1Hit;
     // Set energies normally, but need to zero if below threshold
     if (it->etEm() < HGCalEmTpEtMin)
       l1Hit.setEcalTowerEt(0.);
     else
       l1Hit.setEcalTowerEt(it->etEm());
 
     if (it->etHad() < HGCalHadTpEtMin)
       l1Hit.setHcalTowerEt(0.);
     else
       l1Hit.setHcalTowerEt(it->etHad());
 
     l1Hit.setTowerEta(it->eta());
     l1Hit.setTowerPhi(it->phi());
     l1Hit.setTowerIEta(-98);  // -98 mean HGCal
     l1Hit.setTowerIPhi(-98);
     l1Hit.setIsBarrel(false);
     (*L1CaloTowerCalibratedCollection).push_back(l1Hit);
     if (debug)
       printf("HGCal tower iEta %i iPhi %i eta %f phi %f ecal_et %f hcal_et_sum %f\n",
              (int)l1Hit.towerIEta(),
              (int)l1Hit.towerIPhi(),
              l1Hit.towerEta(),
              l1Hit.towerPhi(),
              l1Hit.ecalTowerEt(),
              l1Hit.hcalTowerEt());
   }
 
   // Loop over Hcal HF tower inputs and create L1CaloTowers and add to
   // L1CaloTowerCalibratedCollection collection
   const auto& decoder = iSetup.getData(decoderTag_);
   for (const auto& hit : *hcalTowerHandle.product()) {
     HcalTrigTowerDetId id = hit.id();
     double et = decoder.hcaletValue(hit.id(), hit.t0());
     if (et < HFTpEtMin)
       continue;
     // Only doing HF so skip outside range
     if (abs(id.ieta()) < l1t::CaloTools::kHFBegin)
       continue;
     if (abs(id.ieta()) > l1t::CaloTools::kHFEnd)
       continue;
 
     l1tp2::CaloTower l1Hit;
     l1Hit.setEcalTowerEt(0.);
     l1Hit.setHcalTowerEt(et);
     l1Hit.setTowerEta(l1t::CaloTools::towerEta(id.ieta()));
     l1Hit.setTowerPhi(l1t::CaloTools::towerPhi(id.ieta(), id.iphi()));
     l1Hit.setTowerIEta(id.ieta());
     l1Hit.setTowerIPhi(id.iphi());
     l1Hit.setIsBarrel(false);
     (*L1CaloTowerCalibratedCollection).push_back(l1Hit);
 
     if (debug)
       printf("HCAL HF tower iEta %i iPhi %i eta %f phi %f ecal_et %f hcal_et_sum %f\n",
              (int)l1Hit.towerIEta(),
              (int)l1Hit.towerIPhi(),
              l1Hit.towerEta(),
              l1Hit.towerPhi(),
              l1Hit.ecalTowerEt(),
              l1Hit.hcalTowerEt());
   }
 
   // N Tower totals
   // For mapping to estimated nvtx in event
   int i_ecal_hits_leq_threshold = 0;
   int i_hgcalEM_hits_leq_threshold = 0;
   int i_hcal_hits_leq_threshold = 0;
   int i_hgcalHad_hits_leq_threshold = 0;
   int i_hf_hits_leq_threshold = 0;
 
   // Loop over the collection containing all hits
   // and calculate the number of hits falling into the
   // "less than or equal" nTowers variable which maps to
   // estimated number of vertices
   for (auto& l1CaloTower : (*L1CaloTowerCalibratedCollection)) {
     // Barrel ECAL
     if (l1CaloTower.ecalTowerEt() > 0. && l1CaloTower.towerIEta() != -98) {
       if (l1CaloTower.ecalTowerEt() <= puThresholdEcalMax && l1CaloTower.ecalTowerEt() >= puThresholdEcalMin) {
         i_ecal_hits_leq_threshold++;
       }
     }
 
     // HGCal EM
     if (l1CaloTower.ecalTowerEt() > 0. && l1CaloTower.towerIEta() == -98) {
       if (l1CaloTower.ecalTowerEt() <= puThresholdHGCalEMMax && l1CaloTower.ecalTowerEt() >= puThresholdHGCalEMMin) {
         i_hgcalEM_hits_leq_threshold++;
       }
     }
 
     // Barrel HCAL
     if (l1CaloTower.hcalTowerEt() > 0. && l1CaloTower.towerIEta() != -98 &&
         abs(l1CaloTower.towerEta()) < 2.0)  // abs(eta) < 2 just keeps us out of HF
     {
       if (l1CaloTower.hcalTowerEt() <= puThresholdHcalMax && l1CaloTower.hcalTowerEt() >= puThresholdHcalMin) {
         i_hcal_hits_leq_threshold++;
       }
     }
 
     // HGCal Had
     if (l1CaloTower.hcalTowerEt() > 0. && l1CaloTower.towerIEta() == -98) {
       if (l1CaloTower.hcalTowerEt() <= puThresholdHGCalHadMax && l1CaloTower.hcalTowerEt() >= puThresholdHGCalHadMin) {
         i_hgcalHad_hits_leq_threshold++;
       }
     }
 
     // HF
     if (l1CaloTower.hcalTowerEt() > 0. && l1CaloTower.towerIEta() != -98 &&
         abs(l1CaloTower.towerEta()) > 2.0)  // abs(eta) > 2 keeps us out of barrel HF
     {
       if (l1CaloTower.hcalTowerEt() <= puThresholdHFMax && l1CaloTower.hcalTowerEt() >= puThresholdHFMin) {
         i_hf_hits_leq_threshold++;
       }
     }
   }
 
   // For each subdetector, map to the estimated number of vertices
   double ecal_nvtx = nHits_to_nvtx_funcs["ecal"].Eval(i_ecal_hits_leq_threshold);
   double hcal_nvtx = nHits_to_nvtx_funcs["hcal"].Eval(i_hcal_hits_leq_threshold);
   double hgcalEM_nvtx = nHits_to_nvtx_funcs["hgcalEM"].Eval(i_hgcalEM_hits_leq_threshold);
   double hgcalHad_nvtx = nHits_to_nvtx_funcs["hgcalHad"].Eval(i_hgcalHad_hits_leq_threshold);
   double hf_nvtx = nHits_to_nvtx_funcs["hf"].Eval(i_hf_hits_leq_threshold);
   // Make sure all values are >= 0
   if (ecal_nvtx < 0)
     ecal_nvtx = 0;
   if (hcal_nvtx < 0)
     hcal_nvtx = 0;
   if (hgcalEM_nvtx < 0)
     hgcalEM_nvtx = 0;
   if (hgcalHad_nvtx < 0)
     hgcalHad_nvtx = 0;
   if (hf_nvtx < 0)
     hf_nvtx = 0;
   // Best estimate is avg of all except HF.
   // This is b/c HF currently has such poor prediction power, it only degrades the avg result
   // NEW, with 10_3_X, hgcal and HF has the best results based on the values I took...
   // skip ECAL and HCAL
   //*EstimatedNvtx = ( ecal_nvtx + hcal_nvtx + hgcalEM_nvtx + hgcalHad_nvtx + hf_nvtx ) / 3.;
   *EstimatedNvtx = (hgcalEM_nvtx + hgcalHad_nvtx + hf_nvtx) / 3.;
 
   if (debug) {
     double lumi = iEvent.eventAuxiliary().luminosityBlock();
     double event = iEvent.eventAuxiliary().event();
 
     printf(
         "L1TowerCalibrater: lumi %.0f evt %.0f nTowers for subdetecters \
             \nECAL:      %i --> nvtx = %.1f \
             \nHGCal EM:  %i --> nvtx = %.1f \
             \nHCAL:      %i --> nvtx = %.1f \
             \nHGCal Had: %i --> nvtx = %.1f \
             \nHCAL HF:   %i --> nvtx = %.1f \
             \nEstimated Nvtx = %.1f\n",
         lumi,
         event,
         i_ecal_hits_leq_threshold,
         ecal_nvtx,
         i_hgcalEM_hits_leq_threshold,
         hgcalEM_nvtx,
         i_hcal_hits_leq_threshold,
         hcal_nvtx,
         i_hgcalHad_hits_leq_threshold,
         hgcalHad_nvtx,
         i_hf_hits_leq_threshold,
         hf_nvtx,
         *EstimatedNvtx);
   }
 
   // Use estimated number of vertices to subtract off PU contributions
   // to each and every hit. In cases where the energy would go negative,
   // limit this to zero.
   if (!skipCalibrations)  // skipCalibrations simply passes the towers through
   {
     for (auto& l1CaloTower : (*L1CaloTowerCalibratedCollection)) {
       // Barrel ECAL eta slices
       if (l1CaloTower.ecalTowerEt() > 0. && l1CaloTower.towerIEta() != -98) {
         if (abs(l1CaloTower.towerIEta()) <= 3) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["ecal"]["er1to3"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 6 && abs(l1CaloTower.towerIEta()) >= 4) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["ecal"]["er4to6"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 9 && abs(l1CaloTower.towerIEta()) >= 7) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["ecal"]["er7to9"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 12 && abs(l1CaloTower.towerIEta()) >= 10) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["ecal"]["er10to12"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 15 && abs(l1CaloTower.towerIEta()) >= 13) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["ecal"]["er13to15"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 18 && abs(l1CaloTower.towerIEta()) >= 16) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["ecal"]["er16to18"].Eval(*EstimatedNvtx) * barrelSF);
         }
       }
 
       // HGCal EM eta slices
       if (l1CaloTower.ecalTowerEt() > 0. && l1CaloTower.towerIEta() == -98) {
         if (abs(l1CaloTower.towerEta()) <= 1.8) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalEM"]["er1p4to1p8"].Eval(*EstimatedNvtx) * hgcalSF);
         }
         if (abs(l1CaloTower.towerEta()) <= 2.1 && abs(l1CaloTower.towerEta()) > 1.8) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalEM"]["er1p8to2p1"].Eval(*EstimatedNvtx) * hgcalSF);
         }
         if (abs(l1CaloTower.towerEta()) <= 2.4 && abs(l1CaloTower.towerEta()) > 2.1) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalEM"]["er2p1to2p4"].Eval(*EstimatedNvtx) * hgcalSF);
         }
         if (abs(l1CaloTower.towerEta()) <= 2.7 && abs(l1CaloTower.towerEta()) > 2.4) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalEM"]["er2p4to2p7"].Eval(*EstimatedNvtx) * hgcalSF);
         }
         if (abs(l1CaloTower.towerEta()) <= 3.1 && abs(l1CaloTower.towerEta()) > 2.7) {
           l1CaloTower.setEcalTowerEt(l1CaloTower.ecalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalEM"]["er2p7to3p1"].Eval(*EstimatedNvtx) * hgcalSF);
         }
       }
 
       // Barrel HCAL eta slices
       if (l1CaloTower.hcalTowerEt() > 0. && l1CaloTower.towerIEta() != -98 &&
           abs(l1CaloTower.towerEta()) < 2.0)  // abs(eta) < 2 just keeps us out of HF
       {
         if (abs(l1CaloTower.towerIEta()) <= 3) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hcal"]["er1to3"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 6 && abs(l1CaloTower.towerIEta()) >= 4) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hcal"]["er4to6"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 9 && abs(l1CaloTower.towerIEta()) >= 7) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hcal"]["er7to9"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 12 && abs(l1CaloTower.towerIEta()) >= 10) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hcal"]["er10to12"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 15 && abs(l1CaloTower.towerIEta()) >= 13) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hcal"]["er13to15"].Eval(*EstimatedNvtx) * barrelSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 18 && abs(l1CaloTower.towerIEta()) >= 16) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hcal"]["er16to18"].Eval(*EstimatedNvtx) * barrelSF);
         }
       }
 
       // HGCal Had eta slices
       if (l1CaloTower.hcalTowerEt() > 0. && l1CaloTower.towerIEta() == -98) {
         if (abs(l1CaloTower.towerEta()) <= 1.8) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalHad"]["er1p4to1p8"].Eval(*EstimatedNvtx) * hgcalSF);
         }
         if (abs(l1CaloTower.towerEta()) <= 2.1 && abs(l1CaloTower.towerEta()) > 1.8) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalHad"]["er1p8to2p1"].Eval(*EstimatedNvtx) * hgcalSF);
         }
         if (abs(l1CaloTower.towerEta()) <= 2.4 && abs(l1CaloTower.towerEta()) > 2.1) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalHad"]["er2p1to2p4"].Eval(*EstimatedNvtx) * hgcalSF);
         }
         if (abs(l1CaloTower.towerEta()) <= 2.7 && abs(l1CaloTower.towerEta()) > 2.4) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalHad"]["er2p4to2p7"].Eval(*EstimatedNvtx) * hgcalSF);
         }
         if (abs(l1CaloTower.towerEta()) <= 3.1 && abs(l1CaloTower.towerEta()) > 2.7) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hgcalHad"]["er2p7to3p1"].Eval(*EstimatedNvtx) * hgcalSF);
         }
       }
 
       // HF eta slices
       if (l1CaloTower.hcalTowerEt() > 0. && l1CaloTower.towerIEta() != -98 &&
           abs(l1CaloTower.towerEta()) > 2.0)  // abs(eta) > 2 keeps us out of barrel HF
       {
         if (abs(l1CaloTower.towerIEta()) <= 33 && abs(l1CaloTower.towerIEta()) >= 29) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hf"]["er29to33"].Eval(*EstimatedNvtx) * hfSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 37 && abs(l1CaloTower.towerIEta()) >= 34) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hf"]["er34to37"].Eval(*EstimatedNvtx) * hfSF);
         }
         if (abs(l1CaloTower.towerIEta()) <= 41 && abs(l1CaloTower.towerIEta()) >= 38) {
           l1CaloTower.setHcalTowerEt(l1CaloTower.hcalTowerEt() -
                                      all_nvtx_to_PU_sub_funcs["hf"]["er38to41"].Eval(*EstimatedNvtx) * hfSF);
         }
       }
 
       // Make sure none are negative
       if (l1CaloTower.ecalTowerEt() < 0.)
         l1CaloTower.setEcalTowerEt(0.);
       if (l1CaloTower.hcalTowerEt() < 0.)
         l1CaloTower.setHcalTowerEt(0.);
     }
   }
 
   iEvent.put(std::move(EstimatedNvtx), "EstimatedNvtx");
   iEvent.put(std::move(L1CaloTowerCalibratedCollection), "L1CaloTowerCalibratedCollection");
 }

Member Data Documentation

◆ all_nvtx_to_PU_sub_funcs

std::map<std::string, std::map<std::string, TF1> > L1TowerCalibrator::all_nvtx_to_PU_sub_funcs

private

Definition at line 115 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator(), and produce().

◆ barrelSF

const double L1TowerCalibrator::barrelSF

private

Definition at line 87 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ debug

const bool L1TowerCalibrator::debug

private

Definition at line 90 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), produce(), and runTauIdMVA.TauIDEmbedder::runTauID().

◆ decoderTag_

const edm::ESGetToken<CaloTPGTranscoder, CaloTPGRecord> L1TowerCalibrator::decoderTag_

private

Definition at line 102 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ ecal_nvtx_to_PU_sub_funcs

std::map<std::string, TF1> L1TowerCalibrator::ecal_nvtx_to_PU_sub_funcs

private

Definition at line 110 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator().

◆ EcalTpEtMin

const double L1TowerCalibrator::EcalTpEtMin

private

Definition at line 71 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hcal_nvtx_to_PU_sub_funcs

std::map<std::string, TF1> L1TowerCalibrator::hcal_nvtx_to_PU_sub_funcs

private

Definition at line 111 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator().

◆ hcalToken_

const edm::EDGetTokenT<HcalTrigPrimDigiCollection> L1TowerCalibrator::hcalToken_

private

Definition at line 100 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hcalTowerHandle

edm::Handle<HcalTrigPrimDigiCollection> L1TowerCalibrator::hcalTowerHandle

private

Definition at line 101 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ HcalTpEtMin

const double L1TowerCalibrator::HcalTpEtMin

private

Definition at line 70 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hf_nvtx_to_PU_sub_funcs

std::map<std::string, TF1> L1TowerCalibrator::hf_nvtx_to_PU_sub_funcs

private

Definition at line 114 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator().

◆ hfSF

const double L1TowerCalibrator::hfSF

private

Definition at line 89 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ HFTpEtMin

const double L1TowerCalibrator::HFTpEtMin

private

Definition at line 74 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hgcalEM_nvtx_to_PU_sub_funcs

std::map<std::string, TF1> L1TowerCalibrator::hgcalEM_nvtx_to_PU_sub_funcs

private

Definition at line 112 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator().

◆ HGCalEmTpEtMin

const double L1TowerCalibrator::HGCalEmTpEtMin

private

Definition at line 73 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hgcalHad_nvtx_to_PU_sub_funcs

std::map<std::string, TF1> L1TowerCalibrator::hgcalHad_nvtx_to_PU_sub_funcs

private

Definition at line 113 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator().

◆ HGCalHadTpEtMin

const double L1TowerCalibrator::HGCalHadTpEtMin

private

Definition at line 72 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hgcalSF

const double L1TowerCalibrator::hgcalSF

private

Definition at line 88 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hgcalTowers

l1t::HGCalTowerBxCollection L1TowerCalibrator::hgcalTowers

private

Definition at line 98 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hgcalTowersHandle

edm::Handle<l1t::HGCalTowerBxCollection> L1TowerCalibrator::hgcalTowersHandle

private

Definition at line 97 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hgcalTowersToken_

const edm::EDGetTokenT<l1t::HGCalTowerBxCollection> L1TowerCalibrator::hgcalTowersToken_

private

Definition at line 96 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ l1CaloTowerHandle

edm::Handle<l1tp2::CaloTowerCollection> L1TowerCalibrator::l1CaloTowerHandle

private

Definition at line 94 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ l1TowerToken_

const edm::EDGetTokenT<l1tp2::CaloTowerCollection> L1TowerCalibrator::l1TowerToken_

private

Definition at line 93 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ nHits_to_nvtx_funcs

std::map<std::string, TF1> L1TowerCalibrator::nHits_to_nvtx_funcs

private

Definition at line 106 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator(), and produce().

◆ nHits_to_nvtx_params

std::vector<edm::ParameterSet> L1TowerCalibrator::nHits_to_nvtx_params

private

Definition at line 105 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator().

◆ nvtx_to_PU_sub_params

std::vector<edm::ParameterSet> L1TowerCalibrator::nvtx_to_PU_sub_params

private

Definition at line 109 of file L1TowerCalibrator.cc.

Referenced by L1TowerCalibrator().

◆ puThreshold

const double L1TowerCalibrator::puThreshold

private

Definition at line 75 of file L1TowerCalibrator.cc.

◆ puThresholdEcalMax

const double L1TowerCalibrator::puThresholdEcalMax

private

Definition at line 80 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdEcalMin

const double L1TowerCalibrator::puThresholdEcalMin

private

Definition at line 79 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHcalMax

const double L1TowerCalibrator::puThresholdHcalMax

private

Definition at line 78 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHcalMin

const double L1TowerCalibrator::puThresholdHcalMin

private

Definition at line 77 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHFMax

const double L1TowerCalibrator::puThresholdHFMax

private

Definition at line 86 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHFMin

const double L1TowerCalibrator::puThresholdHFMin

private

Definition at line 85 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHGCalEMMax

const double L1TowerCalibrator::puThresholdHGCalEMMax

private

Definition at line 82 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHGCalEMMin

const double L1TowerCalibrator::puThresholdHGCalEMMin

private

Definition at line 81 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHGCalHadMax

const double L1TowerCalibrator::puThresholdHGCalHadMax

private

Definition at line 84 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHGCalHadMin

const double L1TowerCalibrator::puThresholdHGCalHadMin

private

Definition at line 83 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdL1eg

const double L1TowerCalibrator::puThresholdL1eg

private

Definition at line 76 of file L1TowerCalibrator.cc.

◆ skipCalibrations

const bool L1TowerCalibrator::skipCalibrations

private

Definition at line 91 of file L1TowerCalibrator.cc.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ L1TowerCalibrator()

Member Function Documentation

◆ produce()

Member Data Documentation

◆ all_nvtx_to_PU_sub_funcs

◆ barrelSF

◆ debug

◆ decoderTag_

◆ ecal_nvtx_to_PU_sub_funcs

◆ EcalTpEtMin

◆ hcal_nvtx_to_PU_sub_funcs

◆ hcalToken_

◆ hcalTowerHandle

◆ HcalTpEtMin

◆ hf_nvtx_to_PU_sub_funcs

◆ hfSF

◆ HFTpEtMin

◆ hgcalEM_nvtx_to_PU_sub_funcs

◆ HGCalEmTpEtMin

◆ hgcalHad_nvtx_to_PU_sub_funcs

◆ HGCalHadTpEtMin

◆ hgcalSF

◆ hgcalTowers

◆ hgcalTowersHandle

◆ hgcalTowersToken_

◆ l1CaloTowerHandle

◆ l1TowerToken_

◆ nHits_to_nvtx_funcs

◆ nHits_to_nvtx_params

◆ nvtx_to_PU_sub_params

◆ puThreshold

◆ puThresholdEcalMax

◆ puThresholdEcalMin

◆ puThresholdHcalMax

◆ puThresholdHcalMin

◆ puThresholdHFMax

◆ puThresholdHFMin

◆ puThresholdHGCalEMMax

◆ puThresholdHGCalEMMin

◆ puThresholdHGCalHadMax

◆ puThresholdHGCalHadMin

◆ puThresholdL1eg

◆ skipCalibrations