#include <L1TRate.h>

Inheritance diagram for L1TRate:

Public Member Functions
	L1TRate (const edm::ParameterSet &ps)

	~L1TRate () override

Public Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMOneEDAnalyzer ()

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	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

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

std::vector< ESResolverIndex > 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

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	beginLuminosityBlock (edm::LuminosityBlock const &lumiBlock, edm::EventSetup const &c) override

void	bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &, const edm::EventSetup &) override

void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

void	endLuminosityBlock (edm::LuminosityBlock const &lumiBlock, edm::EventSetup const &c) override

Protected Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

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

Private Member Functions
bool	getXSexFitsOMDS (const edm::ParameterSet &ps)

bool	getXSexFitsPython (const edm::ParameterSet &ps)

Private Attributes
std::map< TString, int >	m_algoBit

MonitorElement *	m_ErrorMonitor

L1TMenuHelper::Tokens	m_helperTokens

std::map< std::string, bool >	m_inputCategories

edm::EDGetTokenT< L1GlobalTriggerReadoutRecord >	m_l1GtDataDaqInputTag

const edm::ESGetToken< L1GtPrescaleFactors, L1GtPrescaleFactorsAlgoTrigRcd >	m_l1GtPfAlgoToken

L1GtUtils	m_l1GtUtils

const std::vector< std::vector< int > > *	m_listsPrescaleFactors

std::map< int, double >	m_lsLuminosity

std::map< int, int >	m_lsPrescaleIndex

std::map< int, std::map< TString, double > >	m_lsRates

int	m_lsShiftGTRates

int	m_maxNbins

const edm::ESGetToken< L1GtTriggerMenu, L1GtTriggerMenuRcd >	m_menuToken

std::string	m_outputFile

edm::ParameterSet	m_parameters

int	m_refPrescaleSet

edm::EDGetTokenT< LumiScalersCollection >	m_scalersSource_colLScal

edm::EDGetTokenT< Level1TriggerScalersCollection >	m_scalersSource_triggerScalers

std::map< std::string, std::string >	m_selectedTriggers

std::map< TString, TF1 * >	m_templateFunctions

bool	m_verbose

std::map< TString, MonitorElement * >	m_xSecObservedToExpected

std::map< TString, MonitorElement * >	m_xSecVsInstLumi

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

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 Attributes inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 47 of file L1TRate.h.

Constructor & Destructor Documentation

◆ L1TRate()

L1TRate::L1TRate ( const edm::ParameterSet & ps )

Definition at line 32 of file L1TRate.cc.

References L1TSync_cfi::Categories, gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_inputCategories, m_l1GtDataDaqInputTag, m_lsShiftGTRates, m_maxNbins, m_outputFile, m_parameters, m_refPrescaleSet, m_scalersSource_colLScal, m_scalersSource_triggerScalers, and m_verbose.

     : m_menuToken(esConsumes<edm::Transition::BeginRun>()),
       m_l1GtPfAlgoToken(esConsumes<edm::Transition::BeginRun>()),
       m_helperTokens(L1TMenuHelper::consumes<edm::Transition::BeginRun>(consumesCollector())),
       m_l1GtUtils(ps, consumesCollector(), false, *this) {
   m_maxNbins = 2500;  // Maximum LS for each run (for binning purposes)
   m_parameters = ps;
 
   // Mapping parameter input variables
   m_scalersSource_colLScal =
       consumes<LumiScalersCollection>(m_parameters.getParameter<InputTag>("inputTagScalersResults"));
   m_scalersSource_triggerScalers =
       consumes<Level1TriggerScalersCollection>(m_parameters.getParameter<InputTag>("inputTagScalersResults"));
   m_l1GtDataDaqInputTag =
       consumes<L1GlobalTriggerReadoutRecord>(m_parameters.getParameter<InputTag>("inputTagL1GtDataDaq"));
   m_verbose = m_parameters.getUntrackedParameter<bool>("verbose", false);
   m_refPrescaleSet = m_parameters.getParameter<int>("refPrescaleSet");
   m_lsShiftGTRates = m_parameters.getUntrackedParameter<int>("lsShiftGTRates", 0);
 
   // Getting which categories to monitor
   ParameterSet Categories = ps.getParameter<ParameterSet>("categories");
   m_inputCategories["Mu"] = Categories.getUntrackedParameter<bool>("Mu");
   m_inputCategories["EG"] = Categories.getUntrackedParameter<bool>("EG");
   m_inputCategories["IsoEG"] = Categories.getUntrackedParameter<bool>("IsoEG");
   m_inputCategories["Jet"] = Categories.getUntrackedParameter<bool>("Jet");
   m_inputCategories["CenJet"] = Categories.getUntrackedParameter<bool>("CenJet");
   m_inputCategories["ForJet"] = Categories.getUntrackedParameter<bool>("ForJet");
   m_inputCategories["TauJet"] = Categories.getUntrackedParameter<bool>("TauJet");
   m_inputCategories["ETM"] = Categories.getUntrackedParameter<bool>("ETM");
   m_inputCategories["ETT"] = Categories.getUntrackedParameter<bool>("ETT");
   m_inputCategories["HTT"] = Categories.getUntrackedParameter<bool>("HTT");
   m_inputCategories["HTM"] = Categories.getUntrackedParameter<bool>("HTM");
 
   // What to do if we want our output to be saved to a external file
   m_outputFile = ps.getUntrackedParameter<string>("outputFile", "");
 
   if (!m_outputFile.empty()) {
     cout << "L1T Monitoring histograms will be saved to " << m_outputFile.c_str() << endl;
   }
 
   bool disable = ps.getUntrackedParameter<bool>("disableROOToutput", false);
   if (disable) {
     m_outputFile = "";
   }
 }

