CSCDigitizer Class Reference

#include <CSCDigitizer.h>

Inheritance diagram for CSCDigitizer:

Public Types
typedef edm::DetSetVector < StripDigiSimLink >	DigiSimLinks

Public Member Functions
	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 *)
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 41 of file CSCDigitizer.h.

Member Typedef Documentation

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

Definition at line 44 of file CSCDigitizer.h.

Constructor & Destructor Documentation

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

explicit

configurable parameters

Definition at line 16 of file CSCDigitizer.cc.

References edm::ParameterSet::getParameter(), and theNeutronReader.

: 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(0),
  theCSCGeometry(0),
  theLayersNeeded(p.getParameter<unsigned int>("layersNeeded")),
  digitizeBadChambers_(p.getParameter<bool>("digitizeBadChambers"))
{
  if(p.getParameter<bool>("doNeutrons"))
  {
    theNeutronReader = new CSCNeutronReader(p.getParameter<edm::ParameterSet>("neutrons"));
  }
}

CSCDigitizer::~CSCDigitizer

(

)

Definition at line 34 of file CSCDigitizer.cc.

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

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

Member Function Documentation

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

digitize

Definition at line 45 of file CSCDigitizer.cc.

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

Referenced by CSCDigiProducer::produce().

{
  // 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 != 0)
  {
    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);
  }
}

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

private

finds the layer in the geometry associated with this det ID

Definition at line 215 of file CSCDigitizer.cc.

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

Referenced by doAction().

                                                        {
  assert(theCSCGeometry != 0);
  const GeomDetUnit* detUnit = theCSCGeometry->idToDetUnit(CSCDetId(detId));
  if(detUnit == 0)
  {
    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);
}

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

private

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

Definition at line 152 of file CSCDigitizer.cc.

References CSCDetId::chamberId(), i, j, CSCDetId::layer(), python.multivaluedict::map(), query::result, and theLayersNeeded.

Referenced by doAction().

{
  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;
}

void CSCDigitizer::setGeometry ( const CSCGeometry *  geom )

inline

sets geometry

Definition at line 62 of file CSCDigitizer.h.

References relativeConstraints::geom, and theCSCGeometry.

Referenced by CSCDigiProducer::produce().

{theCSCGeometry = geom;}

void CSCDigitizer::setMagneticField

(

const MagneticField *

field

)

sets the magnetic field

Definition at line 197 of file CSCDigitizer.cc.

References CSCDriftSim::setMagneticField(), and theDriftSim.

Referenced by CSCDigiProducer::produce().

                                                               {
  theDriftSim->setMagneticField(field);
}

void CSCDigitizer::setParticleDataTable

(

const ParticleDataTable *

pdt

)

Definition at line 209 of file CSCDigitizer.cc.

References CSCWireHitSim::setParticleDataTable(), and theWireHitSim.

Referenced by CSCDigiProducer::produce().

{
  theWireHitSim->setParticleDataTable(pdt);
}

void CSCDigitizer::setStripConditions

(

CSCStripConditions *

cond

)

Definition at line 202 of file CSCDigitizer.cc.

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

Referenced by CSCDigiProducer::CSCDigiProducer().

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

Member Data Documentation

bool CSCDigitizer::digitizeBadChambers_

private

Definition at line 87 of file CSCDigitizer.h.

Referenced by doAction().

CSCStripConditions* CSCDigitizer::theConditions

private

Definition at line 85 of file CSCDigitizer.h.

Referenced by doAction(), and setStripConditions().

const CSCGeometry* CSCDigitizer::theCSCGeometry

private

Definition at line 84 of file CSCDigitizer.h.

Referenced by findLayer(), and setGeometry().

CSCDriftSim* CSCDigitizer::theDriftSim

private

Definition at line 78 of file CSCDigitizer.h.

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

unsigned int CSCDigitizer::theLayersNeeded

private

Definition at line 86 of file CSCDigitizer.h.

Referenced by doAction(), and layersMissing().

CSCNeutronReader* CSCDigitizer::theNeutronReader

private

Definition at line 83 of file CSCDigitizer.h.

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

CSCStripElectronicsSim* CSCDigitizer::theStripElectronicsSim

private

Definition at line 82 of file CSCDigitizer.h.

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

CSCStripHitSim* CSCDigitizer::theStripHitSim

private

Definition at line 80 of file CSCDigitizer.h.

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

CSCWireElectronicsSim* CSCDigitizer::theWireElectronicsSim

private

Definition at line 81 of file CSCDigitizer.h.

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

CSCWireHitSim* CSCDigitizer::theWireHitSim

private

Definition at line 79 of file CSCDigitizer.h.

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