#include <CSCChamberMasker.cc>

Inheritance diagram for CSCChamberMasker:

Public Member Functions
	CSCChamberMasker (const edm::ParameterSet &)

	~CSCChamberMasker () override

Public Member Functions inherited from edm::stream::EDProducer<>
	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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &)

Private Member Functions
template<typename T , typename C = MuonDigiCollection<CSCDetId, T>>
void	ageDigis (edm::Event &event, edm::EDGetTokenT< C > &digiToken, CLHEP::HepRandomEngine &randGen, std::unique_ptr< C > &filteredDigis)

void	beginRun (edm::Run const &, edm::EventSetup const &) override

template<typename T , typename C = MuonDigiCollection<CSCDetId, T>>
void	copyDigis (edm::Event &event, edm::EDGetTokenT< C > &digiToken, std::unique_ptr< C > &filteredDigis)

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::ESGetToken < MuonSystemAging, MuonSystemAgingRcd >	m_agingObjToken

edm::EDGetTokenT < CSCALCTDigiCollection >	m_alctDigiToken

edm::EDGetTokenT < CSCCLCTDigiCollection >	m_clctDigiToken

std::map< CSCDetId, std::pair < unsigned int, float > >	m_CSCEffs

edm::ESGetToken< CSCGeometry, MuonGeometryRecord >	m_cscGeomToken

edm::EDGetTokenT < CSCStripDigiCollection >	m_stripDigiToken

edm::EDGetTokenT < CSCWireDigiCollection >	m_wireDigiToken

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
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

Detailed Description

SimMuon/CSCDigitizer/plugins/CSCChamberMasker.cc

Description: Class to mask CSC digis on a chamber by chamber basis

Definition at line 69 of file CSCChamberMasker.cc.

Constructor & Destructor Documentation

CSCChamberMasker::CSCChamberMasker ( const edm::ParameterSet & iConfig )

explicit

Definition at line 112 of file CSCChamberMasker.cc.

     : m_stripDigiToken(consumes<CSCStripDigiCollection>(iConfig.getParameter<edm::InputTag>("stripDigiTag"))),
       m_wireDigiToken(consumes<CSCWireDigiCollection>(iConfig.getParameter<edm::InputTag>("wireDigiTag"))),
       m_clctDigiToken(consumes<CSCCLCTDigiCollection>(iConfig.getParameter<edm::InputTag>("clctDigiTag"))),
       m_alctDigiToken(consumes<CSCALCTDigiCollection>(iConfig.getParameter<edm::InputTag>("alctDigiTag"))),
       m_cscGeomToken(esConsumes<CSCGeometry, MuonGeometryRecord>()),
       m_agingObjToken(esConsumes<MuonSystemAging, MuonSystemAgingRcd>()) {
   produces<CSCStripDigiCollection>("MuonCSCStripDigi");
   produces<CSCWireDigiCollection>("MuonCSCWireDigi");
   produces<CSCCLCTDigiCollection>("MuonCSCCLCTDigi");
   produces<CSCALCTDigiCollection>("MuonCSCALCTDigi");
 }

CSCChamberMasker::~CSCChamberMasker ( )

override

Definition at line 125 of file CSCChamberMasker.cc.

125 {}

Member Function Documentation

template<typename T , typename C >

void CSCChamberMasker::ageDigis	(	edm::Event &	event,
		edm::EDGetTokenT< C > &	digiToken,
		CLHEP::HepRandomEngine &	randGen,
		std::unique_ptr< C > &	filteredDigis
	)

private

Definition at line 176 of file CSCChamberMasker.cc.

References CSCDetId::chamber(), CSCDetId, EFF_CHAMBER, EFF_STRIPS, EFF_WIRES, CSCDetId::endcap(), edm::EDGetTokenT< T >::isUninitialized(), dqmiolumiharvest::j, dqmdumpme::last, CSCDetId::layer(), phase1PixelTopology::layer, m_CSCEffs, CSCDetId::ring(), CSCDetId::station(), and relativeConstraints::value.

                                                                  {
   if (!digiToken.isUninitialized()) {
     edm::Handle<C> digis;
     event.getByToken(digiToken, digis);
 
     for (const auto &j : (*digis)) {
       auto digiItr = j.second.first;
       auto last = j.second.second;
 
       CSCDetId const cscDetId = j.first;
 
       // Since lookups are chamber-centric, make new DetId with layer=0
       CSCDetId chId = CSCDetId(cscDetId.endcap(), cscDetId.station(), cscDetId.ring(), cscDetId.chamber(), 0);
 
       for (; digiItr != last; ++digiItr) {
         auto chEffIt = m_CSCEffs.find(chId);
 
         if (chEffIt != m_CSCEffs.end()) {
           std::pair<unsigned int, float> typeEff = chEffIt->second;
           int type = typeEff.first % 10;   // second digit gives type of inefficiency
           int layer = typeEff.first / 10;  // first digit gives layer (0 = chamber level)
 
           bool doRandomize = false;
           if (((std::is_same<T, CSCStripDigi>::value && type == EFF_WIRES) ||
                (std::is_same<T, CSCWireDigi>::value && type == EFF_STRIPS) || type == EFF_CHAMBER) &&
               (layer == 0 || cscDetId.layer() == layer))
             doRandomize = true;
 
           if (!doRandomize || (randGen.flat() <= typeEff.second)) {
             filteredDigis->insertDigi(cscDetId, *digiItr);
           }
         }
       }
     }
   }
 }

void CSCChamberMasker::beginRun	(	edm::Run const &	run,
		edm::EventSetup const &	iSetup
	)

overrideprivate

Definition at line 217 of file CSCChamberMasker.cc.

References CSCDetId::chamber(), chambers, CSCDetId::endcap(), edm::EventSetup::getHandle(), m_agingObjToken, m_CSCEffs, m_cscGeomToken, CSCDetId::rawIdMaker(), CSCDetId::ring(), and CSCDetId::station().

                                                                               {
   m_CSCEffs.clear();
 
   edm::ESHandle<CSCGeometry> cscGeom = iSetup.getHandle(m_cscGeomToken);
 
   edm::ESHandle<MuonSystemAging> agingObj = iSetup.getHandle(m_agingObjToken);
 
   const auto chambers = cscGeom->chambers();
 
   for (const auto *ch : chambers) {
     CSCDetId chId = ch->id();
     unsigned int rawId = chId.rawIdMaker(chId.endcap(), chId.station(), chId.ring(), chId.chamber(), 0);
     float eff = 1.;
     int type = 0;
     for (auto &agingPair : agingObj->m_CSCChambEffs) {
       if (agingPair.first != rawId)
         continue;
 
       type = agingPair.second.first;
       eff = agingPair.second.second;
       m_CSCEffs[chId] = std::make_pair(type, eff);
       break;
     }
   }
 }

template<typename T , typename C >

void CSCChamberMasker::copyDigis	(	edm::Event &	event,
		edm::EDGetTokenT< C > &	digiToken,
		std::unique_ptr< C > &	filteredDigis
	)

private

Definition at line 157 of file CSCChamberMasker.cc.

References edm::EDGetTokenT< T >::isUninitialized(), dqmiolumiharvest::j, and dqmdumpme::last.

                                                                                                                {
   if (!digiToken.isUninitialized()) {
     edm::Handle<C> digis;
     event.getByToken(digiToken, digis);
     for (const auto &j : (*digis)) {
       auto digiItr = j.second.first;
       auto last = j.second.second;
 
       CSCDetId const cscDetId = j.first;
 
       for (; digiItr != last; ++digiItr) {
         filteredDigis->insertDigi(cscDetId, *digiItr);
       }
     }
   }
 }

void CSCChamberMasker::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 245 of file CSCChamberMasker.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

                                                                                   {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("stripDigiTag", edm::InputTag("simMuonCSCDigis:MuonCSCStripDigi"));
   desc.add<edm::InputTag>("wireDigiTag", edm::InputTag("simMuonCSCDigis:MuonCSCWireDigi"));
   desc.add<edm::InputTag>("comparatorDigiTag", edm::InputTag("simMuonCSCDigis:MuonCSCComparatorDigi"));
   desc.add<edm::InputTag>("rpcDigiTag", edm::InputTag("simMuonCSCDigis:MuonCSCRPCDigi"));
   desc.add<edm::InputTag>("alctDigiTag", edm::InputTag("simMuonCSCDigis:MuonCSCALCTDigi"));
   desc.add<edm::InputTag>("clctDigiTag", edm::InputTag("simMuonCSCDigis:MuonCSCCLCTDigi"));
   descriptions.add("cscChamberMasker", desc);
 }

void CSCChamberMasker::produce	(	edm::Event &	event,
		const edm::EventSetup &	conditions
	)

overrideprivate

Definition at line 132 of file CSCChamberMasker.cc.

References edmPickEvents::event, edm::RandomNumberGenerator::getEngine(), m_alctDigiToken, m_clctDigiToken, m_stripDigiToken, m_wireDigiToken, eostools::move(), and edm::Event::streamID().

                                                                                {
   edm::Service<edm::RandomNumberGenerator> randGenService;
   CLHEP::HepRandomEngine &randGen = randGenService->getEngine(event.streamID());
 
   std::unique_ptr<CSCStripDigiCollection> filteredStripDigis(new CSCStripDigiCollection());
   std::unique_ptr<CSCWireDigiCollection> filteredWireDigis(new CSCWireDigiCollection());
   std::unique_ptr<CSCCLCTDigiCollection> filteredCLCTDigis(new CSCCLCTDigiCollection());
   std::unique_ptr<CSCALCTDigiCollection> filteredALCTDigis(new CSCALCTDigiCollection());
 
   // Handle wire and strip digis
   ageDigis<CSCStripDigi>(event, m_stripDigiToken, randGen, filteredStripDigis);
   ageDigis<CSCWireDigi>(event, m_wireDigiToken, randGen, filteredWireDigis);
 
   // Don't touch CLCT or ALCT digis
   copyDigis<CSCCLCTDigi>(event, m_clctDigiToken, filteredCLCTDigis);
   copyDigis<CSCALCTDigi>(event, m_alctDigiToken, filteredALCTDigis);
 
   event.put(std::move(filteredStripDigis), "MuonCSCStripDigi");
   event.put(std::move(filteredWireDigis), "MuonCSCWireDigi");
   event.put(std::move(filteredCLCTDigis), "MuonCSCCLCTDigi");
   event.put(std::move(filteredALCTDigis), "MuonCSCALCTDigi");
 }