◆ ~L1TRate()

L1TRate::~L1TRate ( )

override

Definition at line 79 of file L1TRate.cc.

79 {}

Member Function Documentation

◆ analyze()

void L1TRate::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Reimplemented from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 285 of file L1TRate.cc.

References gather_cfg::cout, edm::ParameterSet::getUntrackedParameter(), Level1TriggerRates::gtAlgoCountsRate(), L1GlobalTriggerReadoutRecord::gtFdlVector(), mps_fire::i, iEvent, muonGEMDigis_cfi::instLumi, edm::HandleBase::isValid(), m_algoBit, m_l1GtDataDaqInputTag, m_lsLuminosity, m_lsPrescaleIndex, m_lsRates, m_lsShiftGTRates, m_parameters, m_scalersSource_colLScal, m_scalersSource_triggerScalers, m_selectedTriggers, and m_verbose.

                                                                        {
   edm::Handle<L1GlobalTriggerReadoutRecord> gtReadoutRecordData;
   edm::Handle<Level1TriggerScalersCollection> triggerScalers;
   edm::Handle<LumiScalersCollection> colLScal;
 
   iEvent.getByToken(m_l1GtDataDaqInputTag, gtReadoutRecordData);
   iEvent.getByToken(m_scalersSource_colLScal, colLScal);
   iEvent.getByToken(m_scalersSource_triggerScalers, triggerScalers);
 
   // Integers
   int EventRun = iEvent.id().run();
   unsigned int eventLS = iEvent.id().luminosityBlock();
 
   // Getting the trigger trigger rates from GT and buffering it
   if (triggerScalers.isValid()) {
     Level1TriggerScalersCollection::const_iterator itL1TScalers = triggerScalers->begin();
     Level1TriggerRates trigRates(*itL1TScalers, EventRun);
 
     int gtLS = (*itL1TScalers).lumiSegmentNr() + m_lsShiftGTRates;
 
     // If we haven't got the data from this LS yet get it
     if (m_lsRates.find(gtLS) == m_lsRates.end()) {
       if (m_verbose) {
         cout << "[L1TRate:] Buffering GT Rates for LS=" << gtLS << endl;
       }
       map<TString, double> bufferRate;
 
       // Buffer the rate informations for all selected bits
       for (map<string, string>::const_iterator i = m_selectedTriggers.begin(); i != m_selectedTriggers.end(); i++) {
         string tTrigger = (*i).second;
 
         // If trigger name is defined we store the rate
         if (tTrigger != "Undefined") {
           unsigned int trigBit = m_algoBit[tTrigger];
           double trigRate = trigRates.gtAlgoCountsRate()[trigBit];
 
           bufferRate[tTrigger] = trigRate;
         }
       }
       m_lsRates[gtLS] = bufferRate;
     }
   }
 
   // Getting from the SCAL the luminosity information and buffering it
   if (colLScal.isValid() && !colLScal->empty()) {
     LumiScalersCollection::const_iterator itLScal = colLScal->begin();
     unsigned int scalLS = itLScal->sectionNumber();
 
     // If we haven't got the data from this SCAL LS yet get it
     if (m_lsLuminosity.find(scalLS) == m_lsLuminosity.end()) {
       if (m_verbose) {
         cout << "[L1TRate:] Buffering SCAL-HF Lumi for LS=" << scalLS << endl;
       }
       double instLumi = itLScal->instantLumi();                  // Getting Instant Lumi from HF (via SCAL)
       double deadTimeNormHF = itLScal->deadTimeNormalization();  // Getting Dead Time Normalization from HF (via SCAL)
 
       // If HF Dead Time Corrections is requested we apply it
       // NOTE: By default this is assumed false since for now WbM fits do NOT assume this correction
       if (m_parameters.getUntrackedParameter<bool>("useHFDeadTimeNormalization", false)) {
         // Protecting for deadtime = 0
         if (deadTimeNormHF == 0) {
           instLumi = 0;
         } else {
           instLumi = instLumi / deadTimeNormHF;
         }
       }
       // Buffering the luminosity information
       m_lsLuminosity[scalLS] = instLumi;
     }
   }
 
   // Getting the prescale index used when this event was triggered
   if (gtReadoutRecordData.isValid()) {
     // If we haven't got the data from this LS yet get it
     if (m_lsPrescaleIndex.find(eventLS) == m_lsPrescaleIndex.end()) {
       if (m_verbose) {
         cout << "[L1TRate:] Buffering Prescale Index for LS=" << eventLS << endl;
       }
 
       // Getting Final Decision Logic (FDL) Data from GT
       const vector<L1GtFdlWord>& gtFdlVectorData = gtReadoutRecordData->gtFdlVector();
 
       // Getting the index for the fdl data for this event
       int indexFDL = 0;
       for (unsigned int i = 0; i < gtFdlVectorData.size(); i++) {
         if (gtFdlVectorData[i].bxInEvent() == 0) {
           indexFDL = i;
           break;
         }
       }
 
       int CurrentPrescalesIndex = gtFdlVectorData[indexFDL].gtPrescaleFactorIndexAlgo();
       m_lsPrescaleIndex[eventLS] = CurrentPrescalesIndex;
     }
   }
 }

