#include <BPhysicsOniaDQM.h>

Inheritance diagram for BPhysicsOniaDQM:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	Get the analysis. More...

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

	BPhysicsOniaDQM (const edm::ParameterSet &)
	Constructor. More...

	~BPhysicsOniaDQM () override
	Destructor. More...

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Private Member Functions
float	computeMass (const math::XYZVector &vec1, const math::XYZVector &vec2)

bool	isMuonInAccept (const reco::Muon &recoMu)

bool	selGlobalMuon (const reco::Muon &recoMu)

bool	selTrackerMuon (const reco::Muon &recoMu)

Detailed Description

DQM offline for quarkonia

Author: S. Bolognesi, Eric - CERN

Definition at line 27 of file BPhysicsOniaDQM.h.

Constructor & Destructor Documentation

BPhysicsOniaDQM::BPhysicsOniaDQM ( const edm::ParameterSet & parameters )

Constructor.

Definition at line 29 of file BPhysicsOniaDQM.cc.

References edm::ParameterSet::getParameter(), and metname.

                                                                {
   // Muon Collection Label
   vertex_ = consumes<reco::VertexCollection>(
       parameters.getParameter<InputTag>("vertex"));
   theMuonCollectionLabel_ = consumes<reco::MuonCollection>(
       parameters.getParameter<InputTag>("MuonCollection"));
   lumiSummaryToken_ = consumes<LumiSummary, edm::InLumi>(
       parameters.getParameter<InputTag>("lumiSummary"));
 
   global_background = nullptr;
   diMuonMass_global = nullptr;
   tracker_background = nullptr;
   diMuonMass_tracker = nullptr;
   standalone_background = nullptr;
   diMuonMass_standalone = nullptr;
 
   glbSigCut = nullptr;
   glbSigNoCut = nullptr;
   glbBkgNoCut = nullptr;
   staSigCut = nullptr;
   staSigNoCut = nullptr;
   staBkgNoCut = nullptr;
   trkSigCut = nullptr;
   trkSigNoCut = nullptr;
   trkBkgNoCut = nullptr;
 
   metname = "oniaAnalyzer";
 }

BPhysicsOniaDQM::~BPhysicsOniaDQM ( )

override

Destructor.

Definition at line 58 of file BPhysicsOniaDQM.cc.

58 {}

Member Function Documentation

void BPhysicsOniaDQM::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Get the analysis.

Definition at line 103 of file BPhysicsOniaDQM.cc.

