RPCRecHitValid Class Reference

#include <RPCRecHitValid.h>

Inheritance diagram for RPCRecHitValid:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	RPCRecHitValid (const edm::ParameterSet &pset)
	~RPCRecHitValid () override
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 Types
typedef MonitorElement *	MEP
typedef RPCRecHitCollection	RecHits
typedef SimHitTPAssociationProducer::SimHitTPAssociationList	SimHitAssoc
typedef edm::PSimHitContainer	SimHits
typedef TrackingParticleCollection	SimParticles

Private Attributes
std::map< int, int >	detIdToIndexMapBarrel_
std::map< int, int >	detIdToIndexMapEndcap_
RPCValidHistograms	h_
MEP	h_eventCount
MEP	h_matchOccupancyBarrel_detId
MEP	h_matchOccupancyEndcap_detId
MEP	h_noiseOccupancyBarrel_detId
MEP	h_noiseOccupancyEndcap_detId
MEP	h_nRPCHitPerRecoMuon
MEP	h_nRPCHitPerRecoMuonBarrel
MEP	h_nRPCHitPerRecoMuonEndcap
MEP	h_nRPCHitPerRecoMuonOverlap
MEP	h_nRPCHitPerSimMuon
MEP	h_nRPCHitPerSimMuonBarrel
MEP	h_nRPCHitPerSimMuonEndcap
MEP	h_nRPCHitPerSimMuonOverlap
MEP	h_recoMuonBarrel_eta
MEP	h_recoMuonBarrel_phi
MEP	h_recoMuonBarrel_pt
MEP	h_recoMuonEndcap_eta
MEP	h_recoMuonEndcap_phi
MEP	h_recoMuonEndcap_pt
MEP	h_recoMuonNoRPC_eta
MEP	h_recoMuonNoRPC_phi
MEP	h_recoMuonNoRPC_pt
MEP	h_recoMuonOverlap_eta
MEP	h_recoMuonOverlap_phi
MEP	h_recoMuonOverlap_pt
MEP	h_recPunchOccupancyBarrel_station
MEP	h_recPunchOccupancyBarrel_wheel
MEP	h_recPunchOccupancyBarrel_wheel_station
MEP	h_recPunchOccupancyEndcap_disk
MEP	h_recPunchOccupancyEndcap_disk_ring
MEP	h_refOccupancyBarrel_detId
MEP	h_refOccupancyEndcap_detId
MEP	h_refPunchOccupancyBarrel_station
MEP	h_refPunchOccupancyBarrel_wheel
MEP	h_refPunchOccupancyBarrel_wheel_station
MEP	h_refPunchOccupancyEndcap_disk
MEP	h_refPunchOccupancyEndcap_disk_ring
MEP	h_rollAreaBarrel_detId
MEP	h_rollAreaEndcap_detId
MEP	h_simMuonBarrel_eta
MEP	h_simMuonBarrel_phi
MEP	h_simMuonBarrel_pt
MEP	h_simMuonEndcap_eta
MEP	h_simMuonEndcap_phi
MEP	h_simMuonEndcap_pt
MEP	h_simMuonNoRPC_eta
MEP	h_simMuonNoRPC_phi
MEP	h_simMuonNoRPC_pt
MEP	h_simMuonOverlap_eta
MEP	h_simMuonOverlap_phi
MEP	h_simMuonOverlap_pt
MEP	h_simParticleType
MEP	h_simParticleTypeBarrel
MEP	h_simParticleTypeEndcap
edm::EDGetTokenT< reco::MuonCollection >	muonToken_
edm::EDGetTokenT< RecHits >	recHitToken_
edm::EDGetTokenT< SimHitAssoc >	simHitAssocToken_
edm::EDGetTokenT< SimHits >	simHitToken_
edm::EDGetTokenT< SimParticles >	simParticleToken_
std::string	subDir_

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

Definition at line 24 of file RPCRecHitValid.h.

Member Typedef Documentation

MEP

typedef MonitorElement* RPCRecHitValid::MEP

private

Definition at line 45 of file RPCRecHitValid.h.

RecHits

typedef RPCRecHitCollection RPCRecHitValid::RecHits

private

Definition at line 34 of file RPCRecHitValid.h.

SimHitAssoc

typedef SimHitTPAssociationProducer::SimHitTPAssociationList RPCRecHitValid::SimHitAssoc

private

Definition at line 36 of file RPCRecHitValid.h.

SimHits

typedef edm::PSimHitContainer RPCRecHitValid::SimHits

private

Definition at line 33 of file RPCRecHitValid.h.

SimParticles

typedef TrackingParticleCollection RPCRecHitValid::SimParticles

private

Definition at line 35 of file RPCRecHitValid.h.

Constructor & Destructor Documentation

RPCRecHitValid()

RPCRecHitValid::RPCRecHitValid

(

const edm::ParameterSet &

pset

)

Definition at line 26 of file RPCRecHitValid.cc.

                                                          {
  simHitToken_ = consumes<SimHits>(pset.getParameter<edm::InputTag>("simHit"));
  recHitToken_ = consumes<RecHits>(pset.getParameter<edm::InputTag>("recHit"));
  simParticleToken_ = consumes<SimParticles>(pset.getParameter<edm::InputTag>("simTrack"));
  simHitAssocToken_ = consumes<SimHitAssoc>(pset.getParameter<edm::InputTag>("simHitAssoc"));
  muonToken_ = consumes<reco::MuonCollection>(pset.getParameter<edm::InputTag>("muon"));
 
  subDir_ = pset.getParameter<std::string>("subDir");
}

References muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

~RPCRecHitValid()

RPCRecHitValid::~RPCRecHitValid ( )

inlineoverride

Definition at line 27 of file RPCRecHitValid.h.

{};

Member Function Documentation

analyze()