◆ beginLuminosityBlock()

void L1TRate::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiBlock,
		edm::EventSetup const &	c
	)

overrideprotected

Definition at line 193 of file L1TRate.cc.

References gather_cfg::cout, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), and m_verbose.

                                                                                         {
   if (m_verbose) {
     cout << "[L1TRate:] Called beginLuminosityBlock at LS=" << lumiBlock.id().luminosityBlock() << endl;
   }
 }

◆ bookHistograms()

void L1TRate::bookHistograms	(	DQMStore::IBooker &	ibooker,
		const edm::Run &	,
		const edm::EventSetup &	iSetup
	)

overrideprotectedvirtual

Implements DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 84 of file L1TRate.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), edm::EventSetup::getData(), L1TMenuHelper::getEtaRangeByAlias(), L1TMenuHelper::getLUSOTrigger(), edm::ParameterSet::getParameter(), L1TMenuHelper::getPrescaleByAlias(), L1TMenuHelper::getQualityAlias(), getXSexFitsOMDS(), getXSexFitsPython(), L1GtPrescaleFactors::gtPrescaleFactors(), mps_fire::i, m_algoBit, m_ErrorMonitor, m_helperTokens, m_inputCategories, m_l1GtPfAlgoToken, m_l1GtUtils, m_listsPrescaleFactors, m_maxNbins, m_menuToken, m_parameters, m_refPrescaleSet, m_selectedTriggers, m_templateFunctions, m_xSecObservedToExpected, m_xSecVsInstLumi, L1TRate_cfi::maxInstantLuminosity, optionsL1T::menu, L1TRate_cfi::minInstantLuminosity, L1GtUtils::retrieveL1EventSetup(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), and L1TRate_cfi::srcAlgoXSecFit.

                                                                                                  {
   const L1GtTriggerMenu& menu = iSetup.getData(m_menuToken);
   const L1GtPrescaleFactors& l1GtPfAlgo = iSetup.getData(m_l1GtPfAlgoToken);
 
   // Initializing DQM Monitor Elements
   ibooker.setCurrentFolder("L1T/L1TRate");
   m_ErrorMonitor = ibooker.book1D("ErrorMonitor", "ErrorMonitor", 5, 0, 5);
   m_ErrorMonitor->setBinLabel(1, "WARNING_DB_CONN_FAILED");         // Errors from L1TOMDSHelper
   m_ErrorMonitor->setBinLabel(2, "WARNING_DB_QUERY_FAILED");        // Errors from L1TOMDSHelper
   m_ErrorMonitor->setBinLabel(3, "WARNING_DB_INCORRECT_NBUNCHES");  // Errors from L1TOMDSHelper
   m_ErrorMonitor->setBinLabel(4, "WARNING_PY_MISSING_FIT");
   m_ErrorMonitor->setBinLabel(5, "UNKNOWN");
 
   // Retriving the list of prescale sets
   m_listsPrescaleFactors = &(l1GtPfAlgo.gtPrescaleFactors());
 
   // Getting Lowest Prescale Single Object Triggers from the menu
   L1TMenuHelper myMenuHelper = L1TMenuHelper(iSetup, m_helperTokens);
   m_l1GtUtils.retrieveL1EventSetup(iSetup);
   m_selectedTriggers = myMenuHelper.getLUSOTrigger(m_inputCategories, m_refPrescaleSet, m_l1GtUtils);
 
   //-> Getting template fits for the algLo cross sections
   int srcAlgoXSecFit = m_parameters.getParameter<int>("srcAlgoXSecFit");
   if (srcAlgoXSecFit == 0) {
     getXSexFitsOMDS(m_parameters);
   } else if (srcAlgoXSecFit == 1) {
     getXSexFitsPython(m_parameters);
   }
 
   for (const auto& algo : menu.gtAlgorithmMap()) {
     m_algoBit[algo.second.algoAlias()] = algo.second.algoBitNumber();
   }
 
   double minInstantLuminosity = m_parameters.getParameter<double>("minInstantLuminosity");
   double maxInstantLuminosity = m_parameters.getParameter<double>("maxInstantLuminosity");
 
   // Initializing DQM Monitor Elements
   for (map<string, string>::const_iterator i = m_selectedTriggers.begin(); i != m_selectedTriggers.end(); i++) {
     TString tCategory = (*i).first;
     TString tTrigger = (*i).second;
 
     TString tErrorMessage = "";
     TF1* tTestFunction;
 
     if (tTrigger != "Undefined" && m_templateFunctions.find(tTrigger) != m_templateFunctions.end()) {
       tTestFunction = m_templateFunctions[tTrigger];
     } else if (tTrigger == "Undefined") {
       TString tFunc = "-1";
       tTestFunction = new TF1("FitParametrization_" + tTrigger, tFunc, 0, double(m_maxNbins) - 0.5);
     } else if (m_templateFunctions.find(tTrigger) == m_templateFunctions.end()) {
       TString tFunc = "-1";
       tTestFunction = new TF1("FitParametrization_" + tTrigger, tFunc, 0, double(m_maxNbins) - 0.5);
       tErrorMessage = " (Undefined Test Function)";
     } else {
       TString tFunc = "-1";
       tTestFunction = new TF1("FitParametrization_" + tTrigger, tFunc, 0, double(m_maxNbins) - 0.5);
     }
 
     if (tTrigger != "Undefined") {
       if (myMenuHelper.getPrescaleByAlias(tCategory, tTrigger) != 1) {
         tErrorMessage += " WARNING: Default Prescale = ";
         tErrorMessage += myMenuHelper.getPrescaleByAlias(tCategory, tTrigger);
       }
 
       if (tCategory == "Mu" && myMenuHelper.getEtaRangeByAlias(tCategory, tTrigger) != 4294967295) {
         tErrorMessage += " WARNING: Eta Range = ";
         tErrorMessage += myMenuHelper.getEtaRangeByAlias(tCategory, tTrigger);
       } else if (tCategory == "EG" && myMenuHelper.getEtaRangeByAlias(tCategory, tTrigger) != 32639) {
         tErrorMessage += " WARNING: Eta Range = ";
         tErrorMessage += myMenuHelper.getEtaRangeByAlias(tCategory, tTrigger);
       } else if (tCategory == "IsoEG" && myMenuHelper.getEtaRangeByAlias(tCategory, tTrigger) != 32639) {
         tErrorMessage += " WARNING: Eta Range = ";
         tErrorMessage += myMenuHelper.getEtaRangeByAlias(tCategory, tTrigger);
       }
 
       if (tCategory == "Mu" && myMenuHelper.getQualityAlias(tCategory, tTrigger) != 240) {
         tErrorMessage += " WARNING: Quality = ";
         tErrorMessage += myMenuHelper.getQualityAlias(tCategory, tTrigger);
       }
     }
 
     ibooker.setCurrentFolder("L1T/L1TRate/TriggerCrossSections");
     m_xSecVsInstLumi[tTrigger] = ibooker.bookProfile(tCategory,
                                                      "Cross Sec. vs Inst. Lumi Algo: " + tTrigger + tErrorMessage,
                                                      m_maxNbins,
                                                      minInstantLuminosity,
                                                      maxInstantLuminosity,
                                                      0,
                                                      500);
     m_xSecVsInstLumi[tTrigger]->setAxisTitle("Instantaneous Luminosity [10^{30}cm^{-2}s^{-1}]", 1);
     m_xSecVsInstLumi[tTrigger]->setAxisTitle("Algorithm #sigma [#mu b]", 2);
     m_xSecVsInstLumi[tTrigger]->getTProfile()->GetListOfFunctions()->Add(tTestFunction);
     m_xSecVsInstLumi[tTrigger]->getTProfile()->SetMarkerStyle(23);
 
     ibooker.setCurrentFolder("L1T/L1TRate/Certification");
     m_xSecObservedToExpected[tTrigger] =
         ibooker.book1D(tCategory, "Algo: " + tTrigger + tErrorMessage, m_maxNbins, -0.5, double(m_maxNbins) - 0.5);
     m_xSecObservedToExpected[tTrigger]->setAxisTitle("Lumi Section", 1);
     m_xSecObservedToExpected[tTrigger]->setAxisTitle("#sigma_{obs} / #sigma_{exp}", 2);
   }
 }

