PSSDigitizerAlgorithm Class Reference

#include <PSSDigitizerAlgorithm.h>

Inheritance diagram for PSSDigitizerAlgorithm:

Public Member Functions
void	accumulateSimHits (const std::vector< PSimHit >::const_iterator inputBegin, const std::vector< PSimHit >::const_iterator inputEnd, const size_t inputBeginGlobalIndex, const unsigned int tofBin, const Phase2TrackerGeomDetUnit *pixdet, const GlobalVector &bfield) override
void	init (const edm::EventSetup &es) override
	PSSDigitizerAlgorithm (const edm::ParameterSet &conf)
	~PSSDigitizerAlgorithm () override
Public Member Functions inherited from Phase2TrackerDigitizerAlgorithm
virtual void	digitize (const Phase2TrackerGeomDetUnit *pixdet, std::map< int, DigitizerUtility::DigiSimInfo > &digi_map, const TrackerTopology *tTopo)
virtual void	initializeEvent (CLHEP::HepRandomEngine &eng)
void	loadAccumulator (unsigned int detId, const std::map< int, float > &accumulator)
	Phase2TrackerDigitizerAlgorithm (const edm::ParameterSet &conf_common, const edm::ParameterSet &conf_specific)
virtual	~Phase2TrackerDigitizerAlgorithm ()

Additional Inherited Members
Protected Types inherited from Phase2TrackerDigitizerAlgorithm
using	Frame = GloballyPositioned< double >
using	Parameters = std::vector< edm::ParameterSet >
using	signal_map_const_iterator = signal_map_type::const_iterator
using	signal_map_iterator = signal_map_type::iterator
using	signal_map_type = std::map< int, DigitizerUtility::Amplitude, std::less< int > >
using	signalMaps = std::map< uint32_t, signal_map_type >
using	simlink_map = std::map< unsigned int, std::vector< float >, std::less< unsigned int > >
Protected Member Functions inherited from Phase2TrackerDigitizerAlgorithm
virtual void	add_noise (const Phase2TrackerGeomDetUnit *pixdet, float thePixelThreshold)
double	calcQ (float x) const
void	drift (const PSimHit &hit, const Phase2TrackerGeomDetUnit *pixdet, const GlobalVector &bfield, const std::vector< DigitizerUtility::EnergyDepositUnit > &ionization_points, std::vector< DigitizerUtility::SignalPoint > &collection_points) const
LocalVector	DriftDirection (const Phase2TrackerGeomDetUnit *pixdet, const GlobalVector &bfield, const DetId &detId) const
void	fluctuateEloss (int particleId, float momentum, float eloss, float length, int NumberOfSegments, std::vector< float > &elossVector) const
void	induce_signal (const PSimHit &hit, const size_t hitIndex, const unsigned int tofBin, const Phase2TrackerGeomDetUnit *pixdet, const std::vector< DigitizerUtility::SignalPoint > &collection_points)
virtual void	module_killing_conf (uint32_t detID)
virtual void	module_killing_DB (uint32_t detID)
virtual void	pixel_inefficiency (const SubdetEfficiencies &eff, const Phase2TrackerGeomDetUnit *pixdet, const TrackerTopology *tTopo)
virtual void	pixel_inefficiency_db (uint32_t detID)
void	primary_ionization (const PSimHit &hit, std::vector< DigitizerUtility::EnergyDepositUnit > &ionization_points) const
Protected Attributes inherited from Phase2TrackerDigitizerAlgorithm
signalMaps	_signal
const bool	addNoise
const bool	addNoisyPixels
const bool	AddPixelInefficiency
const bool	addThresholdSmearing
const bool	addXtalk
const bool	alpha2Order
std::vector< edm::ParameterSet >	badPixels
const float	ClusterWidth
const Parameters	DeadModules
const bool	doDigitalReadout
const std::unique_ptr< SiG4UniversalFluctuation >	fluctuate
const bool	fluctuateCharge
std::unique_ptr< CLHEP::RandGaussQ >	gaussDistribution_
edm::ESHandle< TrackerGeometry >	geom_
const float	GeVperElectron
const float	interstripCoupling
const bool	makeDigiSimLinks_
edm::ESHandle< SiPixelFedCablingMap >	map_
bool	pixelFlag
const double	pseudoRadDamage
const double	pseudoRadDamageRadius
CLHEP::HepRandomEngine *	rengine_
const float	Sigma0
const float	SigmaCoeff
edm::ESHandle< SiPixelQuality >	SiPixelBadModule_
edm::ESHandle< SiPixelLorentzAngle >	SiPixelLorentzAngle_
std::unique_ptr< CLHEP::RandGaussQ >	smearedThreshold_Barrel_
std::unique_ptr< CLHEP::RandGaussQ >	smearedThreshold_Endcap_
const SubdetEfficiencies	subdetEfficiencies_
const float	tanLorentzAnglePerTesla_Barrel
const float	tanLorentzAnglePerTesla_Endcap
const int	theAdcFullScale
const float	theElectronPerADC
const double	theHIPThresholdInE_Barrel
const double	theHIPThresholdInE_Endcap
const float	theNoiseInElectrons
const std::unique_ptr< GaussianTailNoiseGenerator >	theNoiser
const float	theReadoutNoise
const std::unique_ptr< SiPixelGainCalibrationOfflineSimService >	theSiPixelGainCalibrationService_
const float	theThresholdInE_Barrel
const float	theThresholdInE_Endcap
const double	theThresholdSmearing_Barrel
const double	theThresholdSmearing_Endcap
const float	theTofLowerCut
const float	theTofUpperCut
const double	tMax
const bool	use_deadmodule_DB_
const bool	use_ineff_from_db_
const bool	use_LorentzAngle_DB_
const bool	use_module_killing_

Detailed Description

Definition at line 9 of file PSSDigitizerAlgorithm.h.

Constructor & Destructor Documentation

PSSDigitizerAlgorithm::PSSDigitizerAlgorithm

(

const edm::ParameterSet &

conf

)

Definition at line 31 of file PSSDigitizerAlgorithm.cc.

References Phase2TrackerDigitizerAlgorithm::AddPixelInefficiency, Phase2TrackerDigitizerAlgorithm::pixelFlag, Phase2TrackerDigitizerAlgorithm::theAdcFullScale, Phase2TrackerDigitizerAlgorithm::theElectronPerADC, Phase2TrackerDigitizerAlgorithm::theThresholdInE_Barrel, Phase2TrackerDigitizerAlgorithm::theThresholdInE_Endcap, and Phase2TrackerDigitizerAlgorithm::tMax.

                                                                       :
  Phase2TrackerDigitizerAlgorithm(conf.getParameter<ParameterSet>("AlgorithmCommon"),
                  conf.getParameter<ParameterSet>("PSSDigitizerAlgorithm"))
{
  pixelFlag = false;
  LogInfo("PSSDigitizerAlgorithm") << "Algorithm constructed "
                   << "Configuration parameters:"
                   << "Threshold/Gain = "
                   << "threshold in electron Endcap = "
                   << theThresholdInE_Endcap
                   << "threshold in electron Barrel = "
                   << theThresholdInE_Barrel
                   << " " << theElectronPerADC << " " << theAdcFullScale
                   << " The delta cut-off is set to " << tMax
                   << " pix-inefficiency " << AddPixelInefficiency;
}

PSSDigitizerAlgorithm::~PSSDigitizerAlgorithm ( )

override

Definition at line 47 of file PSSDigitizerAlgorithm.cc.

References LogDebug.

                                              {
  LogDebug("PSSDigitizerAlgorithm") << "Algorithm deleted";
}

Member Function Documentation

void PSSDigitizerAlgorithm::accumulateSimHits ( const std::vector< PSimHit >::const_iterator  inputBegin, const std::vector< PSimHit >::const_iterator  inputEnd, const size_t  inputBeginGlobalIndex, const unsigned int  tofBin, const Phase2TrackerGeomDetUnit *  pixdet, const GlobalVector &  bfield  )

overridevirtual

Implements Phase2TrackerDigitizerAlgorithm.

Definition at line 50 of file PSSDigitizerAlgorithm.cc.

References Phase2TrackerDigitizerAlgorithm::drift(), GeomDet::geographicalId(), Phase2TrackerDigitizerAlgorithm::induce_signal(), LogDebug, mag(), Phase2TrackerDigitizerAlgorithm::primary_ionization(), DetId::rawId(), GeomDet::surface(), Phase2TrackerDigitizerAlgorithm::theTofLowerCut, Phase2TrackerDigitizerAlgorithm::theTofUpperCut, and Surface::toGlobal().

                                                      {
  // produce SignalPoint's for all SimHit's in detector
  // Loop over hits
  uint32_t detId = pixdet->geographicalId().rawId();
  size_t simHitGlobalIndex = inputBeginGlobalIndex; // This needs to be stored to create the digi-sim link later
  for (auto it = inputBegin; it != inputEnd; ++it, ++simHitGlobalIndex) {
    // skip hits not in this detector.
    if ((*it).detUnitId() != detId)
      continue;

    LogDebug("PSSDigitizerAlgorithm")
      << (*it).particleType() << " " << (*it).pabs() << " "
      << (*it).energyLoss() << " " << (*it).tof() << " "
      << (*it).trackId() << " " << (*it).processType() << " "
      << (*it).detUnitId()
      << (*it).entryPoint() << " " << (*it).exitPoint();
    
    std::vector<DigitizerUtility::EnergyDepositUnit> ionization_points;
    std::vector<DigitizerUtility::SignalPoint> collection_points;

    // fill collection_points for this SimHit, indpendent of topology
    // Check the TOF cut
    if (((*it).tof() - pixdet->surface().toGlobal((*it).localPosition()).mag()/30.) >= theTofLowerCut &&
    ((*it).tof() - pixdet->surface().toGlobal((*it).localPosition()).mag()/30.) <= theTofUpperCut) {
      primary_ionization(*it, ionization_points); // fills _ionization_points
      drift(*it, pixdet, bfield, ionization_points, collection_points);  // transforms _ionization_points to collection_points

      // compute induced signal on readout elements and add to _signal
      induce_signal(*it, simHitGlobalIndex, tofBin, pixdet, collection_points); // *ihit needed only for SimHit<-->Digi link
    }
  }
}

void PSSDigitizerAlgorithm::init ( const edm::EventSetup &  es )

overridevirtual

Implements Phase2TrackerDigitizerAlgorithm.

Definition at line 28 of file PSSDigitizerAlgorithm.cc.

References edm::EventSetup::get().

                                                        {
  es.get<TrackerDigiGeometryRecord>().get(geom_);
}