void RPCRecHitValid::analyze ( const edm::Event &  event, const edm::EventSetup &  eventSetup  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 354 of file RPCRecHitValid.cc.

                                                                                   {
  h_eventCount->Fill(1);
 
  // Get the RPC Geometry
  edm::ESHandle<RPCGeometry> rpcGeom;
  eventSetup.get<MuonGeometryRecord>().get(rpcGeom);
 
  // Retrieve SimHits from the event
  edm::Handle<edm::PSimHitContainer> simHitHandle;
  if (!event.getByToken(simHitToken_, simHitHandle)) {
    edm::LogInfo("RPCRecHitValid") << "Cannot find simHit collection\n";
    return;
  }
 
  // Retrieve RecHits from the event
  edm::Handle<RPCRecHitCollection> recHitHandle;
  if (!event.getByToken(recHitToken_, recHitHandle)) {
    edm::LogInfo("RPCRecHitValid") << "Cannot find recHit collection\n";
    return;
  }
 
  // Get SimParticles
  edm::Handle<TrackingParticleCollection> simParticleHandle;
  if (!event.getByToken(simParticleToken_, simParticleHandle)) {
    edm::LogInfo("RPCRecHitValid") << "Cannot find TrackingParticle collection\n";
    return;
  }
 
  // Get SimParticle to SimHit association map
  edm::Handle<SimHitTPAssociationProducer::SimHitTPAssociationList> simHitsTPAssoc;
  if (!event.getByToken(simHitAssocToken_, simHitsTPAssoc)) {
    edm::LogInfo("RPCRecHitValid") << "Cannot find TrackingParticle to SimHit association map\n";
    return;
  }
 
  // Get RecoMuons
  edm::Handle<reco::MuonCollection> muonHandle;
  if (!event.getByToken(muonToken_, muonHandle)) {
    edm::LogInfo("RPCRecHitValid") << "Cannot find muon collection\n";
    return;
  }
 
  typedef edm::PSimHitContainer::const_iterator SimHitIter;
  typedef RPCRecHitCollection::const_iterator RecHitIter;
  typedef std::vector<TrackPSimHitRef> SimHitRefs;
 
  // TrackingParticles with (and without) RPC simHits
  SimHitRefs muonSimHits, pthrSimHits;
 
  for (int i = 0, n = simParticleHandle->size(); i < n; ++i) {
    TrackingParticleRef simParticle(simParticleHandle, i);
    if (simParticle->pt() < 1.0 or simParticle->p() < 2.5)
      continue;  // globalMuon acceptance
 
    // Collect SimHits from this Tracking Particle
    SimHitRefs simHitsFromParticle;
    auto range = std::equal_range(simHitsTPAssoc->begin(),
                                  simHitsTPAssoc->end(),
                                  std::make_pair(simParticle, TrackPSimHitRef()),
                                  SimHitTPAssociationProducer::simHitTPAssociationListGreater);
    for (auto simParticleToHit = range.first; simParticleToHit != range.second; ++simParticleToHit) {
      auto simHit = simParticleToHit->second;
      const DetId detId(simHit->detUnitId());
      if (detId.det() != DetId::Muon or detId.subdetId() != MuonSubdetId::RPC)
        continue;
 
      simHitsFromParticle.push_back(simParticleToHit->second);
    }
    const int nRPCHit = simHitsFromParticle.size();
    const bool hasRPCHit = nRPCHit > 0;
 
    if (abs(simParticle->pdgId()) == 13) {
      muonSimHits.insert(muonSimHits.end(), simHitsFromParticle.begin(), simHitsFromParticle.end());
 
      // Count number of Barrel hits and Endcap hits
      int nRPCHitBarrel = 0;
      int nRPCHitEndcap = 0;
      for (auto simHit : simHitsFromParticle) {
        const RPCDetId rpcDetId = static_cast<const RPCDetId>(simHit->detUnitId());
        const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(rpcDetId));
        if (!roll)
          continue;
 
        if (rpcDetId.region() == 0)
          ++nRPCHitBarrel;
        else
          ++nRPCHitEndcap;
      }
 
      // Fill TrackingParticle related histograms
      h_nRPCHitPerSimMuon->Fill(nRPCHit);
      if (nRPCHitBarrel and nRPCHitEndcap) {
        h_nRPCHitPerSimMuonOverlap->Fill(nRPCHit);
        h_simMuonOverlap_pt->Fill(simParticle->pt());
        h_simMuonOverlap_eta->Fill(simParticle->eta());
        h_simMuonOverlap_phi->Fill(simParticle->phi());
      } else if (nRPCHitBarrel) {
        h_nRPCHitPerSimMuonBarrel->Fill(nRPCHit);
        h_simMuonBarrel_pt->Fill(simParticle->pt());
        h_simMuonBarrel_eta->Fill(simParticle->eta());
        h_simMuonBarrel_phi->Fill(simParticle->phi());
      } else if (nRPCHitEndcap) {
        h_nRPCHitPerSimMuonEndcap->Fill(nRPCHit);
        h_simMuonEndcap_pt->Fill(simParticle->pt());
        h_simMuonEndcap_eta->Fill(simParticle->eta());
        h_simMuonEndcap_phi->Fill(simParticle->phi());
      } else {
        h_simMuonNoRPC_pt->Fill(simParticle->pt());
        h_simMuonNoRPC_eta->Fill(simParticle->eta());
        h_simMuonNoRPC_phi->Fill(simParticle->phi());
      }
    } else {
      pthrSimHits.insert(pthrSimHits.end(), simHitsFromParticle.begin(), simHitsFromParticle.end());
    }
 
    if (hasRPCHit) {
      switch (simParticle->pdgId()) {
        case 13:
          h_simParticleType->Fill(0);
          break;
        case -13:
          h_simParticleType->Fill(1);
          break;
        case 11:
          h_simParticleType->Fill(2);
          break;
        case -11:
          h_simParticleType->Fill(3);
          break;
        case 211:
          h_simParticleType->Fill(4);
          break;
        case -211:
          h_simParticleType->Fill(5);
          break;
        case 321:
          h_simParticleType->Fill(6);
          break;
        case -321:
          h_simParticleType->Fill(7);
          break;
        case 2212:
          h_simParticleType->Fill(8);
          break;
        case -2212:
          h_simParticleType->Fill(9);
          break;
        default:
          h_simParticleType->Fill(10);
          break;
      }
    }
  }
 
  // Loop over muon simHits, fill histograms which does not need associations
  int nRefHitBarrel = 0, nRefHitEndcap = 0;
  for (auto simHit : muonSimHits) {
    const RPCDetId detId = static_cast<const RPCDetId>(simHit->detUnitId());
    const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(detId));
 
    const int region = roll->id().region();
    const int ring = roll->id().ring();
    // const int sector = roll->id().sector();
    const int station = roll->id().station();
    // const int layer = roll->id().layer();
    // const int subSector = roll->id().subsector();
 
    if (region == 0) {
      ++nRefHitBarrel;
      h_.refHitOccupancyBarrel_wheel->Fill(ring);
      h_.refHitOccupancyBarrel_station->Fill(station);
      h_.refHitOccupancyBarrel_wheel_station->Fill(ring, station);
 
      h_refOccupancyBarrel_detId->Fill(detIdToIndexMapBarrel_[simHit->detUnitId()]);
    } else {
      ++nRefHitEndcap;
      h_.refHitOccupancyEndcap_disk->Fill(region * station);
      h_.refHitOccupancyEndcap_disk_ring->Fill(region * station, ring);
 
      h_refOccupancyEndcap_detId->Fill(detIdToIndexMapEndcap_[simHit->detUnitId()]);
    }
  }
 
  // Loop over punch-through simHits, fill histograms which does not need
  // associations
  for (auto simHit : pthrSimHits) {
    const RPCDetId detId = static_cast<const RPCDetId>(simHit->detUnitId());
    const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(detId()));
 
    const int region = roll->id().region();
    const int ring = roll->id().ring();
    // const int sector = roll->id().sector();
    const int station = roll->id().station();
    // const int layer = roll->id().layer();
    // const int subSector = roll->id().subsector();
 
    if (region == 0) {
      ++nRefHitBarrel;
      h_refPunchOccupancyBarrel_wheel->Fill(ring);
      h_refPunchOccupancyBarrel_station->Fill(station);
      h_refPunchOccupancyBarrel_wheel_station->Fill(ring, station);
 
      h_refOccupancyBarrel_detId->Fill(detIdToIndexMapBarrel_[simHit->detUnitId()]);
    } else {
      ++nRefHitEndcap;
      h_refPunchOccupancyEndcap_disk->Fill(region * station);
      h_refPunchOccupancyEndcap_disk_ring->Fill(region * station, ring);
 
      h_refOccupancyEndcap_detId->Fill(detIdToIndexMapEndcap_[simHit->detUnitId()]);
    }
  }
  h_.nRefHitBarrel->Fill(nRefHitBarrel);
  h_.nRefHitEndcap->Fill(nRefHitEndcap);
 
  // Loop over recHits, fill histograms which does not need associations
  int sumClusterSizeBarrel = 0, sumClusterSizeEndcap = 0;
  int nRecHitBarrel = 0, nRecHitEndcap = 0;
  for (RecHitIter recHitIter = recHitHandle->begin(); recHitIter != recHitHandle->end(); ++recHitIter) {
    const RPCDetId detId = static_cast<const RPCDetId>(recHitIter->rpcId());
    const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(detId()));
    if (!roll)
      continue;
 
    const int region = roll->id().region();
    const int ring = roll->id().ring();
    // const int sector = roll->id().sector();
    const int station = roll->id().station();
    // const int layer = roll->id().layer();
    // const int subSector = roll->id().subsector();
 
    const double time = recHitIter->timeError() >= 0 ? recHitIter->time() : recHitIter->BunchX() * 25;
 
    h_.clusterSize->Fill(recHitIter->clusterSize());
 
    if (region == 0) {
      ++nRecHitBarrel;
      sumClusterSizeBarrel += recHitIter->clusterSize();
      h_.clusterSizeBarrel->Fill(recHitIter->clusterSize());
      h_.recHitOccupancyBarrel_wheel->Fill(ring);
      h_.recHitOccupancyBarrel_station->Fill(station);
      h_.recHitOccupancyBarrel_wheel_station->Fill(ring, station);
 
      h_.timeBarrel->Fill(time);
    } else {
      ++nRecHitEndcap;
      sumClusterSizeEndcap += recHitIter->clusterSize();
      h_.clusterSizeEndcap->Fill(recHitIter->clusterSize());
      h_.recHitOccupancyEndcap_disk->Fill(region * station);
      h_.recHitOccupancyEndcap_disk_ring->Fill(region * station, ring);
 
      h_.timeEndcap->Fill(time);
    }
 
    if (roll->isIRPC()) {
      h_.timeIRPC->Fill(time);
    } else {
      h_.timeCRPC->Fill(time);
    }
  }
  const double nRecHit = nRecHitBarrel + nRecHitEndcap;
  h_.nRecHitBarrel->Fill(nRecHitBarrel);
  h_.nRecHitEndcap->Fill(nRecHitEndcap);
  if (nRecHit > 0) {
    const int sumClusterSize = sumClusterSizeBarrel + sumClusterSizeEndcap;
    h_.avgClusterSize->Fill(double(sumClusterSize) / nRecHit);
 
    if (nRecHitBarrel > 0) {
      h_.avgClusterSizeBarrel->Fill(double(sumClusterSizeBarrel) / nRecHitBarrel);
    }
    if (nRecHitEndcap > 0) {
      h_.avgClusterSizeEndcap->Fill(double(sumClusterSizeEndcap) / nRecHitEndcap);
    }
  }
 
  // Start matching SimHits to RecHits
  typedef std::map<TrackPSimHitRef, RecHitIter> SimToRecHitMap;
  SimToRecHitMap simToRecHitMap;
 
  for (auto simHit : muonSimHits) {
    const RPCDetId simDetId = static_cast<const RPCDetId>(simHit->detUnitId());
    // const RPCRoll* simRoll = dynamic_cast<const
    // RPCRoll*>(rpcGeom->roll(simDetId));
 
    const double simX = simHit->localPosition().x();
 
    for (RecHitIter recHitIter = recHitHandle->begin(); recHitIter != recHitHandle->end(); ++recHitIter) {
      const RPCDetId recDetId = static_cast<const RPCDetId>(recHitIter->rpcId());
      const RPCRoll *recRoll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(recDetId));
      if (!recRoll)
        continue;
 
      if (simDetId != recDetId)
        continue;
 
      const double recX = recHitIter->localPosition().x();
      const double newDx = fabs(recX - simX);
 
      // Associate SimHit to RecHit
      SimToRecHitMap::const_iterator prevSimToReco = simToRecHitMap.find(simHit);
      if (prevSimToReco == simToRecHitMap.end()) {
        simToRecHitMap.insert(std::make_pair(simHit, recHitIter));
      } else {
        const double oldDx = fabs(prevSimToReco->second->localPosition().x() - simX);
 
        if (newDx < oldDx) {
          simToRecHitMap[simHit] = recHitIter;
        }
      }
    }
  }
 
  // Now we have simHit-recHit mapping
  // So we can fill up relavant histograms
  int nMatchHitBarrel = 0, nMatchHitEndcap = 0;
  for (auto match : simToRecHitMap) {
    TrackPSimHitRef simHit = match.first;
    RecHitIter recHitIter = match.second;
 
    const RPCDetId detId = static_cast<const RPCDetId>(simHit->detUnitId());
    const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(detId));
 
    const int region = roll->id().region();
    const int ring = roll->id().ring();
    // const int sector = roll->id().sector();
    const int station = roll->id().station();
    // const int layer = roll->id().layer();
    // const int subsector = roll->id().subsector();
 
    const double simX = simHit->localPosition().x();
    const double recX = recHitIter->localPosition().x();
    const double errX = sqrt(recHitIter->localPositionError().xx());
    const double dX = recX - simX;
    const double pull = errX == 0 ? -999 : dX / errX;
 
    // const GlobalPoint simPos = roll->toGlobal(simHitIter->localPosition());
    // const GlobalPoint recPos = roll->toGlobal(recHitIter->localPosition());
 
    if (region == 0) {
      ++nMatchHitBarrel;
      h_.resBarrel->Fill(dX);
      h_.pullBarrel->Fill(pull);
      h_.matchOccupancyBarrel_wheel->Fill(ring);
      h_.matchOccupancyBarrel_station->Fill(station);
      h_.matchOccupancyBarrel_wheel_station->Fill(ring, station);
 
      h_.res_wheel_res->Fill(ring, dX);
      h_.res_station_res->Fill(station, dX);
      h_.pull_wheel_pull->Fill(ring, pull);
      h_.pull_station_pull->Fill(station, pull);
 
      h_matchOccupancyBarrel_detId->Fill(detIdToIndexMapBarrel_[detId.rawId()]);
    } else {
      ++nMatchHitEndcap;
      h_.resEndcap->Fill(dX);
      h_.pullEndcap->Fill(pull);
      h_.matchOccupancyEndcap_disk->Fill(region * station);
      h_.matchOccupancyEndcap_disk_ring->Fill(region * station, ring);
 
      h_.res_disk_res->Fill(region * station, dX);
      h_.res_ring_res->Fill(ring, dX);
      h_.pull_disk_pull->Fill(region * station, pull);
      h_.pull_ring_pull->Fill(ring, pull);
 
      h_matchOccupancyEndcap_detId->Fill(detIdToIndexMapEndcap_[detId.rawId()]);
    }
  }
  h_.nMatchHitBarrel->Fill(nMatchHitBarrel);
  h_.nMatchHitEndcap->Fill(nMatchHitEndcap);
 
  // Reco Muon hits
  for (reco::MuonCollection::const_iterator muon = muonHandle->begin(); muon != muonHandle->end(); ++muon) {
    if (!muon->isGlobalMuon())
      continue;
 
    int nRPCHitBarrel = 0;
    int nRPCHitEndcap = 0;
 
    const reco::TrackRef glbTrack = muon->globalTrack();
    for (trackingRecHit_iterator recHit = glbTrack->recHitsBegin(); recHit != glbTrack->recHitsEnd(); ++recHit) {
      if (!(*recHit)->isValid())
        continue;
      const DetId detId = (*recHit)->geographicalId();
      if (detId.det() != DetId::Muon or detId.subdetId() != MuonSubdetId::RPC)
        continue;
      const RPCDetId rpcDetId = static_cast<const RPCDetId>(detId);
 
      if (rpcDetId.region() == 0)
        ++nRPCHitBarrel;
      else
        ++nRPCHitEndcap;
    }
 
    const int nRPCHit = nRPCHitBarrel + nRPCHitEndcap;
    h_nRPCHitPerRecoMuon->Fill(nRPCHit);
    if (nRPCHitBarrel and nRPCHitEndcap) {
      h_nRPCHitPerRecoMuonOverlap->Fill(nRPCHit);
      h_recoMuonOverlap_pt->Fill(muon->pt());
      h_recoMuonOverlap_eta->Fill(muon->eta());
      h_recoMuonOverlap_phi->Fill(muon->phi());
    } else if (nRPCHitBarrel) {
      h_nRPCHitPerRecoMuonBarrel->Fill(nRPCHit);
      h_recoMuonBarrel_pt->Fill(muon->pt());
      h_recoMuonBarrel_eta->Fill(muon->eta());
      h_recoMuonBarrel_phi->Fill(muon->phi());
    } else if (nRPCHitEndcap) {
      h_nRPCHitPerRecoMuonEndcap->Fill(nRPCHit);
      h_recoMuonEndcap_pt->Fill(muon->pt());
      h_recoMuonEndcap_eta->Fill(muon->eta());
      h_recoMuonEndcap_phi->Fill(muon->phi());
    } else {
      h_recoMuonNoRPC_pt->Fill(muon->pt());
      h_recoMuonNoRPC_eta->Fill(muon->eta());
      h_recoMuonNoRPC_phi->Fill(muon->phi());
    }
  }
 
  // Find Non-muon hits
  for (RecHitIter recHitIter = recHitHandle->begin(); recHitIter != recHitHandle->end(); ++recHitIter) {
    const RPCDetId detId = static_cast<const RPCDetId>(recHitIter->rpcId());
    const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(detId));
 
    const int region = roll->id().region();
    const int ring = roll->id().ring();
    // const int sector = roll->id().sector();
    const int station = roll->id().station();
    // const int layer = roll->id().layer();
    // const int subsector = roll->id().subsector();
 
    bool matched = false;
    for (auto match : simToRecHitMap) {
      if (recHitIter == match.second) {
        matched = true;
        break;
      }
    }
 
    if (!matched) {
      int nPunchMatched = 0;
      // Check if this recHit came from non-muon simHit
      for (auto simHit : pthrSimHits) {
        const int absSimHitPType = abs(simHit->particleType());
        if (absSimHitPType == 13)
          continue;
 
        const RPCDetId simDetId = static_cast<const RPCDetId>(simHit->detUnitId());
        if (simDetId == detId)
          ++nPunchMatched;
      }
 
      if (nPunchMatched > 0) {
        if (region == 0) {
          h_recPunchOccupancyBarrel_wheel->Fill(ring);
          h_recPunchOccupancyBarrel_station->Fill(station);
          h_recPunchOccupancyBarrel_wheel_station->Fill(ring, station);
        } else {
          h_recPunchOccupancyEndcap_disk->Fill(region * station);
          h_recPunchOccupancyEndcap_disk_ring->Fill(region * station, ring);
        }
      }
    }
  }
 
  // Find noise recHits : RecHits without SimHit match
  for (RecHitIter recHitIter = recHitHandle->begin(); recHitIter != recHitHandle->end(); ++recHitIter) {
    const RPCDetId recDetId = static_cast<const RPCDetId>(recHitIter->rpcId());
    const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(recDetId));
 
    const int region = roll->id().region();
    //    const int ring = roll->id().ring(); // UNUSED VARIABLE
    // const int sector = roll->id().sector();
    //    const int station = roll->id().station(); // UNUSED VARIABLE
    // const int layer = roll->id().layer();
    // const int subsector = roll->id().subsector();
 
    const double recX = recHitIter->localPosition().x();
    const double recErrX = sqrt(recHitIter->localPositionError().xx());
 
    bool matched = false;
    for (SimHitIter simHitIter = simHitHandle->begin(); simHitIter != simHitHandle->end(); ++simHitIter) {
      const RPCDetId simDetId = static_cast<const RPCDetId>(simHitIter->detUnitId());
      const RPCRoll *simRoll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(simDetId));
      if (!simRoll)
        continue;
 
      if (simDetId != recDetId)
        continue;
 
      const double simX = simHitIter->localPosition().x();
      const double dX = fabs(recX - simX);
 
      if (dX / recErrX < 5) {
        matched = true;
        break;
      }
    }
 
    if (!matched) {
      if (region == 0) {
        h_noiseOccupancyBarrel_detId->Fill(detIdToIndexMapBarrel_[recDetId.rawId()]);
      } else {
        h_noiseOccupancyEndcap_detId->Fill(detIdToIndexMapEndcap_[recDetId.rawId()]);
      }
    }
  }
 
  h_eventCount->Fill(2);
}

