L1TRate_Offline Class Reference

#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 21 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 63 of file L1TRate_Offline.cc.

{}

Member Function Documentation

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

protectedvirtual

Implements edm::EDAnalyzer.

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

      if(gtFdlVectorData.size() != 0)
        {
      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 66 of file L1TRate_Offline.cc.

References gather_cfg::cout, AlcaSiStripGainsHarvester_cff::DQMStore, and cppFunctionSkipper::operator.

                                  {

  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 254 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 96 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");

  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_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 82 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 261 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 248 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 532 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.