CMS 3D CMS Logo

List of all members | Public Member Functions | Public Attributes
PixelDigitizerAlgorithm Class Reference

#include <PixelDigitizerAlgorithm.h>

Inheritance diagram for PixelDigitizerAlgorithm:
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 add_cross_talk (const Phase2TrackerGeomDetUnit *pixdet) override
 
void init (const edm::EventSetup &es) override
 
 PixelDigitizerAlgorithm (const edm::ParameterSet &conf)
 
 ~PixelDigitizerAlgorithm () 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 ()
 

Public Attributes

const double even_column_interchannelCoupling_next_column
 
const double even_row_interchannelCoupling_next_row
 
const double odd_column_interchannelCoupling_next_column
 
const double odd_row_interchannelCoupling_next_row
 

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)
 
virtual void add_noisy_cells (const Phase2TrackerGeomDetUnit *pixdet, float thePixelThreshold)
 
double calcQ (float x) const
 
int convertSignalToAdc (uint32_t detID, float signal_in_elec, float threshold)
 
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 std::unique_ptr< SiG4UniversalFluctuationfluctuate
 
const bool fluctuateCharge
 
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_
 
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< GaussianTailNoiseGeneratortheNoiser
 
const int thePhase2ReadoutMode
 
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 PixelDigitizerAlgorithm.h.

Constructor & Destructor Documentation

PixelDigitizerAlgorithm::PixelDigitizerAlgorithm ( const edm::ParameterSet conf)

Definition at line 64 of file PixelDigitizerAlgorithm.cc.

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

64  :
66  conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm")),
67  odd_row_interchannelCoupling_next_row(conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm").getParameter<double>("Odd_row_interchannelCoupling_next_row")),
68  even_row_interchannelCoupling_next_row(conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm").getParameter<double>("Even_row_interchannelCoupling_next_row")),
69  odd_column_interchannelCoupling_next_column(conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm").getParameter<double>("Odd_column_interchannelCoupling_next_column")),
70  even_column_interchannelCoupling_next_column(conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm").getParameter<double>("Even_column_interchannelCoupling_next_column"))
71 {
72  pixelFlag = true;
73  LogInfo("PixelDigitizerAlgorithm") << "Algorithm constructed "
74  << "Configuration parameters:"
75  << "Threshold/Gain = "
76  << "threshold in electron Endcap = "
78  << "threshold in electron Barrel = "
80  << " " << theElectronPerADC << " " << theAdcFullScale
81  << " The delta cut-off is set to " << tMax
82  << " pix-inefficiency " << AddPixelInefficiency;
83 }
const double even_row_interchannelCoupling_next_row
T getParameter(std::string const &) const
const double odd_column_interchannelCoupling_next_column
const double even_column_interchannelCoupling_next_column
Phase2TrackerDigitizerAlgorithm(const edm::ParameterSet &conf_common, const edm::ParameterSet &conf_specific)
const double odd_row_interchannelCoupling_next_row
PixelDigitizerAlgorithm::~PixelDigitizerAlgorithm ( )
override

Definition at line 84 of file PixelDigitizerAlgorithm.cc.

References LogDebug.

84  {
85  LogDebug("PixelDigitizerAlgorithm") << "Algorithm deleted";
86 }
#define LogDebug(id)

Member Function Documentation

void PixelDigitizerAlgorithm::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 87 of file PixelDigitizerAlgorithm.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().

92  {
93  // produce SignalPoint's for all SimHit's in detector
94  // Loop over hits
95  uint32_t detId = pixdet->geographicalId().rawId();
96  size_t simHitGlobalIndex = inputBeginGlobalIndex; // This needs to be stored to create the digi-sim link later
97  for (auto it = inputBegin; it != inputEnd; ++it, ++simHitGlobalIndex) {
98  // skip hits not in this detector.
99  if ((*it).detUnitId() != detId)
100  continue;
101 
102  LogDebug ("PixelDigitizerAlgorithm")
103  << (*it).particleType() << " " << (*it).pabs() << " "
104  << (*it).energyLoss() << " " << (*it).tof() << " "
105  << (*it).trackId() << " " << (*it).processType() << " "
106  << (*it).detUnitId()
107  << (*it).entryPoint() << " " << (*it).exitPoint();
108 
109  std::vector<DigitizerUtility::EnergyDepositUnit> ionization_points;
110  std::vector<DigitizerUtility::SignalPoint> collection_points;
111 
112  // fill collection_points for this SimHit, indpendent of topology
113  // Check the TOF cut
114  if (((*it).tof() - pixdet->surface().toGlobal((*it).localPosition()).mag()/30.) >= theTofLowerCut &&
115  ((*it).tof() - pixdet->surface().toGlobal((*it).localPosition()).mag()/30.) <= theTofUpperCut) {
116  primary_ionization(*it, ionization_points); // fills _ionization_points
117  drift (*it, pixdet, bfield, ionization_points, collection_points); // transforms _ionization_points to collection_points
118 
119  // compute induced signal on readout elements and add to _signal
120  induce_signal(*it, simHitGlobalIndex, tofBin, pixdet, collection_points); // *ihit needed only for SimHit<-->Digi link
121  }
122  }
123 }
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
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:50
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:42
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 PixelDigitizerAlgorithm::add_cross_talk ( const Phase2TrackerGeomDetUnit pixdet)
overridevirtual

Reimplemented from Phase2TrackerDigitizerAlgorithm.

Definition at line 129 of file PixelDigitizerAlgorithm.cc.

References Phase2TrackerDigitizerAlgorithm::_signal, officialStyle::chan, PixelDigi::channelToPixel(), even_column_interchannelCoupling_next_column, even_row_interchannelCoupling_next_row, GeomDet::geographicalId(), checklumidiff::l, PixelTopology::ncolumns(), PixelTopology::nrows(), odd_column_interchannelCoupling_next_column, odd_row_interchannelCoupling_next_row, PixelTopology::pitch(), Phase2TrackerDigitizerAlgorithm::pixelFlag, Phase2TrackerDigi::pixelToChannel(), PixelDigi::pixelToChannel(), DetId::rawId(), alignCSCRings::s, and PixelGeomDetUnit::specificTopology().

129  {
130  if (!pixelFlag) return;
131 
132  const Phase2TrackerTopology* topol = &pixdet->specificTopology();
133 
134  // cross-talk calculation valid for the case of 25x100 pixels
135  const float pitch_first = 0.0025;
136  const float pitch_second = 0.0100;
137 
138  if ( topol->pitch().first != pitch_first || topol->pitch().second != pitch_second ) return; // please check that units are really cm!
139 
140  uint32_t detID = pixdet->geographicalId().rawId();
141  signal_map_type& theSignal = _signal[detID];
142  signal_map_type signalNew;
143  int numRows = topol->nrows();
144  int numColumns = topol->ncolumns();
145 
146  for (auto & s : theSignal) {
147  float signalInElectrons = s.second.ampl(); // signal in electrons
148 
149  auto hitChan = PixelDigi::channelToPixel(s.first);
150 
151  float signalInElectrons_odd_row_Xtalk_next_row = signalInElectrons * odd_row_interchannelCoupling_next_row;
152  float signalInElectrons_even_row_Xtalk_next_row = signalInElectrons * even_row_interchannelCoupling_next_row;
153  float signalInElectrons_odd_column_Xtalk_next_column = signalInElectrons * odd_column_interchannelCoupling_next_column;
154  float signalInElectrons_even_column_Xtalk_next_column = signalInElectrons * even_column_interchannelCoupling_next_column;
155 
156  //subtract the charge which will be shared
157  s.second.set(signalInElectrons-signalInElectrons_odd_row_Xtalk_next_row-signalInElectrons_even_row_Xtalk_next_row-signalInElectrons_odd_column_Xtalk_next_column-signalInElectrons_even_column_Xtalk_next_column);
158 
159  if (hitChan.first != 0) {
160  auto XtalkPrev = std::make_pair(hitChan.first-1, hitChan.second);
161  int chanXtalkPrev = (pixelFlag) ? PixelDigi::pixelToChannel(XtalkPrev.first, XtalkPrev.second)
162  : Phase2TrackerDigi::pixelToChannel(XtalkPrev.first, XtalkPrev.second);
163  if (hitChan.first % 2 == 1) signalNew.emplace(chanXtalkPrev, DigitizerUtility::Amplitude(signalInElectrons_even_row_Xtalk_next_row, nullptr, -1.0));
164  else signalNew.emplace(chanXtalkPrev, DigitizerUtility::Amplitude(signalInElectrons_odd_row_Xtalk_next_row, nullptr, -1.0));
165  }
166  if (hitChan.first < (numRows-1)) {
167  auto XtalkNext = std::make_pair(hitChan.first+1, hitChan.second);
168  int chanXtalkNext = (pixelFlag) ? PixelDigi::pixelToChannel(XtalkNext.first, XtalkNext.second)
169  : Phase2TrackerDigi::pixelToChannel(XtalkNext.first, XtalkNext.second);
170  if (hitChan.first % 2 == 1) signalNew.emplace(chanXtalkNext, DigitizerUtility::Amplitude(signalInElectrons_odd_row_Xtalk_next_row, nullptr, -1.0));
171  else signalNew.emplace(chanXtalkNext, DigitizerUtility::Amplitude(signalInElectrons_even_row_Xtalk_next_row, nullptr, -1.0));
172  }
173 
174  if (hitChan.second != 0) {
175  auto XtalkPrev = std::make_pair(hitChan.first, hitChan.second-1);
176  int chanXtalkPrev = (pixelFlag) ? PixelDigi::pixelToChannel(XtalkPrev.first, XtalkPrev.second)
177  : Phase2TrackerDigi::pixelToChannel(XtalkPrev.first, XtalkPrev.second);
178  if (hitChan.second % 2 == 1) signalNew.emplace(chanXtalkPrev, DigitizerUtility::Amplitude(signalInElectrons_even_column_Xtalk_next_column, nullptr, -1.0));
179  else signalNew.emplace(chanXtalkPrev, DigitizerUtility::Amplitude(signalInElectrons_odd_column_Xtalk_next_column, nullptr, -1.0));
180  }
181  if (hitChan.second < (numColumns-1)) {
182  auto XtalkNext = std::make_pair(hitChan.first, hitChan.second+1);
183  int chanXtalkNext = (pixelFlag) ? PixelDigi::pixelToChannel(XtalkNext.first, XtalkNext.second)
184  : Phase2TrackerDigi::pixelToChannel(XtalkNext.first, XtalkNext.second);
185  if (hitChan.second % 2 == 1) signalNew.emplace(chanXtalkNext, DigitizerUtility::Amplitude(signalInElectrons_odd_column_Xtalk_next_column, nullptr, -1.0));
186  else signalNew.emplace(chanXtalkNext, DigitizerUtility::Amplitude(signalInElectrons_even_column_Xtalk_next_column, nullptr, -1.0));
187  }
188  }
189  for (auto const & l : signalNew) {
190  int chan = l.first;
191  auto iter = theSignal.find(chan);
192  if (iter != theSignal.end()) {
193  theSignal[chan] += l.second.ampl();
194  } else {
195  theSignal.emplace(chan, DigitizerUtility::Amplitude(l.second.ampl(), nullptr, -1.0));
196  }
197  }
198 }
const double even_row_interchannelCoupling_next_row
virtual int nrows() const =0
const double odd_column_interchannelCoupling_next_column
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:50
virtual std::pair< float, float > pitch() const =0
const double even_column_interchannelCoupling_next_column
static int pixelToChannel(int row, int col)
Definition: PixelDigi.h:68
static PackedDigiType pixelToChannel(unsigned int row, unsigned int col)
DetId geographicalId() const
The label of this GeomDet.
Definition: GeomDet.h:79
std::map< int, DigitizerUtility::Amplitude, std::less< int > > signal_map_type
const double odd_row_interchannelCoupling_next_row
static std::pair< int, int > channelToPixel(int ch)
Definition: PixelDigi.h:62
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
chan
lumi = TPaveText(lowX+0.38, lowY+0.061, lowX+0.45, lowY+0.161, "NDC") lumi.SetBorderSize( 0 ) lumi...
virtual int ncolumns() const =0
void PixelDigitizerAlgorithm::init ( const edm::EventSetup es)
overridevirtual

Implements Phase2TrackerDigitizerAlgorithm.

Definition at line 49 of file PixelDigitizerAlgorithm.cc.

References edm::EventSetup::get().

49  {
50  if (use_ineff_from_db_) // load gain calibration service fromdb...
51  theSiPixelGainCalibrationService_->setESObjects(es);
52 
55 
56  if (use_LorentzAngle_DB_) // Get Lorentz angle from DB record
58 
59  // gets the map and geometry from the DB (to kill ROCs)
60  es.get<SiPixelFedCablingMapRcd>().get(map_);
62 }
edm::ESHandle< SiPixelFedCablingMap > map_
const std::unique_ptr< SiPixelGainCalibrationOfflineSimService > theSiPixelGainCalibrationService_
edm::ESHandle< SiPixelLorentzAngle > SiPixelLorentzAngle_
edm::ESHandle< TrackerGeometry > geom_
T get() const
Definition: EventSetup.h:71
edm::ESHandle< SiPixelQuality > SiPixelBadModule_

Member Data Documentation

const double PixelDigitizerAlgorithm::even_column_interchannelCoupling_next_column

Definition at line 30 of file PixelDigitizerAlgorithm.h.

Referenced by add_cross_talk().

const double PixelDigitizerAlgorithm::even_row_interchannelCoupling_next_row

Definition at line 28 of file PixelDigitizerAlgorithm.h.

Referenced by add_cross_talk().

const double PixelDigitizerAlgorithm::odd_column_interchannelCoupling_next_column

Definition at line 29 of file PixelDigitizerAlgorithm.h.

Referenced by add_cross_talk().

const double PixelDigitizerAlgorithm::odd_row_interchannelCoupling_next_row

Definition at line 27 of file PixelDigitizerAlgorithm.h.

Referenced by add_cross_talk().