References funct::abs(), DetId::det(), fftjetvertexadder_cfi::errX, options_cfi::eventSetup, get, mps_fire::i, match(), muonTagProbeFilters_cff::matched, DetId::Muon, HLT_2018_cff::muon, dqmiodumpmetadata::n, or, FastTimerService_cff::range, DetId::rawId(), rpcPointValidation_cfi::recHit, RPCDetId::region(), HLT_2018_cff::region, relativeConstraints::ring, RPCGeometry::roll(), MuonSubdetId::RPC, rpcPointValidation_cfi::simHit, SimHitTPAssociationProducer::simHitTPAssociationListGreater(), mathSSE::sqrt(), relativeConstraints::station, DetId::subdetId(), and ntuplemaker::time.

bookHistograms()

void RPCRecHitValid::bookHistograms ( DQMStore::IBooker &  booker, edm::Run const &  run, edm::EventSetup const &  eventSetup  )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 36 of file RPCRecHitValid.cc.

                                                                                                               {
  // Book MonitorElements
  h_.bookHistograms(booker, subDir_);
 
  // SimHit plots, not compatible to RPCPoint-RPCRecHit comparison
  booker.setCurrentFolder(subDir_ + "/HitProperty");
  h_simParticleType = booker.book1D("SimHitPType", "SimHit particle type", 11, 0, 11);
  h_simParticleType->getTH1()->SetMinimum(0);
  if (TH1 *h = h_simParticleType->getTH1()) {
    h->GetXaxis()->SetBinLabel(1, "#mu^{-}");
    h->GetXaxis()->SetBinLabel(2, "#mu^{+}");
    h->GetXaxis()->SetBinLabel(3, "e^{-}");
    h->GetXaxis()->SetBinLabel(4, "e^{+}");
    h->GetXaxis()->SetBinLabel(5, "#pi^{+}");
    h->GetXaxis()->SetBinLabel(6, "#pi^{-}");
    h->GetXaxis()->SetBinLabel(7, "K^{+}");
    h->GetXaxis()->SetBinLabel(8, "K^{-}");
    h->GetXaxis()->SetBinLabel(9, "p^{+}");
    h->GetXaxis()->SetBinLabel(10, "p^{-}");
    h->GetXaxis()->SetBinLabel(11, "Other");
  }
 
  booker.setCurrentFolder(subDir_ + "/Track");
 
  h_nRPCHitPerSimMuon = booker.book1D("NRPCHitPerSimMuon", "Number of RPC SimHit per SimMuon", 11, -0.5, 10.5);
  h_nRPCHitPerSimMuonBarrel =
      booker.book1D("NRPCHitPerSimMuonBarrel", "Number of RPC SimHit per SimMuon", 11, -0.5, 10.5);
  h_nRPCHitPerSimMuonOverlap =
      booker.book1D("NRPCHitPerSimMuonOverlap", "Number of RPC SimHit per SimMuon", 11, -0.5, 10.5);
  h_nRPCHitPerSimMuonEndcap =
      booker.book1D("NRPCHitPerSimMuonEndcap", "Number of RPC SimHit per SimMuon", 11, -0.5, 10.5);
 
  h_nRPCHitPerRecoMuon = booker.book1D("NRPCHitPerRecoMuon", "Number of RPC RecHit per RecoMuon", 11, -0.5, 10.5);
  h_nRPCHitPerRecoMuonBarrel =
      booker.book1D("NRPCHitPerRecoMuonBarrel", "Number of RPC RecHit per RecoMuon", 11, -0.5, 10.5);
  h_nRPCHitPerRecoMuonOverlap =
      booker.book1D("NRPCHitPerRecoMuonOverlap", "Number of RPC RecHit per RecoMuon", 11, -0.5, 10.5);
  h_nRPCHitPerRecoMuonEndcap =
      booker.book1D("NRPCHitPerRecoMuonEndcap", "Number of RPC RecHit per RecoMuon", 11, -0.5, 10.5);
 
  h_nRPCHitPerSimMuon->getTH1()->SetMinimum(0);
  h_nRPCHitPerSimMuonBarrel->getTH1()->SetMinimum(0);
  h_nRPCHitPerSimMuonOverlap->getTH1()->SetMinimum(0);
  h_nRPCHitPerSimMuonEndcap->getTH1()->SetMinimum(0);
 
  h_nRPCHitPerRecoMuon->getTH1()->SetMinimum(0);
  h_nRPCHitPerRecoMuonBarrel->getTH1()->SetMinimum(0);
  h_nRPCHitPerRecoMuonOverlap->getTH1()->SetMinimum(0);
  h_nRPCHitPerRecoMuonEndcap->getTH1()->SetMinimum(0);
 
  float ptBins[] = {0, 1, 2, 5, 10, 20, 30, 50, 100, 200, 300, 500};
  const int nPtBins = sizeof(ptBins) / sizeof(float) - 1;
  h_simMuonBarrel_pt =
      booker.book1D("SimMuonBarrel_pt", "SimMuon RPCHit in Barrel  p_{T};p_{T} [GeV/c^{2}]", nPtBins, ptBins);
  h_simMuonOverlap_pt =
      booker.book1D("SimMuonOverlap_pt", "SimMuon RPCHit in Overlap p_{T};p_{T} [GeV/c^{2}]", nPtBins, ptBins);
  h_simMuonEndcap_pt =
      booker.book1D("SimMuonEndcap_pt", "SimMuon RPCHit in Endcap  p_{T};p_{T} [GeV/c^{2}]", nPtBins, ptBins);
  h_simMuonNoRPC_pt =
      booker.book1D("SimMuonNoRPC_pt", "SimMuon without RPCHit p_{T};p_{T} [GeV/c^{2}]", nPtBins, ptBins);
  h_simMuonBarrel_eta = booker.book1D("SimMuonBarrel_eta", "SimMuon RPCHit in Barrel  #eta;#eta", 50, -2.5, 2.5);
  h_simMuonOverlap_eta = booker.book1D("SimMuonOverlap_eta", "SimMuon RPCHit in Overlap #eta;#eta", 50, -2.5, 2.5);
  h_simMuonEndcap_eta = booker.book1D("SimMuonEndcap_eta", "SimMuon RPCHit in Endcap  #eta;#eta", 50, -2.5, 2.5);
  h_simMuonNoRPC_eta = booker.book1D("SimMuonNoRPC_eta", "SimMuon without RPCHit #eta;#eta", 50, -2.5, 2.5);
  h_simMuonBarrel_phi =
      booker.book1D("SimMuonBarrel_phi", "SimMuon RPCHit in Barrel  #phi;#phi", 36, -TMath::Pi(), TMath::Pi());
  h_simMuonOverlap_phi =
      booker.book1D("SimMuonOverlap_phi", "SimMuon RPCHit in Overlap #phi;#phi", 36, -TMath::Pi(), TMath::Pi());
  h_simMuonEndcap_phi =
      booker.book1D("SimMuonEndcap_phi", "SimMuon RPCHit in Endcap  #phi;#phi", 36, -TMath::Pi(), TMath::Pi());
  h_simMuonNoRPC_phi =
      booker.book1D("SimMuonNoRPC_phi", "SimMuon without RPCHit #phi;#phi", 36, -TMath::Pi(), TMath::Pi());
 
  h_recoMuonBarrel_pt =
      booker.book1D("RecoMuonBarrel_pt", "RecoMuon RPCHit in Barrel  p_{T};p_{T} [GeV/c^{2}]", nPtBins, ptBins);
  h_recoMuonOverlap_pt =
      booker.book1D("RecoMuonOverlap_pt", "RecoMuon RPCHit in Overlap p_{T};p_{T} [GeV/c^{2}]", nPtBins, ptBins);
  h_recoMuonEndcap_pt =
      booker.book1D("RecoMuonEndcap_pt", "RecoMuon RPCHit in Endcap  p_{T};p_{T} [GeV/c^{2}]", nPtBins, ptBins);
  h_recoMuonNoRPC_pt =
      booker.book1D("RecoMuonNoRPC_pt", "RecoMuon without RPCHit p_{T};p_{T} [GeV/c^{2}]", nPtBins, ptBins);
  h_recoMuonBarrel_eta = booker.book1D("RecoMuonBarrel_eta", "RecoMuon RPCHit in Barrel  #eta;#eta", 50, -2.5, 2.5);
  h_recoMuonOverlap_eta = booker.book1D("RecoMuonOverlap_eta", "RecoMuon RPCHit in Overlap #eta;#eta", 50, -2.5, 2.5);
  h_recoMuonEndcap_eta = booker.book1D("RecoMuonEndcap_eta", "RecoMuon RPCHit in Endcap  #eta;#eta", 50, -2.5, 2.5);
  h_recoMuonNoRPC_eta = booker.book1D("RecoMuonNoRPC_eta", "RecoMuon without RPCHit #eta;#eta", 50, -2.5, 2.5);
  h_recoMuonBarrel_phi =
      booker.book1D("RecoMuonBarrel_phi", "RecoMuon RPCHit in Barrel  #phi;#phi", 36, -TMath::Pi(), TMath::Pi());
  h_recoMuonOverlap_phi =
      booker.book1D("RecoMuonOverlap_phi", "RecoMuon RPCHit in Overlap #phi;#phi", 36, -TMath::Pi(), TMath::Pi());
  h_recoMuonEndcap_phi =
      booker.book1D("RecoMuonEndcap_phi", "RecoMuon RPCHit in Endcap  #phi;#phi", 36, -TMath::Pi(), TMath::Pi());
  h_recoMuonNoRPC_phi =
      booker.book1D("RecoMuonNoRPC_phi", "RecoMuon without RPCHit #phi;#phi", 36, -TMath::Pi(), TMath::Pi());
 
  h_simMuonBarrel_pt->getTH1()->SetMinimum(0);
  h_simMuonOverlap_pt->getTH1()->SetMinimum(0);
  h_simMuonEndcap_pt->getTH1()->SetMinimum(0);
  h_simMuonNoRPC_pt->getTH1()->SetMinimum(0);
  h_simMuonBarrel_eta->getTH1()->SetMinimum(0);
  h_simMuonOverlap_eta->getTH1()->SetMinimum(0);
  h_simMuonEndcap_eta->getTH1()->SetMinimum(0);
  h_simMuonNoRPC_eta->getTH1()->SetMinimum(0);
  h_simMuonBarrel_phi->getTH1()->SetMinimum(0);
  h_simMuonOverlap_phi->getTH1()->SetMinimum(0);
  h_simMuonEndcap_phi->getTH1()->SetMinimum(0);
  h_simMuonNoRPC_phi->getTH1()->SetMinimum(0);
 
  h_recoMuonBarrel_pt->getTH1()->SetMinimum(0);
  h_recoMuonOverlap_pt->getTH1()->SetMinimum(0);
  h_recoMuonEndcap_pt->getTH1()->SetMinimum(0);
  h_recoMuonNoRPC_pt->getTH1()->SetMinimum(0);
  h_recoMuonBarrel_eta->getTH1()->SetMinimum(0);
  h_recoMuonOverlap_eta->getTH1()->SetMinimum(0);
  h_recoMuonEndcap_eta->getTH1()->SetMinimum(0);
  h_recoMuonNoRPC_eta->getTH1()->SetMinimum(0);
  h_recoMuonBarrel_phi->getTH1()->SetMinimum(0);
  h_recoMuonOverlap_phi->getTH1()->SetMinimum(0);
  h_recoMuonEndcap_phi->getTH1()->SetMinimum(0);
  h_recoMuonNoRPC_phi->getTH1()->SetMinimum(0);
 
  booker.setCurrentFolder(subDir_ + "/Occupancy");
 
  h_eventCount = booker.book1D("EventCount", "Event count", 3, 1, 4);
  h_eventCount->getTH1()->SetMinimum(0);
  if (h_eventCount) {
    TH1 *h = h_eventCount->getTH1();
    h->GetXaxis()->SetBinLabel(1, "eventBegin");
    h->GetXaxis()->SetBinLabel(2, "eventEnd");
    h->GetXaxis()->SetBinLabel(3, "run");
  }
  h_eventCount->Fill(3);
 
  h_refPunchOccupancyBarrel_wheel =
      booker.book1D("RefPunchOccupancyBarrel_wheel", "RefPunchthrough occupancy", 5, -2.5, 2.5);
  h_refPunchOccupancyEndcap_disk =
      booker.book1D("RefPunchOccupancyEndcap_disk", "RefPunchthrough occupancy", 9, -4.5, 4.5);
  h_refPunchOccupancyBarrel_station =
      booker.book1D("RefPunchOccupancyBarrel_station", "RefPunchthrough occupancy", 4, 0.5, 4.5);
  h_recPunchOccupancyBarrel_wheel =
      booker.book1D("RecPunchOccupancyBarrel_wheel", "Punchthrough recHit occupancy", 5, -2.5, 2.5);
  h_recPunchOccupancyEndcap_disk =
      booker.book1D("RecPunchOccupancyEndcap_disk", "Punchthrough recHit occupancy", 9, -4.5, 4.5);
  h_recPunchOccupancyBarrel_station =
      booker.book1D("RecPunchOccupancyBarrel_station", "Punchthrough recHit occupancy", 4, 0.5, 4.5);
 
  h_refPunchOccupancyBarrel_wheel->getTH1()->SetMinimum(0);
  h_refPunchOccupancyEndcap_disk->getTH1()->SetMinimum(0);
  h_refPunchOccupancyBarrel_station->getTH1()->SetMinimum(0);
  h_recPunchOccupancyBarrel_wheel->getTH1()->SetMinimum(0);
  h_recPunchOccupancyEndcap_disk->getTH1()->SetMinimum(0);
  h_recPunchOccupancyBarrel_station->getTH1()->SetMinimum(0);
 
  h_refPunchOccupancyBarrel_wheel_station =
      booker.book2D("RefPunchOccupancyBarrel_wheel_station", "RefPunchthrough occupancy", 5, -2.5, 2.5, 4, 0.5, 4.5);
  h_refPunchOccupancyEndcap_disk_ring =
      booker.book2D("RefPunchOccupancyEndcap_disk_ring", "RefPunchthrough occupancy", 9, -4.5, 4.5, 4, 0.5, 4.5);
  h_recPunchOccupancyBarrel_wheel_station = booker.book2D(
      "RecPunchOccupancyBarrel_wheel_station", "Punchthrough recHit occupancy", 5, -2.5, 2.5, 4, 0.5, 4.5);
  h_recPunchOccupancyEndcap_disk_ring =
      booker.book2D("RecPunchOccupancyEndcap_disk_ring", "Punchthrough recHit occupancy", 9, -4.5, 4.5, 4, 0.5, 4.5);
 
  h_refPunchOccupancyBarrel_wheel_station->getTH2F()->SetOption("COLZ");
  h_refPunchOccupancyEndcap_disk_ring->getTH2F()->SetOption("COLZ");
  h_recPunchOccupancyBarrel_wheel_station->getTH2F()->SetOption("COLZ");
  h_recPunchOccupancyEndcap_disk_ring->getTH2F()->SetOption("COLZ");
 
  h_refPunchOccupancyBarrel_wheel_station->getTH2F()->SetContour(10);
  h_refPunchOccupancyEndcap_disk_ring->getTH2F()->SetContour(10);
  h_recPunchOccupancyBarrel_wheel_station->getTH2F()->SetContour(10);
  h_recPunchOccupancyEndcap_disk_ring->getTH2F()->SetContour(10);
 
  h_refPunchOccupancyBarrel_wheel_station->getTH2F()->SetStats(false);
  h_refPunchOccupancyEndcap_disk_ring->getTH2F()->SetStats(false);
  h_recPunchOccupancyBarrel_wheel_station->getTH2F()->SetStats(false);
  h_recPunchOccupancyEndcap_disk_ring->getTH2F()->SetStats(false);
 
  h_refPunchOccupancyBarrel_wheel_station->getTH2F()->SetMinimum(0);
  h_refPunchOccupancyEndcap_disk_ring->getTH2F()->SetMinimum(0);
  h_recPunchOccupancyBarrel_wheel_station->getTH2F()->SetMinimum(0);
  h_recPunchOccupancyEndcap_disk_ring->getTH2F()->SetMinimum(0);
 
  for (int i = 1; i <= 5; ++i) {
    TString binLabel = Form("Wheel %d", i - 3);
    h_refPunchOccupancyBarrel_wheel->getTH1()->GetXaxis()->SetBinLabel(i, binLabel);
    h_refPunchOccupancyBarrel_wheel_station->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
    h_recPunchOccupancyBarrel_wheel->getTH1()->GetXaxis()->SetBinLabel(i, binLabel);
    h_recPunchOccupancyBarrel_wheel_station->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
  }
 
  for (int i = 1; i <= 9; ++i) {
    TString binLabel = Form("Disk %d", i - 5);
    h_refPunchOccupancyEndcap_disk->getTH1()->GetXaxis()->SetBinLabel(i, binLabel);
    h_refPunchOccupancyEndcap_disk_ring->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
    h_recPunchOccupancyEndcap_disk->getTH1()->GetXaxis()->SetBinLabel(i, binLabel);
    h_recPunchOccupancyEndcap_disk_ring->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
  }
 
  for (int i = 1; i <= 4; ++i) {
    TString binLabel = Form("Station %d", i);
    h_refPunchOccupancyBarrel_station->getTH1()->GetXaxis()->SetBinLabel(i, binLabel);
    h_refPunchOccupancyBarrel_wheel_station->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
    h_recPunchOccupancyBarrel_station->getTH1()->GetXaxis()->SetBinLabel(i, binLabel);
    h_recPunchOccupancyBarrel_wheel_station->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
  }
 
  for (int i = 1; i <= 4; ++i) {
    TString binLabel = Form("Ring %d", i);
    h_refPunchOccupancyEndcap_disk_ring->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
    h_recPunchOccupancyEndcap_disk_ring->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
  }
 
  // Book roll-by-roll histograms
  edm::ESHandle<RPCGeometry> rpcGeom;
  eventSetup.get<MuonGeometryRecord>().get(rpcGeom);
 
  int nRPCRollBarrel = 0, nRPCRollEndcap = 0;
 
  TrackingGeometry::DetContainer rpcDets = rpcGeom->dets();
  for (auto det : rpcDets) {
    auto rpcCh = dynamic_cast<const RPCChamber *>(det);
    if (!rpcCh)
      continue;
 
    std::vector<const RPCRoll *> rolls = rpcCh->rolls();
    for (auto roll : rolls) {
      if (!roll)
        continue;
 
      // RPCGeomServ rpcSrv(roll->id());
      const int rawId = roll->geographicalId().rawId();
      // if ( !roll->specs()->isRPC() ) { cout << "\nNoRPC : " << rpcSrv.name()
      // << ' ' << rawId << endl; continue; }
 
      if (roll->isBarrel()) {
        detIdToIndexMapBarrel_[rawId] = nRPCRollBarrel;
        // rollIdToNameMapBarrel_[rawId] = rpcSrv.name();
        ++nRPCRollBarrel;
      } else {
        detIdToIndexMapEndcap_[rawId] = nRPCRollEndcap;
        // rollIdToNameMapEndcap_[rawId] = rpcSrv.name();
        ++nRPCRollEndcap;
      }
    }
  }
 
  booker.setCurrentFolder(subDir_ + "/Occupancy");
  h_matchOccupancyBarrel_detId = booker.book1D("MatchOccupancyBarrel_detId",
                                               "Matched hit occupancy;roll index (can be arbitrary)",
                                               nRPCRollBarrel,
                                               0,
                                               nRPCRollBarrel);
  h_matchOccupancyEndcap_detId = booker.book1D("MatchOccupancyEndcap_detId",
                                               "Matched hit occupancy;roll index (can be arbitrary)",
                                               nRPCRollEndcap,
                                               0,
                                               nRPCRollEndcap);
  h_refOccupancyBarrel_detId = booker.book1D("RefOccupancyBarrel_detId",
                                             "Reference hit occupancy;roll index (can be arbitrary)",
                                             nRPCRollBarrel,
                                             0,
                                             nRPCRollBarrel);
  h_refOccupancyEndcap_detId = booker.book1D("RefOccupancyEndcap_detId",
                                             "Reference hit occupancy;roll index (can be arbitrary)",
                                             nRPCRollEndcap,
                                             0,
                                             nRPCRollEndcap);
  h_noiseOccupancyBarrel_detId = booker.book1D(
      "NoiseOccupancyBarrel_detId", "Noise occupancy;roll index (can be arbitrary)", nRPCRollBarrel, 0, nRPCRollBarrel);
  h_noiseOccupancyEndcap_detId = booker.book1D(
      "NoiseOccupancyEndcap_detId", "Noise occupancy;roll index (can be arbitrary)", nRPCRollEndcap, 0, nRPCRollEndcap);
 
  h_matchOccupancyBarrel_detId->getTH1()->SetMinimum(0);
  h_matchOccupancyEndcap_detId->getTH1()->SetMinimum(0);
  h_refOccupancyBarrel_detId->getTH1()->SetMinimum(0);
  h_refOccupancyEndcap_detId->getTH1()->SetMinimum(0);
  h_noiseOccupancyBarrel_detId->getTH1()->SetMinimum(0);
  h_noiseOccupancyEndcap_detId->getTH1()->SetMinimum(0);
 
  h_rollAreaBarrel_detId = booker.bookProfile(
      "RollAreaBarrel_detId", "Roll area;roll index;Area", nRPCRollBarrel, 0., 1. * nRPCRollBarrel, 0., 1e5);
  h_rollAreaEndcap_detId = booker.bookProfile(
      "RollAreaEndcap_detId", "Roll area;roll index;Area", nRPCRollEndcap, 0., 1. * nRPCRollEndcap, 0., 1e5);
 
  for (auto detIdToIndex : detIdToIndexMapBarrel_) {
    const int rawId = detIdToIndex.first;
    const int index = detIdToIndex.second;
 
    const RPCDetId rpcDetId = static_cast<const RPCDetId>(rawId);
    const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(rpcDetId));
 
    // RPCGeomServ rpcSrv(roll->id());
    // if ( !roll->specs()->isRPC() ) { cout << "\nNoRPC : " << rpcSrv.name() <<
    // ' ' << rawId << endl; continue; }
 
    const StripTopology &topol = roll->specificTopology();
    const double area = topol.stripLength() * topol.nstrips() * topol.pitch();
 
    h_rollAreaBarrel_detId->Fill(index, area);
  }
 
  for (auto detIdToIndex : detIdToIndexMapEndcap_) {
    const int rawId = detIdToIndex.first;
    const int index = detIdToIndex.second;
 
    const RPCDetId rpcDetId = static_cast<const RPCDetId>(rawId);
    const RPCRoll *roll = dynamic_cast<const RPCRoll *>(rpcGeom->roll(rpcDetId));
 
    // RPCGeomServ rpcSrv(roll->id());
    // if ( !roll->specs()->isRPC() ) { cout << "\nNoRPC : " << rpcSrv.name() <<
    // ' ' << rawId << endl; continue; }
 
    const StripTopology &topol = roll->specificTopology();
    const double area = topol.stripLength() * topol.nstrips() * topol.pitch();
 
    h_rollAreaEndcap_detId->Fill(index, area);
  }
}

