#include <L1TowerCalibrator.cc>

Inheritance diagram for L1TowerCalibrator:

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

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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 &	operator= (EDConsumerBase &&)=default

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

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

double	barrelSF

bool	debug

edm::ESGetToken< CaloTPGTranscoder, CaloTPGRecord >	decoderTag_

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

double	EcalTpEtMin

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

edm::EDGetTokenT< HcalTrigPrimDigiCollection >	hcalToken_

edm::Handle< HcalTrigPrimDigiCollection >	hcalTowerHandle

double	HcalTpEtMin

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

double	hfSF

double	HFTpEtMin

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

double	HGCalEmTpEtMin

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

double	HGCalHadTpEtMin

double	hgcalSF

l1t::HGCalTowerBxCollection	hgcalTowers

edm::Handle< l1t::HGCalTowerBxCollection >	hgcalTowersHandle

edm::EDGetTokenT< l1t::HGCalTowerBxCollection >	hgcalTowersToken_

edm::Handle< l1tp2::CaloTowerCollection >	l1CaloTowerHandle

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

double	puThreshold

double	puThresholdEcalMax

double	puThresholdEcalMin

double	puThresholdHcalMax

double	puThresholdHcalMin

double	puThresholdHFMax

double	puThresholdHFMin

double	puThresholdHGCalEMMax

double	puThresholdHGCalEMMin

double	puThresholdHGCalHadMax

double	puThresholdHGCalHadMin

double	puThresholdL1eg

bool	skipCalibrations

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
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::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

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

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

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)

template<Transition B>
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 (std::string instanceName)

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

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

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

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

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

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

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

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

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

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.

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

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, L1TowerCalibrationProducer_cfi::iEta, nHits_to_nvtx_funcs, nHits_to_nvtx_params, nvtx_to_PU_sub_params, p1, p2, muonDTDigis_cfi::pset, AlCaHLTBitMon_QueryRunRegistry::string, and parallelization::uint.

Member Function Documentation

◆ produce()

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

overrideprivatevirtual

Implements edm::EDProducer.

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

double L1TowerCalibrator::barrelSF

private

Definition at line 87 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ debug

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_

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

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_

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

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

double L1TowerCalibrator::hfSF

private

Definition at line 89 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ HFTpEtMin

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

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

double L1TowerCalibrator::HGCalHadTpEtMin

private

Definition at line 72 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ hgcalSF

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_

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_

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

double L1TowerCalibrator::puThreshold

private

Definition at line 75 of file L1TowerCalibrator.cc.

◆ puThresholdEcalMax

double L1TowerCalibrator::puThresholdEcalMax

private

Definition at line 80 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdEcalMin

double L1TowerCalibrator::puThresholdEcalMin

private

Definition at line 79 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHcalMax

double L1TowerCalibrator::puThresholdHcalMax

private

Definition at line 78 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHcalMin

double L1TowerCalibrator::puThresholdHcalMin

private

Definition at line 77 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHFMax

double L1TowerCalibrator::puThresholdHFMax

private

Definition at line 86 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHFMin

double L1TowerCalibrator::puThresholdHFMin

private

Definition at line 85 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHGCalEMMax

double L1TowerCalibrator::puThresholdHGCalEMMax

private

Definition at line 82 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHGCalEMMin

double L1TowerCalibrator::puThresholdHGCalEMMin

private

Definition at line 81 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHGCalHadMax

double L1TowerCalibrator::puThresholdHGCalHadMax

private

Definition at line 84 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdHGCalHadMin

double L1TowerCalibrator::puThresholdHGCalHadMin

private

Definition at line 83 of file L1TowerCalibrator.cc.

Referenced by produce().

◆ puThresholdL1eg

double L1TowerCalibrator::puThresholdL1eg

private

Definition at line 76 of file L1TowerCalibrator.cc.

◆ skipCalibrations

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