#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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 >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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 ( const edm::ParameterSet & parameters )

Constructor.

Definition at line 28 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 54 of file BPhysicsOniaDQM.cc.

54 {}

Member Function Documentation

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.

References edm::Event::getByToken(), muon::isGoodMuon(), edm::HandleBase::isValid(), LogTrace, metname, patZpeak::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 &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 56 of file BPhysicsOniaDQM.cc.

References dqm::implementation::IBooker::book1D(), and dqm::implementation::NavigatorBase::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 198 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 220 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 226 of file BPhysicsOniaDQM.cc.

References reco::Muon::globalTrack(), reco::Muon::innerTrack(), reco::HitPattern::numberOfValidMuonHits(), AlCaHLTBitMon_ParallelJobs::p, reco::HitPattern::pixelLayersWithMeasurement(), and submitPVResolutionJobs::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 240 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 53 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::diMuonMass_standalone

private

Definition at line 55 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::diMuonMass_tracker

private

Definition at line 54 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::glbBkgNoCut

private

Definition at line 62 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::glbSigCut

private

Definition at line 60 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::glbSigNoCut

private

Definition at line 61 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::global_background

private

Definition at line 56 of file BPhysicsOniaDQM.h.

edm::EDGetTokenT<LumiSummary> BPhysicsOniaDQM::lumiSummaryToken_

private

Definition at line 47 of file BPhysicsOniaDQM.h.

std::string BPhysicsOniaDQM::metname

private

Definition at line 50 of file BPhysicsOniaDQM.h.

math::XYZPoint BPhysicsOniaDQM::RefVtx

private

Definition at line 70 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::staBkgNoCut

private

Definition at line 65 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::standalone_background

private

Definition at line 58 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::staSigCut

private

Definition at line 63 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::staSigNoCut

private

Definition at line 64 of file BPhysicsOniaDQM.h.

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

private

Definition at line 46 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::tracker_background

private

Definition at line 57 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::trkBkgNoCut

private

Definition at line 68 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::trkSigCut

private

Definition at line 66 of file BPhysicsOniaDQM.h.

MonitorElement* BPhysicsOniaDQM::trkSigNoCut

private

Definition at line 67 of file BPhysicsOniaDQM.h.

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

Member Function Documentation

Member Data Documentation