BeamMonitor.h

Go to the documentation of this file.

#ifndef BeamMonitor_H
#define BeamMonitor_H

// C++
#include <string>
// CMS
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "RecoVertex/BeamSpotProducer/interface/BSTrkParameters.h"
#include "RecoVertex/BeamSpotProducer/interface/BeamFitter.h"
#include <fstream>


//
// class declaration
//

class BeamMonitor : public edm::EDAnalyzer {
  public:

    BeamMonitor( const edm::ParameterSet& );
    ~BeamMonitor();

  protected:

    // BeginJob
    void beginJob();

    // BeginRun
    void beginRun(const edm::Run& r, const edm::EventSetup& c);

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

    void beginLuminosityBlock(const edm::LuminosityBlock& lumiSeg,
        const edm::EventSetup& context) ;

    void endLuminosityBlock(const edm::LuminosityBlock& lumiSeg,
        const edm::EventSetup& c);
    // EndRun
    void endRun(const edm::Run& r, const edm::EventSetup& c);
    // Endjob
    void endJob(const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& c);

  private:

    void FitAndFill(const edm::LuminosityBlock& lumiSeg,int&,int&,int&);
    void RestartFitting();
    void scrollTH1(TH1 *, std::time_t);
    bool testScroll(std::time_t &, std::time_t &);
    const char * formatFitTime( const std::time_t &);
    edm::ParameterSet parameters_;
    std::string monitorName_;
    edm::InputTag bsSrc_; // beam spot
    edm::InputTag tracksLabel_;
    edm::InputTag pvSrc_; // primary vertex
    edm::InputTag hltSrc_;//hlt collection

    int fitNLumi_;
    int fitPVNLumi_;
    int resetFitNLumi_;
    int resetPVNLumi_;
    int intervalInSec_;
    bool debug_;
    bool onlineMode_;
    std::vector<std::string> jetTrigger_;

    DQMStore* dbe_;
    BeamFitter * theBeamFitter;

    int countEvt_;       //counter
    int countLumi_;      //counter
    int beginLumiOfBSFit_;
    int endLumiOfBSFit_;
    int beginLumiOfPVFit_;
    int endLumiOfPVFit_;
    int lastlumi_; // previous LS processed
    int nextlumi_; // next LS of Fit
    std::time_t refBStime[2];
    std::time_t refPVtime[2];
    unsigned int nthBSTrk_;
    int nFitElements_;
    int nFits_;
    double deltaSigCut_;
    unsigned int min_Ntrks_;
    double maxZ_;
    unsigned int minNrVertices_;
    double minVtxNdf_;
    double minVtxWgt_;

    bool resetHistos_;
    bool StartAverage_;
    int firstAverageFit_;
    int countGapLumi_;

    bool processed_;
   
    // ----------member data ---------------------------

    //   std::vector<BSTrkParameters> fBSvector;
    reco::BeamSpot refBS;
    reco::BeamSpot preBS;

    // MonitorElements:
    MonitorElement * h_nTrk_lumi;
    MonitorElement * h_d0_phi0;
    MonitorElement * h_trk_z0;
    MonitorElement * h_vx_vy;
    MonitorElement * h_vx_dz;
    MonitorElement * h_vy_dz;
    MonitorElement * h_trkPt;
    MonitorElement * h_trkVz;
    MonitorElement * fitResults;
    MonitorElement * h_x0;
    MonitorElement * h_y0;
    MonitorElement * h_z0;
    MonitorElement * h_sigmaX0;
    MonitorElement * h_sigmaY0;
    MonitorElement * h_sigmaZ0;
    MonitorElement * h_nVtx;
    MonitorElement * h_nVtx_st;
    MonitorElement * h_PVx[2];
    MonitorElement * h_PVy[2];
    MonitorElement * h_PVz[2];
    MonitorElement * h_PVxz;
    MonitorElement * h_PVyz;
    MonitorElement * pvResults;
    std::map<TString, MonitorElement*> hs;

    // The histo of the primary vertex for  DQM gui
    std::map<int, std::vector<float> > mapPVx,mapPVy,mapPVz;
    std::map<int, std::vector<int> > mapNPV;
    //keep track of beginLuminosity block and time
    std::map<int, int> mapBeginBSLS,mapBeginPVLS;
    std::map<int, std::time_t> mapBeginBSTime, mapBeginPVTime;
    //these maps are needed to keep track of size of theBSvector and pvStore
    std::map<int, std::size_t> mapLSBSTrkSize;
    std::map<int, size_t>mapLSPVStoreSize;
    //to correct the cutFlot Table
    std::map<int, TH1F> mapLSCF;

    // Summary:
    Float_t reportSummary_;
    Float_t summarySum_;
    Float_t summaryContent_[3];
    MonitorElement * reportSummary;
    MonitorElement * reportSummaryContents[3];
    MonitorElement * reportSummaryMap;
    MonitorElement * cutFlowTable;
    // variables for beam fit

    //
    std::time_t tmpTime;
    std::time_t startTime;
    std::time_t refTime;
    edm::TimeValue_t ftimestamp;
    // Current run
    int frun;
    int lastNZbin; // last non zero bin of time histos
};

#endif