MuonIsolationDQM.h

Go to the documentation of this file.

// -*- C++ -*-
// MuonIsolationDQM.h
// Package:    Muon Isolation DQM
// Class:      MuonIsolationDQM
// 
/*
  
Description: Muon Isolation DQM class

NOTE: The static member variable declarations *should* include the key word "static", but 
I haven't found an elegant way to initalize the vectors.  Static primatives (e.g. int, 
float, ...) and simple static objects are easy to initialze.  Outside of the class 
decleration, you would write
    
int MuonIsolationDQM::CONST_INT = 5;
FooType MuonIsolationDQM::CONST_FOOT = Foo(constructor_argument);
    
but you can't do this if you want to, say, initalize a std::vector with a bunch of 
different values.  So, you can't make them static and you have to initialize them using 
a member method.  To keep it consistent, I've just initialized them all in the same 
method, even the simple types.
*/

//
// Original Author:  "C. Jess Riedel", UC Santa Barbara
//         Created:  Tue Jul 17 15:58:24 CDT 2007
//

//Base class
#include "FWCore/Framework/interface/EDAnalyzer.h"

//Member types
#include "FWCore/Utilities/interface/InputTag.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

//Other include files
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/RecoCandidate/interface/IsoDeposit.h"
#include "DataFormats/RecoCandidate/interface/IsoDepositFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"

#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"

//----------------------------------------

//Forward declarations
class TH1;
class TH1I;
class TH1D;
class TH2;
class TProfile;

//------------------------------------------
//  Class Declaration: MuonIsolationDQM
//--------------------------------------
class MuonIsolationDQM : public edm::EDAnalyzer {
  //---------namespace and typedefs--------------
  typedef edm::View<reco::Muon>::const_iterator MuonIterator;
  typedef edm::RefToBase<reco::Muon> MuonBaseRef;
  typedef edm::Handle<reco::IsoDepositMap> MuIsoDepHandle;
  typedef const reco::IsoDeposit MuIsoDepRef;
  
public:
  //---------methods----------------------------
  explicit MuonIsolationDQM(const edm::ParameterSet&);
  ~MuonIsolationDQM();
  
  
private:
  //---------methods----------------------------
  virtual void beginJob(void) ;
  virtual void beginRun(void) ;
  virtual void analyze(const edm::Event&, const edm::EventSetup&);
  virtual void endJob() ;
  void InitStatics();
  void RecordData(const reco::Muon& muon);//Fills Histograms with info from single muo
  //  void doPFIsoPlots(MuonIterator muon); //Fills Histograms with PF info from single muo (only for GLB)
  void InitHistos();//adds title, bin information to member histograms
  void FillHistos(int);//Fills histograms with data
  void FillNVtxHistos(int);
  void NormalizeHistos(); //Normalize to number of muons

  //----- helper methods
  int  GetNVtxBin(int); 
  TH1* GetTH1FromMonitorElement(MonitorElement* me);

  //----------Static Variables---------------
  
  //Collection labels
  edm::EDGetTokenT<reco::VertexCollection> theVertexCollectionLabel_;
  edm::EDGetTokenT<reco::MuonCollection>   theMuonCollectionLabel_;

  //root file name
  std::string rootfilename;
  // Directories within the rootfile
  std::string dirName;

  //Histogram parameters
  static const int NUM_VARS      = 48; // looking at R03 and R05.  Total of 54 histos.
  static const int NUM_VARS_2D   = 10; // looking only at R03.  Total of 8 TH2F. 
  static const int NUM_VARS_NVTX = 6 ;
  
  double L_BIN_WIDTH;//large bins
  double S_BIN_WIDTH;//small bins
  int LOG_BINNING_ENABLED;//pseudo log binning for profile plots
  int NUM_LOG_BINS;
  double LOG_BINNING_RATIO;
  bool requireGLBMuon;
  bool requireSTAMuon;
  bool requireTRKMuon;

  std::string title_sam;
  std::string title_cone;
  //  std::string title_cd;
  
  std::vector<std::string> main_titles;//[NUM_VARS]
  std::vector<std::string> axis_titles;//[NUM_VARS]
  std::vector<std::string> names;//[NUM_VARS]
  std::vector< std::vector<double> > param;//[NUM_VARS][3]
  std::vector<int> isContinuous;//[NUM_VARS]
  
  std::vector<std::string> titles_2D;     //[NUM_VARS]
  std::vector<std::string> names_2D;      //[NUM_VARS]

  std::vector<std::string> main_titles_NVtxs;
  std::vector<std::string> names_NVtxs;
  std::vector<std::string> axis_titles_NVtxs;
  //---------------Dynamic Variables---------------------
  
  //MonitorElement
  DQMStore* dbe;
  
  //The Data
  double theData[NUM_VARS];
  double theData2D[NUM_VARS_2D];
  double theDataNVtx[NUM_VARS_NVTX];

  //Histograms
  MonitorElement* h_nMuons;
  std::vector<MonitorElement*> h_1D;     //[NUM_VARS]
  std::vector<MonitorElement*> h_2D;     //[NUM_VARS_2D]
  std::vector<MonitorElement*> h_1D_NVTX;//[NUM_VARS_NVTX]
  
  //  std::vector<MonitorElement*> cd_plots;//[NUM_VARS]
  
  //Counters
  int nEvents;
  int nSTAMuons;
  int nGLBMuons;
  int nTRKMuons;
  
  //enums for monitorElement
  enum {NOAXIS,XAXIS,YAXIS,ZAXIS};
};