MuonShowerDigiFiller Class Reference

#include <RecoMuon/MuonIdentification/src/MuonShowerDigiFiller.cc>

Public Member Functions
void	fill (reco::MuonChamberMatch &muChMatch) const
void	fillDefault (reco::MuonChamberMatch &muChMatch) const
void	getDigis (edm::Event &iEvent)
void	getES (const edm::EventSetup &iSetup)
	MuonShowerDigiFiller ()
	MuonShowerDigiFiller (const edm::ParameterSet &, edm::ConsumesCollector &&iC)

Private Attributes
edm::Handle< CSCStripDigiCollection >	m_cscDigis
edm::EDGetTokenT< CSCStripDigiCollection >	m_cscDigisToken
edm::ESHandle< CSCGeometry >	m_cscGeometry
edm::ESGetToken< CSCGeometry, MuonGeometryRecord >	m_cscGeometryToken
double	m_digiMaxDistanceX
edm::Handle< DTDigiCollection >	m_dtDigis
edm::EDGetTokenT< DTDigiCollection >	m_dtDigisToken
edm::ESHandle< DTGeometry >	m_dtGeometry
edm::ESGetToken< DTGeometry, MuonGeometryRecord >	m_dtGeometryToken

Detailed Description

Description: Class filling shower information using DT and CSC digis

Implementation: <Notes on="" implementation>="">

Definition at line 47 of file MuonShowerDigiFiller.h.

Constructor & Destructor Documentation

MuonShowerDigiFiller() [1/2]

MuonShowerDigiFiller::MuonShowerDigiFiller ( )

inline

Definition at line 49 of file MuonShowerDigiFiller.h.

{}

MuonShowerDigiFiller() [2/2]

MuonShowerDigiFiller::MuonShowerDigiFiller	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 27 of file MuonShowerDigiFiller.cc.

    : m_digiMaxDistanceX(iConfig.getParameter<double>("digiMaxDistanceX")),
      m_dtDigisToken(iC.consumes<DTDigiCollection>(iConfig.getParameter<edm::InputTag>("dtDigiCollectionLabel"))),
      m_cscDigisToken(
          iC.consumes<CSCStripDigiCollection>(iConfig.getParameter<edm::InputTag>("cscDigiCollectionLabel"))),
      m_dtGeometryToken(iC.esConsumes<edm::Transition::BeginRun>()),
      m_cscGeometryToken(iC.esConsumes<edm::Transition::BeginRun>()) {}

Member Function Documentation

fill()

void MuonShowerDigiFiller::fill

(

reco::MuonChamberMatch &

muChMatch

)

const

Definition at line 49 of file MuonShowerDigiFiller.cc.

References funct::abs(), MuonSubdetId::CSC, reco::MuonChamberMatch::detector(), change_name::diff, MuonSubdetId::DT, dqmMemoryStats::float, CSCLayer::geometry(), reco::MuonChamberMatch::id, CSCGeometry::layer(), DTGeometry::layer(), phase1PixelTopology::layer, m_cscDigis, m_cscGeometry, m_digiMaxDistanceX, m_dtDigis, m_dtGeometry, DTChamberId::maxLayerId, CSCDetId::maxLayerId(), DTChamberId::maxSuperLayerId, reco::MuonChamberMatch::nDigisInRange, EcalCondDBWriter_cfi::pedestal, FastTimerService_cff::range, DetId::rawId(), DTLayer::specificTopology(), remoteMonitoring_LASER_era2018_cfg::threshold, DTTopology::wirePosition(), reco::MuonChamberMatch::x, CSCLayerGeometry::xOfStrip(), and reco::MuonChamberMatch::y.

                                                                     {
  int nDigisInRange = 0;

  // DT chamber
  if (muChMatch.detector() == MuonSubdetId::DT) {
    double xTrack = muChMatch.x;

    for (int sl = 1; sl <= DTChamberId::maxSuperLayerId; sl += 2) {
      for (int layer = 1; layer <= DTChamberId::maxLayerId; ++layer) {
        const DTLayerId layerId(DTChamberId(muChMatch.id.rawId()), sl, layer);

        auto range = m_dtDigis->get(layerId);

        for (auto digiIt = range.first; digiIt != range.second; ++digiIt) {
          const auto topo = m_dtGeometry->layer(layerId)->specificTopology();

          double xWire = topo.wirePosition((*digiIt).wire());
          double dX = std::abs(xWire - xTrack);

          if (dX < m_digiMaxDistanceX)
            nDigisInRange++;
        }
      }
    }
  }

  else if (muChMatch.detector() == MuonSubdetId::CSC) {
    double xTrack = muChMatch.x;
    double yTrack = muChMatch.y;

    for (int iLayer = 1; iLayer <= CSCDetId::maxLayerId(); ++iLayer) {
      const CSCDetId chId(muChMatch.id.rawId());
      const CSCDetId layerId(chId.endcap(), chId.station(), chId.ring(), chId.chamber(), iLayer);

      auto range = m_cscDigis->get(layerId);

      for (auto digiIt = range.first; digiIt != range.second; ++digiIt) {
        std::vector<int> adcVals = digiIt->getADCCounts();
        bool hasFired = false;
        float pedestal = 0.5 * (float)(adcVals[0] + adcVals[1]);
        float threshold = 13.3;
        float diff = 0.;
        for (const auto& adcVal : adcVals) {
          diff = (float)adcVal - pedestal;
          if (diff > threshold) {
            hasFired = true;
            break;
          }
        }

        if (!hasFired)
          continue;

        const CSCLayerGeometry* layerGeom = m_cscGeometry->layer(layerId)->geometry();

        Float_t xStrip = layerGeom->xOfStrip(digiIt->getStrip(), yTrack);
        float dX = std::abs(xStrip - xTrack);

        if (dX < m_digiMaxDistanceX)
          nDigisInRange++;
      }
    }
  }

  muChMatch.nDigisInRange = nDigisInRange;
}

