#include <DTOccupancyEfficiency.h>

Inheritance diagram for DTOccupancyEfficiency:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

	DTOccupancyEfficiency (const edm::ParameterSet &pset)
	Constructor. More...

	~DTOccupancyEfficiency () 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

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

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
bool	debug

MonitorElement *	digisPerEvent

std::map< int, std::map< int, MonitorElement * > >	digisPerRing

edm::EDGetTokenT< DTDigiCollection >	dtDigiToken_

edm::ParameterSet	parameters

edm::EDGetTokenT< DTRecSegment4DCollection >	recHits4DToken_

MonitorElement *	recHitsPer4DSegment

MonitorElement *	recHitsPerEvent

edm::EDGetTokenT< DTRecHitCollection >	recHitToken_

MonitorElement *	segments4DPerEvent

MonitorElement *	t0From4DPhiSegment

MonitorElement *	t0From4DZSegment

MonitorElement *	timeBoxesPerEvent

std::map< int, std::map< int, MonitorElement * > >	timeBoxesPerRing

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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

DQM Analysis of 4D DT segments, DTdigis and DTRecHits :

Definition at line 25 of file DTOccupancyEfficiency.h.

Constructor & Destructor Documentation

◆ DTOccupancyEfficiency()

DTOccupancyEfficiency::DTOccupancyEfficiency ( const edm::ParameterSet & pset )

Constructor.

Definition at line 22 of file DTOccupancyEfficiency.cc.

                                                                      {
   debug = pset.getUntrackedParameter<bool>("debug", false);
   // label for dtDigis
   dtDigiToken_ = consumes<DTDigiCollection>(edm::InputTag(pset.getParameter<string>("digiLabel")));
   // the name of the 4D rec hits collection
   recHits4DToken_ = consumes<DTRecSegment4DCollection>(edm::InputTag(pset.getParameter<string>("recHits4DLabel")));
   // the name of the rechits collection
   recHitToken_ = consumes<DTRecHitCollection>(edm::InputTag(pset.getParameter<string>("recHitLabel")));
  
   parameters = pset;
 }

References debug, and muonDTDigis_cfi::pset.

◆ ~DTOccupancyEfficiency()

DTOccupancyEfficiency::~DTOccupancyEfficiency ( )

override

Destructor.

Definition at line 34 of file DTOccupancyEfficiency.cc.

34 {}

Member Function Documentation

◆ analyze()

void DTOccupancyEfficiency::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 88 of file DTOccupancyEfficiency.cc.

                                                                                      {
   if (debug)
     cout << "[DTOccupancyEfficiency] Analyze #Run: " << event.id().run() << " #Event: " << event.id().event() << endl;
  
   // Digi collection
   edm::Handle<DTDigiCollection> dtdigis;
   event.getByToken(dtDigiToken_, dtdigis);
  
   int numberOfDigis = 0;
   std::map<int, std::map<int, int>> numberOfDigisPerRing;
  
   DTDigiCollection::DigiRangeIterator dtLayerId_It;
   for (dtLayerId_It = dtdigis->begin(); dtLayerId_It != dtdigis->end(); ++dtLayerId_It) {  // Loop over layers
     for (DTDigiCollection::const_iterator digiIt = ((*dtLayerId_It).second).first;
          digiIt != ((*dtLayerId_It).second).second;
          ++digiIt) {
       numberOfDigis++;
       int wheelId = ((*dtLayerId_It).first).wheel();
       int stationId = ((*dtLayerId_It).first).station();
       (numberOfDigisPerRing[stationId])[wheelId] += 1;
  
       timeBoxesPerEvent->Fill((*digiIt).countsTDC());
       (timeBoxesPerRing[stationId])[wheelId]->Fill((*digiIt).countsTDC());
     }
   }
  
   // Total number of Digis per Event
   digisPerEvent->Fill(numberOfDigis);
  
   // Digis per Ring in every wheel
   for (int station = 1; station < 5; station++)
     for (int wheel = -2; wheel < 3; wheel++)
       (digisPerRing[station])[wheel]->Fill((numberOfDigisPerRing[station])[wheel]);
  
   // Get the 4D segment collection from the event
   edm::Handle<DTRecSegment4DCollection> all4DSegments;
   event.getByToken(recHits4DToken_, all4DSegments);
  
   segments4DPerEvent->Fill(all4DSegments->size());
  
   // Get the rechit collection from the event
   Handle<DTRecHitCollection> dtRecHits;
   event.getByToken(recHitToken_, dtRecHits);
  
   recHitsPerEvent->Fill(dtRecHits->size());
  
   // Number of RecHits per segment and t0 from Segment
   // Loop on all segments
   for (DTRecSegment4DCollection::const_iterator segment = all4DSegments->begin(); segment != all4DSegments->end();
        ++segment) {
     unsigned int nHits = (segment->hasPhi() ? (segment->phiSegment())->recHits().size() : 0);
     nHits += (segment->hasZed() ? (segment->zSegment())->recHits().size() : 0);
     recHitsPer4DSegment->Fill(nHits);
  
     if (segment->hasPhi()) {
       double segmentPhiT0 = segment->phiSegment()->t0();
       if (segment->phiSegment()->ist0Valid())
         t0From4DPhiSegment->Fill(segmentPhiT0);
     }
     if (segment->hasZed()) {
       double segmentZT0 = segment->zSegment()->t0();
       if (segment->zSegment()->ist0Valid())
         t0From4DZSegment->Fill(segmentZT0);
     }
   }
 }

