#include <L1TRate_Offline.h>

Inheritance diagram for L1TRate_Offline:

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

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

virtual	~L1TRate_Offline ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

void	beginJob ()

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

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

void	endJob ()

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

void	endRun (const edm::Run &run, const edm::EventSetup &iSetup)

bool	getXSexFitsPython (const edm::ParameterSet &ps)

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 Attributes
DQMStore *	dbe

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

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

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 43 of file L1TRate_Offline.h.

Member Enumeration Documentation

enum L1TRate_Offline::Errors

Enumerator
UNKNOWN
WARNING_PY_MISSING_FIT

Definition at line 47 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 22 of file L1TRate_Offline.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and cppFunctionSkipper::operator.

                                                        {
 
   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;
   }
 
   dbe = Service < DQMStore > ().operator->();
   dbe->setVerbose(0);
   dbe->setCurrentFolder("L1T/L1TRate");
 
   // Inicializing Variables
   if (m_verbose) {cout << "[L1TRate_Offline:] Pointer for DQM Store: " << dbe << endl;}
 }

L1TRate_Offline::~L1TRate_Offline ( )

virtual

Definition at line 64 of file L1TRate_Offline.cc.

64 {}

Member Function Documentation

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

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 391 of file L1TRate_Offline.cc.

References gather_cfg::cout, Level1TriggerRates::deadtimePercent(), edm::Event::getByToken(), i, edm::EventBase::id(), edm::HandleBase::isValid(), edm::EventID::luminosityBlock(), 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->size()){
 
     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->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_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;}
       }
 
       int CurrentPrescalesIndex  = gtFdlVectorData[indexFDL].gtPrescaleFactorIndexAlgo(); // <###### WE NEED TO STORE THIS
       m_lsPrescaleIndex[eventLS] = CurrentPrescalesIndex;
 
     }
 
   }
 
 
 }

void L1TRate_Offline::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 67 of file L1TRate_Offline.cc.

References gather_cfg::cout, cppFunctionSkipper::operator, DQMStore::rmdir(), and DQMStore::setCurrentFolder().

                                   {
 
   if (m_verbose) {cout << "[L1TRate_Offline:] Called beginJob." << endl;}
 
   // get hold of back-end interface
   DQMStore *dbe = 0;
   dbe = Service < DQMStore > ().operator->();
 
   if (dbe) {
     dbe->setCurrentFolder("L1T/L1TRate");
     dbe->rmdir("L1T/L1TRate");
   }
 
 }

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 255 of file L1TRate_Offline.cc.

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

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

void L1TRate_Offline::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 97 of file L1TRate_Offline.cc.

References gather_cfg::cout, edm::EventSetup::get(), L1TMenuHelper::getEtaRangeByAlias(), L1TMenuHelper::getLUSOTrigger(), L1TMenuHelper::getPrescaleByAlias(), L1TMenuHelper::getQualityAlias(), L1GtTriggerMenu::gtAlgorithmMap(), L1GtPrescaleFactors::gtPrescaleFactors(), i, relval_steps::menu, edm::ESHandle< class >::product(), and pat::UNKNOWN.

                                                                             {
 
   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
   dbe->setCurrentFolder("L1T/L1TRate");
   m_ErrorMonitor = dbe->book1D("ErrorMonitor", "ErrorMonitor",2,0,2);
   m_ErrorMonitor->setBinLabel(UNKNOWN               ,"UNKNOWN");
   m_ErrorMonitor->setBinLabel(WARNING_PY_MISSING_FIT,"WARNING_PY_MISSING_FIT");
 
   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_selectedTriggers = myMenuHelper.getLUSOTrigger(m_inputCategories,m_refPrescaleSet);
 
   //-> 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);
       }
 
     }
 
     dbe->setCurrentFolder("L1T/L1TRate/xSecDelivLumi"); // trigger counts...
     m_xSecObservedVsDelivLumi[tTrigger] = dbe->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);
 
     dbe->setCurrentFolder("L1T/L1TRate/xSecRecorLumi"); // trigger counts...
     m_xSecObservedVsRecorLumi[tTrigger] = dbe->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);
 
     dbe->setCurrentFolder("L1T/L1TRate/TriggerCounts"); // trigger counts...
     m_CountsVsLS[tTrigger] = dbe->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
 
     dbe->setCurrentFolder("L1T/L1TRate/xSecObs");
     m_xSecObservedVsLS[tTrigger] = dbe->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);
 
     dbe->setCurrentFolder("L1T/L1TRate/Delivered");
     m_DelivLumiVsLS[tTrigger] = dbe->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);
 
     dbe->setCurrentFolder("L1T/L1TRate/Recorded");
     m_RecorLumiVsLS[tTrigger] = dbe->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);
 
     dbe->setCurrentFolder("L1T/L1TRate/Ratio");
     m_xSecObservedToExpected[tTrigger] = dbe->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);
 
 //     dbe->setCurrentFolder("L1T/L1TRate/DeadTime");
 //     m_DeadTimeVsLS[tTrigger] = dbe->book1D(tCategory, "Dead Time: "+tTrigger+tErrorMessage,m_maxNbins,-0.5,double(m_maxNbins)-0.5);
 //     m_DeadTimeVsLS[tTrigger] ->setAxisTitle("Lumi Section" ,1);
 //     m_DeadTimeVsLS[tTrigger] ->setAxisTitle("Dead Time" ,2);
 
   }
 
 }