References custom_jme_cff::area, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), RPCGeometry::dets(), options_cfi::eventSetup, get, dqm::legacy::MonitorElement::getTH1(), dqm::legacy::MonitorElement::getTH2F(), mps_fire::i, HLTTauDQMOffline_cfi::nPtBins, Pi, L1EGammaClusterEmuProducer_cfi::ptBins, RPCGeometry::roll(), dqm::implementation::NavigatorBase::setCurrentFolder(), and RPCRoll::specificTopology().

Member Data Documentation

detIdToIndexMapBarrel_

std::map<int, int> RPCRecHitValid::detIdToIndexMapBarrel_

private

Definition at line 74 of file RPCRecHitValid.h.

detIdToIndexMapEndcap_

std::map<int, int> RPCRecHitValid::detIdToIndexMapEndcap_

private

Definition at line 74 of file RPCRecHitValid.h.

h_

RPCValidHistograms RPCRecHitValid::h_

private

Definition at line 46 of file RPCRecHitValid.h.

h_eventCount

MEP RPCRecHitValid::h_eventCount

private

Definition at line 48 of file RPCRecHitValid.h.

h_matchOccupancyBarrel_detId

MEP RPCRecHitValid::h_matchOccupancyBarrel_detId

private

Definition at line 65 of file RPCRecHitValid.h.