References gather_cfg::cout, debug, HcalObjRepresent::Fill(), nHits, FastTrackerRecHitMaskProducer_cfi::recHits, edm::second(), relativeConstraints::station, and makeMuonMisalignmentScenario::wheel.

◆ bookHistograms()

void DTOccupancyEfficiency::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	context
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 36 of file DTOccupancyEfficiency.cc.

                                                                          {
   ibooker.setCurrentFolder("DT/04-OccupancyEfficiency/digisPerRing");
   // Digis per ring
   for (int station = 1; station < 5; station++) {
     string station_s = to_string(station);
     for (int wheel = -2; wheel < 3; wheel++) {
       string wheel_s = to_string(wheel);
       if (wheel > 0)
         wheel_s = "+" + wheel_s;
       string histoName = "digisPerMB" + station_s + "W" + wheel_s;
       string histoTitle = "Number of digis in MB" + station_s + "YB" + wheel_s;
       (digisPerRing[station])[wheel] = ibooker.book1D(histoName, histoTitle, 100, 0, 150);
     }
   }
  
   ibooker.setCurrentFolder("DT/04-OccupancyEfficiency/timeBoxesPerRing");
   // TimeBoxes per ring
   for (int station = 1; station < 5; station++) {
     string station_s = to_string(station);
     for (int wheel = -2; wheel < 3; wheel++) {
       string wheel_s = to_string(wheel);
       if (wheel > 0)
         wheel_s = "+" + wheel_s;
       string histoName = "timeBoxesPerMB" + station_s + "W" + wheel_s;
       string histoTitle = "Number of TDC counts in MB" + station_s + "YB" + wheel_s;
       (timeBoxesPerRing[station])[wheel] = ibooker.book1D(histoName, histoTitle, 400, 0, 1600);
     }
   }
  
   ibooker.setCurrentFolder("DT/04-OccupancyEfficiency");
  
   // TimeBoxes
   timeBoxesPerEvent = ibooker.book1D("timeBoxesPerEvent", "TDC counts per event", 400, 0, 1600);
  
   // Digis
   digisPerEvent = ibooker.book1D("digisPerEvent", "Number of digis per event", 100, 0, 900);
  
   // RecHits
   recHitsPerEvent = ibooker.book1D("recHitsPerEvent", "Number of RecHits per event", 100, 0, 250);
  
   // 4D segments
   segments4DPerEvent = ibooker.book1D("segments4DPerEvent", "Number of 4D Segments per event", 50, 0, 50);
  
   recHitsPer4DSegment = ibooker.book1D("recHitsPer4DSegment", "Number of RecHits per segment", 16, 0.5, 16.5);
  
   // T0 from segements
   t0From4DPhiSegment = ibooker.book1D("t0From4DPhiSegment", "T0 from 4D Phi segments", 100, -150, 150);
   t0From4DZSegment = ibooker.book1D("t0From4DZSegment", "T0 from 4D Z segments", 100, -150, 150);
 }

