RecoMuonValidator.h

Go to the documentation of this file.

#ifndef Validation_RecoMuon_RecoMuonValidator_H
#define Validation_RecoMuon_RecoMuonValidator_H

#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/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include "SimTracker/Common/interface/TrackingParticleSelector.h"

#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "SimDataFormats/TrackingAnalysis/interface/TrackingParticle.h"

#include "SimMuon/MCTruth/interface/MuonAssociatorByHits.h"
#include <DQMServices/Core/interface/DQMEDAnalyzer.h>

// for selection cut
#include "CommonTools/Utils/interface/StringCutObjectSelector.h"

class DQMStore;
class MonitorElement;
class MuonServiceProxy;
class TrackAssociatorBase;

class RecoMuonValidator : public thread_unsafe::DQMEDAnalyzer
{
 public:
  RecoMuonValidator(const edm::ParameterSet& pset);
  ~RecoMuonValidator();

  virtual void dqmBeginRun(const edm::Run&, const edm::EventSetup& eventSetup);
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  virtual void endRun();
  virtual void analyze(const edm::Event& event, const edm::EventSetup& eventSetup);
  virtual int countMuonHits(const reco::Track& track) const;
  virtual int countTrackerHits(const reco::Track& track) const;

 protected:
  unsigned int verbose_;

  edm::InputTag simLabel_;
  edm::InputTag muonLabel_;
  std::string muonSelection_;
  edm::EDGetTokenT<TrackingParticleCollection> simToken_;
  edm::EDGetTokenT<edm::View<reco::Muon> > muonToken_;

  edm::InputTag muAssocLabel_;
  const MuonAssociatorByHits * assoByHits;
  //  edm::EDGetTokenT<> muAssocToken_;
  
  edm::InputTag beamspotLabel_;
  edm::InputTag primvertexLabel_;
  edm::EDGetTokenT<reco::BeamSpot> beamspotToken_;
  edm::EDGetTokenT<reco::VertexCollection> primvertexToken_;

  std::string outputFileName_;
  std::string subDir_;
  std::string subsystemname_;
  edm::ParameterSet pset;

  MuonServiceProxy * theMuonService;
  DQMStore * dbe_;
  
  bool doAbsEta_;
  bool doAssoc_;
  bool usePFMuon_;


  TrackingParticleSelector tpSelector_;

  // Track to use
  MuonAssociatorByHits::MuonTrackType trackType_;

  struct MuonME;
  MuonME * muonME_;

  struct CommonME;
  CommonME * commonME_;

//
//struct for histogram dimensions
//
struct HistoDimensions {
  //p
  unsigned int nBinP;
  double minP, maxP;
  //pt
  unsigned int nBinPt;
  double minPt, maxPt;
  //if abs eta
  bool doAbsEta;
  //eta
  unsigned int nBinEta;
  double minEta, maxEta;
  //phi
  unsigned int nBinPhi;
  double minPhi, maxPhi;
  //dxy
  unsigned int nBinDxy;
  double minDxy, maxDxy;
  //dz
  unsigned int nBinDz; 
  double minDz, maxDz;
  //pulls
  unsigned int nBinPull;
  double wPull;
  //resolustions
  unsigned int nBinErr;
  double minErrP, maxErrP;
  double minErrPt, maxErrPt;
  double minErrQPt, maxErrQPt;
  double minErrEta, maxErrEta;
  double minErrPhi, maxErrPhi;
  double minErrDxy, maxErrDxy;
  double minErrDz, maxErrDz;
  //track multiplicities
  unsigned int nTrks, nAssoc;
  unsigned int nDof;
  // for PF muons
  bool usePFMuon;
};

  HistoDimensions hDim;
 private:
  StringCutObjectSelector<reco::Muon> selector_;
  bool wantTightMuon_;

};

#endif
/* vim:set ts=2 sts=2 sw=2 expandtab: */