◆ dqmBeginRun()

void L1TRate::dqmBeginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Reimplemented from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 186 of file L1TRate.cc.

References gather_cfg::cout, and m_verbose.

                                                              {
   //
   if (m_verbose) {
     cout << "[L1TRate:] Called beginRun." << endl;
   }
 }

◆ endLuminosityBlock()

void L1TRate::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiBlock,
		edm::EventSetup const &	c
	)

overrideprotected

Definition at line 200 of file L1TRate.cc.

References gather_cfg::cout, mps_fire::i, edm::LuminosityBlockBase::id(), eostools::ls(), edm::LuminosityBlockID::luminosityBlock(), m_algoBit, m_lsLuminosity, m_lsPrescaleIndex, m_lsRates, m_selectedTriggers, m_verbose, m_xSecObservedToExpected, m_xSecVsInstLumi, and genParticles_cff::map.

                                                                                       {
   int eventLS = lumiBlock.id().luminosityBlock();
   if (m_verbose) {
     cout << "[L1TRate:] Called endLuminosityBlock at LS=" << eventLS << endl;
   }
 
   // We can certify LS -1 since we should have available:
   // gt rates: (current LS)-1
   // prescale: current LS
   // lumi    : current LS
   //eventLS--;
 
   // Checking if all necessary quantities are defined for our calculations
   bool isDefRate, isDefLumi, isDefPrescaleIndex;
   map<TString, double>* rates = nullptr;
   double lumi = 0;
   int prescalesIndex = 0;
 
   // Reseting MonitorElements so we can refill them
   for (map<string, string>::const_iterator i = m_selectedTriggers.begin(); i != m_selectedTriggers.end(); i++) {
     string tTrigger = (*i).second;
     m_xSecObservedToExpected[tTrigger]->getTH1()->Reset("ICE");
     m_xSecVsInstLumi[tTrigger]->getTH1()->Reset("ICE");
   }
 
   for (map<int, map<TString, double> >::iterator i = m_lsRates.begin(); i != m_lsRates.end(); i++) {
     unsigned int ls = (*i).first;
     rates = &(*i).second;
     isDefRate = true;
 
     if (m_lsLuminosity.find(ls) == m_lsLuminosity.end()) {
       isDefLumi = false;
     } else {
       isDefLumi = true;
       lumi = m_lsLuminosity[ls];
     }
 
     if (m_lsPrescaleIndex.find(ls) == m_lsPrescaleIndex.end()) {
       isDefPrescaleIndex = false;
     } else {
       isDefPrescaleIndex = true;
       prescalesIndex = m_lsPrescaleIndex[ls];
     }
 
     if (isDefRate && isDefLumi && isDefPrescaleIndex) {
       const vector<int>& currentPrescaleFactors = (*m_listsPrescaleFactors).at(prescalesIndex);
 
       for (map<string, string>::const_iterator i = m_selectedTriggers.begin(); i != m_selectedTriggers.end(); i++) {
         string tTrigger = (*i).second;
         TF1* tTestFunction = (TF1*)m_xSecVsInstLumi[tTrigger]->getTProfile()->GetListOfFunctions()->First();
 
         // If trigger name is defined we get the rate fit parameters
         if (tTrigger != "Undefined") {
           unsigned int trigBit = m_algoBit[tTrigger];
           double trigPrescale = currentPrescaleFactors[trigBit];
           double trigRate = (*rates)[tTrigger];
 
           if (lumi != 0 && trigPrescale != 0 && trigRate != 0) {
             double AlgoXSec = (trigPrescale * trigRate) / lumi;
             double TemplateFunctionValue = tTestFunction->Eval(lumi);
 
             // Checking against Template function
             int ibin = m_xSecObservedToExpected[tTrigger]->getTH1()->FindBin(ls);
             m_xSecObservedToExpected[tTrigger]->setBinContent(ibin, AlgoXSec / TemplateFunctionValue);
             m_xSecVsInstLumi[tTrigger]->Fill(lumi, AlgoXSec);
 
             if (m_verbose) {
               cout << "[L1TRate:] ls=" << ls << " Algo=" << tTrigger << " XSec=" << AlgoXSec
                    << " Test=" << AlgoXSec / TemplateFunctionValue << endl;
             }
 
           } else {
             int ibin = m_xSecObservedToExpected[tTrigger]->getTH1()->FindBin(ls);
             m_xSecObservedToExpected[tTrigger]->setBinContent(ibin, 0.000001);
             if (m_verbose) {
               cout << "[L1TRate:] Algo=" << tTrigger << " XSec=Failed" << endl;
             }
           }
         }
       }
     }
   }
 }

