#include <L1TRate_Offline.h>

Inheritance diagram for L1TRate_Offline:

Public Types
enum	Errors { UNKNOWN = 1, WARNING_PY_MISSING_FIT = 2 }

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

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

	~L1TRate_Offline () 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 &run, const edm::EventSetup &iSetup) override

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

bool	getXSexFitsPython (const edm::ParameterSet &ps)

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

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 Attributes
std::map< TString, int >	m_algoBit

std::map< TString, TF1 * >	m_algoFit

std::map< TString, MonitorElement * >	m_CountsVsLS

std::map< TString, MonitorElement * >	m_DelivLumiVsLS

MonitorElement *	m_ErrorMonitor

L1TMenuHelper::Tokens	m_helperTokens

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

std::map< TString, MonitorElement * >	m_InstLumiVsLS

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, std::map< TString, double > >	m_lsCounts

std::map< int, double >	m_lsDeadTime

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

std::map< TString, MonitorElement * >	m_PrescIndexVsLS

std::map< TString, MonitorElement * >	m_RecorLumiVsLS

int	m_refPrescaleSet

edm::EDGetTokenT< Level1TriggerScalersCollection >	m_scalersSource_L1TSCollection

edm::EDGetTokenT< LumiScalersCollection >	m_scalersSource_LSCollection

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_xSecObservedVsDelivLumi

std::map< TString, MonitorElement * >	m_xSecObservedVsLS

std::map< TString, MonitorElement * >	m_xSecObservedVsRecorLumi

std::map< TString, MonitorElement * >	m_xSecVsInstLumi

Additional Inherited Members
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 43 of file L1TRate_Offline.h.

Member Enumeration Documentation

◆ Errors

enum L1TRate_Offline::Errors

Enumerator
UNKNOWN
WARNING_PY_MISSING_FIT

Definition at line 45 of file L1TRate_Offline.h.

45 { UNKNOWN = 1, WARNING_PY_MISSING_FIT = 2 };

L1TRate_Offline::UNKNOWN

Definition: L1TRate_Offline.h:45

L1TRate_Offline::WARNING_PY_MISSING_FIT

Definition: L1TRate_Offline.h:45

Constructor & Destructor Documentation

◆ L1TRate_Offline()

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

Definition at line 18 of file L1TRate_Offline.cc.

References L1TSync_cfi::Categories, gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_inputCategories, m_l1GtDataDaqInputTag, m_lsShiftGTRates, m_maxNbins, m_parameters, m_refPrescaleSet, m_scalersSource_L1TSCollection, m_scalersSource_LSCollection, 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_LSCollection =
       consumes<LumiScalersCollection>(m_parameters.getParameter<InputTag>("inputTagScalersResults"));
   m_scalersSource_L1TSCollection =
       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");
 
   // Initializing Variables
   if (m_verbose) {
     cout << "[L1TRate_Offline:] ____________ Storage initialization ____________ " << endl;
     cout << "[L1TRate_Offline:] Setting up dbe folder: L1T/L1TRate" << endl;
   }
 }

◆ ~L1TRate_Offline()

L1TRate_Offline::~L1TRate_Offline ( )

override

Definition at line 59 of file L1TRate_Offline.cc.

59 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

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

Definition at line 338 of file L1TRate_Offline.cc.

