#include <CSCDigitizer.h>

Public Types
typedef edm::DetSetVector< StripDigiSimLink >	DigiSimLinks

Public Member Functions
	CSCDigitizer (const CSCDigitizer &)=delete

	CSCDigitizer (const edm::ParameterSet &p)
	configurable parameters More...

void	doAction (MixCollection< PSimHit > &simHits, CSCWireDigiCollection &wireDigis, CSCStripDigiCollection &stripDigis, CSCComparatorDigiCollection &comparators, DigiSimLinks &wireDigiSimLinks, DigiSimLinks &stripDigiSimLinks, CLHEP::HepRandomEngine *)

CSCDigitizer &	operator= (const CSCDigitizer &)=delete

void	setGeometry (const CSCGeometry *geom)
	sets geometry More...

void	setMagneticField (const MagneticField *field)
	sets the magnetic field More...

void	setParticleDataTable (const ParticleDataTable *pdt)

void	setStripConditions (CSCStripConditions *cond)

	~CSCDigitizer ()

Private Member Functions
const CSCLayer *	findLayer (int detId) const
	finds the layer in the geometry associated with this det ID More...

std::list< int >	layersMissing (const CSCStripDigiCollection &stripDigis) const
	finds which layers, 1-6, aren't in the current list More...

Private Attributes
bool	digitizeBadChambers_

CSCStripConditions *	theConditions

const CSCGeometry *	theCSCGeometry

CSCDriftSim *	theDriftSim

unsigned int	theLayersNeeded

CSCNeutronReader *	theNeutronReader

CSCStripElectronicsSim *	theStripElectronicsSim

CSCStripHitSim *	theStripHitSim

CSCWireElectronicsSim *	theWireElectronicsSim

CSCWireHitSim *	theWireHitSim

Detailed Description

Digitizer class for endcap muon CSCs.

Author: Rick Wilkinson

Definition at line 40 of file CSCDigitizer.h.

Member Typedef Documentation

◆ DigiSimLinks

typedef edm::DetSetVector<StripDigiSimLink> CSCDigitizer::DigiSimLinks

Definition at line 42 of file CSCDigitizer.h.

Constructor & Destructor Documentation

◆ ~CSCDigitizer()

CSCDigitizer::~CSCDigitizer ( )

Definition at line 30 of file CSCDigitizer.cc.

                             {
   delete theNeutronReader;
   delete theStripElectronicsSim;
   delete theWireElectronicsSim;
   delete theStripHitSim;
   delete theWireHitSim;
   delete theDriftSim;
 }

References theDriftSim, theNeutronReader, theStripElectronicsSim, theStripHitSim, theWireElectronicsSim, and theWireHitSim.

◆ CSCDigitizer() [1/2]

CSCDigitizer::CSCDigitizer ( const CSCDigitizer & )

delete

◆ CSCDigitizer() [2/2]

CSCDigitizer::CSCDigitizer ( const edm::ParameterSet & p )

explicit

configurable parameters

Definition at line 15 of file CSCDigitizer.cc.

     : theDriftSim(new CSCDriftSim()),
       theWireHitSim(new CSCWireHitSim(theDriftSim, p)),
       theStripHitSim(new CSCStripHitSim()),
       theWireElectronicsSim(new CSCWireElectronicsSim(p.getParameter<edm::ParameterSet>("wires"))),
       theStripElectronicsSim(new CSCStripElectronicsSim(p.getParameter<edm::ParameterSet>("strips"))),
       theNeutronReader(nullptr),
       theCSCGeometry(nullptr),
       theLayersNeeded(p.getParameter<unsigned int>("layersNeeded")),
       digitizeBadChambers_(p.getParameter<bool>("digitizeBadChambers")) {
   if (p.getParameter<bool>("doNeutrons")) {
     theNeutronReader = new CSCNeutronReader(p.getParameter<edm::ParameterSet>("neutrons"));
   }
 }

References AlCaHLTBitMon_ParallelJobs::p, and theNeutronReader.

Member Function Documentation

◆ doAction()

void CSCDigitizer::doAction	(	MixCollection< PSimHit > &	simHits,
		CSCWireDigiCollection &	wireDigis,
		CSCStripDigiCollection &	stripDigis,
		CSCComparatorDigiCollection &	comparators,
		DigiSimLinks &	wireDigiSimLinks,
		DigiSimLinks &	stripDigiSimLinks,
		CLHEP::HepRandomEngine *	engine
	)

digitize

