MuonRecoAnalyzer.h

Go to the documentation of this file.

#ifndef MuonRecoAnalyzer_H
#define MuonRecoAnalyzer_H

#include <memory>
#include <fstream>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"

#include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"

#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/Scalers/interface/DcsStatus.h"

class MuonRecoAnalyzer : public DQMEDAnalyzer {
public:
  MuonRecoAnalyzer(const edm::ParameterSet&);

  ~MuonRecoAnalyzer() override;

  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;

  //calculate residual & pull:
  void GetRes(reco::TrackRef t1, reco::TrackRef t2, std::string par, float& res, float& pull);

  //Functions needed by the SoftMuon MVA monitoring
  double getDeltaR(reco::Track track1, reco::Track track2);

  int getPv(int tidx, const reco::VertexCollection* vc);

private:
  // ----------member data ---------------------------
  MuonServiceProxy* theService;
  edm::ParameterSet parameters;

  edm::EDGetTokenT<edm::View<reco::Muon> > theMuonCollectionLabel_;
  edm::EDGetTokenT<reco::VertexCollection> theVertexLabel_;
  edm::EDGetTokenT<reco::BeamSpot> theBeamSpotLabel_;
  edm::EDGetTokenT<DcsStatusCollection> dcsStatusCollection_;

  // Switch for verbosity
  std::string metname;
  bool doMVA;

  //histo binning parameters
  int etaBin;
  double etaMin;
  double etaMax;

  int thetaBin;
  double thetaMin;
  double thetaMax;

  int phiBin;
  double phiMin;
  double phiMax;

  int chi2Bin;
  double chi2Min;
  double chi2Max;

  int pBin;
  double pMin;
  double pMax;

  int ptBin;
  double ptMin;
  double ptMax;

  int pResBin;
  double pResMin;
  double pResMax;

  int rhBin;
  double rhMin;
  double rhMax;

  int tunePBin;
  double tunePMin;
  double tunePMax;

  //the histos
  MonitorElement* muReco;
  // global muon
  std::vector<MonitorElement*> etaGlbTrack;
  std::vector<MonitorElement*> etaResolution;
  std::vector<MonitorElement*> thetaGlbTrack;
  std::vector<MonitorElement*> thetaResolution;
  std::vector<MonitorElement*> phiGlbTrack;
  std::vector<MonitorElement*> phiResolution;
  std::vector<MonitorElement*> chi2OvDFGlbTrack;
  std::vector<MonitorElement*> probchi2GlbTrack;
  std::vector<MonitorElement*> pGlbTrack;
  std::vector<MonitorElement*> ptGlbTrack;
  std::vector<MonitorElement*> qGlbTrack;
  std::vector<MonitorElement*> qOverpResolution;
  std::vector<MonitorElement*> qOverptResolution;
  std::vector<MonitorElement*> oneOverpResolution;
  std::vector<MonitorElement*> oneOverptResolution;
  std::vector<MonitorElement*> rhAnalysis;
  std::vector<MonitorElement*> muVStkSytemRotation;
  std::vector<MonitorElement*> phiVsetaGlbTrack;
  std::vector<MonitorElement*> phiVsetaGlbTrack_badlumi;

  //Soft MVA Muon
  MonitorElement* ptSoftMuonMVA;
  MonitorElement* deltaRSoftMuonMVA;
  MonitorElement* gNchi2SoftMuonMVA;
  MonitorElement* vMuHitsSoftMuonMVA;
  MonitorElement* mNuStationsSoftMuonMVA;
  MonitorElement* dxyRefSoftMuonMVA;
  MonitorElement* dzRefSoftMuonMVA;
  MonitorElement* LWHSoftMuonMVA;
  MonitorElement* valPixHitsSoftMuonMVA;
  MonitorElement* innerChi2SoftMuonMVA;
  MonitorElement* outerChi2SoftMuonMVA;
  MonitorElement* iValFracSoftMuonMVA;
  MonitorElement* segCompSoftMuonMVA;
  MonitorElement* chi2LocMomSoftMuonMVA;
  MonitorElement* chi2LocPosSoftMuonMVA;
  MonitorElement* glbTrackTailProbSoftMuonMVA;
  MonitorElement* NTrkVHitsSoftMuonMVA;
  MonitorElement* kinkFinderSoftMuonMVA;
  MonitorElement* vRPChitsSoftMuonMVA;
  MonitorElement* glbKinkFinderSoftMuonMVA;
  MonitorElement* glbKinkFinderLogSoftMuonMVA;
  MonitorElement* staRelChi2SoftMuonMVA;
  MonitorElement* glbDeltaEtaPhiSoftMuonMVA;
  MonitorElement* trkRelChi2SoftMuonMVA;
  MonitorElement* vDThitsSoftMuonMVA;
  MonitorElement* vCSChitsSoftMuonMVA;
  MonitorElement* timeAtIpInOutSoftMuonMVA;
  MonitorElement* timeAtIpInOutErrSoftMuonMVA;
  MonitorElement* getMuonHitsPerStationSoftMuonMVA;
  MonitorElement* QprodSoftMuonMVA;

  MonitorElement* tunePResolution;

  MonitorElement* etaPull;
  MonitorElement* thetaPull;
  MonitorElement* phiPull;
  MonitorElement* qOverpPull;
  MonitorElement* qOverptPull;
  MonitorElement* oneOverpPull;
  MonitorElement* oneOverptPull;

  // tracker muon
  MonitorElement* etaTrack;
  MonitorElement* thetaTrack;
  MonitorElement* phiTrack;
  MonitorElement* chi2OvDFTrack;
  MonitorElement* probchi2Track;
  MonitorElement* pTrack;
  MonitorElement* ptTrack;
  MonitorElement* qTrack;
  // sta muon
  MonitorElement* etaStaTrack;
  MonitorElement* thetaStaTrack;
  MonitorElement* phiStaTrack;
  MonitorElement* chi2OvDFStaTrack;
  MonitorElement* probchi2StaTrack;
  MonitorElement* pStaTrack;
  MonitorElement* ptStaTrack;
  MonitorElement* qStaTrack;
  // efficiency
  std::vector<MonitorElement*> etaEfficiency;
  std::vector<MonitorElement*> phiEfficiency;

  bool IsminiAOD;
  std::string theFolder;
};
#endif