CMS 3D CMS Logo

List of all members | Public Member Functions
PSSDigitizerAlgorithm Class Reference

#include <PSSDigitizerAlgorithm.h>

Inheritance diagram for PSSDigitizerAlgorithm:
Phase2TrackerDigitizerAlgorithm

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, CLHEP::HepRandomEngine &)
 
 ~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 ()
 
 Phase2TrackerDigitizerAlgorithm (const edm::ParameterSet &conf_common, const edm::ParameterSet &conf_specific, CLHEP::HepRandomEngine &)
 
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::ParameterSetbadPixels
 
const float ClusterWidth
 
const Parameters DeadModules
 
const bool doDigitalReadout
 
const std::unique_ptr< CLHEP::RandFlat > flatDistribution_
 
const std::unique_ptr< SiG4UniversalFluctuationfluctuate
 
const bool fluctuateCharge
 
const std::unique_ptr< CLHEP::RandGaussQ > gaussDistribution_
 
edm::ESHandle< TrackerGeometrygeom_
 
const float GeVperElectron
 
const float interstripCoupling
 
const bool makeDigiSimLinks_
 
edm::ESHandle< SiPixelFedCablingMapmap_
 
bool pixelFlag
 
const double pseudoRadDamage
 
const double pseudoRadDamageRadius
 
CLHEP::HepRandomEngine * rengine_
 
const float Sigma0
 
const float SigmaCoeff
 
edm::ESHandle< SiPixelQualitySiPixelBadModule_
 
edm::ESHandle< SiPixelLorentzAngleSiPixelLorentzAngle_
 
const std::unique_ptr< CLHEP::RandGaussQ > smearedThreshold_Barrel_
 
const 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< GaussianTailNoiseGeneratortheNoiser
 
const float theReadoutNoise
 
const std::unique_ptr< SiPixelGainCalibrationOfflineSimServicetheSiPixelGainCalibrationService_
 
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,
CLHEP::HepRandomEngine &  eng 
)

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.

31  :
33  conf.getParameter<ParameterSet>("PSSDigitizerAlgorithm"),
34  eng)
35 {
36  pixelFlag = false;
37  LogInfo("PSSDigitizerAlgorithm") << "Algorithm constructed "
38  << "Configuration parameters:"
39  << "Threshold/Gain = "
40  << "threshold in electron Endcap = "
42  << "threshold in electron Barrel = "
44  << " " << theElectronPerADC << " " << theAdcFullScale
45  << " The delta cut-off is set to " << tMax
46  << " pix-inefficiency " << AddPixelInefficiency;
47 }
T getParameter(std::string const &) const
Phase2TrackerDigitizerAlgorithm(const edm::ParameterSet &conf_common, const edm::ParameterSet &conf_specific, CLHEP::HepRandomEngine &)
PSSDigitizerAlgorithm::~PSSDigitizerAlgorithm ( )
override

Definition at line 48 of file PSSDigitizerAlgorithm.cc.

References LogDebug.

48  {
49  LogDebug("PSSDigitizerAlgorithm") << "Algorithm deleted";
50 }
#define LogDebug(id)

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 51 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().

56  {
57  // produce SignalPoint's for all SimHit's in detector
58  // Loop over hits
59  uint32_t detId = pixdet->geographicalId().rawId();
60  size_t simHitGlobalIndex = inputBeginGlobalIndex; // This needs to be stored to create the digi-sim link later
61  for (auto it = inputBegin; it != inputEnd; ++it, ++simHitGlobalIndex) {
62  // skip hits not in this detector.
63  if ((*it).detUnitId() != detId)
64  continue;
65 
66  LogDebug("PSSDigitizerAlgorithm")
67  << (*it).particleType() << " " << (*it).pabs() << " "
68  << (*it).energyLoss() << " " << (*it).tof() << " "
69  << (*it).trackId() << " " << (*it).processType() << " "
70  << (*it).detUnitId()
71  << (*it).entryPoint() << " " << (*it).exitPoint();
72 
73  std::vector<DigitizerUtility::EnergyDepositUnit> ionization_points;
74  std::vector<DigitizerUtility::SignalPoint> collection_points;
75 
76  // fill collection_points for this SimHit, indpendent of topology
77  // Check the TOF cut
78  if (((*it).tof() - pixdet->surface().toGlobal((*it).localPosition()).mag()/30.) >= theTofLowerCut &&
79  ((*it).tof() - pixdet->surface().toGlobal((*it).localPosition()).mag()/30.) <= theTofUpperCut) {
80  primary_ionization(*it, ionization_points); // fills _ionization_points
81  drift(*it, pixdet, bfield, ionization_points, collection_points); // transforms _ionization_points to collection_points
82 
83  // compute induced signal on readout elements and add to _signal
84  induce_signal(*it, simHitGlobalIndex, tofBin, pixdet, collection_points); // *ihit needed only for SimHit<-->Digi link
85  }
86  }
87 }
void primary_ionization(const PSimHit &hit, std::vector< DigitizerUtility::EnergyDepositUnit > &ionization_points) const
#define LogDebug(id)
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:106
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
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
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:42
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
DetId geographicalId() const
The label of this GeomDet.
Definition: GeomDet.h:79
void induce_signal(const PSimHit &hit, const size_t hitIndex, const unsigned int tofBin, const Phase2TrackerGeomDetUnit *pixdet, const std::vector< DigitizerUtility::SignalPoint > &collection_points)
void PSSDigitizerAlgorithm::init ( const edm::EventSetup es)
overridevirtual

Implements Phase2TrackerDigitizerAlgorithm.

Definition at line 28 of file PSSDigitizerAlgorithm.cc.

References edm::EventSetup::get().

28  {
30 }
edm::ESHandle< TrackerGeometry > geom_
const T & get() const
Definition: EventSetup.h:55