Definition at line 39 of file CSCDigitizer.cc.

                                                           {
   // arrange the hits by layer
   std::map<int, edm::PSimHitContainer> hitMap;
   for (MixCollection<PSimHit>::MixItr hitItr = simHits.begin(); hitItr != simHits.end(); ++hitItr) {
     hitMap[hitItr->detUnitId()].push_back(*hitItr);
   }
  
   // count how many layers on each chamber are hit
   std::map<int, std::set<int>> layersInChamberHit;
   for (std::map<int, edm::PSimHitContainer>::const_iterator hitMapItr = hitMap.begin(); hitMapItr != hitMap.end();
        ++hitMapItr) {
     CSCDetId cscDetId(hitMapItr->first);
     int chamberId = cscDetId.chamberId();
     layersInChamberHit[chamberId].insert(cscDetId.layer());
   }
  
   // add neutron background, if needed
   if (theNeutronReader != nullptr) {
     theNeutronReader->addHits(hitMap, engine);
   }
  
   // now loop over layers and run the simulation for each one
   for (std::map<int, edm::PSimHitContainer>::const_iterator hitMapItr = hitMap.begin(); hitMapItr != hitMap.end();
        ++hitMapItr) {
     CSCDetId detId = CSCDetId(hitMapItr->first);
     int chamberId = detId.chamberId();
     int endc = detId.endcap();
     int stat = detId.station();
     int ring = detId.ring();
     int cham = detId.chamber();
  
     unsigned int nLayersInChamberHitForWireDigis = 0;
     if (stat == 1 && ring == 1) {  // ME1b
       std::set<int> layersInME1a = layersInChamberHit[CSCDetId(endc, stat, 4, cham, 0)];
       std::set<int> layersInME11 = layersInChamberHit[chamberId];
       layersInME11.insert(layersInME1a.begin(), layersInME1a.end());
       nLayersInChamberHitForWireDigis = layersInME11.size();
     } else if (stat == 1 && ring == 4) {  // ME1a
       std::set<int> layersInME1b = layersInChamberHit[CSCDetId(endc, stat, 1, cham, 0)];
       std::set<int> layersInME11 = layersInChamberHit[chamberId];
       layersInME11.insert(layersInME1b.begin(), layersInME1b.end());
       nLayersInChamberHitForWireDigis = layersInME11.size();
     } else
       nLayersInChamberHitForWireDigis = layersInChamberHit[chamberId].size();
  
     unsigned int nLayersInChamberHitForStripDigis = layersInChamberHit[chamberId].size();
  
     if (nLayersInChamberHitForWireDigis < theLayersNeeded && nLayersInChamberHitForStripDigis < theLayersNeeded)
       continue;
     // skip bad chambers
     if (!digitizeBadChambers_ && theConditions->isInBadChamber(CSCDetId(hitMapItr->first)))
       continue;
  
     const CSCLayer *layer = findLayer(hitMapItr->first);
     const edm::PSimHitContainer &layerSimHits = hitMapItr->second;
  
     std::vector<CSCDetectorHit> newWireHits, newStripHits;
  
     //    LogTrace("CSCDigitizer") << "CSCDigitizer: found " <<
     //    layerSimHits.size() <<" hit(s) in layer"
     edm::LogVerbatim("CSCDigitizer") << "[CSCDigitizer] found " << layerSimHits.size() << " hit(s) in layer"
                                      << layer->id();
     //       << " E" << layer->id().endcap() << " S" << layer->id().station() <<
     //       " R" << layer->id().ring()
     //       << " C" << layer->id().chamber() << " L" << layer->id().layer();
  
     // turn the edm::PSimHits into WireHits, using the WireHitSim
     { newWireHits.swap(theWireHitSim->simulate(layer, layerSimHits, engine)); }
     if (!newWireHits.empty()) {
       newStripHits.swap(theStripHitSim->simulate(layer, newWireHits));
     }
  
     // turn the hits into wire digis, using the electronicsSim
     if (nLayersInChamberHitForWireDigis >= theLayersNeeded) {
       theWireElectronicsSim->simulate(layer, newWireHits, engine);
       theWireElectronicsSim->fillDigis(wireDigis, engine);
       wireDigiSimLinks.insert(theWireElectronicsSim->digiSimLinks());
     }
     if (nLayersInChamberHitForStripDigis >= theLayersNeeded) {
       theStripElectronicsSim->simulate(layer, newStripHits, engine);
       theStripElectronicsSim->fillDigis(stripDigis, comparators, engine);
       stripDigiSimLinks.insert(theStripElectronicsSim->digiSimLinks());
     }
   }
  
   // fill in the layers were missing from this chamber
   std::list<int> missingLayers = layersMissing(stripDigis);
   for (std::list<int>::const_iterator missingLayerItr = missingLayers.begin(); missingLayerItr != missingLayers.end();
        ++missingLayerItr) {
     const CSCLayer *layer = findLayer(*missingLayerItr);
     theStripElectronicsSim->fillMissingLayer(layer, comparators, stripDigis, engine);
   }
 }

References CSCNeutronReader::addHits(), relativeConstraints::cham, CSCDetId::chamber(), CSCDetId::chamberId(), CSCBaseElectronicsSim::digiSimLinks(), digitizeBadChambers_, CSCDetId::endcap(), CSCWireElectronicsSim::fillDigis(), CSCStripElectronicsSim::fillDigis(), CSCStripElectronicsSim::fillMissingLayer(), findLayer(), edm::DetSetVector< T >::insert(), CSCStripConditions::isInBadChamber(), CSCDetId::layer(), phase1PixelTopology::layer, layersMissing(), relativeConstraints::ring, CSCDetId::ring(), FastTrackerRecHitCombiner_cfi::simHits, CSCStripHitSim::simulate(), CSCWireHitSim::simulate(), CSCBaseElectronicsSim::simulate(), edm_modernize_messagelogger::stat, CSCDetId::station(), theConditions, theLayersNeeded, theNeutronReader, theStripElectronicsSim, theStripHitSim, theWireElectronicsSim, and theWireHitSim.

Referenced by CSCDigiProducer::produce().

◆ findLayer()

const CSCLayer * CSCDigitizer::findLayer ( int detId ) const

private

finds the layer in the geometry associated with this det ID

Definition at line 185 of file CSCDigitizer.cc.

                                                        {
   assert(theCSCGeometry != nullptr);
   const GeomDetUnit *detUnit = theCSCGeometry->idToDetUnit(CSCDetId(detId));
   if (detUnit == nullptr) {
     throw cms::Exception("CSCDigiProducer") << "Invalid DetUnit: " << CSCDetId(detId)
                                             << "\nPerhaps your signal or pileup dataset are not compatible with "
                                                "the current release?";
   }
   return dynamic_cast<const CSCLayer *>(detUnit);
 }

References cms::cuda::assert(), Exception, CSCGeometry::idToDetUnit(), and theCSCGeometry.

Referenced by doAction().

◆ layersMissing()

std::list< int > CSCDigitizer::layersMissing ( const CSCStripDigiCollection & stripDigis ) const

private

finds which layers, 1-6, aren't in the current list

Definition at line 139 of file CSCDigitizer.cc.

                                                                                        {
   std::list<int> result;
  
   std::map<int, std::list<int>> layersInChamberWithDigi;
   for (CSCStripDigiCollection::DigiRangeIterator j = stripDigis.begin(); j != stripDigis.end(); j++) {
     CSCDetId layerId((*j).first);
     // make sure the vector of digis isn't empty
     if ((*j).second.first != (*j).second.second) {
       int chamberId = layerId.chamberId();
       layersInChamberWithDigi[chamberId].push_back(layerId.layer());
     }
   }
  
   std::list<int> oneThruSix;
   for (int i = 1; i <= 6; ++i)
     oneThruSix.push_back(i);
  
   for (std::map<int, std::list<int>>::iterator layersInChamberWithDigiItr = layersInChamberWithDigi.begin();
        layersInChamberWithDigiItr != layersInChamberWithDigi.end();
        ++layersInChamberWithDigiItr) {
     std::list<int> &layersHit = layersInChamberWithDigiItr->second;
     if (layersHit.size() < 6 && layersHit.size() >= theLayersNeeded) {
       layersHit.sort();
       std::list<int> missingLayers(6);
       std::list<int>::iterator lastLayerMissing = set_difference(
           oneThruSix.begin(), oneThruSix.end(), layersHit.begin(), layersHit.end(), missingLayers.begin());
       int chamberId = layersInChamberWithDigiItr->first;
       for (std::list<int>::iterator layerMissingItr = missingLayers.begin(); layerMissingItr != lastLayerMissing;
            ++layerMissingItr) {
         // got from layer 1-6 to layer ID
         result.push_back(chamberId + *layerMissingItr);
       }
     }
   }
   return result;
 }

References CSCDetId::chamberId(), mps_fire::i, dqmiolumiharvest::j, CSCDetId::layer(), genParticles_cff::map, mps_fire::result, and theLayersNeeded.

Referenced by doAction().

◆ operator=()

CSCDigitizer& CSCDigitizer::operator= ( const CSCDigitizer & )

delete

◆ setGeometry()

void CSCDigitizer::setGeometry ( const CSCGeometry * geom )

inline

sets geometry

Definition at line 62 of file CSCDigitizer.h.

62 { theCSCGeometry = geom; }

References relativeConstraints::geom, and theCSCGeometry.

Referenced by CSCDigiProducer::produce().

◆ setMagneticField()

void CSCDigitizer::setMagneticField ( const MagneticField * field )

sets the magnetic field

Definition at line 176 of file CSCDigitizer.cc.

176 { theDriftSim->setMagneticField(field); }

References CSCDriftSim::setMagneticField(), and theDriftSim.

Referenced by CSCDigiProducer::produce().

◆ setParticleDataTable()

void CSCDigitizer::setParticleDataTable ( const ParticleDataTable * pdt )

Definition at line 183 of file CSCDigitizer.cc.

183 { theWireHitSim->setParticleDataTable(pdt); }

References CSCWireHitSim::setParticleDataTable(), and theWireHitSim.

Referenced by CSCDigiProducer::produce().

◆ setStripConditions()

void CSCDigitizer::setStripConditions ( CSCStripConditions * cond )

Definition at line 178 of file CSCDigitizer.cc.

                                                               {
   theConditions = cond;                              // cache here
   theStripElectronicsSim->setStripConditions(cond);  // percolate downwards
 }

References CSCStripElectronicsSim::setStripConditions(), theConditions, and theStripElectronicsSim.

Referenced by CSCDigiProducer::CSCDigiProducer().

Member Data Documentation

◆ digitizeBadChambers_

bool CSCDigitizer::digitizeBadChambers_

private

Definition at line 87 of file CSCDigitizer.h.

Referenced by doAction().

◆ theConditions

CSCStripConditions* CSCDigitizer::theConditions

private

Definition at line 85 of file CSCDigitizer.h.

Referenced by doAction(), and setStripConditions().

◆ theCSCGeometry

const CSCGeometry* CSCDigitizer::theCSCGeometry

private

Definition at line 84 of file CSCDigitizer.h.

Referenced by findLayer(), and setGeometry().

◆ theDriftSim

CSCDriftSim* CSCDigitizer::theDriftSim

private

Definition at line 78 of file CSCDigitizer.h.

Referenced by setMagneticField(), and ~CSCDigitizer().

◆ theLayersNeeded

unsigned int CSCDigitizer::theLayersNeeded

private

Definition at line 86 of file CSCDigitizer.h.

Referenced by doAction(), and layersMissing().

◆ theNeutronReader

CSCNeutronReader* CSCDigitizer::theNeutronReader

private

Definition at line 83 of file CSCDigitizer.h.

Referenced by CSCDigitizer(), doAction(), and ~CSCDigitizer().

◆ theStripElectronicsSim

CSCStripElectronicsSim* CSCDigitizer::theStripElectronicsSim

private

Definition at line 82 of file CSCDigitizer.h.

Referenced by doAction(), setStripConditions(), and ~CSCDigitizer().

◆ theStripHitSim

CSCStripHitSim* CSCDigitizer::theStripHitSim

private

Definition at line 80 of file CSCDigitizer.h.

Referenced by doAction(), and ~CSCDigitizer().

◆ theWireElectronicsSim

CSCWireElectronicsSim* CSCDigitizer::theWireElectronicsSim

private

Definition at line 81 of file CSCDigitizer.h.

Referenced by doAction(), and ~CSCDigitizer().

◆ theWireHitSim

CSCWireHitSim* CSCDigitizer::theWireHitSim

private

Definition at line 79 of file CSCDigitizer.h.

Referenced by doAction(), setParticleDataTable(), and ~CSCDigitizer().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ DigiSimLinks

Constructor & Destructor Documentation

◆ ~CSCDigitizer()

◆ CSCDigitizer() [1/2]

◆ CSCDigitizer() [2/2]

Member Function Documentation

◆ doAction()

◆ findLayer()

◆ layersMissing()

◆ operator=()

◆ setGeometry()

◆ setMagneticField()

◆ setParticleDataTable()

◆ setStripConditions()

Member Data Documentation

◆ digitizeBadChambers_

◆ theConditions

◆ theCSCGeometry

◆ theDriftSim

◆ theLayersNeeded

◆ theNeutronReader

◆ theStripElectronicsSim

◆ theStripHitSim

◆ theWireElectronicsSim

◆ theWireHitSim