#include <L1TRate_Harvest.h>

Inheritance diagram for L1TRate_Harvest:

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

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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

virtual	~L1TRate_Harvest ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) 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)

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

CurrentProcessingContext const *	currentContext () const

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

std::map< TString, int >	m_algoBit

MonitorElement *	m_ErrorMonitor

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

edm::InputTag	m_l1GtDataDaqInputTag

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::InputTag	m_scalersSource

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
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 42 of file L1TRate_Harvest.h.

Member Enumeration Documentation

enum L1TRate_Harvest::Errors

Enumerator
UNKNOWN
WARNING_PY_MISSING_FIT

Definition at line 46 of file L1TRate_Harvest.h.

              {
     UNKNOWN                = 1,
     WARNING_PY_MISSING_FIT = 2
   };

Constructor & Destructor Documentation

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

Definition at line 41 of file L1TRate_Harvest.cc.

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

                                                        {
 
   m_maxNbins   = 2500; // Maximum LS for each run (for binning purposes)
   m_parameters = ps;
 
   // Mapping parameter input variables
   m_scalersSource       = m_parameters.getParameter         <InputTag>("inputTagScalersResults");
   m_l1GtDataDaqInputTag = 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"); 
 
   // Inicializing Variables
   dbe = NULL;
 
   if (ps.getUntrackedParameter < bool > ("dqmStore", false)) {
     dbe = Service < DQMStore > ().operator->();
     dbe->setVerbose(0);
   }
   
   // 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 = "";}
   
   if (dbe != NULL) {dbe->setCurrentFolder("L1T/L1TRate_Harvest");}
   
 }

L1TRate_Harvest::~L1TRate_Harvest ( )

virtual

Definition at line 90 of file L1TRate_Harvest.cc.

90 {}