◆ getXSexFitsOMDS()

bool L1TRate::getXSexFitsOMDS ( const edm::ParameterSet & ps )

private

Definition at line 391 of file L1TRate.cc.

References a, WbMTriggerXSecFit::bitName, L1TOMDSHelper::connect(), dqm::impl::MonitorElement::Fill(), WbMTriggerXSecFit::fitFunction, edm::ParameterSet::getParameter(), m_ErrorMonitor, m_parameters, m_selectedTriggers, m_templateFunctions, L1TRate_cfi::maxInstantLuminosity, L1TRate_cfi::minInstantLuminosity, L1TOMDSHelper::NO_ERROR, L1TBPTX_cfi::oracleDB, WbMTriggerXSecFit::p0, WbMTriggerXSecFit::p1, WbMTriggerXSecFit::p2, L1TBPTX_cfi::pathCondDB, and WbMTriggerXSecFit::pm1.

Referenced by bookHistograms().

                                                        {
   bool noError = true;
 
   string oracleDB = ps.getParameter<string>("oracleDB");
   string pathCondDB = ps.getParameter<string>("pathCondDB");
 
   L1TOMDSHelper myOMDSHelper;
   int conError;
   myOMDSHelper.connect(oracleDB, pathCondDB, conError);
 
   map<string, WbMTriggerXSecFit> wbmFits;
 
   if (conError == L1TOMDSHelper::NO_ERROR) {
     int errorRetrive;
     wbmFits = myOMDSHelper.getWbMAlgoXsecFits(errorRetrive);
 
     // Filling errors if they exist
     if (errorRetrive != L1TOMDSHelper::NO_ERROR) {
       noError = false;
       string eName = myOMDSHelper.enumToStringError(errorRetrive);
       m_ErrorMonitor->Fill(eName);
     }
   } else {
     noError = false;
     string eName = myOMDSHelper.enumToStringError(conError);
     m_ErrorMonitor->Fill(eName);
   }
 
   double minInstantLuminosity = m_parameters.getParameter<double>("minInstantLuminosity");
   double maxInstantLuminosity = m_parameters.getParameter<double>("maxInstantLuminosity");
 
   // Getting rate fit parameters for all input triggers
   for (map<string, string>::const_iterator a = m_selectedTriggers.begin(); a != m_selectedTriggers.end(); a++) {
     string tTrigger = (*a).second;
 
     // If trigger name is defined we get the rate fit parameters
     if (tTrigger != "Undefined") {
       if (wbmFits.find(tTrigger) != wbmFits.end()) {
         WbMTriggerXSecFit tWbMParameters = wbmFits[tTrigger];
 
         vector<double> tParameters;
         tParameters.push_back(tWbMParameters.pm1);
         tParameters.push_back(tWbMParameters.p0);
         tParameters.push_back(tWbMParameters.p1);
         tParameters.push_back(tWbMParameters.p2);
 
         // Retriving and populating the m_templateFunctions array
         m_templateFunctions[tTrigger] = new TF1("FitParametrization_" + tWbMParameters.bitName,
                                                 tWbMParameters.fitFunction,
                                                 minInstantLuminosity,
                                                 maxInstantLuminosity);
         m_templateFunctions[tTrigger]->SetParameters(&tParameters[0]);
         m_templateFunctions[tTrigger]->SetLineWidth(1);
         m_templateFunctions[tTrigger]->SetLineColor(kRed);
 
       } else {
         noError = false;
       }
     }
   }
 
   return noError;
 }