References gather_cfg::cout, Level1TriggerRates::deadtimePercent(), edm::ParameterSet::getUntrackedParameter(), L1GlobalTriggerReadoutRecord::gtFdlVector(), mps_fire::i, iEvent, muonGEMDigis_cfi::instLumi, edm::HandleBase::isValid(), m_algoBit, m_l1GtDataDaqInputTag, m_lsDeadTime, m_lsLuminosity, m_lsPrescaleIndex, m_lsRates, m_lsShiftGTRates, m_parameters, m_scalersSource_L1TSCollection, m_scalersSource_LSCollection, 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_LSCollection, colLScal);
   iEvent.getByToken(m_scalersSource_L1TSCollection, 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() && !triggerScalers->empty()) {
     Level1TriggerScalersCollection::const_iterator itL1TScalers = triggerScalers->begin();
     Level1TriggerRates trigRates(*itL1TScalers, EventRun);
 
     //Trying to retrieve GT DeadTime
     //unsigned long long deadtime = itL1TScalers->deadtime();
     //double deadtime = itL1TScalers->deadtime();
     double deadtime = trigRates.deadtimePercent();  //correct DeadTime % to retrieve
 
     // Trying to get the trigger counts
     const std::vector<unsigned int> gtAlgoCounts = itL1TScalers->gtAlgoCounts();
 
     //    int lumisegment = (*itL1TScalers).lumiSegmentNr();
 
     //     cout << "deadtime =" << deadtime << " --
     //    cout << "lumisegment = " << lumisegment << endl;
 
     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()) {
       map<TString, double> bufferCount;
 
       // 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 trigCount = gtAlgoCounts[trigBit];
 
           bufferCount[tTrigger] = trigCount;
         }
       }
       m_lsRates[gtLS] = bufferCount;
       m_lsDeadTime[gtLS] = deadtime;
     }
   }
 
   // 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_Offline:] Buffering SCAL-HF Lumi for LS=" << scalLS << endl;
       }
       double instLumi =
           itLScal->instantLumi();  // Getting Instant Lumi from HF (via SCAL) // <###### WE NEED TO STORE THIS
       double deadTimeNormHF = itLScal->deadTimeNormalization();  // Getting Dead Time Normalization from HF (via SCAL)
 
       //       double mylumiFill      = itLScal->lumiFill();      // Integrated lumi since beginning of fill, delivered
       //       double mylumiRun       = itLScal->lumiRun();       // Integrated lumi since beginning of run, delivered
       //       double myliveLumiFill  = itLScal->liveLumiFill();  // Integrated lumi since beginning of fill, live
       //       double myliveLumiRun   = itLScal->liveLumiRun();   // Integrated lumi since beginning of run, live
       //       double mysectionNumber = itLScal->sectionNumber(); // Lumi section  number for this info
       //       double mynumOrbits     = itLScal->numOrbits();     // Number of orbits that have passed this run
 
       //       cout << "instantLumi = " << instLumi << endl;
       //       cout << "lumiFill = " << mylumiFill << endl;
       //       cout << "lumiRun = " << mylumiRun << endl;
       //       cout << "livelumiFill = " << myliveLumiFill << endl;
       //       cout << "livelumiRun = " << myliveLumiRun << endl;
       //       cout << "sectionNumber = " << mysectionNumber << endl;
       //       cout << "numOrbits = " << mynumOrbits << endl;
 
       // 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_Offline:] 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;
         }
       }
 
       if (!gtFdlVectorData.empty()) {
         int CurrentPrescalesIndex =
             gtFdlVectorData[indexFDL].gtPrescaleFactorIndexAlgo();  // <###### WE NEED TO STORE THIS
         m_lsPrescaleIndex[eventLS] = CurrentPrescalesIndex;
       }
     }
   }
 }

◆ beginLuminosityBlock()

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

overrideprotected

Definition at line 210 of file L1TRate_Offline.cc.

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

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

◆ bookHistograms()

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

overrideprotectedvirtual

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