References edm::Event::getByToken(), muon::isGoodMuon(), edm::HandleBase::isValid(), LogTrace, metname, extraflags_cff::muons, and muon::TrackerMuonArbitrated.

                                                                            {
   LogTrace(metname) << "[BPhysicsOniaDQM] Analysis of event # ";
 
   // Take the STA muon container
   Handle<MuonCollection> muons;
   iEvent.getByToken(theMuonCollectionLabel_, muons);
 
   Handle<reco::VertexCollection> privtxs;
   iEvent.getByToken(vertex_, privtxs);
   VertexCollection::const_iterator privtx;
 
   if (privtxs->begin() != privtxs->end()) {
     privtx = privtxs->begin();
     RefVtx = privtx->position();
   } else {
     RefVtx.SetXYZ(0., 0., 0.);
   }
 
   if (muons.isValid()) {
     for (MuonCollection::const_iterator recoMu1 = muons->begin();
          recoMu1 != muons->end(); ++recoMu1) {
       // only loop over the remaining muons if recoMu1 is one of the following
       if (recoMu1->isGlobalMuon() || recoMu1->isTrackerMuon() ||
           recoMu1->isStandAloneMuon()) {
         for (MuonCollection::const_iterator recoMu2 = recoMu1 + 1;
              recoMu2 != muons->end(); ++recoMu2) {
           // fill the relevant histograms if recoMu2 satisfies one of the
           // following
           if (recoMu1->isGlobalMuon() && recoMu2->isGlobalMuon()) {
             math::XYZVector vec1 = recoMu1->globalTrack()->momentum();
             math::XYZVector vec2 = recoMu2->globalTrack()->momentum();
             float massJPsi = computeMass(vec1, vec2);
 
             // if opposite charges, fill glbSig, else fill glbBkg
             if (((*recoMu1).charge() * (*recoMu2).charge()) < 0) {
               if (diMuonMass_global != nullptr) {  // BPhysicsOniaDQM original one
                 diMuonMass_global->Fill(massJPsi);
               }
 
               if (glbSigNoCut != nullptr) {
                 glbSigNoCut->Fill(massJPsi);
                 if (selGlobalMuon(*recoMu1) && selGlobalMuon(*recoMu2)) {
                   if (glbSigCut != nullptr) glbSigCut->Fill(massJPsi);
                 }
               }
             } else {
               if (global_background != nullptr) {  // BPhysicsOniaDQM original one
                 global_background->Fill(massJPsi);
               }
 
               if (glbBkgNoCut != nullptr) {
                 glbBkgNoCut->Fill(massJPsi);
               }
             }
           }
 
           if (recoMu1->isStandAloneMuon() && recoMu2->isStandAloneMuon() &&
               fabs(recoMu1->outerTrack()->d0()) < 5 &&
               fabs(recoMu1->outerTrack()->dz()) < 30 &&
               fabs(recoMu2->outerTrack()->d0()) < 5 &&
               fabs(recoMu2->outerTrack()->dz()) < 30) {
             math::XYZVector vec1 = recoMu1->outerTrack()->momentum();
             math::XYZVector vec2 = recoMu2->outerTrack()->momentum();
             float massJPsi = computeMass(vec1, vec2);
 
             // if opposite charges, fill staSig, else fill staBkg
             if (((*recoMu1).charge() * (*recoMu2).charge()) < 0) {
               if (diMuonMass_standalone != nullptr) {
                 diMuonMass_standalone->Fill(massJPsi);
               }
 
               if (staSigNoCut != nullptr) {
                 staSigNoCut->Fill(massJPsi);
               }
             } else {
               if (standalone_background != nullptr) {
                 standalone_background->Fill(massJPsi);
               }
 
               if (staBkgNoCut != nullptr) {
                 staBkgNoCut->Fill(massJPsi);
               }
             }
           }
 
           if (recoMu1->isTrackerMuon() && recoMu2->isTrackerMuon() &&
               muon::isGoodMuon(*recoMu1, muon::TrackerMuonArbitrated) &&
               muon::isGoodMuon(*recoMu2, muon::TrackerMuonArbitrated)) {
             math::XYZVector vec1 = recoMu1->innerTrack()->momentum();
             math::XYZVector vec2 = recoMu2->innerTrack()->momentum();
             float massJPsi = computeMass(vec1, vec2);
 
             // if opposite charges, fill trkSig, else fill trkBkg
             if (((*recoMu1).charge() * (*recoMu2).charge()) < 0) {
               if (diMuonMass_tracker != nullptr) {
                 diMuonMass_tracker->Fill(massJPsi);
               }
 
               if (trkSigNoCut != nullptr) {
                 trkSigNoCut->Fill(massJPsi);
                 if (selTrackerMuon(*recoMu1) && selTrackerMuon(*recoMu2)) {
                   if (trkSigCut != nullptr) trkSigCut->Fill(massJPsi);
                 }
               }
             } else {
               if (tracker_background != nullptr) {
                 tracker_background->Fill(massJPsi);
               }
 
               if (trkBkgNoCut != nullptr) {
                 trkBkgNoCut->Fill(massJPsi);
               }
             }
           }
         }  // end of 2nd MuonCollection
       }    // end of GLB,STA,TRK muon check
     }      // end of 1st MuonCollection
   }        // Is this MuonCollection vaild?
 }

void BPhysicsOniaDQM::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

override

Definition at line 60 of file BPhysicsOniaDQM.cc.