fillDefault()

void MuonShowerDigiFiller::fillDefault

(

reco::MuonChamberMatch &

muChMatch

)

const

Definition at line 116 of file MuonShowerDigiFiller.cc.

References reco::MuonChamberMatch::nDigisInRange.

{ muChMatch.nDigisInRange = 0; }

getDigis()

void MuonShowerDigiFiller::getDigis

(

edm::Event &

iEvent

)

Definition at line 44 of file MuonShowerDigiFiller.cc.

References iEvent, m_cscDigis, m_cscDigisToken, m_dtDigis, and m_dtDigisToken.

                                                    {
  iEvent.getByToken(m_dtDigisToken, m_dtDigis);
  iEvent.getByToken(m_cscDigisToken, m_cscDigis);
}

getES()

void MuonShowerDigiFiller::getES

(

const edm::EventSetup &

iSetup

)

Definition at line 39 of file MuonShowerDigiFiller.cc.

References edm::EventSetup::getHandle(), m_cscGeometry, m_cscGeometryToken, m_dtGeometry, and m_dtGeometryToken.

                                                            {
  m_dtGeometry = iSetup.getHandle(m_dtGeometryToken);
  m_cscGeometry = iSetup.getHandle(m_cscGeometryToken);
}

Member Data Documentation

m_cscDigis

edm::Handle<CSCStripDigiCollection> MuonShowerDigiFiller::m_cscDigis

private

Definition at line 71 of file MuonShowerDigiFiller.h.

Referenced by fill(), and getDigis().

m_cscDigisToken

edm::EDGetTokenT<CSCStripDigiCollection> MuonShowerDigiFiller::m_cscDigisToken

private

Definition at line 62 of file MuonShowerDigiFiller.h.

Referenced by getDigis().

m_cscGeometry

edm::ESHandle<CSCGeometry> MuonShowerDigiFiller::m_cscGeometry

private

Definition at line 68 of file MuonShowerDigiFiller.h.

Referenced by fill(), and getES().

m_cscGeometryToken

edm::ESGetToken<CSCGeometry, MuonGeometryRecord> MuonShowerDigiFiller::m_cscGeometryToken

private

Definition at line 65 of file MuonShowerDigiFiller.h.

Referenced by getES().

m_digiMaxDistanceX

double MuonShowerDigiFiller::m_digiMaxDistanceX

private

Definition at line 59 of file MuonShowerDigiFiller.h.

Referenced by fill().

m_dtDigis

edm::Handle<DTDigiCollection> MuonShowerDigiFiller::m_dtDigis

private

Definition at line 70 of file MuonShowerDigiFiller.h.

Referenced by fill(), and getDigis().

m_dtDigisToken

edm::EDGetTokenT<DTDigiCollection> MuonShowerDigiFiller::m_dtDigisToken

private

Definition at line 61 of file MuonShowerDigiFiller.h.

Referenced by getDigis().

m_dtGeometry

edm::ESHandle<DTGeometry> MuonShowerDigiFiller::m_dtGeometry

private

Definition at line 67 of file MuonShowerDigiFiller.h.

Referenced by fill(), and getES().

m_dtGeometryToken

edm::ESGetToken<DTGeometry, MuonGeometryRecord> MuonShowerDigiFiller::m_dtGeometryToken

private

Definition at line 64 of file MuonShowerDigiFiller.h.

Referenced by getES().