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
	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 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::ESGetToken< RPCGeometry, MuonGeometryRecord >	rpcGeomToken_
edm::ESGetToken< RPCGeometry, MuonGeometryRecord >	rpcGeomTokenInRun_
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 >
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

Definition at line 26 of file RPCRecHitValid.h.

Member Typedef Documentation

MEP

typedef MonitorElement* RPCRecHitValid::MEP

private

Definition at line 50 of file RPCRecHitValid.h.

RecHits

typedef RPCRecHitCollection RPCRecHitValid::RecHits

private

Definition at line 36 of file RPCRecHitValid.h.

SimHitAssoc

typedef SimHitTPAssociationProducer::SimHitTPAssociationList RPCRecHitValid::SimHitAssoc

private

Definition at line 38 of file RPCRecHitValid.h.

SimHits

typedef edm::PSimHitContainer RPCRecHitValid::SimHits

private

Definition at line 35 of file RPCRecHitValid.h.

SimParticles

typedef TrackingParticleCollection RPCRecHitValid::SimParticles

private

Definition at line 37 of file RPCRecHitValid.h.

Constructor & Destructor Documentation

RPCRecHitValid()

RPCRecHitValid::RPCRecHitValid

(

const edm::ParameterSet &

pset

)

Definition at line 24 of file RPCRecHitValid.cc.

References deDxTools::esConsumes(), muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

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

  rpcGeomToken_ = esConsumes();
  rpcGeomTokenInRun_ = esConsumes<edm::Transition::BeginRun>();
}

~RPCRecHitValid()

RPCRecHitValid::~RPCRecHitValid ( )

inlineoverride

Definition at line 29 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.

References funct::abs(), DetId::det(), fftjetvertexadder_cfi::errX, options_cfi::eventSetup, mps_fire::i, RPCRoll::id(), RPCRoll::isIRPC(), muonTagProbeFilters_cff::matched, DetId::Muon, HLT_2023v12_cff::muon, dqmiodumpmetadata::n, or, FastTimerService_cff::range, DetId::rawId(), rpcPointValidation_cfi::recHit, RPCDetId::region(), HLT_2023v12_cff::region, RPCDetId::ring(), relativeConstraints::ring, MuonSubdetId::RPC, rpcPointValidation_cfi::simHit, SimHitTPAssociationProducer::simHitTPAssociationListGreater(), mathSSE::sqrt(), relativeConstraints::station, RPCDetId::station(), DetId::subdetId(), and protons_cff::time.

                                                                                   {
  h_eventCount->Fill(1);

  // Get the RPC Geometry
  auto rpcGeom = eventSetup.getHandle(rpcGeomToken_);

  // 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 (const 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 (const 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 (const 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 (const 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 (const 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 (const 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 (const 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);
}

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 37 of file RPCRecHitValid.cc.

