#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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)=0

virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

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

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)

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

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

edm::EDGetTokenT< L1GlobalTriggerReadoutRecord >	m_l1GtDataDaqInputTag

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

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
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 48 of file L1TRate.h.

Constructor & Destructor Documentation

◆ L1TRate()

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

Definition at line 33 of file L1TRate.cc.

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

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.

◆ ~L1TRate()

L1TRate::~L1TRate ( )

override

Definition at line 76 of file L1TRate.cc.

76 {}

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 288 of file L1TRate.cc.

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

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.

◆ beginLuminosityBlock()

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

overrideprotected

Definition at line 196 of file L1TRate.cc.

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

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

◆ bookHistograms()

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

overrideprotected

Definition at line 81 of file L1TRate.cc.

                                                                                                  {
   ESHandle<L1GtTriggerMenu> menuRcd;
   ESHandle<L1GtPrescaleFactors> l1GtPfAlgo;
  
   iSetup.get<L1GtTriggerMenuRcd>().get(menuRcd);
   iSetup.get<L1GtPrescaleFactorsAlgoTrigRcd>().get(l1GtPfAlgo);
  
   const L1GtTriggerMenu* menu = menuRcd.product();
   const L1GtPrescaleFactors* m_l1GtPfAlgo = l1GtPfAlgo.product();
  
   // 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 = &(m_l1GtPfAlgo->gtPrescaleFactors());
  
   // Getting Lowest Prescale Single Object Triggers from the menu
   L1TMenuHelper myMenuHelper = L1TMenuHelper(iSetup);
   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 (CItAlgo algo = menu->gtAlgorithmMap().begin(); algo != menu->gtAlgorithmMap().end(); ++algo) {
     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);
   }
 }

References edm::EventSetup::get(), get, L1TMenuHelper::getEtaRangeByAlias(), L1TMenuHelper::getLUSOTrigger(), edm::ParameterSet::getParameter(), L1TMenuHelper::getPrescaleByAlias(), L1TMenuHelper::getQualityAlias(), getXSexFitsOMDS(), getXSexFitsPython(), L1GtPrescaleFactors::gtPrescaleFactors(), mps_fire::i, m_algoBit, m_ErrorMonitor, m_inputCategories, m_l1GtUtils, m_listsPrescaleFactors, m_maxNbins, m_parameters, m_refPrescaleSet, m_selectedTriggers, m_templateFunctions, m_xSecObservedToExpected, m_xSecVsInstLumi, L1TRate_cfi::maxInstantLuminosity, optionsL1T::menu, L1TRate_cfi::minInstantLuminosity, edm::ESHandle< T >::product(), L1GtUtils::retrieveL1EventSetup(), dqm::impl::MonitorElement::setBinLabel(), L1TRate_cfi::srcAlgoXSecFit, and tools::TF1.

◆ dqmBeginRun()

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

overrideprotectedvirtual

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

Definition at line 189 of file L1TRate.cc.

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

References gather_cfg::cout, and m_verbose.

◆ endLuminosityBlock()

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

overrideprotected

Definition at line 203 of file L1TRate.cc.

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

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, genParticles_cff::map, and tools::TF1.

◆ getXSexFitsOMDS()

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

private

Definition at line 394 of file L1TRate.cc.

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

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, WbMTriggerXSecFit::pm1, and tools::TF1.

Referenced by bookHistograms().

◆ getXSexFitsPython()

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

private

Definition at line 467 of file L1TRate.cc.

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

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

Referenced by bookHistograms().

Member Data Documentation

◆ m_algoBit

std::map<TString, int> L1TRate::m_algoBit

private

Definition at line 89 of file L1TRate.h.

Referenced by analyze(), bookHistograms(), and endLuminosityBlock().

◆ m_ErrorMonitor

MonitorElement* L1TRate::m_ErrorMonitor

private

Definition at line 105 of file L1TRate.h.

Referenced by bookHistograms(), getXSexFitsOMDS(), and getXSexFitsPython().

◆ m_inputCategories

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

private

Definition at line 90 of file L1TRate.h.

Referenced by bookHistograms(), and L1TRate().

◆ m_l1GtDataDaqInputTag

edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> L1TRate::m_l1GtDataDaqInputTag

private

Definition at line 99 of file L1TRate.h.

Referenced by analyze(), and L1TRate().

◆ m_l1GtUtils

L1GtUtils L1TRate::m_l1GtUtils

private