h_matchOccupancyEndcap_detId

MEP RPCRecHitValid::h_matchOccupancyEndcap_detId

private

Definition at line 66 of file RPCRecHitValid.h.

h_noiseOccupancyBarrel_detId

MEP RPCRecHitValid::h_noiseOccupancyBarrel_detId

private

Definition at line 69 of file RPCRecHitValid.h.

h_noiseOccupancyEndcap_detId

MEP RPCRecHitValid::h_noiseOccupancyEndcap_detId

private

Definition at line 70 of file RPCRecHitValid.h.

h_nRPCHitPerRecoMuon

MEP RPCRecHitValid::h_nRPCHitPerRecoMuon

private

Definition at line 51 of file RPCRecHitValid.h.

h_nRPCHitPerRecoMuonBarrel

MEP RPCRecHitValid::h_nRPCHitPerRecoMuonBarrel

private

Definition at line 51 of file RPCRecHitValid.h.

h_nRPCHitPerRecoMuonEndcap

MEP RPCRecHitValid::h_nRPCHitPerRecoMuonEndcap

private

Definition at line 51 of file RPCRecHitValid.h.

h_nRPCHitPerRecoMuonOverlap

MEP RPCRecHitValid::h_nRPCHitPerRecoMuonOverlap

private

Definition at line 51 of file RPCRecHitValid.h.

h_nRPCHitPerSimMuon

MEP RPCRecHitValid::h_nRPCHitPerSimMuon

private

Definition at line 50 of file RPCRecHitValid.h.

h_nRPCHitPerSimMuonBarrel

MEP RPCRecHitValid::h_nRPCHitPerSimMuonBarrel

private

Definition at line 50 of file RPCRecHitValid.h.

h_nRPCHitPerSimMuonEndcap

MEP RPCRecHitValid::h_nRPCHitPerSimMuonEndcap

private

Definition at line 50 of file RPCRecHitValid.h.

h_nRPCHitPerSimMuonOverlap

MEP RPCRecHitValid::h_nRPCHitPerSimMuonOverlap

private

Definition at line 50 of file RPCRecHitValid.h.

h_recoMuonBarrel_eta

MEP RPCRecHitValid::h_recoMuonBarrel_eta

private

Definition at line 56 of file RPCRecHitValid.h.

h_recoMuonBarrel_phi

MEP RPCRecHitValid::h_recoMuonBarrel_phi

private

Definition at line 57 of file RPCRecHitValid.h.

h_recoMuonBarrel_pt

MEP RPCRecHitValid::h_recoMuonBarrel_pt

private

Definition at line 55 of file RPCRecHitValid.h.

h_recoMuonEndcap_eta