References dqm::implementation::IBooker::book1D(), HltBtagPostValidation_cff::histoName, dqm::implementation::NavigatorBase::setCurrentFolder(), relativeConstraints::station, and makeMuonMisalignmentScenario::wheel.

Member Data Documentation

◆ debug

bool DTOccupancyEfficiency::debug

private

Definition at line 42 of file DTOccupancyEfficiency.h.

Referenced by util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), util.rrapi.RRApi::get(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), runTauIdMVA.TauIDEmbedder::runTauID(), and pkg.AbstractPkg::write().

◆ digisPerEvent

MonitorElement* DTOccupancyEfficiency::digisPerEvent

private

Definition at line 56 of file DTOccupancyEfficiency.h.

◆ digisPerRing

std::map<int, std::map<int, MonitorElement*> > DTOccupancyEfficiency::digisPerRing

private

Definition at line 64 of file DTOccupancyEfficiency.h.

◆ dtDigiToken_

edm::EDGetTokenT<DTDigiCollection> DTOccupancyEfficiency::dtDigiToken_

private

Definition at line 45 of file DTOccupancyEfficiency.h.

◆ parameters

edm::ParameterSet DTOccupancyEfficiency::parameters

private

Definition at line 53 of file DTOccupancyEfficiency.h.

◆ recHits4DToken_

edm::EDGetTokenT<DTRecSegment4DCollection> DTOccupancyEfficiency::recHits4DToken_

private

Definition at line 48 of file DTOccupancyEfficiency.h.

◆ recHitsPer4DSegment

MonitorElement* DTOccupancyEfficiency::recHitsPer4DSegment

private

Definition at line 59 of file DTOccupancyEfficiency.h.

◆ recHitsPerEvent

MonitorElement* DTOccupancyEfficiency::recHitsPerEvent

private

Definition at line 58 of file DTOccupancyEfficiency.h.

◆ recHitToken_

edm::EDGetTokenT<DTRecHitCollection> DTOccupancyEfficiency::recHitToken_

private

Definition at line 51 of file DTOccupancyEfficiency.h.

◆ segments4DPerEvent

MonitorElement* DTOccupancyEfficiency::segments4DPerEvent

private

Definition at line 57 of file DTOccupancyEfficiency.h.

◆ t0From4DPhiSegment

MonitorElement* DTOccupancyEfficiency::t0From4DPhiSegment

private

Definition at line 60 of file DTOccupancyEfficiency.h.

◆ t0From4DZSegment

MonitorElement* DTOccupancyEfficiency::t0From4DZSegment

private

Definition at line 61 of file DTOccupancyEfficiency.h.

◆ timeBoxesPerEvent

MonitorElement* DTOccupancyEfficiency::timeBoxesPerEvent

private

Definition at line 55 of file DTOccupancyEfficiency.h.

◆ timeBoxesPerRing

std::map<int, std::map<int, MonitorElement*> > DTOccupancyEfficiency::timeBoxesPerRing

private

Definition at line 63 of file DTOccupancyEfficiency.h.

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTOccupancyEfficiency()

◆ ~DTOccupancyEfficiency()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

Member Data Documentation

◆ debug

◆ digisPerEvent

◆ digisPerRing

◆ dtDigiToken_

◆ parameters

◆ recHits4DToken_

◆ recHitsPer4DSegment

◆ recHitsPerEvent

◆ recHitToken_

◆ segments4DPerEvent

◆ t0From4DPhiSegment

◆ t0From4DZSegment

◆ timeBoxesPerEvent

◆ timeBoxesPerRing