void L1TRate_Offline::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 83 of file L1TRate_Offline.cc.

References gather_cfg::cout.

                                 {
 
   if (m_verbose) {cout << "[L1TRate_Offline:] Called endJob." << endl;}
 
   if (m_outputFile.size() != 0 && dbe)
     dbe->save(m_outputFile);
 
   return;
 
 }

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 262 of file L1TRate_Offline.cc.

References gather_cfg::cout, i, edm::LuminosityBlockBase::id(), j, fjr2json::lumi, edm::LuminosityBlockID::luminosityBlock(), and python.multivaluedict::map().

                                                                                               {
 
   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;
   int                  prescalesIndex=0;
 
   bool isDefCount;
   map<TString,double>* counts=0;
 
   // 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){
 
       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);
          }
         }
       }
     }
   }
 }

void L1TRate_Offline::endRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 249 of file L1TRate_Offline.cc.

References gather_cfg::cout.

                                                                           {
 
   if (m_verbose) {cout << "[L1TRate_Offline:] Called endRun." << endl;}
 
 }

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

protected

Definition at line 530 of file L1TRate_Offline.cc.

References a, b, and edm::ParameterSet::getParameter().

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

Member Data Documentation

DQMStore* L1TRate_Offline::dbe

private

Definition at line 132 of file L1TRate_Offline.h.

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

private

Definition at line 96 of file L1TRate_Offline.h.

std::map<TString,TF1*> L1TRate_Offline::m_algoFit

private

Definition at line 97 of file L1TRate_Offline.h.

std::map<TString,MonitorElement*> L1TRate_Offline::m_CountsVsLS

private

Definition at line 107 of file L1TRate_Offline.h.

std::map<TString,MonitorElement*> L1TRate_Offline::m_DelivLumiVsLS

private

Definition at line 113 of file L1TRate_Offline.h.

MonitorElement* L1TRate_Offline::m_ErrorMonitor

private

Definition at line 129 of file L1TRate_Offline.h.

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

private

Definition at line 98 of file L1TRate_Offline.h.

std::map<TString,MonitorElement*> L1TRate_Offline::m_InstLumiVsLS

private

Definition at line 108 of file L1TRate_Offline.h.

edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> L1TRate_Offline::m_l1GtDataDaqInputTag

private

Definition at line 123 of file L1TRate_Offline.h.

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

private

Definition at line 89 of file L1TRate_Offline.h.

std::map<int,std::map<TString,double> > L1TRate_Offline::m_lsCounts

private

Definition at line 118 of file L1TRate_Offline.h.

std::map<int,double> L1TRate_Offline::m_lsDeadTime

private

Definition at line 92 of file L1TRate_Offline.h.

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

private

Definition at line 94 of file L1TRate_Offline.h.

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

private

Definition at line 93 of file L1TRate_Offline.h.

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

private

Definition at line 95 of file L1TRate_Offline.h.

int L1TRate_Offline::m_lsShiftGTRates

private

Definition at line 83 of file L1TRate_Offline.h.

int L1TRate_Offline::m_maxNbins

private

Definition at line 82 of file L1TRate_Offline.h.

std::string L1TRate_Offline::m_outputFile

private

Definition at line 86 of file L1TRate_Offline.h.

edm::ParameterSet L1TRate_Offline::m_parameters

private

Definition at line 126 of file L1TRate_Offline.h.

std::map<TString,MonitorElement*> L1TRate_Offline::m_PrescIndexVsLS

private

Definition at line 109 of file L1TRate_Offline.h.

std::map<TString,MonitorElement*> L1TRate_Offline::m_RecorLumiVsLS

private

Definition at line 114 of file L1TRate_Offline.h.

int L1TRate_Offline::m_refPrescaleSet

private

Definition at line 81 of file L1TRate_Offline.h.

edm::EDGetTokenT<Level1TriggerScalersCollection> L1TRate_Offline::m_scalersSource_L1TSCollection

private

Definition at line 122 of file L1TRate_Offline.h.

edm::EDGetTokenT<LumiScalersCollection> L1TRate_Offline::m_scalersSource_LSCollection

private

Definition at line 121 of file L1TRate_Offline.h.

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

private

Definition at line 99 of file L1TRate_Offline.h.

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

private

Definition at line 116 of file L1TRate_Offline.h.

bool L1TRate_Offline::m_verbose

private

Definition at line 78 of file L1TRate_Offline.h.

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

private

Definition at line 100 of file L1TRate_Offline.h.

std::map<TString,MonitorElement*> L1TRate_Offline::m_xSecObservedVsDelivLumi

private

Definition at line 103 of file L1TRate_Offline.h.

std::map<TString,MonitorElement*> L1TRate_Offline::m_xSecObservedVsLS

private

Definition at line 112 of file L1TRate_Offline.h.

std::map<TString,MonitorElement*> L1TRate_Offline::m_xSecObservedVsRecorLumi

private

Definition at line 104 of file L1TRate_Offline.h.

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

private

Definition at line 101 of file L1TRate_Offline.h.

Public Types

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation