L1TBPTX.h

Go to the documentation of this file.

#ifndef L1TBPTX_H
#define L1TBPTX_H

/*
 * \file L1TBPTX.h
 *
 * $Date: 2012/08/20 12:10:13 $
 * $Revision: 1.1 $
 * \author J. Pela
 *
*/

// system include files
#include <memory>
#include <unistd.h>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DQM/L1TMonitor/interface/L1TOMDSHelper.h"

#include <TString.h>

#include <iostream>
#include <fstream>
#include <string>
#include <vector>

class RateBuffer{
  
  public:
  
    void   fill         (int ls,bool isAlgo,int bit,double rate){    
      if(isAlgo){m_lsAlgoRate[std::pair<int,int>(ls,bit)]=rate;}
      else      {m_lsTechRate[std::pair<int,int>(ls,bit)]=rate;}
    }
    double getLSRate    (int ls,bool isAlgo,int bit,double rate){
      if(isAlgo){return m_lsAlgoRate[std::pair<int,int>(ls,bit)];}
      else      {return m_lsTechRate[std::pair<int,int>(ls,bit)];}
    }
    double getLSAlgoRate(int ls,int bit,double rate){
      return m_lsAlgoRate[std::pair<int,int>(ls,bit)];
    }
    double getLSTechRate(int ls,int bit,double rate){
      return m_lsTechRate[std::pair<int,int>(ls,bit)];
    }

  private:
  
    std::map<std::pair<int,int>,double> m_lsAlgoRate;
    std::map<std::pair<int,int>,double> m_lsTechRate;
  
};

class L1TBPTX : public edm::EDAnalyzer {

  public:

    enum BeamMode{
      NOMODE=1,
      SETUP=2,
      INJPILOT=3,
      INJINTR=4,
      INJNOMN=5,
      PRERAMP=6,
      RAMP=7,
      FLATTOP=8,
      SQUEEZE=9,
      ADJUST=10,
      STABLE=11,
      UNSTABLE=12,
      BEAMDUMP=13,
      RAMPDOWN=14,
      RECOVERY=15,
      INJDUMP=16,
      CIRCDUMP=17,
      ABORT=18,
      CYCLING=19,
      WBDUMP=20,
      NOBEAM=21
    };

    enum Errors{
      UNKNOWN                       = 1,
      WARNING_DB_CONN_FAILED        = 2,
      WARNING_DB_QUERY_FAILED       = 3,
      WARNING_DB_INCORRECT_NBUNCHES = 4,
      ERROR_UNABLE_RETRIVE_PRODUCT  = 5,
      ERROR_TRIGGERALIAS_NOTVALID   = 6,
      ERROR_LSBLOCK_NOTVALID        = 7
    };

  public:

    L1TBPTX(const edm::ParameterSet& ps);   // Constructor
    virtual ~L1TBPTX();                     // Destructor

  protected:

    void analyze (const edm::Event& e, const edm::EventSetup& c);  // Analyze
    void beginJob();                                               // BeginJob
    void endJob  (void);                                           // EndJob
    void beginRun(const edm::Run& run, const edm::EventSetup& iSetup);
    void endRun  (const edm::Run& run, const edm::EventSetup& iSetup);

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


  // Private Methods
  private:

    void getBeamConfOMDS();
    void doFractionInSync(bool iForce=false, bool iBad=false);
    void certifyLSBlock  (std::string iTrigger, int iInitLs, int iEndLs ,float iValue);

  // Variables
  private:

    DQMStore * dbe; // The DQM Service Handle

    edm::ParameterSet                      m_parameters;
    std::vector<edm::ParameterSet>         m_monitorBits;
    std::vector<edm::ParameterSet>         m_monitorRates;
    std::string                            m_outputFile;          // file name for ROOT ouput

    // bool
    bool                                   m_verbose;
    bool                                   m_currentLSValid;
    bool*                                  m_processedLS;

    // Int
    std::map<TString,int>                  m_effNumerator;
    std::map<TString,int>                  m_effDenominator;
    std::map<TString,int>                  m_missFireNumerator;
    std::map<TString,int>                  m_missFireDenominator;

    int                                    m_refPrescaleSet;
    int                                    m_currentPrescalesIndex;
    unsigned int                           m_currentLS;            // Current LS
    unsigned int                           m_currentGTLS;
    //unsigned int                         m_eventLS;
    unsigned int                           m_lhcFill;             //

    // Vectors
    BeamConfiguration                      m_beamConfig;          // Current Bunch Structure
    std::vector< std::pair<int,int> >      m_selAlgoBit;
    std::vector< std::pair<int,int> >      m_selTechBit;
    
    // Const Vectors                                          
    const std::vector<std::vector<int> >* m_prescaleFactorsAlgoTrig;
    const std::vector<std::vector<int> >* m_prescaleFactorsTechTrig;
                                           
    // Maps
    std::map<int,TString>                  m_algoBit_Alias;
    std::map<int,TString>                  m_techBit_Alias;


    std::map<TString,MonitorElement*>            m_meAlgoEfficiency;
    std::map<TString,MonitorElement*>            m_meAlgoMissFire;
    std::map<TString,MonitorElement*>            m_meTechEfficiency;
    std::map<TString,MonitorElement*>            m_meTechMissFire;

    
    std::map<std::pair<bool,int>,MonitorElement*> m_meRate;
    std::map<std::pair<bool,int>,double>          m_l1Rate;
    
    // MonitorElement
    MonitorElement*                        m_ErrorMonitor;

    // Input tags
    edm::InputTag                          m_scalersSource;       // Where to get L1 Scalers
    edm::InputTag                          m_l1GtEvmSource;
    edm::InputTag                          m_l1GtDataDaqInputTag;

};

#endif