Member Function Documentation

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

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 321 of file L1TRate_Harvest.cc.

                                                                                 {
 
 //   edm::Handle<L1GlobalTriggerReadoutRecord>   gtReadoutRecordData;
 //   edm::Handle<Level1TriggerScalersCollection> triggerScalers;
 //   edm::Handle<LumiScalersCollection>          colLScal;
  
 //   iEvent.getByLabel(m_l1GtDataDaqInputTag, gtReadoutRecordData);
 //   iEvent.getByLabel(m_scalersSource      , colLScal);
 //   iEvent.getByLabel(m_scalersSource      , 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_Harvest:] 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_Harvest:] 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_Harvest:] 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_Harvest::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 93 of file L1TRate_Harvest.cc.

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

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

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 231 of file L1TRate_Harvest.cc.

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

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 123 of file L1TRate_Harvest.cc.

                                                                             {
 
 //   if (m_verbose) {cout << "[L1TRate_Harvest:] 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_Harvest");
 
 //   // 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
 //   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 minRate = m_parameters.getParameter<double>("minRate");
 // //   double maxRate = m_parameters.getParameter<double>("maxRate");
 
 //   // 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_Harvest/TriggerRates"); // rate of the trigger...
 //     m_xSecVsInstLumi[tTrigger] = dbe->book1D(tCategory,
 //                       "Rate "+tTrigger+tErrorMessage,
 //                       m_maxNbins,0.5,m_maxNbins+0.5); 
 //     m_xSecVsInstLumi[tTrigger] ->setAxisTitle("Luminosity Section" ,1);
 //     m_xSecVsInstLumi[tTrigger] ->setAxisTitle("Rate [Hz]" ,2);
 //     //    m_xSecVsInstLumi[tTrigger] ->getTProfile()->GetListOfFunctions()->Add(tTestFunction);
 //     m_xSecVsInstLumi[tTrigger] ->getTH1()->SetMarkerStyle(23);
 
 //   }  
 
 }

void L1TRate_Harvest::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 109 of file L1TRate_Harvest.cc.

References gather_cfg::cout.

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

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 238 of file L1TRate_Harvest.cc.

                                                                                               {
 
 //   int eventLS = lumiBlock.id().luminosityBlock();  
 //   if (m_verbose) {cout << "[L1TRate_Harvest:] 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_xSecVsInstLumi[tTrigger]->getTH1()->FindBin(ls);
 //      //  m_xSecObservedToExpected[tTrigger]->setBinContent(ibin,AlgoXSec/TemplateFunctionValue);
 //             m_xSecVsInstLumi[tTrigger]->setBinContent(trigPrescale*trigRate,ibin);
   
 //      //            if(m_verbose){cout<<"[L1TRate_Harvest:] 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_Harvest:] Algo="<< tTrigger<< " XSec=Failed" << endl;}
 //    //}
 //         }
 //       }  
 //     }    
 //   }
 }

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 226 of file L1TRate_Harvest.cc.

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

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

private

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

private

Definition at line 425 of file L1TRate_Harvest.cc.

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

                                                                 {
 
   // error meaning
   bool noError = true;
   cout<<"WP01"<<endl;
   // Getting fit parameters
   std::vector<edm::ParameterSet>  m_fitParameters = ps.getParameter< vector<ParameterSet> >("fitParameters");
   cout<<"WP02"<<endl;
   double minInstantLuminosity = m_parameters.getParameter<double>("minInstantLuminosity");
   double maxInstantLuminosity = m_parameters.getParameter<double>("maxInstantLuminosity");
   cout<<"WP03"<<endl;
   // Getting rate fit parameters for all input triggers
   for(map<string,string>::const_iterator a=m_selectedTriggers.begin() ; a!=m_selectedTriggers.end() ; a++){
   cout<<"WP04-1"<<endl;
     string tTrigger = (*a).second;
 
     // If trigger name is defined we get the rate fit parameters 
     if(tTrigger != "Undefined"){
       cout<<"WP04-2"<<endl;
       bool foundFit = false;
 
       for(unsigned int b=0 ; b<m_fitParameters.size() ; b++){
     cout<<"WP05-1"<<endl;
         if(tTrigger == m_fitParameters[b].getParameter<string>("AlgoName")){
       cout<<"WP05-2"<<endl;
           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){
           cout<<"WP05-3"<<endl;
           noError = false;
 
           int eCount = m_ErrorMonitor->getTH1()->GetBinContent(WARNING_PY_MISSING_FIT);
           eCount++;
           m_ErrorMonitor->getTH1()->SetBinContent(WARNING_PY_MISSING_FIT,eCount);
 
           cout<<"WP05-4"<<endl;
         }
       }
     }
   }
 
    return noError;
 
 }

Member Data Documentation

DQMStore* L1TRate_Harvest::dbe

private

Definition at line 112 of file L1TRate_Harvest.h.

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

private

Definition at line 94 of file L1TRate_Harvest.h.

MonitorElement* L1TRate_Harvest::m_ErrorMonitor

private

Definition at line 109 of file L1TRate_Harvest.h.

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

private

Definition at line 95 of file L1TRate_Harvest.h.

edm::InputTag L1TRate_Harvest::m_l1GtDataDaqInputTag

private

Definition at line 103 of file L1TRate_Harvest.h.

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

private

Definition at line 88 of file L1TRate_Harvest.h.

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

private

Definition at line 92 of file L1TRate_Harvest.h.

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

private

Definition at line 91 of file L1TRate_Harvest.h.

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

private

Definition at line 93 of file L1TRate_Harvest.h.

int L1TRate_Harvest::m_lsShiftGTRates

private

Definition at line 82 of file L1TRate_Harvest.h.

int L1TRate_Harvest::m_maxNbins

private

Definition at line 81 of file L1TRate_Harvest.h.

std::string L1TRate_Harvest::m_outputFile

private

Definition at line 85 of file L1TRate_Harvest.h.

edm::ParameterSet L1TRate_Harvest::m_parameters

private

Definition at line 106 of file L1TRate_Harvest.h.

int L1TRate_Harvest::m_refPrescaleSet

private

Definition at line 80 of file L1TRate_Harvest.h.

edm::InputTag L1TRate_Harvest::m_scalersSource

private

Definition at line 102 of file L1TRate_Harvest.h.

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

private

Definition at line 96 of file L1TRate_Harvest.h.

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

private

Definition at line 99 of file L1TRate_Harvest.h.

bool L1TRate_Harvest::m_verbose

private

Definition at line 77 of file L1TRate_Harvest.h.

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

private

Definition at line 97 of file L1TRate_Harvest.h.

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

private

Definition at line 98 of file L1TRate_Harvest.h.

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation