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 (const edm::ParameterSet &, edm::ConsumesCollector &&iC)

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

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::MuonShowerDigiFiller	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 30 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")))
{
   
}

Member Function Documentation

void MuonShowerDigiFiller::fill

(

reco::MuonChamberMatch &

muChMatch

)

const

Definition at line 62 of file MuonShowerDigiFiller.cc.

References funct::abs(), MuonSubdetId::CSC, reco::MuonChamberMatch::detector(), diffTreeTool::diff, MuonSubdetId::DT, objects.autophobj::float, CSCLayer::geometry(), reco::MuonChamberMatch::id, CSCGeometry::layer(), DTGeometry::layer(), m_cscDigis, m_cscGeometry, m_digiMaxDistanceX, m_dtDigis, m_dtGeometry, DTChamberId::maxLayerId, CSCDetId::maxLayerId(), DTChamberId::maxSuperLayerId, reco::MuonChamberMatch::nDigisInRange, muonCSCDigis_cfi::pedestal, DetId::rawId(), DTLayer::specificTopology(), electronIdCutBased_cfi::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;
  
}

void MuonShowerDigiFiller::fillDefault

(

reco::MuonChamberMatch &

muChMatch

)

const

Definition at line 145 of file MuonShowerDigiFiller.cc.

References reco::MuonChamberMatch::nDigisInRange.

{

  muChMatch.nDigisInRange = 0;

}

void MuonShowerDigiFiller::getDigis

(

edm::Event &

iEvent

)

Definition at line 53 of file MuonShowerDigiFiller.cc.

References edm::Event::getByToken(), m_cscDigis, m_cscDigisToken, m_dtDigis, and m_dtDigisToken.

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

}

void MuonShowerDigiFiller::getES

(

const edm::EventSetup &

iSetup

)

Definition at line 44 of file MuonShowerDigiFiller.cc.

References edm::EventSetup::get(), m_cscGeometry, and m_dtGeometry.

{

  iSetup.get<MuonGeometryRecord>().get(m_dtGeometry);
  iSetup.get<MuonGeometryRecord>().get(m_cscGeometry);

}

Member Data Documentation

edm::Handle<CSCStripDigiCollection> MuonShowerDigiFiller::m_cscDigis

private

Definition at line 71 of file MuonShowerDigiFiller.h.

Referenced by fill(), and getDigis().

edm::EDGetTokenT<CSCStripDigiCollection> MuonShowerDigiFiller::m_cscDigisToken

private

Definition at line 65 of file MuonShowerDigiFiller.h.

Referenced by getDigis().

edm::ESHandle<CSCGeometry> MuonShowerDigiFiller::m_cscGeometry

private

Definition at line 68 of file MuonShowerDigiFiller.h.

Referenced by fill(), and getES().

double MuonShowerDigiFiller::m_digiMaxDistanceX

private

Definition at line 62 of file MuonShowerDigiFiller.h.

Referenced by fill().

edm::Handle<DTDigiCollection> MuonShowerDigiFiller::m_dtDigis

private

Definition at line 70 of file MuonShowerDigiFiller.h.

Referenced by fill(), and getDigis().

edm::EDGetTokenT<DTDigiCollection> MuonShowerDigiFiller::m_dtDigisToken

private

Definition at line 64 of file MuonShowerDigiFiller.h.

Referenced by getDigis().

edm::ESHandle<DTGeometry> MuonShowerDigiFiller::m_dtGeometry

private

Definition at line 67 of file MuonShowerDigiFiller.h.

Referenced by fill(), and getES().