References custom_jme_cff::area, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), TrackValidation_cff::e5, options_cfi::eventSetup, dqm::legacy::MonitorElement::getTH1(), dqm::legacy::MonitorElement::getTH2F(), h, mps_fire::i, HLTTauDQMOffline_cfi::nPtBins, StripTopology::nstrips(), Pi, StripTopology::pitch(), L1EGammaClusterEmuProducer_cfi::ptBins, RPCChamber::rolls(), dqm::implementation::NavigatorBase::setCurrentFolder(), RPCRoll::specificTopology(), and StripTopology::stripLength().

                                                                                                               {
  // 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
  auto rpcGeom = eventSetup.getHandle(rpcGeomTokenInRun_);

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

Member Data Documentation

detIdToIndexMapBarrel_

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

private

Definition at line 79 of file RPCRecHitValid.h.

detIdToIndexMapEndcap_

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

private

Definition at line 79 of file RPCRecHitValid.h.

h_

RPCValidHistograms RPCRecHitValid::h_

private

Definition at line 51 of file RPCRecHitValid.h.

h_eventCount

MEP RPCRecHitValid::h_eventCount

private

Definition at line 53 of file RPCRecHitValid.h.

h_matchOccupancyBarrel_detId

MEP RPCRecHitValid::h_matchOccupancyBarrel_detId

private

Definition at line 70 of file RPCRecHitValid.h.

h_matchOccupancyEndcap_detId

MEP RPCRecHitValid::h_matchOccupancyEndcap_detId

private

Definition at line 71 of file RPCRecHitValid.h.

h_noiseOccupancyBarrel_detId

MEP RPCRecHitValid::h_noiseOccupancyBarrel_detId

private

Definition at line 74 of file RPCRecHitValid.h.

h_noiseOccupancyEndcap_detId

MEP RPCRecHitValid::h_noiseOccupancyEndcap_detId

private

Definition at line 75 of file RPCRecHitValid.h.

h_nRPCHitPerRecoMuon

MEP RPCRecHitValid::h_nRPCHitPerRecoMuon

private

Definition at line 56 of file RPCRecHitValid.h.

h_nRPCHitPerRecoMuonBarrel

MEP RPCRecHitValid::h_nRPCHitPerRecoMuonBarrel

private

Definition at line 56 of file RPCRecHitValid.h.

h_nRPCHitPerRecoMuonEndcap

MEP RPCRecHitValid::h_nRPCHitPerRecoMuonEndcap

private

Definition at line 56 of file RPCRecHitValid.h.

h_nRPCHitPerRecoMuonOverlap

MEP RPCRecHitValid::h_nRPCHitPerRecoMuonOverlap

private

Definition at line 56 of file RPCRecHitValid.h.

h_nRPCHitPerSimMuon

MEP RPCRecHitValid::h_nRPCHitPerSimMuon

private

Definition at line 55 of file RPCRecHitValid.h.

h_nRPCHitPerSimMuonBarrel

MEP RPCRecHitValid::h_nRPCHitPerSimMuonBarrel

private

Definition at line 55 of file RPCRecHitValid.h.

h_nRPCHitPerSimMuonEndcap

MEP RPCRecHitValid::h_nRPCHitPerSimMuonEndcap

private

Definition at line 55 of file RPCRecHitValid.h.

h_nRPCHitPerSimMuonOverlap

MEP RPCRecHitValid::h_nRPCHitPerSimMuonOverlap

private

Definition at line 55 of file RPCRecHitValid.h.

h_recoMuonBarrel_eta

MEP RPCRecHitValid::h_recoMuonBarrel_eta

private

Definition at line 61 of file RPCRecHitValid.h.

h_recoMuonBarrel_phi

MEP RPCRecHitValid::h_recoMuonBarrel_phi

private

Definition at line 62 of file RPCRecHitValid.h.

h_recoMuonBarrel_pt

MEP RPCRecHitValid::h_recoMuonBarrel_pt

private

Definition at line 60 of file RPCRecHitValid.h.

h_recoMuonEndcap_eta

MEP RPCRecHitValid::h_recoMuonEndcap_eta

private

Definition at line 61 of file RPCRecHitValid.h.

h_recoMuonEndcap_phi

MEP RPCRecHitValid::h_recoMuonEndcap_phi

private

Definition at line 62 of file RPCRecHitValid.h.

h_recoMuonEndcap_pt

MEP RPCRecHitValid::h_recoMuonEndcap_pt

private

Definition at line 60 of file RPCRecHitValid.h.

h_recoMuonNoRPC_eta

MEP RPCRecHitValid::h_recoMuonNoRPC_eta

private

Definition at line 61 of file RPCRecHitValid.h.

h_recoMuonNoRPC_phi

MEP RPCRecHitValid::h_recoMuonNoRPC_phi

private

Definition at line 62 of file RPCRecHitValid.h.

h_recoMuonNoRPC_pt

MEP RPCRecHitValid::h_recoMuonNoRPC_pt

private

Definition at line 60 of file RPCRecHitValid.h.

h_recoMuonOverlap_eta

MEP RPCRecHitValid::h_recoMuonOverlap_eta

private

Definition at line 61 of file RPCRecHitValid.h.

h_recoMuonOverlap_phi

MEP RPCRecHitValid::h_recoMuonOverlap_phi

private

Definition at line 62 of file RPCRecHitValid.h.

h_recoMuonOverlap_pt

MEP RPCRecHitValid::h_recoMuonOverlap_pt

private

Definition at line 60 of file RPCRecHitValid.h.

h_recPunchOccupancyBarrel_station

MEP RPCRecHitValid::h_recPunchOccupancyBarrel_station

private

Definition at line 67 of file RPCRecHitValid.h.

h_recPunchOccupancyBarrel_wheel

MEP RPCRecHitValid::h_recPunchOccupancyBarrel_wheel

private

Definition at line 67 of file RPCRecHitValid.h.

h_recPunchOccupancyBarrel_wheel_station

MEP RPCRecHitValid::h_recPunchOccupancyBarrel_wheel_station

private

Definition at line 68 of file RPCRecHitValid.h.

h_recPunchOccupancyEndcap_disk

MEP RPCRecHitValid::h_recPunchOccupancyEndcap_disk

private

Definition at line 67 of file RPCRecHitValid.h.

h_recPunchOccupancyEndcap_disk_ring

MEP RPCRecHitValid::h_recPunchOccupancyEndcap_disk_ring

private

Definition at line 68 of file RPCRecHitValid.h.

h_refOccupancyBarrel_detId

MEP RPCRecHitValid::h_refOccupancyBarrel_detId

private

Definition at line 72 of file RPCRecHitValid.h.

h_refOccupancyEndcap_detId

MEP RPCRecHitValid::h_refOccupancyEndcap_detId

private

Definition at line 73 of file RPCRecHitValid.h.

h_refPunchOccupancyBarrel_station

MEP RPCRecHitValid::h_refPunchOccupancyBarrel_station

private

Definition at line 65 of file RPCRecHitValid.h.

h_refPunchOccupancyBarrel_wheel

MEP RPCRecHitValid::h_refPunchOccupancyBarrel_wheel

private

Definition at line 65 of file RPCRecHitValid.h.

h_refPunchOccupancyBarrel_wheel_station

MEP RPCRecHitValid::h_refPunchOccupancyBarrel_wheel_station

private

Definition at line 66 of file RPCRecHitValid.h.

h_refPunchOccupancyEndcap_disk

MEP RPCRecHitValid::h_refPunchOccupancyEndcap_disk

private

Definition at line 65 of file RPCRecHitValid.h.

h_refPunchOccupancyEndcap_disk_ring

MEP RPCRecHitValid::h_refPunchOccupancyEndcap_disk_ring

private

Definition at line 66 of file RPCRecHitValid.h.

h_rollAreaBarrel_detId

MEP RPCRecHitValid::h_rollAreaBarrel_detId

private

Definition at line 76 of file RPCRecHitValid.h.

h_rollAreaEndcap_detId

MEP RPCRecHitValid::h_rollAreaEndcap_detId

private

Definition at line 77 of file RPCRecHitValid.h.

h_simMuonBarrel_eta

MEP RPCRecHitValid::h_simMuonBarrel_eta

private

Definition at line 58 of file RPCRecHitValid.h.

h_simMuonBarrel_phi

MEP RPCRecHitValid::h_simMuonBarrel_phi

private

Definition at line 59 of file RPCRecHitValid.h.

h_simMuonBarrel_pt

MEP RPCRecHitValid::h_simMuonBarrel_pt

private

Definition at line 57 of file RPCRecHitValid.h.

h_simMuonEndcap_eta

MEP RPCRecHitValid::h_simMuonEndcap_eta

private

Definition at line 58 of file RPCRecHitValid.h.

h_simMuonEndcap_phi

MEP RPCRecHitValid::h_simMuonEndcap_phi

private

Definition at line 59 of file RPCRecHitValid.h.

h_simMuonEndcap_pt

MEP RPCRecHitValid::h_simMuonEndcap_pt

private

Definition at line 57 of file RPCRecHitValid.h.

h_simMuonNoRPC_eta

MEP RPCRecHitValid::h_simMuonNoRPC_eta

private

Definition at line 58 of file RPCRecHitValid.h.

h_simMuonNoRPC_phi

MEP RPCRecHitValid::h_simMuonNoRPC_phi

private

Definition at line 59 of file RPCRecHitValid.h.

h_simMuonNoRPC_pt

MEP RPCRecHitValid::h_simMuonNoRPC_pt

private

Definition at line 57 of file RPCRecHitValid.h.

h_simMuonOverlap_eta

MEP RPCRecHitValid::h_simMuonOverlap_eta

private

Definition at line 58 of file RPCRecHitValid.h.

h_simMuonOverlap_phi

MEP RPCRecHitValid::h_simMuonOverlap_phi

private

Definition at line 59 of file RPCRecHitValid.h.

h_simMuonOverlap_pt

MEP RPCRecHitValid::h_simMuonOverlap_pt

private

Definition at line 57 of file RPCRecHitValid.h.

h_simParticleType

MEP RPCRecHitValid::h_simParticleType

private

Definition at line 63 of file RPCRecHitValid.h.

h_simParticleTypeBarrel

MEP RPCRecHitValid::h_simParticleTypeBarrel

private

Definition at line 63 of file RPCRecHitValid.h.

h_simParticleTypeEndcap

MEP RPCRecHitValid::h_simParticleTypeEndcap

private

Definition at line 63 of file RPCRecHitValid.h.

muonToken_

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

private

Definition at line 45 of file RPCRecHitValid.h.

recHitToken_

edm::EDGetTokenT<RecHits> RPCRecHitValid::recHitToken_

private

Definition at line 42 of file RPCRecHitValid.h.

rpcGeomToken_

edm::ESGetToken<RPCGeometry, MuonGeometryRecord> RPCRecHitValid::rpcGeomToken_

private

Definition at line 47 of file RPCRecHitValid.h.

rpcGeomTokenInRun_

edm::ESGetToken<RPCGeometry, MuonGeometryRecord> RPCRecHitValid::rpcGeomTokenInRun_

private

Definition at line 48 of file RPCRecHitValid.h.

simHitAssocToken_

edm::EDGetTokenT<SimHitAssoc> RPCRecHitValid::simHitAssocToken_

private

Definition at line 44 of file RPCRecHitValid.h.

simHitToken_

edm::EDGetTokenT<SimHits> RPCRecHitValid::simHitToken_

private

Definition at line 41 of file RPCRecHitValid.h.

simParticleToken_

edm::EDGetTokenT<SimParticles> RPCRecHitValid::simParticleToken_

private

Definition at line 43 of file RPCRecHitValid.h.

subDir_

std::string RPCRecHitValid::subDir_

private

Definition at line 40 of file RPCRecHitValid.h.