MEP RPCRecHitValid::h_recoMuonEndcap_eta

private

Definition at line 56 of file RPCRecHitValid.h.

h_recoMuonEndcap_phi

MEP RPCRecHitValid::h_recoMuonEndcap_phi

private

Definition at line 57 of file RPCRecHitValid.h.

h_recoMuonEndcap_pt

MEP RPCRecHitValid::h_recoMuonEndcap_pt

private

Definition at line 55 of file RPCRecHitValid.h.

h_recoMuonNoRPC_eta

MEP RPCRecHitValid::h_recoMuonNoRPC_eta

private

Definition at line 56 of file RPCRecHitValid.h.

h_recoMuonNoRPC_phi

MEP RPCRecHitValid::h_recoMuonNoRPC_phi

private

Definition at line 57 of file RPCRecHitValid.h.

h_recoMuonNoRPC_pt

MEP RPCRecHitValid::h_recoMuonNoRPC_pt

private

Definition at line 55 of file RPCRecHitValid.h.

h_recoMuonOverlap_eta

MEP RPCRecHitValid::h_recoMuonOverlap_eta

private

Definition at line 56 of file RPCRecHitValid.h.

h_recoMuonOverlap_phi

MEP RPCRecHitValid::h_recoMuonOverlap_phi

private

