BPhysicsOniaDQM.h

Go to the documentation of this file.

#ifndef BPhysicsOniaDQM_H
#define BPhysicsOniaDQM_H

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/Luminosity/interface/LumiSummary.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"

#include <string>
#include <cmath>
#include <map>

class DQMStore;
class MonitorElement;

class BPhysicsOniaDQM : public edm::EDAnalyzer {
 public:
  BPhysicsOniaDQM(const edm::ParameterSet&);

  virtual ~BPhysicsOniaDQM();

  void beginJob();

  void analyze(const edm::Event&, const edm::EventSetup&);
  void beginLuminosityBlock(const edm::LuminosityBlock& lumiBlock,
                            const edm::EventSetup& iSetup);
  void endLuminosityBlock(const edm::LuminosityBlock& lumiBlock,
                          const edm::EventSetup& iSetup);
  void beginRun(const edm::Run& iRun, const edm::EventSetup& iSetup);
  void endRun(const edm::Run& iRun, const edm::EventSetup& iSetup);
  void endJob(void);

 private:
  float computeMass(const math::XYZVector& vec1, const math::XYZVector& vec2);
  bool isMuonInAccept(const reco::Muon& recoMu);
  bool selGlobalMuon(const reco::Muon& recoMu);
  bool selTrackerMuon(const reco::Muon& recoMu);

  // ----------member data ---------------------------

  DQMStore* theDbe;

  edm::EDGetTokenT<reco::VertexCollection> vertex_;
  // Muon Label
  edm::EDGetTokenT<reco::MuonCollection> theMuonCollectionLabel_;
  edm::EDGetTokenT<LumiSummary> lumiSummaryToken_;

  // Switch for verbosity
  std::string metname;

  // The histos
  MonitorElement* diMuonMass_global;
  MonitorElement* diMuonMass_tracker;
  MonitorElement* diMuonMass_standalone;
  MonitorElement* global_background;
  MonitorElement* tracker_background;
  MonitorElement* standalone_background;

  MonitorElement* glbSigCut;
  MonitorElement* glbSigNoCut;
  MonitorElement* glbBkgNoCut;
  MonitorElement* staSigCut;
  MonitorElement* staSigNoCut;
  MonitorElement* staBkgNoCut;
  MonitorElement* trkSigCut;
  MonitorElement* trkSigNoCut;
  MonitorElement* trkBkgNoCut;

  //   MonitorElement* JPsiGlbYdLumi;
  //   MonitorElement* JPsiStaYdLumi;
  //   MonitorElement* JPsiTrkYdLumi;

  // Yield of dimuon objects
  int jpsiGlbSigPerLS;
  int jpsiStaSigPerLS;
  int jpsiTrkSigPerLS;
  std::map<int, int> jpsiGlbSig;
  std::map<int, int> jpsiStaSig;
  std::map<int, int> jpsiTrkSig;

  math::XYZPoint RefVtx;
};
#endif

// Local Variables:
// show-trailing-whitespace: t
// truncate-lines: t
// End: