#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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () 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 data-class-funcs::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_