◆ getXSexFitsPython()

bool L1TRate::getXSexFitsPython ( const edm::ParameterSet & ps )

private

Definition at line 464 of file L1TRate.cc.

References a, b, dqm::impl::MonitorElement::Fill(), edm::ParameterSet::getParameter(), m_ErrorMonitor, m_parameters, m_selectedTriggers, m_templateFunctions, L1TRate_cfi::maxInstantLuminosity, and L1TRate_cfi::minInstantLuminosity.

Referenced by bookHistograms().

                                                          {
   // error meaning
   bool noError = true;
 
   // Getting fit parameters
   std::vector<edm::ParameterSet> m_fitParameters = ps.getParameter<vector<ParameterSet> >("fitParameters");
 
   double minInstantLuminosity = m_parameters.getParameter<double>("minInstantLuminosity");
   double maxInstantLuminosity = m_parameters.getParameter<double>("maxInstantLuminosity");
 
   // Getting rate fit parameters for all input triggers
   for (map<string, string>::const_iterator a = m_selectedTriggers.begin(); a != m_selectedTriggers.end(); a++) {
     string tTrigger = (*a).second;
 
     // If trigger name is defined we get the rate fit parameters
     if (tTrigger != "Undefined") {
       bool foundFit = false;
 
       for (unsigned int b = 0; b < m_fitParameters.size(); b++) {
         if (tTrigger == m_fitParameters[b].getParameter<string>("AlgoName")) {
           TString tAlgoName = m_fitParameters[b].getParameter<string>("AlgoName");
           TString tTemplateFunction = m_fitParameters[b].getParameter<string>("TemplateFunction");
           vector<double> tParameters = m_fitParameters[b].getParameter<vector<double> >("Parameters");
 
           // Retriving and populating the m_templateFunctions array
           m_templateFunctions[tTrigger] =
               new TF1("FitParametrization_" + tAlgoName, tTemplateFunction, minInstantLuminosity, maxInstantLuminosity);
           m_templateFunctions[tTrigger]->SetParameters(&tParameters[0]);
           m_templateFunctions[tTrigger]->SetLineWidth(1);
           m_templateFunctions[tTrigger]->SetLineColor(kRed);
 
           foundFit = true;
           break;
         }
       }
       if (!foundFit) {
         noError = false;
         string eName = "WARNING_PY_MISSING_FIT";
         m_ErrorMonitor->Fill(eName);
       }
     }
   }
 
   return noError;
 }