Definition at line 110 of file L1TRate.h.

Referenced by bookHistograms().

◆ m_listsPrescaleFactors

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

private

Definition at line 83 of file L1TRate.h.

Referenced by bookHistograms().

◆ m_lsLuminosity

std::map<int, double> L1TRate::m_lsLuminosity

private

Definition at line 87 of file L1TRate.h.

Referenced by analyze(), and endLuminosityBlock().

◆ m_lsPrescaleIndex

std::map<int, int> L1TRate::m_lsPrescaleIndex

private

Definition at line 86 of file L1TRate.h.

Referenced by analyze(), and endLuminosityBlock().

◆ m_lsRates

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

private

Definition at line 88 of file L1TRate.h.

Referenced by analyze(), and endLuminosityBlock().

◆ m_lsShiftGTRates

int L1TRate::m_lsShiftGTRates

private

Definition at line 77 of file L1TRate.h.

Referenced by analyze(), and L1TRate().

◆ m_maxNbins

int L1TRate::m_maxNbins

private

Definition at line 76 of file L1TRate.h.

Referenced by bookHistograms(), and L1TRate().

◆ m_outputFile

std::string L1TRate::m_outputFile

private

Definition at line 80 of file L1TRate.h.

Referenced by L1TRate().

◆ m_parameters

edm::ParameterSet L1TRate::m_parameters

private

Definition at line 102 of file L1TRate.h.

Referenced by analyze(), bookHistograms(), getXSexFitsOMDS(), getXSexFitsPython(), and L1TRate().

◆ m_refPrescaleSet

int L1TRate::m_refPrescaleSet

private

Definition at line 75 of file L1TRate.h.

Referenced by bookHistograms(), and L1TRate().

◆ m_scalersSource_colLScal

edm::EDGetTokenT<LumiScalersCollection> L1TRate::m_scalersSource_colLScal

private

Definition at line 97 of file L1TRate.h.

Referenced by analyze(), and L1TRate().

◆ m_scalersSource_triggerScalers

edm::EDGetTokenT<Level1TriggerScalersCollection> L1TRate::m_scalersSource_triggerScalers

private

Definition at line 98 of file L1TRate.h.

Referenced by analyze(), and L1TRate().

◆ m_selectedTriggers

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

private

Definition at line 91 of file L1TRate.h.

Referenced by analyze(), bookHistograms(), endLuminosityBlock(), getXSexFitsOMDS(), and getXSexFitsPython().

◆ m_templateFunctions

std::map<TString, TF1*> L1TRate::m_templateFunctions

private

Definition at line 94 of file L1TRate.h.

Referenced by bookHistograms(), getXSexFitsOMDS(), and getXSexFitsPython().

◆ m_verbose

bool L1TRate::m_verbose

private

Definition at line 72 of file L1TRate.h.

Referenced by analyze(), beginLuminosityBlock(), dqmBeginRun(), endLuminosityBlock(), and L1TRate().

◆ m_xSecObservedToExpected

std::map<TString, MonitorElement*> L1TRate::m_xSecObservedToExpected

private

Definition at line 92 of file L1TRate.h.

Referenced by bookHistograms(), and endLuminosityBlock().

◆ m_xSecVsInstLumi

std::map<TString, MonitorElement*> L1TRate::m_xSecVsInstLumi

private

Definition at line 93 of file L1TRate.h.

Referenced by bookHistograms(), and endLuminosityBlock().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ L1TRate()

◆ ~L1TRate()

Member Function Documentation

◆ analyze()

◆ beginLuminosityBlock()

◆ bookHistograms()

◆ dqmBeginRun()

◆ endLuminosityBlock()

◆ getXSexFitsOMDS()

◆ getXSexFitsPython()

Member Data Documentation

◆ m_algoBit

◆ m_ErrorMonitor

◆ m_inputCategories

◆ m_l1GtDataDaqInputTag

◆ m_l1GtUtils

◆ m_listsPrescaleFactors

◆ m_lsLuminosity

◆ m_lsPrescaleIndex

◆ m_lsRates

◆ m_lsShiftGTRates

◆ m_maxNbins

◆ m_outputFile

◆ m_parameters

◆ m_refPrescaleSet

◆ m_scalersSource_colLScal

◆ m_scalersSource_triggerScalers

◆ m_selectedTriggers

◆ m_templateFunctions

◆ m_verbose

◆ m_xSecObservedToExpected

◆ m_xSecVsInstLumi