References DQMStore::IBooker::book1D(), and DQMStore::IBooker::setCurrentFolder().

                                                             {
   iBooker.setCurrentFolder("Physics/BPhysics");  // Use folder with name of PAG
 
   global_background = iBooker.book1D(
       "global_background", "Same-sign global-global dimuon mass", 750, 0, 15);
   diMuonMass_global =
       iBooker.book1D("diMuonMass_global",
                      "Opposite-sign global-global dimuon mass", 750, 0, 15);
   tracker_background = iBooker.book1D(
       "tracker_background",
       "Same-sign tracker-tracker (arbitrated) dimuon mass", 750, 0, 15);
   diMuonMass_tracker = iBooker.book1D(
       "diMuonMass_tracker",
       "Opposite-sign tracker-tracker (arbitrated) dimuon mass", 750, 0, 15);
   standalone_background =
       iBooker.book1D("standalone_background",
                      "Same-sign standalone-standalone dimuon mass", 500, 0, 15);
   diMuonMass_standalone = iBooker.book1D(
       "diMuonMass_standalone",
       "Opposite-sign standalone-standalone dimuon mass", 500, 0, 15);
 
   glbSigCut = iBooker.book1D("glbSigCut", "Opposite-sign glb-glb dimuon mass",
                              650, 0, 130);
   glbSigNoCut = iBooker.book1D(
       "glbSigNoCut", "Opposite-sign glb-glb dimuon mass (no cut)", 650, 0, 130);
   glbBkgNoCut = iBooker.book1D(
       "glbBkgNoCut", "Same-sign glb-glb dimuon mass (no cut)", 650, 0, 130);
   staSigCut = iBooker.book1D("staSigCut", "Opposite-sign sta-sta dimuon mass",
                              430, 0, 129);
   staSigNoCut = iBooker.book1D(
       "staSigNoCut", "Opposite-sign sta-sta dimuon mass (no cut)", 430, 0, 129);
   staBkgNoCut = iBooker.book1D(
       "staBkgNoCut", "Same-sign sta-sta dimuon mass (no cut)", 430, 0, 129);
   trkSigCut = iBooker.book1D("trkSigCut", "Opposite-sign trk-trk dimuon mass",
                              650, 0, 130);
   trkSigNoCut = iBooker.book1D(
       "trkSigNoCut", "Opposite-sign trk-trk dimuon mass (no cut)", 650, 0, 130);
   trkBkgNoCut = iBooker.book1D(
       "trkBkgNoCutt", "Same-sign trk-trk dimuon mass (no cut)", 650, 0, 130);
 }

float BPhysicsOniaDQM::computeMass	(	const math::XYZVector &	vec1,
		const math::XYZVector &	vec2
	)

private

Definition at line 223 of file BPhysicsOniaDQM.cc.

References mathSSE::sqrt().

                                                               {
   // mass of muon
   float massMu = 0.10566;
   float eMu1 = -999;
   if (massMu * massMu + vec1.Mag2() > 0)
     eMu1 = sqrt(massMu * massMu + vec1.Mag2());
   float eMu2 = -999;
   if (massMu * massMu + vec2.Mag2() > 0)
     eMu2 = sqrt(massMu * massMu + vec2.Mag2());
 
   float pJPsi = -999;
   if ((vec1 + vec2).Mag2() > 0) pJPsi = sqrt((vec1 + vec2).Mag2());
   float eJPsi = eMu1 + eMu2;
 
   float massJPsi = -999;
   if ((eJPsi * eJPsi - pJPsi * pJPsi) > 0)
     massJPsi = sqrt(eJPsi * eJPsi - pJPsi * pJPsi);
 
   return massJPsi;
 }

bool BPhysicsOniaDQM::isMuonInAccept ( const reco::Muon & recoMu )

private

Definition at line 245 of file BPhysicsOniaDQM.cc.

References reco::LeafCandidate::eta(), reco::LeafCandidate::p(), and reco::LeafCandidate::pt().

                                                            {
   return (fabs(recoMu.eta()) < 2.4 &&
           ((fabs(recoMu.eta()) < 1.3 && recoMu.pt() > 3.3) ||
            (fabs(recoMu.eta()) > 1.3 && fabs(recoMu.eta()) < 2.2 &&
             recoMu.p() > 2.9) ||
            (fabs(recoMu.eta()) > 2.2 && recoMu.pt() > 0.8)));
 }

bool BPhysicsOniaDQM::selGlobalMuon ( const reco::Muon & recoMu )

private

Definition at line 253 of file BPhysicsOniaDQM.cc.

