#include <L1TRate_Offline.h>

Inheritance diagram for L1TRate_Offline:

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

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

	~L1TRate_Offline () override

Public Member Functions inherited from one::DQMEDAnalyzer< edm::one::WatchLuminosityBlocks >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

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

bool	getXSexFitsPython (const edm::ParameterSet &ps)

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

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

std::map< TString, MonitorElement * >	m_InstLumiVsLS

edm::EDGetTokenT< L1GlobalTriggerReadoutRecord >	m_l1GtDataDaqInputTag

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

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

Detailed Description

Definition at line 46 of file L1TRate_Offline.h.

Member Enumeration Documentation

enum L1TRate_Offline::Errors

Enumerator
UNKNOWN
WARNING_PY_MISSING_FIT

Definition at line 50 of file L1TRate_Offline.h.

              {
     UNKNOWN                = 1,
     WARNING_PY_MISSING_FIT = 2
   };

Constructor & Destructor Documentation

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

Definition at line 19 of file L1TRate_Offline.cc.

References TagProbeFitTreeAnalyzer_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_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 ( )

override

Definition at line 55 of file L1TRate_Offline.cc.

55 {}

Member Function Documentation

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

overrideprotected

Definition at line 347 of file L1TRate_Offline.cc.

References gather_cfg::cout, Level1TriggerRates::deadtimePercent(), edm::Event::getByToken(), edm::ParameterSet::getUntrackedParameter(), L1GlobalTriggerReadoutRecord::gtFdlVector(), mps_fire::i, edm::EventBase::id(), muonGEMDigis_cfi::instLumi, edm::HandleBase::isValid(), edm::EventID::luminosityBlock(), m_algoBit, m_l1GtDataDaqInputTag, m_lsDeadTime, m_lsLuminosity, m_lsPrescaleIndex, m_lsRates, m_lsShiftGTRates, m_parameters, m_scalersSource_L1TSCollection, m_scalersSource_LSCollection, 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_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;
     }
 
     }
 
   }
 }

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

overrideprotected

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

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

overrideprotected

Definition at line 60 of file L1TRate_Offline.cc.

References patPFMETCorrections_cff::algo, DQMStore::IBooker::book1D(), DQMStore::IBooker::bookProfile(), gather_cfg::cout, edm::EventSetup::get(), L1TMenuHelper::getEtaRangeByAlias(), L1TMenuHelper::getLUSOTrigger(), edm::ParameterSet::getParameter(), L1TMenuHelper::getPrescaleByAlias(), L1TMenuHelper::getQualityAlias(), getXSexFitsPython(), L1GtTriggerMenu::gtAlgorithmMap(), L1GtPrescaleFactors::gtPrescaleFactors(), mps_fire::i, m_algoBit, m_algoFit, m_CountsVsLS, m_DelivLumiVsLS, m_ErrorMonitor, m_inputCategories, m_l1GtUtils, m_listsPrescaleFactors, m_maxNbins, m_parameters, m_RecorLumiVsLS, m_refPrescaleSet, m_selectedTriggers, m_templateFunctions, m_verbose, m_xSecObservedToExpected, m_xSecObservedVsDelivLumi, m_xSecObservedVsLS, m_xSecObservedVsRecorLumi, edm::ESHandle< T >::product(), L1GtUtils::retrieveL1EventSetup(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), UNKNOWN, and WARNING_PY_MISSING_FIT.

                                                                                                             {
 
   if (m_verbose) {cout << "[L1TRate_Offline:] Called beginRun." << endl;}
 
   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",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 = &(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
   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){ //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);
   }
 }

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

overrideprotected

Definition at line 207 of file L1TRate_Offline.cc.

                                                                     {
 
 }

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

overrideprotected

Definition at line 218 of file L1TRate_Offline.cc.

References gather_cfg::cout, mps_fire::i, edm::LuminosityBlockBase::id(), 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);
          }
         }
       }
     }
   }
 }

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

protected

Definition at line 487 of file L1TRate_Offline.cc.

References a, b, DEFINE_FWK_MODULE, edm::ParameterSet::getParameter(), MonitorElement::getTH1(), m_ErrorMonitor, m_parameters, m_selectedTriggers, m_templateFunctions, 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;
 
 }