#include <L1TRate.h>

Inheritance diagram for L1TRate:

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

virtual	~L1TRate ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

virtual void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

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

virtual	~EDConsumerBase () noexcept(false)

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

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

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

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

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

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

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

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

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 edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 49 of file L1TRate.h.

Constructor & Destructor Documentation

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

Definition at line 35 of file L1TRate.cc.

References TagProbeFitTreeAnalyzer_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_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.size() != 0) {
     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 ( )

virtual

Definition at line 77 of file L1TRate.cc.

77 {}

Member Function Documentation

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

overrideprotected

Definition at line 293 of file L1TRate.cc.

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

                                                                         {
 
   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->size()){
     
     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;   
     }    
   }
 }

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

overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 203 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;}
 
 }

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

overrideprotectedvirtual

Definition at line 82 of file L1TRate.cc.

References patPFMETCorrections_cff::algo, DQMStore::IBooker::book1D(), DQMStore::IBooker::bookProfile(), edm::EventSetup::get(), L1TMenuHelper::getEtaRangeByAlias(), L1TMenuHelper::getLUSOTrigger(), edm::ParameterSet::getParameter(), L1TMenuHelper::getPrescaleByAlias(), L1TMenuHelper::getQualityAlias(), getXSexFitsOMDS(), getXSexFitsPython(), L1GtTriggerMenu::gtAlgorithmMap(), 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, edm::ESHandle< T >::product(), L1GtUtils::retrieveL1EventSetup(), MonitorElement::setBinLabel(), and DQMStore::IBooker::setCurrentFolder().

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

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 198 of file L1TRate.cc.

References gather_cfg::cout, and m_verbose.

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

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

overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 210 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=0;
   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;}
           }
         }
       }  
     }    
   }
 }

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

private

Definition at line 397 of file L1TRate.cc.

References a, WbMTriggerXSecFit::bitName, L1TOMDSHelper::connect(), MonitorElement::Fill(), WbMTriggerXSecFit::fitFunction, edm::ParameterSet::getParameter(), m_ErrorMonitor, m_parameters, m_selectedTriggers, m_templateFunctions, L1TOMDSHelper::NO_ERROR, WbMTriggerXSecFit::p0, WbMTriggerXSecFit::p1, WbMTriggerXSecFit::p2, 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;
 
 }

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

private

Definition at line 472 of file L1TRate.cc.

References a, b, MonitorElement::Fill(), edm::ParameterSet::getParameter(), m_ErrorMonitor, m_parameters, m_selectedTriggers, and m_templateFunctions.

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