References reco::Muon::globalTrack(), reco::Muon::innerTrack(), reco::HitPattern::numberOfValidMuonHits(), AlCaHLTBitMon_ParallelJobs::p, reco::HitPattern::pixelLayersWithMeasurement(), and lumiQueryAPI::q.

                                                           {
   TrackRef iTrack = recoMu.innerTrack();
   const reco::HitPattern &p = iTrack->hitPattern();
 
   TrackRef gTrack = recoMu.globalTrack();
   const reco::HitPattern &q = gTrack->hitPattern();
 
   return (isMuonInAccept(recoMu) && iTrack->found() > 11 &&
           gTrack->chi2() / gTrack->ndof() < 20.0 &&
           q.numberOfValidMuonHits() > 0 &&
           iTrack->chi2() / iTrack->ndof() < 4.0 &&
           // recoMu.muonID("TrackerMuonArbitrated") &&
           // recoMu.muonID("TMLastStationAngTight") &&
           p.pixelLayersWithMeasurement() > 1 &&
           fabs(iTrack->dxy(RefVtx)) < 3.0 && fabs(iTrack->dz(RefVtx)) < 15.0);
 }

bool BPhysicsOniaDQM::selTrackerMuon ( const reco::Muon & recoMu )

private

Definition at line 270 of file BPhysicsOniaDQM.cc.

References reco::Muon::innerTrack(), AlCaHLTBitMon_ParallelJobs::p, and reco::HitPattern::pixelLayersWithMeasurement().

                                                            {
   TrackRef iTrack = recoMu.innerTrack();
   const reco::HitPattern &p = iTrack->hitPattern();
 
   return (isMuonInAccept(recoMu) && iTrack->found() > 11 &&
           iTrack->chi2() / iTrack->ndof() < 4.0 &&
           // recoMu.muonID("TrackerMuonArbitrated") &&
           // recoMu.muonID("TMLastStationAngTight") &&
           p.pixelLayersWithMeasurement() > 1 &&
           fabs(iTrack->dxy(RefVtx)) < 3.0 && fabs(iTrack->dz(RefVtx)) < 15.0);
 }

Member Data Documentation

MonitorElement* BPhysicsOniaDQM::diMuonMass_global

private

Definition at line 56 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::diMuonMass_standalone

private

Definition at line 58 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::diMuonMass_tracker

private

Definition at line 57 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::glbBkgNoCut

private

Definition at line 65 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::glbSigCut

private

Definition at line 63 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::glbSigNoCut

private

Definition at line 64 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::global_background

private

Definition at line 59 of file BPhysicsOniaDQM.h.

edm::EDGetTokenT<LumiSummary> BPhysicsOniaDQM::lumiSummaryToken_

private

Definition at line 50 of file BPhysicsOniaDQM.h.

std::string BPhysicsOniaDQM::metname

private

Definition at line 53 of file BPhysicsOniaDQM.h.

math::XYZPoint BPhysicsOniaDQM::RefVtx

private

Definition at line 73 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::staBkgNoCut

private

Definition at line 68 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::standalone_background

private

Definition at line 61 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::staSigCut

private

Definition at line 66 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::staSigNoCut

private

Definition at line 67 of file BPhysicsOniaDQM.h.

edm::EDGetTokenT<reco::MuonCollection> BPhysicsOniaDQM::theMuonCollectionLabel_

private

Definition at line 49 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::tracker_background

private

Definition at line 60 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::trkBkgNoCut

private

Definition at line 71 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::trkSigCut

private

Definition at line 69 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::trkSigNoCut

private

Definition at line 70 of file BPhysicsOniaDQM.h.

edm::EDGetTokenT<reco::VertexCollection> BPhysicsOniaDQM::vertex_

private

Definition at line 47 of file BPhysicsOniaDQM.h.

Private Attributes
MonitorElement *	diMuonMass_global

MonitorElement *	diMuonMass_standalone

MonitorElement *	diMuonMass_tracker

MonitorElement *	glbBkgNoCut

MonitorElement *	glbSigCut

MonitorElement *	glbSigNoCut

MonitorElement *	global_background

edm::EDGetTokenT< LumiSummary >	lumiSummaryToken_

std::string	metname

math::XYZPoint	RefVtx

MonitorElement *	staBkgNoCut

MonitorElement *	standalone_background

MonitorElement *	staSigCut

MonitorElement *	staSigNoCut

edm::EDGetTokenT< reco::MuonCollection >	theMuonCollectionLabel_

MonitorElement *	tracker_background

MonitorElement *	trkBkgNoCut

MonitorElement *	trkSigCut

MonitorElement *	trkSigNoCut

edm::EDGetTokenT< reco::VertexCollection >	vertex_

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation