CMS 3D CMS Logo
List of all members | Classes | Public Member Functions | Public Attributes
PixelDigitizerAlgorithm Class Reference

#include <PixelDigitizerAlgorithm.h>

Inheritance diagram for PixelDigitizerAlgorithm:
Phase2TrackerDigitizerAlgorithm Pixel3DDigitizerAlgorithm PixelBrickedDigitizerAlgorithm

Classes

class  TimewalkCurve
 
class  TimewalkModel
 

Public Member Functions

void add_cross_talk (const Phase2TrackerGeomDetUnit *pixdet) override
 
void init (const edm::EventSetup &es) override
 
bool isAboveThreshold (const DigitizerUtility::SimHitInfo *hitInfo, float charge, float thr) const override
 
void module_killing_DB (const Phase2TrackerGeomDetUnit *pixdet) override
 
 PixelDigitizerAlgorithm (const edm::ParameterSet &conf, edm::ConsumesCollector iC)
 
bool select_hit (const PSimHit &hit, double tCorr, double &sigScale) const override
 
 ~PixelDigitizerAlgorithm () override
 
- Public Member Functions inherited from Phase2TrackerDigitizerAlgorithm
virtual void accumulateSimHits (const std::vector< PSimHit >::const_iterator inputBegin, const std::vector< PSimHit >::const_iterator inputEnd, const size_t inputBeginGlobalIndex, const uint32_t tofBin, const Phase2TrackerGeomDetUnit *pixdet, const GlobalVector &bfield)
 
virtual void digitize (const Phase2TrackerGeomDetUnit *pixdet, std::map< int, DigitizerUtility::DigiSimInfo > &digi_map, const TrackerTopology *tTopo)
 
virtual void initializeEvent (CLHEP::HepRandomEngine &eng)
 
void loadAccumulator (uint32_t detId, const std::map< int, float > &accumulator)
 
 Phase2TrackerDigitizerAlgorithm (const edm::ParameterSet &conf_common, const edm::ParameterSet &conf_specific, edm::ConsumesCollector iC)
 
virtual ~Phase2TrackerDigitizerAlgorithm ()
 

Public Attributes

bool apply_timewalk_
 
double even_column_interchannelCoupling_next_column_
 
double even_row_interchannelCoupling_next_row_
 
const edm::ESGetToken< SiPixelFedCablingMap, SiPixelFedCablingMapRcd > fedCablingMapToken_
 
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
 
double odd_column_interchannelCoupling_next_column_
 
double odd_row_interchannelCoupling_next_row_
 
edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd > siPixelBadModuleToken_
 
edm::ESGetToken< SiPixelLorentzAngle, SiPixelLorentzAngleSimRcd > siPixelLorentzAngleToken_
 
const TimewalkModel timewalk_model_
 

Additional Inherited Members

- Protected Types inherited from Phase2TrackerDigitizerAlgorithm
using Frame = GloballyPositioned< double >
 
using Parameters = std::vector< edm::ParameterSet >
 
using signal_map_type = std::map< int, DigitizerUtility::Amplitude, std::less< int > >
 
using signalMaps = std::map< uint32_t, signal_map_type >
 
- Protected Member Functions inherited from Phase2TrackerDigitizerAlgorithm
virtual void add_noise (const Phase2TrackerGeomDetUnit *pixdet)
 
virtual void add_noisy_cells (const Phase2TrackerGeomDetUnit *pixdet, float thePixelThreshold)
 
float calcQ (float x)
 
int convertSignalToAdc (uint32_t detID, float signal_in_elec, float threshold)
 
virtual std::vector< DigitizerUtility::SignalPoint > drift (const PSimHit &hit, const Phase2TrackerGeomDetUnit *pixdet, const GlobalVector &bfield, const std::vector< DigitizerUtility::EnergyDepositUnit > &ionization_points) const
 
LocalVector DriftDirection (const Phase2TrackerGeomDetUnit *pixdet, const GlobalVector &bfield, const DetId &detId) const
 
virtual std::vector< float > fluctuateEloss (int particleId, float momentum, float eloss, float length, int NumberOfSegments) const
 
virtual void induce_signal (const PSimHit &hit, const size_t hitIndex, const uint32_t tofBin, const Phase2TrackerGeomDetUnit *pixdet, const std::vector< DigitizerUtility::SignalPoint > &collection_points)
 
virtual void module_killing_conf (uint32_t detID)
 
virtual void pixel_inefficiency (const SubdetEfficiencies &eff, const Phase2TrackerGeomDetUnit *pixdet, const TrackerTopology *tTopo)
 
virtual void pixel_inefficiency_db (uint32_t detID)
 
virtual std::vector< DigitizerUtility::EnergyDepositUnit > primary_ionization (const PSimHit &hit) 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_
 
Parameters badPixels_
 
const float clusterWidth_
 
const Parameters deadModules_
 
const SiPixelFedCablingMap * fedCablingMap_
 
const std::unique_ptr< SiG4UniversalFluctuation > fluctuate_
 
const bool fluctuateCharge_
 
std::unique_ptr< CLHEP::RandGaussQ > gaussDistribution_
 
const TrackerGeometry * geom_
 
const float GeVperElectron_
 
const float interstripCoupling_
 
const bool makeDigiSimLinks_
 
bool pixelFlag_
 
const double pseudoRadDamage_
 
const double pseudoRadDamageRadius_
 
CLHEP::HepRandomEngine * rengine_
 
const float Sigma0_
 
const float SigmaCoeff_
 
const SiPhase2OuterTrackerLorentzAngle * siPhase2OTLorentzAngle_
 
const SiPixelQuality * siPixelBadModule_
 
const 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 int thePhase2ReadoutMode_
 
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 12 of file PixelDigitizerAlgorithm.h.

Constructor & Destructor Documentation

◆ PixelDigitizerAlgorithm()

PixelDigitizerAlgorithm::PixelDigitizerAlgorithm ( const edm::ParameterSet &  conf,
edm::ConsumesCollector  iC 
)

Definition at line 40 of file PixelDigitizerAlgorithm.cc.

References Phase2TrackerDigitizerAlgorithm::addPixelInefficiency_, edm::ConsumesCollector::esConsumes(), LogDebug, Phase2TrackerDigitizerAlgorithm::pixelFlag_, siPixelBadModuleToken_, siPixelLorentzAngleToken_, Phase2TrackerDigitizerAlgorithm::theAdcFullScale_, Phase2TrackerDigitizerAlgorithm::theElectronPerADC_, Phase2TrackerDigitizerAlgorithm::theThresholdInE_Barrel_, Phase2TrackerDigitizerAlgorithm::theThresholdInE_Endcap_, Phase2TrackerDigitizerAlgorithm::tMax_, Phase2TrackerDigitizerAlgorithm::use_deadmodule_DB_, and Phase2TrackerDigitizerAlgorithm::use_LorentzAngle_DB_.

41  : Phase2TrackerDigitizerAlgorithm(conf.getParameter<ParameterSet>("AlgorithmCommon"),
42  conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm"),
43  iC),
44  odd_row_interchannelCoupling_next_row_(conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm")
45  .getParameter<double>("Odd_row_interchannelCoupling_next_row")),
46  even_row_interchannelCoupling_next_row_(conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm")
47  .getParameter<double>("Even_row_interchannelCoupling_next_row")),
48  odd_column_interchannelCoupling_next_column_(
49  conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm")
50  .getParameter<double>("Odd_column_interchannelCoupling_next_column")),
51  even_column_interchannelCoupling_next_column_(
52  conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm")
53  .getParameter<double>("Even_column_interchannelCoupling_next_column")),
54  apply_timewalk_(conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm").getParameter<bool>("ApplyTimewalk")),
55  timewalk_model_(
56  conf.getParameter<ParameterSet>("PixelDigitizerAlgorithm").getParameter<edm::ParameterSet>("TimewalkModel")),
57  fedCablingMapToken_(iC.esConsumes()),
58  geomToken_(iC.esConsumes()) {
59  if (use_deadmodule_DB_)
60  siPixelBadModuleToken_ = iC.esConsumes();
61  if (use_LorentzAngle_DB_)
62  siPixelLorentzAngleToken_ = iC.esConsumes();
63  pixelFlag_ = true;
64  LogDebug("PixelDigitizerAlgorithm") << "Algorithm constructed "
65  << "Configuration parameters:"
66  << "Threshold/Gain = "
67  << "threshold in electron Endcap = " << theThresholdInE_Endcap_
68  << "threshold in electron Barrel = " << theThresholdInE_Barrel_ << " "
69  << theElectronPerADC_ << " " << theAdcFullScale_
70  << " The delta cut-off is set to " << tMax_ << " pix-inefficiency "
71  << addPixelInefficiency_;
72 }
PixelDigitizerAlgorithm::fedCablingMapToken_
const edm::ESGetToken< SiPixelFedCablingMap, SiPixelFedCablingMapRcd > fedCablingMapToken_
Definition: PixelDigitizerAlgorithm.h:67
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
PixelDigitizerAlgorithm::siPixelBadModuleToken_
edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd > siPixelBadModuleToken_
Definition: PixelDigitizerAlgorithm.h:65
PixelDigitizerAlgorithm::timewalk_model_
const TimewalkModel timewalk_model_
Definition: PixelDigitizerAlgorithm.h:63
PixelDigitizerAlgorithm::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: PixelDigitizerAlgorithm.h:68
Phase2TrackerDigitizerAlgorithm::Phase2TrackerDigitizerAlgorithm
Phase2TrackerDigitizerAlgorithm(const edm::ParameterSet &conf_common, const edm::ParameterSet &conf_specific, edm::ConsumesCollector iC)
Definition: Phase2TrackerDigitizerAlgorithm.cc:44
PixelDigitizerAlgorithm::siPixelLorentzAngleToken_
edm::ESGetToken< SiPixelLorentzAngle, SiPixelLorentzAngleSimRcd > siPixelLorentzAngleToken_
Definition: PixelDigitizerAlgorithm.h:66
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233

◆ ~PixelDigitizerAlgorithm()

PixelDigitizerAlgorithm::~PixelDigitizerAlgorithm ( )
override

Definition at line 73 of file PixelDigitizerAlgorithm.cc.

References LogDebug.

73 { LogDebug("PixelDigitizerAlgorithm") << "Algorithm deleted"; }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233

Member Function Documentation

◆ add_cross_talk()

void PixelDigitizerAlgorithm::add_cross_talk ( const Phase2TrackerGeomDetUnit *  pixdet)
overridevirtual

Reimplemented from Phase2TrackerDigitizerAlgorithm.

Definition at line 87 of file PixelDigitizerAlgorithm.cc.

References Phase2TrackerDigitizerAlgorithm::_signal, funct::abs(), officialStyle::chan, PixelDigi::channelToPixel(), even_column_interchannelCoupling_next_column_, even_row_interchannelCoupling_next_row_, GeomDet::geographicalId(), cmsLHEtoEOSManager::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().

87  {
88  if (!pixelFlag_)
89  return;
90 
91  const Phase2TrackerTopology* topol = &pixdet->specificTopology();
92 
93  // cross-talk calculation valid for the case of 25x100 pixels
94  const float pitch_first = 0.0025;
95  const float pitch_second = 0.0100;
96 
97  // 0.5 um tolerance when comparing the pitch to accommodate the small changes in different TK geometrie (temporary fix)
98  const double pitch_tolerance(0.0005);
99 
100  if (std::abs(topol->pitch().first - pitch_first) > pitch_tolerance ||
101  std::abs(topol->pitch().second - pitch_second) > pitch_tolerance)
102  return;
103 
104  uint32_t detID = pixdet->geographicalId().rawId();
105  signal_map_type& theSignal = _signal[detID];
106  signal_map_type signalNew;
107 
108  int numRows = topol->nrows();
109  int numColumns = topol->ncolumns();
110 
111  for (auto& s : theSignal) {
112  float signalInElectrons = s.second.ampl(); // signal in electrons
113 
114  auto hitChan = PixelDigi::channelToPixel(s.first);
115 
116  float signalInElectrons_odd_row_Xtalk_next_row = signalInElectrons * odd_row_interchannelCoupling_next_row_;
117  float signalInElectrons_even_row_Xtalk_next_row = signalInElectrons * even_row_interchannelCoupling_next_row_;
118  float signalInElectrons_odd_column_Xtalk_next_column =
119  signalInElectrons * odd_column_interchannelCoupling_next_column_;
120  float signalInElectrons_even_column_Xtalk_next_column =
121  signalInElectrons * even_column_interchannelCoupling_next_column_;
122 
123  // subtract the charge which will be shared
124  s.second.set(signalInElectrons - signalInElectrons_odd_row_Xtalk_next_row -
125  signalInElectrons_even_row_Xtalk_next_row - signalInElectrons_odd_column_Xtalk_next_column -
126  signalInElectrons_even_column_Xtalk_next_column);
127 
128  if (hitChan.first != 0) {
129  auto XtalkPrev = std::make_pair(hitChan.first - 1, hitChan.second);
130  int chanXtalkPrev = pixelFlag_ ? PixelDigi::pixelToChannel(XtalkPrev.first, XtalkPrev.second)
131  : Phase2TrackerDigi::pixelToChannel(XtalkPrev.first, XtalkPrev.second);
132  if (hitChan.first % 2 == 1)
133  signalNew.emplace(chanXtalkPrev,
134  DigitizerUtility::Amplitude(signalInElectrons_even_row_Xtalk_next_row, nullptr, -1.0));
135  else
136  signalNew.emplace(chanXtalkPrev,
137  DigitizerUtility::Amplitude(signalInElectrons_odd_row_Xtalk_next_row, nullptr, -1.0));
138  }
139  if (hitChan.first < numRows - 1) {
140  auto XtalkNext = std::make_pair(hitChan.first + 1, hitChan.second);
141  int chanXtalkNext = pixelFlag_ ? PixelDigi::pixelToChannel(XtalkNext.first, XtalkNext.second)
142  : Phase2TrackerDigi::pixelToChannel(XtalkNext.first, XtalkNext.second);
143  if (hitChan.first % 2 == 1)
144  signalNew.emplace(chanXtalkNext,
145  DigitizerUtility::Amplitude(signalInElectrons_odd_row_Xtalk_next_row, nullptr, -1.0));
146  else
147  signalNew.emplace(chanXtalkNext,
148  DigitizerUtility::Amplitude(signalInElectrons_even_row_Xtalk_next_row, nullptr, -1.0));
149  }
150 
151  if (hitChan.second != 0) {
152  auto XtalkPrev = std::make_pair(hitChan.first, hitChan.second - 1);
153  int chanXtalkPrev = pixelFlag_ ? PixelDigi::pixelToChannel(XtalkPrev.first, XtalkPrev.second)
154  : Phase2TrackerDigi::pixelToChannel(XtalkPrev.first, XtalkPrev.second);
155  if (hitChan.second % 2 == 1)
156  signalNew.emplace(chanXtalkPrev,
157  DigitizerUtility::Amplitude(signalInElectrons_even_column_Xtalk_next_column, nullptr, -1.0));
158  else
159  signalNew.emplace(chanXtalkPrev,
160  DigitizerUtility::Amplitude(signalInElectrons_odd_column_Xtalk_next_column, nullptr, -1.0));
161  }
162  if (hitChan.second < numColumns - 1) {
163  auto XtalkNext = std::make_pair(hitChan.first, hitChan.second + 1);
164  int chanXtalkNext = pixelFlag_ ? PixelDigi::pixelToChannel(XtalkNext.first, XtalkNext.second)
165  : Phase2TrackerDigi::pixelToChannel(XtalkNext.first, XtalkNext.second);
166  if (hitChan.second % 2 == 1)
167  signalNew.emplace(chanXtalkNext,
168  DigitizerUtility::Amplitude(signalInElectrons_odd_column_Xtalk_next_column, nullptr, -1.0));
169  else
170  signalNew.emplace(chanXtalkNext,
171  DigitizerUtility::Amplitude(signalInElectrons_even_column_Xtalk_next_column, nullptr, -1.0));
172  }
173  }
174  for (auto const& l : signalNew) {
175  int chan = l.first;
176  auto iter = theSignal.find(chan);
177  if (iter != theSignal.end()) {
178  iter->second += l.second.ampl();
179  } else {
180  theSignal.emplace(chan, DigitizerUtility::Amplitude(l.second.ampl(), nullptr, -1.0));
181  }
182  }
183 }
PixelDigi::channelToPixel
static std::pair< int, int > channelToPixel(int ch)
Definition: PixelDigi.h:69
PixelTopology::ncolumns
virtual int ncolumns() const =0
PixelTopology::nrows
virtual int nrows() const =0
PixelDigi::pixelToChannel
static int pixelToChannel(int row, int col)
Definition: PixelDigi.h:75
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
GeomDet::geographicalId
DetId geographicalId() const
The label of this GeomDet.
Definition: GeomDet.h:64
Phase2TrackerDigitizerAlgorithm::signal_map_type
std::map< int, DigitizerUtility::Amplitude, std::less< int > > signal_map_type
Definition: Phase2TrackerDigitizerAlgorithm.h:103
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
officialStyle.chan
chan
lumi = TPaveText(lowX+0.38, lowY+0.061, lowX+0.45, lowY+0.161, "NDC") lumi.SetBorderSize( 0 ) lumi...
Definition: officialStyle.py:106
PixelGeomDetUnit::specificTopology
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
Definition: PixelGeomDetUnit.cc:17
PixelTopology::pitch
virtual std::pair< float, float > pitch() const =0
pixelgpudetails::pixelToChannel
constexpr uint32_t pixelToChannel(int row, int col)
Definition: SiPixelRawToClusterGPUKernel.h:69

◆ init()

void PixelDigitizerAlgorithm::init ( const edm::EventSetup &  es)
overridevirtual

Implements Phase2TrackerDigitizerAlgorithm.

Definition at line 25 of file PixelDigitizerAlgorithm.cc.

References edm::EventSetup::getData().

25  {
26  if (use_ineff_from_db_) // load gain calibration service fromdb...
27  theSiPixelGainCalibrationService_->setESObjects(es);
28 
29  if (use_deadmodule_DB_)
30  siPixelBadModule_ = &es.getData(siPixelBadModuleToken_);
31 
32  if (use_LorentzAngle_DB_) // Get Lorentz angle from DB record
33  siPixelLorentzAngle_ = &es.getData(siPixelLorentzAngleToken_);
34 
35  // gets the map and geometry from the DB (to kill ROCs)
36  fedCablingMap_ = &es.getData(fedCablingMapToken_);
37  geom_ = &es.getData(geomToken_);
38 }
PixelDigitizerAlgorithm::fedCablingMapToken_
const edm::ESGetToken< SiPixelFedCablingMap, SiPixelFedCablingMapRcd > fedCablingMapToken_
Definition: PixelDigitizerAlgorithm.h:67
PixelDigitizerAlgorithm::siPixelBadModuleToken_
edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd > siPixelBadModuleToken_
Definition: PixelDigitizerAlgorithm.h:65
Phase2TrackerDigitizerAlgorithm::siPixelLorentzAngle_
const SiPixelLorentzAngle * siPixelLorentzAngle_
Definition: Phase2TrackerDigitizerAlgorithm.h:85
Phase2TrackerDigitizerAlgorithm::theSiPixelGainCalibrationService_
const std::unique_ptr< SiPixelGainCalibrationOfflineSimService > theSiPixelGainCalibrationService_
Definition: Phase2TrackerDigitizerAlgorithm.h:210
PixelDigitizerAlgorithm::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: PixelDigitizerAlgorithm.h:68
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:122
PixelDigitizerAlgorithm::siPixelLorentzAngleToken_
edm::ESGetToken< SiPixelLorentzAngle, SiPixelLorentzAngleSimRcd > siPixelLorentzAngleToken_
Definition: PixelDigitizerAlgorithm.h:66

◆ isAboveThreshold()

bool PixelDigitizerAlgorithm::isAboveThreshold ( const DigitizerUtility::SimHitInfo *  hitInfo,
float  charge,
float  thr 
) const
overridevirtual

Reimplemented from Phase2TrackerDigitizerAlgorithm.

Definition at line 237 of file PixelDigitizerAlgorithm.cc.

References apply_timewalk_, ALCARECOTkAlJpsiMuMu_cff::charge, Phase2TrackerDigitizerAlgorithm::theTofLowerCut_, Phase2TrackerDigitizerAlgorithm::theTofUpperCut_, DigitizerUtility::SimHitInfo::time(), protons_cff::time, and timewalk_model_.

239  {
240  if (charge < thr)
241  return false;
242  if (apply_timewalk_ && hitInfo) {
243  float corrected_time = hitInfo->time();
244  double time = corrected_time + timewalk_model_(charge, thr);
245  return (time >= theTofLowerCut_ && time < theTofUpperCut_);
246  } else
247  return true;
248 }
PixelDigitizerAlgorithm::timewalk_model_
const TimewalkModel timewalk_model_
Definition: PixelDigitizerAlgorithm.h:63
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
protons_cff.time
time
Definition: protons_cff.py:36

◆ module_killing_DB()

void PixelDigitizerAlgorithm::module_killing_DB ( const Phase2TrackerGeomDetUnit *  pixdet)
overridevirtual

Implements Phase2TrackerDigitizerAlgorithm.

Definition at line 252 of file PixelDigitizerAlgorithm.cc.

References Phase2TrackerDigitizerAlgorithm::_signal, funct::abs(), Phase2TrackerDigitizerAlgorithm::badPixels_, Phase2TrackerDigi::channelToPixel(), PixelDigi::channelToPixel(), SiPixelQuality::disabledModuleType::errorType, Phase2TrackerDigitizerAlgorithm::fedCablingMap_, SiPixelFedCablingMap::findItem(), GeomDet::geographicalId(), SiPixelQuality::getBadComponentList(), sipixelobjects::PixelROC::idInDetUnit(), SiPixelQuality::IsRocBad(), dqmiolumiharvest::j, dqmdumpme::k, DTRecHitClients_cfi::local, mod(), PixelTopology::ncolumns(), AlCaHLTBitMon_ParallelJobs::p, castor_dqm_sourceclient_file_cfg::path, SiPixelFedCablingMap::pathToDetUnit(), Phase2TrackerDigitizerAlgorithm::pixelFlag_, DetId::rawId(), alignCSCRings::s, Phase2TrackerDigitizerAlgorithm::siPixelBadModule_, PixelGeomDetUnit::specificTopology(), and sipixelobjects::PixelROC::toGlobal().

252  {
253  bool isbad = false;
254  uint32_t detID = pixdet->geographicalId().rawId();
255  int ncol = pixdet->specificTopology().ncolumns();
256  if (ncol < 0)
257  return;
258  std::vector<SiPixelQuality::disabledModuleType> disabledModules = siPixelBadModule_->getBadComponentList();
259 
260  SiPixelQuality::disabledModuleType badmodule;
261  for (const auto& mod : disabledModules) {
262  if (detID == mod.DetID) {
263  isbad = true;
264  badmodule = mod;
265  break;
266  }
267  }
268 
269  if (!isbad)
270  return;
271 
272  signal_map_type& theSignal = _signal[detID]; // check validity
273  if (badmodule.errorType == 0) { // this is a whole dead module.
274  for (auto& s : theSignal)
275  s.second.set(0.); // reset amplitude
276  } else { // all other module types: half-modules and single ROCs.
277  // Get Bad ROC position:
278  // follow the example of getBadRocPositions in CondFormats/SiPixelObjects/src/SiPixelQuality.cc
279  std::vector<GlobalPixel> badrocpositions;
280  for (size_t j = 0; j < static_cast<size_t>(ncol); j++) {
281  if (siPixelBadModule_->IsRocBad(detID, j)) {
282  std::vector<CablingPathToDetUnit> path = fedCablingMap_->pathToDetUnit(detID);
283  for (auto const& p : path) {
284  const PixelROC* myroc = fedCablingMap_->findItem(p);
285  if (myroc->idInDetUnit() == j) {
286  LocalPixel::RocRowCol local = {39, 25}; //corresponding to center of ROC row, col
287  GlobalPixel global = myroc->toGlobal(LocalPixel(local));
288  badrocpositions.push_back(global);
289  break;
290  }
291  }
292  }
293  }
294 
295  for (auto& s : theSignal) {
296  std::pair<int, int> ip;
297  if (pixelFlag_)
298  ip = PixelDigi::channelToPixel(s.first);
299  else
300  ip = Phase2TrackerDigi::channelToPixel(s.first);
301 
302  for (auto const& p : badrocpositions) {
303  for (auto& k : badPixels_) {
304  if (p.row == k.getParameter<int>("row") && ip.first == k.getParameter<int>("row") &&
305  std::abs(ip.second - p.col) < k.getParameter<int>("col")) {
306  s.second.set(0.);
307  }
308  }
309  }
310  }
311  }
312 }
SiPixelQuality::IsRocBad
bool IsRocBad(const uint32_t &detid, const short &rocNb) const
Definition: SiPixelQuality.cc:107
PixelDigi::channelToPixel
static std::pair< int, int > channelToPixel(int ch)
Definition: PixelDigi.h:69
PixelTopology::ncolumns
virtual int ncolumns() const =0
castor_dqm_sourceclient_file_cfg.path
path
Definition: castor_dqm_sourceclient_file_cfg.py:37
SiPixelFedCablingMap::pathToDetUnit
std::vector< sipixelobjects::CablingPathToDetUnit > pathToDetUnit(uint32_t rawDetId) const final
Definition: SiPixelFedCablingMap.cc:150
sipixelobjects::LocalPixel
identify pixel inside single ROC
Definition: LocalPixel.h:7
sipixelobjects::GlobalPixel
global coordinates (row and column in DetUnit, as in PixelDigi)
Definition: GlobalPixel.h:6
SiPixelQuality::getBadComponentList
const std::vector< disabledModuleType > getBadComponentList() const
Definition: SiPixelQuality.h:100
sipixelobjects::PixelROC::toGlobal
GlobalPixel toGlobal(const LocalPixel &loc) const
Definition: PixelROC.h:55
SiPixelFedCablingMap::findItem
const sipixelobjects::PixelROC * findItem(const sipixelobjects::CablingPathToDetUnit &path) const final
Definition: SiPixelFedCablingMap.cc:124
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
GeomDet::geographicalId
DetId geographicalId() const
The label of this GeomDet.
Definition: GeomDet.h:64
Phase2TrackerDigitizerAlgorithm::signal_map_type
std::map< int, DigitizerUtility::Amplitude, std::less< int > > signal_map_type
Definition: Phase2TrackerDigitizerAlgorithm.h:103
DTRecHitClients_cfi.local
local
Definition: DTRecHitClients_cfi.py:10
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
sipixelobjects::LocalPixel::RocRowCol
row and collumn in ROC representation
Definition: LocalPixel.h:13
Phase2TrackerDigi::channelToPixel
static std::pair< unsigned int, unsigned int > channelToPixel(unsigned int ch)
Definition: Phase2TrackerDigi.h:39
PixelGeomDetUnit::specificTopology
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
Definition: PixelGeomDetUnit.cc:17
mod
T mod(const T &a, const T &b)
Definition: ecalDccMap.h:4
sipixelobjects::PixelROC::idInDetUnit
unsigned int idInDetUnit() const
id of this ROC in DetUnit etermined by token path
Definition: PixelROC.h:37
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153

◆ select_hit()

bool PixelDigitizerAlgorithm::select_hit ( const PSimHit &  hit,
double  tCorr,
double &  sigScale 
) const
overridevirtual

Reimplemented from Phase2TrackerDigitizerAlgorithm.

Definition at line 77 of file PixelDigitizerAlgorithm.cc.

References Phase2TrackerDigitizerAlgorithm::theTofLowerCut_, Phase2TrackerDigitizerAlgorithm::theTofUpperCut_, and protons_cff::time.

77  {
78  double time = hit.tof() - tCorr;
79  return (time >= theTofLowerCut_ && time < theTofUpperCut_);
80 }
protons_cff.time
time
Definition: protons_cff.py:36

Member Data Documentation

◆ apply_timewalk_

bool PixelDigitizerAlgorithm::apply_timewalk_

Definition at line 62 of file PixelDigitizerAlgorithm.h.

Referenced by isAboveThreshold().

◆ even_column_interchannelCoupling_next_column_

double PixelDigitizerAlgorithm::even_column_interchannelCoupling_next_column_

Definition at line 59 of file PixelDigitizerAlgorithm.h.

Referenced by add_cross_talk().

◆ even_row_interchannelCoupling_next_row_

double PixelDigitizerAlgorithm::even_row_interchannelCoupling_next_row_

Definition at line 57 of file PixelDigitizerAlgorithm.h.

Referenced by add_cross_talk(), and PixelBrickedDigitizerAlgorithm::PixelBrickedDigitizerAlgorithm().

◆ fedCablingMapToken_

const edm::ESGetToken<SiPixelFedCablingMap, SiPixelFedCablingMapRcd> PixelDigitizerAlgorithm::fedCablingMapToken_

Definition at line 67 of file PixelDigitizerAlgorithm.h.

◆ geomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> PixelDigitizerAlgorithm::geomToken_

Definition at line 68 of file PixelDigitizerAlgorithm.h.

◆ odd_column_interchannelCoupling_next_column_

double PixelDigitizerAlgorithm::odd_column_interchannelCoupling_next_column_

Definition at line 58 of file PixelDigitizerAlgorithm.h.

Referenced by add_cross_talk().

◆ odd_row_interchannelCoupling_next_row_

double PixelDigitizerAlgorithm::odd_row_interchannelCoupling_next_row_

Definition at line 56 of file PixelDigitizerAlgorithm.h.

Referenced by add_cross_talk().

◆ siPixelBadModuleToken_

edm::ESGetToken<SiPixelQuality, SiPixelQualityRcd> PixelDigitizerAlgorithm::siPixelBadModuleToken_

Definition at line 65 of file PixelDigitizerAlgorithm.h.

Referenced by PixelDigitizerAlgorithm().

◆ siPixelLorentzAngleToken_

edm::ESGetToken<SiPixelLorentzAngle, SiPixelLorentzAngleSimRcd> PixelDigitizerAlgorithm::siPixelLorentzAngleToken_

Definition at line 66 of file PixelDigitizerAlgorithm.h.

Referenced by PixelDigitizerAlgorithm().

◆ timewalk_model_

const TimewalkModel PixelDigitizerAlgorithm::timewalk_model_

Definition at line 63 of file PixelDigitizerAlgorithm.h.

Referenced by isAboveThreshold().

Generated for CMSSW Reference Manual by doxygen 1.8.14