diff --git a/DQMOffline/Alignment/interface/DiMuonVertexMonitor.h b/DQMOffline/Alignment/interface/DiMuonVertexMonitor.h index 567f538a80db9..f79003a9f051c 100644 --- a/DQMOffline/Alignment/interface/DiMuonVertexMonitor.h +++ b/DQMOffline/Alignment/interface/DiMuonVertexMonitor.h @@ -76,6 +76,8 @@ class DiMuonVertexMonitor : public DQMEDAnalyzer { MonitorElement *hCosPhi3D_; MonitorElement *hCosPhiInv_; MonitorElement *hCosPhiInv3D_; + MonitorElement *hCosPhiUnbalance_; + MonitorElement *hCosPhi3DUnbalance_; MonitorElement *hInvMass_; MonitorElement *hCutFlow_; diff --git a/DQMOffline/Alignment/src/DiMuonVertexMonitor.cc b/DQMOffline/Alignment/src/DiMuonVertexMonitor.cc index b55a78a39eadb..ad9979c813e8f 100644 --- a/DQMOffline/Alignment/src/DiMuonVertexMonitor.cc +++ b/DQMOffline/Alignment/src/DiMuonVertexMonitor.cc @@ -91,6 +91,8 @@ void DiMuonVertexMonitor::bookHistograms(DQMStore::IBooker& iBooker, edm::Run co hCosPhi3D_ = iBooker.book1D("CosPhi3D", fmt::sprintf("%s;cos(#phi_{3D});%s", histTit, ps), 50, -1., 1.); hCosPhiInv_ = iBooker.book1D("CosPhiInv", fmt::sprintf("%s;inverted cos(#phi_{xy});%s", histTit, ps), 50, -1., 1.); hCosPhiInv3D_ = iBooker.book1D("CosPhiInv3D", fmt::sprintf("%s;inverted cos(#phi_{3D});%s", histTit, ps), 50, -1., 1.); + hCosPhiUnbalance_ = iBooker.book1D("CosPhiUnbalance", fmt::sprintf("%s;cos(#phi_{xy}) unbalance;#Delta%s", histTit, ps), 50, -1.,1.); + hCosPhi3DUnbalance_ = iBooker.book1D("CosPhi3DUnbalance", fmt::sprintf("%s;cos(#phi_{3D}) unbalance;#Delta%s", histTit, ps), 50, -1., 1.); hdxy_ = iBooker.book1D("dxy", fmt::sprintf("%s;muon track d_{xy}(PV) [#mum];muon tracks", histTit), 150, -300, 300); hdz_ = iBooker.book1D("dz", fmt::sprintf("%s;muon track d_{z}(PV) [#mum];muon tracks", histTit), 150, -300, 300); @@ -252,6 +254,12 @@ void DiMuonVertexMonitor::analyze(const edm::Event& iEvent, const edm::EventSetu // inverted hCosPhiInv_->Fill(-cosphi); hCosPhiInv3D_->Fill(-cosphi3D); + + // unbalance + hCosPhiUnbalance_->Fill(cosphi, 1.); + hCosPhiUnbalance_->Fill(-cosphi, -1.); + hCosPhi3DUnbalance_->Fill(cosphi3D, 1.); + hCosPhi3DUnbalance_->Fill(-cosphi3D, -1.); } } else { edm::LogWarning("DiMuonVertexMonitor") << "hardest primary vertex in the event is not valid!";