MuonShowerDigiFiller.cc

Go to the documentation of this file.

//
// Package:    MuonShowerDigiFiller
// Class:      MuonShowerDigiFiller
// 
//
// Original Author:  Carlo Battilana, INFN BO 
//         Created:  Sat Mar 23 14:36:22 CET 2019
//
//


// system include files

// user include files
#include "RecoMuon/MuonIdentification/interface/MuonShowerDigiFiller.h"

#include "Geometry/Records/interface/MuonGeometryRecord.h"

//
// constructors and destructor
//

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

void 
MuonShowerDigiFiller::getES( const edm::EventSetup& iSetup )
{

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

}

void 
MuonShowerDigiFiller::getDigis( edm::Event& iEvent )
{
  
  iEvent.getByToken(m_dtDigisToken,  m_dtDigis);
  iEvent.getByToken(m_cscDigisToken, m_cscDigis);

}

void 
MuonShowerDigiFiller::fill( reco::MuonChamberMatch & muChMatch ) const 
{

  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 
{

  muChMatch.nDigisInRange = 0;

}