Definition at line 64 of file L1TRate_Offline.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), gather_cfg::cout, edm::EventSetup::getData(), L1TMenuHelper::getEtaRangeByAlias(), L1TMenuHelper::getLUSOTrigger(), edm::ParameterSet::getParameter(), L1TMenuHelper::getPrescaleByAlias(), L1TMenuHelper::getQualityAlias(), getXSexFitsPython(), L1GtPrescaleFactors::gtPrescaleFactors(), mps_fire::i, m_algoBit, m_algoFit, m_CountsVsLS, m_DelivLumiVsLS, m_ErrorMonitor, m_helperTokens, m_inputCategories, m_l1GtPfAlgoToken, m_l1GtUtils, m_listsPrescaleFactors, m_maxNbins, m_menuToken, m_parameters, m_RecorLumiVsLS, m_refPrescaleSet, m_selectedTriggers, m_templateFunctions, m_verbose, m_xSecObservedToExpected, m_xSecObservedVsDelivLumi, m_xSecObservedVsLS, m_xSecObservedVsRecorLumi, L1TRate_cfi::maxInstantLuminosity, optionsL1T::menu, L1TRate_cfi::minInstantLuminosity, L1GtUtils::retrieveL1EventSetup(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), UNKNOWN, and WARNING_PY_MISSING_FIT.

                                                                                                              {
   if (m_verbose) {
     cout << "[L1TRate_Offline:] Called beginRun." << endl;
   }
 
   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", 2, 0, 2);
   m_ErrorMonitor->setBinLabel(UNKNOWN, "UNKNOWN");
   m_ErrorMonitor->setBinLabel(WARNING_PY_MISSING_FIT, "WARNING_PY_MISSING_FIT");
 
   if (m_verbose) {
     cout << "[L1TRate_Offline:] m_ErrorMonitor: " << m_ErrorMonitor << endl;
   }
 
   // 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
   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) {  //hexadecimal of the whole range
         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/xSecDelivLumi");  // trigger counts...
     m_xSecObservedVsDelivLumi[tTrigger] = ibooker.bookProfile(tCategory,
                                                               "Cross Sec. vs Deliv. Lumi: " + tTrigger + tErrorMessage,
                                                               m_maxNbins,
                                                               minInstantLuminosity,
                                                               maxInstantLuminosity,
                                                               0,
                                                               500);
     m_xSecObservedVsDelivLumi[tTrigger]->setAxisTitle("Delivered Luminosity [10^{30}cm^{-2}s^{-1}]", 1);
     m_xSecObservedVsDelivLumi[tTrigger]->setAxisTitle("Algorithm #sigma [#mu b]", 2);
 
     ibooker.setCurrentFolder("L1T/L1TRate/xSecRecorLumi");  // trigger counts...
     m_xSecObservedVsRecorLumi[tTrigger] = ibooker.bookProfile(tCategory,
                                                               "Cross Sec. vs Recor. Lumi: " + tTrigger + tErrorMessage,
                                                               m_maxNbins,
                                                               minInstantLuminosity,
                                                               maxInstantLuminosity,
                                                               0,
                                                               500);
     m_xSecObservedVsRecorLumi[tTrigger]->setAxisTitle("Recorded Luminosity [10^{30}cm^{-2}s^{-1}]", 1);
     m_xSecObservedVsRecorLumi[tTrigger]->setAxisTitle("Algorithm #sigma [#mu b]", 2);
 
     ibooker.setCurrentFolder("L1T/L1TRate/TriggerCounts");  // trigger counts...
     m_CountsVsLS[tTrigger] = ibooker.bookProfile(tCategory,
                                                  "Cross Sec. vs Inst. Lumi Algo: " + tTrigger + tErrorMessage,
                                                  m_maxNbins,
                                                  minInstantLuminosity,
                                                  maxInstantLuminosity,
                                                  0,
                                                  500);
     m_CountsVsLS[tTrigger]->setAxisTitle("Instantaneous Luminosity [10^{30}cm^{-2}s^{-1}]", 1);
     m_CountsVsLS[tTrigger]->setAxisTitle("Algorithm #sigma [#mu b]", 2);
     m_CountsVsLS[tTrigger]->getTProfile()->GetListOfFunctions()->Add(tTestFunction);
     m_CountsVsLS[tTrigger]->getTProfile()->SetMarkerStyle(23);
 
     m_algoFit[tTrigger] = (TF1*)tTestFunction->Clone("Fit_" + tTrigger);  // NOTE: Workaround
 
     ibooker.setCurrentFolder("L1T/L1TRate/xSecObs");
     m_xSecObservedVsLS[tTrigger] =
         ibooker.book1D(tCategory, "Algo: " + tTrigger + tErrorMessage, m_maxNbins, -0.5, double(m_maxNbins) - 0.5);
     m_xSecObservedVsLS[tTrigger]->setAxisTitle("Lumi Section", 1);
     m_xSecObservedVsLS[tTrigger]->setAxisTitle("#sigma_{obs}", 2);
 
     ibooker.setCurrentFolder("L1T/L1TRate/Delivered");
     m_DelivLumiVsLS[tTrigger] =
         ibooker.book1D(tCategory, "Algo: " + tTrigger + tErrorMessage, m_maxNbins, -0.5, double(m_maxNbins) - 0.5);
     m_DelivLumiVsLS[tTrigger]->setAxisTitle("Lumi Section", 1);
     m_DelivLumiVsLS[tTrigger]->setAxisTitle("Deliv. Lumi", 2);
 
     ibooker.setCurrentFolder("L1T/L1TRate/Recorded");
     m_RecorLumiVsLS[tTrigger] =
         ibooker.book1D(tCategory, "Algo: " + tTrigger + tErrorMessage, m_maxNbins, -0.5, double(m_maxNbins) - 0.5);
     m_RecorLumiVsLS[tTrigger]->setAxisTitle("Lumi Section", 1);
     m_RecorLumiVsLS[tTrigger]->setAxisTitle("Recor. Lumi", 2);
 
     ibooker.setCurrentFolder("L1T/L1TRate/Ratio");
     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);
   }
 }

