CSCWireHitSim Class Reference

#include <CSCWireHitSim.h>

Public Member Functions
	CSCWireHitSim (CSCDriftSim *driftSim, const edm::ParameterSet &p)
void	setParticleDataTable (const ParticleDataTable *pdt)
std::vector< CSCDetectorHit > &	simulate (const CSCLayer *layer, const edm::PSimHitContainer &simHits, CLHEP::HepRandomEngine *)
	~CSCWireHitSim ()

Private Member Functions
CSCDetectorHit	driftElectronsToWire ()
std::vector< Local3DPoint >	getIonizationClusters (const PSimHit &hit, const CSCLayer *, CLHEP::HepRandomEngine *)

Private Attributes
CSCDriftSim *	theDriftSim
CSCGasCollisions *	theGasIonizer
std::vector< CSCDetectorHit >	theNewWireHits

Detailed Description

Class used to simulate hit on wire in Endcap Muon CSC.

Author

Rick Wilkinson

Tim Cox

Definition at line 27 of file CSCWireHitSim.h.

Constructor & Destructor Documentation

CSCWireHitSim()

CSCWireHitSim::CSCWireHitSim ( CSCDriftSim *  driftSim, const edm::ParameterSet &  p  )

explicit

Definition at line 10 of file CSCWireHitSim.cc.

    : theDriftSim(driftSim), theGasIonizer(new CSCGasCollisions(p)), theNewWireHits() {}

~CSCWireHitSim()

CSCWireHitSim::~CSCWireHitSim

(

)

Definition at line 13 of file CSCWireHitSim.cc.

{ delete theGasIonizer; }

References theGasIonizer.

Member Function Documentation

driftElectronsToWire()

CSCDetectorHit CSCWireHitSim::driftElectronsToWire ( )

private

getIonizationClusters()

std::vector< LocalPoint > CSCWireHitSim::getIonizationClusters ( const PSimHit &  hit, const CSCLayer *  layer, CLHEP::HepRandomEngine *  engine  )

private

Definition at line 40 of file CSCWireHitSim.cc.

                                                                                           {
  const LocalPoint &entryPoint = simHit.entryPoint();
  const LocalPoint &exitPoint = simHit.exitPoint();
 
  LogTrace("CSCWireHitSim") << "CSCWireHitSim:"
                            << " type=" << simHit.particleType() << " mom=" << simHit.pabs() << "\n Local entry "
                            << entryPoint << " exit " << exitPoint;
 
  std::vector<LocalPoint> positions;
  std::vector<int> electrons;
  theGasIonizer->simulate(simHit, positions, electrons, engine);
 
  std::vector<LocalPoint> results;  // start empty
 
  int j = 0;
  for (std::vector<LocalPoint>::const_iterator pointItr = positions.begin(); pointItr != positions.end(); ++pointItr) {
    ++j;
    // some verification
    if (layer->geometry()->inside(*pointItr)) {
      // push the point for each electron at this point
 
      for (int ie = 1; ie <= electrons[j - 1]; ++ie) {
        // probability of getting attached
        float f_att = 0.5;
        if (CLHEP::RandFlat::shoot(engine) > f_att) {
          results.push_back(*pointItr);
        }
      }
    }
  }
  LogTrace("CSCWireHitSim") << "CSCWireHitSim: there are " << results.size()
                            << " clusters identified with each electron.";
  return results;
}

References pwdgSkimBPark_cfi::electrons, ecal::raw::entryPoint(), dqmiolumiharvest::j, phase1PixelTopology::layer, LogTrace, bookConverter::results, rpcPointValidation_cfi::simHit, CSCGasCollisions::simulate(), and theGasIonizer.

Referenced by simulate().

setParticleDataTable()

void CSCWireHitSim::setParticleDataTable

(

const ParticleDataTable *

pdt

)

Definition at line 77 of file CSCWireHitSim.cc.

{ theGasIonizer->setParticleDataTable(pdt); }

References CSCGasCollisions::setParticleDataTable(), and theGasIonizer.

Referenced by CSCDigitizer::setParticleDataTable().

simulate()

std::vector< CSCDetectorHit > & CSCWireHitSim::simulate	(	const CSCLayer *	layer,
		const edm::PSimHitContainer &	simHits,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 15 of file CSCWireHitSim.cc.

                                                                                   {
  const CSCLayerGeometry *geom = layer->geometry();
 
  theNewWireHits.clear();
  for (edm::PSimHitContainer::const_iterator hitItr = simHits.begin(); hitItr != simHits.end(); ++hitItr) {
    std::vector<LocalPoint> ionClusters = getIonizationClusters(*hitItr, layer, engine);
 
    unsigned nClusters = ionClusters.size();
    theNewWireHits.reserve(theNewWireHits.size() + nClusters);
 
    for (unsigned icl = 0; icl < nClusters; ++icl) {
      // Drift the electrons in the cluster to the nearest wire...
      int nearestWire = geom->nearestWire(ionClusters[icl]);
 
      // The wire hit contains wire# and position measured _along the wire_
      // from where it intersects local y axis.
 
      theNewWireHits.push_back(theDriftSim->getWireHit(ionClusters[icl], layer, nearestWire, *hitItr, engine));
    }
  }
  return theNewWireHits;
}

References relativeConstraints::geom, getIonizationClusters(), CSCDriftSim::getWireHit(), phase1PixelTopology::layer, heavyionUCCDQM_cfi::nClusters, FastTrackerRecHitCombiner_cfi::simHits, theDriftSim, and theNewWireHits.

Referenced by CSCDigitizer::doAction().

Member Data Documentation

theDriftSim

CSCDriftSim* CSCWireHitSim::theDriftSim

private

Definition at line 45 of file CSCWireHitSim.h.

Referenced by simulate().

theGasIonizer

CSCGasCollisions* CSCWireHitSim::theGasIonizer

private

Definition at line 46 of file CSCWireHitSim.h.

Referenced by getIonizationClusters(), setParticleDataTable(), and ~CSCWireHitSim().

theNewWireHits

std::vector<CSCDetectorHit> CSCWireHitSim::theNewWireHits

private

Definition at line 47 of file CSCWireHitSim.h.

Referenced by simulate().