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/MonitorElement.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