◆ endLuminosityBlock()

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

overrideprotected

Definition at line 217 of file L1TRate_Offline.cc.

References dqmiodumpmetadata::counts, gather_cfg::cout, mps_fire::i, edm::LuminosityBlockBase::id(), dqmiolumiharvest::j, edm::LuminosityBlockID::luminosityBlock(), m_algoBit, m_algoFit, m_CountsVsLS, m_DelivLumiVsLS, m_lsDeadTime, m_lsLuminosity, m_lsPrescaleIndex, m_lsRates, m_RecorLumiVsLS, m_selectedTriggers, m_verbose, m_xSecObservedToExpected, m_xSecObservedVsDelivLumi, m_xSecObservedVsLS, m_xSecObservedVsRecorLumi, genParticles_cff::map, and findQualityFiles::size.

                                                                                               {
   int eventLS = lumiBlock.id().luminosityBlock();
   if (m_verbose) {
     cout << "[L1TRate_Offline:] 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;
   bool isDefLumi, isDefPrescaleIndex;
   //map<TString,double>* rates=0;
   double lumi = 0;
   double deadtime = 0;
   unsigned int prescalesIndex = 0;
 
   bool isDefCount;
   map<TString, double>* counts = nullptr;
 
   // Resetting 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_DeadTimeVsLS          [tTrigger]->getTH1()->Reset("ICE");
     m_CountsVsLS[tTrigger]->getTH1()->Reset("ICE");
     m_xSecObservedToExpected[tTrigger]->getTH1()->Reset("ICE");
 
     m_xSecObservedVsLS[tTrigger]->getTH1()->Reset("ICE");
     m_DelivLumiVsLS[tTrigger]->getTH1()->Reset("ICE");
     m_RecorLumiVsLS[tTrigger]->getTH1()->Reset("ICE");
 
     m_xSecObservedVsDelivLumi[tTrigger]->getTH1()->Reset("ICE");
     m_xSecObservedVsRecorLumi[tTrigger]->getTH1()->Reset("ICE");
   }
 
   //Trying to do the same with Counts....
   for (map<int, map<TString, double> >::iterator j = m_lsRates.begin(); j != m_lsRates.end(); j++) {
     unsigned int lsOffline = (*j).first;
     counts = &(*j).second;
     isDefCount = true;
 
     unsigned int lsPreInd;
 
     if (m_lsLuminosity.find(lsOffline) == m_lsLuminosity.end()) {
       isDefLumi = false;
     } else {
       isDefLumi = true;
       lumi = m_lsLuminosity[lsOffline];
       deadtime = m_lsDeadTime[lsOffline];
     }
 
     lsPreInd = lsOffline + 1;  // NOTE: Workaround
 
     if (m_lsPrescaleIndex.find(lsPreInd) == m_lsPrescaleIndex.end()) {
       isDefPrescaleIndex = false;
     } else {
       isDefPrescaleIndex = true;
       prescalesIndex = m_lsPrescaleIndex[lsPreInd];
     }
 
     if (isDefCount && isDefLumi && isDefPrescaleIndex && (prescalesIndex < m_listsPrescaleFactors->size())) {
       const vector<int>& currentPrescaleFactors = (*m_listsPrescaleFactors).at(prescalesIndex);
 
       for (map<string, string>::const_iterator j = m_selectedTriggers.begin(); j != m_selectedTriggers.end(); j++) {
         string tTrigger = (*j).second;
         double trigCount = (*counts)[tTrigger];
 
         //   TF1*   tTestFunction = (TF1*) m_CountsVsLS[tTrigger]->getTProfile()->GetListOfFunctions()->First();
         TF1* tTestFunction = m_algoFit[tTrigger];  // NOTE: Workaround....
 
         // If trigger name is defined we get the rate fit parameters
         if (tTrigger != "Undefined") {
           unsigned int trigBit = m_algoBit[tTrigger];
           double trigPrescale = currentPrescaleFactors[trigBit];
 
           if (lumi != 0 && trigCount != 0 && trigPrescale != 0) {
             double RecLumi = lumi * (1. - deadtime / 100.);
 
             double AlgoXSec = (trigPrescale * trigCount) / RecLumi;
             double TemplateFunctionValue = tTestFunction->Eval(lumi);
 
             // Checking against Template function
             //            m_DeadTimeVsLS [tTrigger]->Fill(lumi,deadtime);
 
             m_xSecObservedVsRecorLumi[tTrigger]->Fill(RecLumi, AlgoXSec);
 
             m_CountsVsLS[tTrigger]->Fill(lumi, AlgoXSec);
 
             int ibin = m_xSecObservedToExpected[tTrigger]->getTH1()->FindBin(lsOffline);
             m_xSecObservedToExpected[tTrigger]->setBinContent(ibin, AlgoXSec / TemplateFunctionValue);
 
             m_DelivLumiVsLS[tTrigger]->setBinContent(ibin, lumi);
             m_RecorLumiVsLS[tTrigger]->setBinContent(ibin, RecLumi);
 
             m_xSecObservedVsLS[tTrigger]->setBinContent(ibin, AlgoXSec);
 
           } else {
             //            m_DeadTimeVsLS [tTrigger]->Fill(0.000001,0.000001);
 
             m_xSecObservedVsRecorLumi[tTrigger]->Fill(0.000001, 0.000001);
 
             m_CountsVsLS[tTrigger]->Fill(0.000001, 0.000001);
 
             int ibin = m_xSecObservedToExpected[tTrigger]->getTH1()->FindBin(lsOffline);
             m_xSecObservedToExpected[tTrigger]->setBinContent(ibin, 0.000001);
 
             m_DelivLumiVsLS[tTrigger]->setBinContent(ibin, 0.000001);
             m_RecorLumiVsLS[tTrigger]->setBinContent(ibin, 0.000001);
 
             m_xSecObservedVsLS[tTrigger]->setBinContent(ibin, 0.000001);
           }
         }
       }
     }
   }
 }

◆ getXSexFitsPython()

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

protected

Definition at line 474 of file L1TRate_Offline.cc.

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

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;
 
         int eCount = m_ErrorMonitor->getTH1()->GetBinContent(WARNING_PY_MISSING_FIT);
         eCount++;
         m_ErrorMonitor->getTH1()->SetBinContent(WARNING_PY_MISSING_FIT, eCount);
       }
     }
   }
 
   return noError;
 }