Definition at line 57 of file RPCRecHitValid.h.

h_recoMuonOverlap_pt

MEP RPCRecHitValid::h_recoMuonOverlap_pt

private

Definition at line 55 of file RPCRecHitValid.h.

h_recPunchOccupancyBarrel_station

MEP RPCRecHitValid::h_recPunchOccupancyBarrel_station

private

Definition at line 62 of file RPCRecHitValid.h.

h_recPunchOccupancyBarrel_wheel

MEP RPCRecHitValid::h_recPunchOccupancyBarrel_wheel

private

Definition at line 62 of file RPCRecHitValid.h.

h_recPunchOccupancyBarrel_wheel_station

MEP RPCRecHitValid::h_recPunchOccupancyBarrel_wheel_station

private

Definition at line 63 of file RPCRecHitValid.h.

h_recPunchOccupancyEndcap_disk

MEP RPCRecHitValid::h_recPunchOccupancyEndcap_disk

private

Definition at line 62 of file RPCRecHitValid.h.

h_recPunchOccupancyEndcap_disk_ring

MEP RPCRecHitValid::h_recPunchOccupancyEndcap_disk_ring

private

Definition at line 63 of file RPCRecHitValid.h.

h_refOccupancyBarrel_detId

MEP RPCRecHitValid::h_refOccupancyBarrel_detId

private

Definition at line 67 of file RPCRecHitValid.h.

h_refOccupancyEndcap_detId

MEP RPCRecHitValid::h_refOccupancyEndcap_detId

private

Definition at line 68 of file RPCRecHitValid.h.

h_refPunchOccupancyBarrel_station

MEP RPCRecHitValid::h_refPunchOccupancyBarrel_station

private

Definition at line 60 of file RPCRecHitValid.h.

h_refPunchOccupancyBarrel_wheel

MEP RPCRecHitValid::h_refPunchOccupancyBarrel_wheel

private

Definition at line 60 of file RPCRecHitValid.h.

h_refPunchOccupancyBarrel_wheel_station

MEP RPCRecHitValid::h_refPunchOccupancyBarrel_wheel_station

private

Definition at line 61 of file RPCRecHitValid.h.

h_refPunchOccupancyEndcap_disk

MEP RPCRecHitValid::h_refPunchOccupancyEndcap_disk

private

Definition at line 60 of file RPCRecHitValid.h.

h_refPunchOccupancyEndcap_disk_ring

MEP RPCRecHitValid::h_refPunchOccupancyEndcap_disk_ring

private

Definition at line 61 of file RPCRecHitValid.h.

h_rollAreaBarrel_detId

MEP RPCRecHitValid::h_rollAreaBarrel_detId

private

Definition at line 71 of file RPCRecHitValid.h.

h_rollAreaEndcap_detId

MEP RPCRecHitValid::h_rollAreaEndcap_detId

private

Definition at line 72 of file RPCRecHitValid.h.

h_simMuonBarrel_eta

MEP RPCRecHitValid::h_simMuonBarrel_eta

private

Definition at line 53 of file RPCRecHitValid.h.

h_simMuonBarrel_phi

MEP RPCRecHitValid::h_simMuonBarrel_phi

private

Definition at line 54 of file RPCRecHitValid.h.

h_simMuonBarrel_pt

MEP RPCRecHitValid::h_simMuonBarrel_pt

private

Definition at line 52 of file RPCRecHitValid.h.

h_simMuonEndcap_eta

MEP RPCRecHitValid::h_simMuonEndcap_eta

private

Definition at line 53 of file RPCRecHitValid.h.

h_simMuonEndcap_phi

MEP RPCRecHitValid::h_simMuonEndcap_phi

private

Definition at line 54 of file RPCRecHitValid.h.

h_simMuonEndcap_pt

MEP RPCRecHitValid::h_simMuonEndcap_pt

private

Definition at line 52 of file RPCRecHitValid.h.

h_simMuonNoRPC_eta

MEP RPCRecHitValid::h_simMuonNoRPC_eta

private

Definition at line 53 of file RPCRecHitValid.h.

h_simMuonNoRPC_phi

MEP RPCRecHitValid::h_simMuonNoRPC_phi

private

Definition at line 54 of file RPCRecHitValid.h.

h_simMuonNoRPC_pt

MEP RPCRecHitValid::h_simMuonNoRPC_pt

private

Definition at line 52 of file RPCRecHitValid.h.

h_simMuonOverlap_eta

MEP RPCRecHitValid::h_simMuonOverlap_eta

private

Definition at line 53 of file RPCRecHitValid.h.

h_simMuonOverlap_phi

MEP RPCRecHitValid::h_simMuonOverlap_phi

private

Definition at line 54 of file RPCRecHitValid.h.

h_simMuonOverlap_pt

MEP RPCRecHitValid::h_simMuonOverlap_pt

private

Definition at line 52 of file RPCRecHitValid.h.

h_simParticleType

MEP RPCRecHitValid::h_simParticleType

private

Definition at line 58 of file RPCRecHitValid.h.

h_simParticleTypeBarrel

MEP RPCRecHitValid::h_simParticleTypeBarrel

private

Definition at line 58 of file RPCRecHitValid.h.

h_simParticleTypeEndcap

MEP RPCRecHitValid::h_simParticleTypeEndcap

private

Definition at line 58 of file RPCRecHitValid.h.

muonToken_

edm::EDGetTokenT<reco::MuonCollection> RPCRecHitValid::muonToken_

private

Definition at line 43 of file RPCRecHitValid.h.

recHitToken_

edm::EDGetTokenT<RecHits> RPCRecHitValid::recHitToken_

private

Definition at line 40 of file RPCRecHitValid.h.

simHitAssocToken_

edm::EDGetTokenT<SimHitAssoc> RPCRecHitValid::simHitAssocToken_

private

Definition at line 42 of file RPCRecHitValid.h.

simHitToken_

edm::EDGetTokenT<SimHits> RPCRecHitValid::simHitToken_

private

Definition at line 39 of file RPCRecHitValid.h.

simParticleToken_

edm::EDGetTokenT<SimParticles> RPCRecHitValid::simParticleToken_

private

Definition at line 41 of file RPCRecHitValid.h.

subDir_

std::string RPCRecHitValid::subDir_

private

Definition at line 38 of file RPCRecHitValid.h.