#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 DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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)

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

DQM offline for quarkonia

Author: S. Bolognesi, Eric - CERN

Definition at line 25 of file BPhysicsOniaDQM.h.

Constructor & Destructor Documentation

◆ BPhysicsOniaDQM()

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

Constructor.

Definition at line 28 of file BPhysicsOniaDQM.cc.

                                                                {
   // 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";
 }

References metname.

◆ ~BPhysicsOniaDQM()

BPhysicsOniaDQM::~BPhysicsOniaDQM ( )

override

Destructor.

Definition at line 54 of file BPhysicsOniaDQM.cc.

54 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Get the analysis.

Reimplemented from DQMEDAnalyzer.

Definition at line 81 of file BPhysicsOniaDQM.cc.

                                                                            {
   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?
 }

References iEvent, muon::isGoodMuon(), LogTrace, metname, PDWG_BPHSkim_cff::muons, and muon::TrackerMuonArbitrated.

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 56 of file BPhysicsOniaDQM.cc.

                                                                                                     {
   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);
 }

References dqm::implementation::IBooker::book1D(), and dqm::implementation::NavigatorBase::setCurrentFolder().

◆ computeMass()

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

private

Definition at line 198 of file BPhysicsOniaDQM.cc.

                                                                                        {
   // 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;
 }

References mathSSE::sqrt().

◆ isMuonInAccept()

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

private

Definition at line 220 of file BPhysicsOniaDQM.cc.

                                                            {
   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)));
 }

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

◆ selGlobalMuon()

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

private

Definition at line 226 of file BPhysicsOniaDQM.cc.

                                                           {
   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);
 }

References reco::Muon::globalTrack(), reco::Muon::innerTrack(), AlCaHLTBitMon_ParallelJobs::p, and submitPVResolutionJobs::q.

◆ selTrackerMuon()

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

private

Definition at line 240 of file BPhysicsOniaDQM.cc.

                                                            {
   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);
 }

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

Member Data Documentation

◆ diMuonMass_global

MonitorElement* BPhysicsOniaDQM::diMuonMass_global

private

Definition at line 53 of file BPhysicsOniaDQM.h.

◆ diMuonMass_standalone

MonitorElement* BPhysicsOniaDQM::diMuonMass_standalone

private

Definition at line 55 of file BPhysicsOniaDQM.h.

◆ diMuonMass_tracker

MonitorElement* BPhysicsOniaDQM::diMuonMass_tracker

private

Definition at line 54 of file BPhysicsOniaDQM.h.

◆ glbBkgNoCut

MonitorElement* BPhysicsOniaDQM::glbBkgNoCut

private

Definition at line 62 of file BPhysicsOniaDQM.h.

◆ glbSigCut

MonitorElement* BPhysicsOniaDQM::glbSigCut

private

Definition at line 60 of file BPhysicsOniaDQM.h.

◆ glbSigNoCut

MonitorElement* BPhysicsOniaDQM::glbSigNoCut

private

Definition at line 61 of file BPhysicsOniaDQM.h.

◆ global_background

MonitorElement* BPhysicsOniaDQM::global_background

private

Definition at line 56 of file BPhysicsOniaDQM.h.

◆ lumiSummaryToken_

edm::EDGetTokenT<LumiSummary> BPhysicsOniaDQM::lumiSummaryToken_

private

Definition at line 47 of file BPhysicsOniaDQM.h.

◆ metname

std::string BPhysicsOniaDQM::metname

private

Definition at line 50 of file BPhysicsOniaDQM.h.

◆ RefVtx

math::XYZPoint BPhysicsOniaDQM::RefVtx

private

Definition at line 70 of file BPhysicsOniaDQM.h.

◆ staBkgNoCut

MonitorElement* BPhysicsOniaDQM::staBkgNoCut

private

Definition at line 65 of file BPhysicsOniaDQM.h.

◆ standalone_background

MonitorElement* BPhysicsOniaDQM::standalone_background

private

Definition at line 58 of file BPhysicsOniaDQM.h.

◆ staSigCut

MonitorElement* BPhysicsOniaDQM::staSigCut

private

Definition at line 63 of file BPhysicsOniaDQM.h.

◆ staSigNoCut

MonitorElement* BPhysicsOniaDQM::staSigNoCut

private

Definition at line 64 of file BPhysicsOniaDQM.h.

◆ theMuonCollectionLabel_

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

private

Definition at line 46 of file BPhysicsOniaDQM.h.

◆ tracker_background

MonitorElement* BPhysicsOniaDQM::tracker_background

private

Definition at line 57 of file BPhysicsOniaDQM.h.

◆ trkBkgNoCut

MonitorElement* BPhysicsOniaDQM::trkBkgNoCut

private

Definition at line 68 of file BPhysicsOniaDQM.h.

◆ trkSigCut

MonitorElement* BPhysicsOniaDQM::trkSigCut

private

Definition at line 66 of file BPhysicsOniaDQM.h.

◆ trkSigNoCut

MonitorElement* BPhysicsOniaDQM::trkSigNoCut

private

Definition at line 67 of file BPhysicsOniaDQM.h.

◆ vertex_

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

private

Definition at line 44 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

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ BPhysicsOniaDQM()

◆ ~BPhysicsOniaDQM()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ computeMass()

◆ isMuonInAccept()

◆ selGlobalMuon()

◆ selTrackerMuon()

Member Data Documentation

◆ diMuonMass_global

◆ diMuonMass_standalone

◆ diMuonMass_tracker

◆ glbBkgNoCut

◆ glbSigCut

◆ glbSigNoCut

◆ global_background

◆ lumiSummaryToken_

◆ metname

◆ RefVtx

◆ staBkgNoCut

◆ standalone_background

◆ staSigCut

◆ staSigNoCut

◆ theMuonCollectionLabel_

◆ tracker_background

◆ trkBkgNoCut

◆ trkSigCut

◆ trkSigNoCut

◆ vertex_