1 #ifndef SiPixelDigitizerAlgorithm_h 2 #define SiPixelDigitizerAlgorithm_h 20 #include "boost/multi_array.hpp" 29 class HepRandomEngine;
64 void accumulateSimHits(
const std::vector<PSimHit>::const_iterator inputBegin,
65 const std::vector<PSimHit>::const_iterator inputEnd,
66 const size_t inputBeginGlobalIndex,
67 const unsigned int tofBin,
71 CLHEP::HepRandomEngine*);
73 std::vector<PixelDigi>& digis,
74 std::vector<PixelDigiSimLink>& simlinks,
76 CLHEP::HepRandomEngine*);
81 void calculateInstlumiFactor(
const std::vector<PileupSummaryInfo> &ps,
int bunchSpacing);
82 void setSimAccumulator(
const std::map<uint32_t, std::map<int, int> >& signalMap);
107 _amp(amp), _frac(1, frac) {
116 _amp(amp), _frac(1, frac) {
124 _hitInfos.emplace_back(hitp, hitIndex, tofBin);
129 operator float()
const {
return _amp;}
130 float ampl()
const {
return _amp;}
132 const std::vector<SimHitInfoForLinks>&
hitInfos()
const {
return _hitInfos; }
138 if (other.
_frac[0]>-0.5){
142 _frac.insert(_frac.end(), other.
_frac.begin(), other.
_frac.end());
180 _energy(energy),_position(x,y,z){}
182 _energy(energy),_position(position){}
183 float x()
const{
return _position.x();}
184 float y()
const{
return _position.y();}
185 float z()
const{
return _position.z();}
202 _sigma_x(1.), _sigma_y(1.), _hitp(
nullptr) {}
205 float t,
float a=1.0) :
206 _pos(x,y), _time(t), _amplitude(
a), _sigma_x(sigma_x),
207 _sigma_y(sigma_y), _hitp(
nullptr) {}
211 _pos(x,y), _time(t), _amplitude(
a), _sigma_x(sigma_x),
212 _sigma_y(sigma_y),_hitp(&hit) {}
215 float x()
const {
return _pos.x();}
216 float y()
const {
return _pos.y();}
219 float time()
const {
return _time;}
250 double thePixelEfficiency[20];
251 double thePixelColEfficiency[20];
252 double thePixelChipEfficiency[20];
253 std::vector<double> theLadderEfficiency_BPix[20];
254 std::vector<double> theModuleEfficiency_BPix[20];
255 double theInnerEfficiency_FPix[20];
256 double theOuterEfficiency_FPix[20];
265 std::map<uint32_t, size_t >
iPU;
269 static const int rocIdMaskBits = 0x1F;
282 float thePixelPseudoRadDamage[20];
408 const std::unique_ptr<SiG4UniversalFluctuation>
fluctuate;
409 const std::unique_ptr<GaussianTailNoiseGenerator>
theNoiser;
412 const std::map<int,CalParameters,std::less<int> >
calmap;
417 std::map<int,CalParameters,std::less<int> > initCal()
const;
418 void primary_ionization(
const PSimHit&
hit, std::vector<EnergyDepositUnit>& ionization_points, CLHEP::HepRandomEngine*)
const;
423 const std::vector<EnergyDepositUnit>& ionization_points,
424 std::vector<SignalPoint>& collection_points)
const;
425 void induce_signal(std::vector<PSimHit>::const_iterator inputBegin,
426 std::vector<PSimHit>::const_iterator inputEnd,
428 const size_t hitIndex,
429 const unsigned int tofBin,
431 const std::vector<SignalPoint>& collection_points);
432 void fluctuateEloss(
int particleId,
float momentum,
float eloss,
433 float length,
int NumberOfSegments,
435 CLHEP::HepRandomEngine*)
const;
437 float thePixelThreshold,
438 CLHEP::HepRandomEngine*);
439 void make_digis(
float thePixelThresholdInE,
442 std::vector<PixelDigi>& digis,
443 std::vector<PixelDigiSimLink>& simlinks,
448 CLHEP::HepRandomEngine*);
450 void pixel_inefficiency_db(uint32_t detID);
461 const DetId& detId)
const;
463 void module_killing_conf(uint32_t detID);
464 void module_killing_DB(uint32_t detID);
467 int PixelTempRewgt2D(
int id_gen,
int id_rewgt,
470 std::map<
int,
float, std::less<int> >& hit_signal,
471 const size_t hitIndex,
472 const unsigned int tofBin,
475 signal_map_type& theSignal,
476 unsigned short int processType);
477 void printCluster(array_2d& cluster);
478 void printCluster(
float arr[
BXM2][
BYM2]);
494 return 0.5*(1.0-std::copysign(
std::sqrt(1.
f- unsafe_expf<4>(-
xx*(1.
f+0.2733
f/(1.
f+0.147
f*
xx)) )),x));
const SiPixel2DTemplateDBObject * dbobject_den
void operator+=(const Amplitude &other)
Amplitude(float amp, const PSimHit *hitp, size_t hitIndex, unsigned int tofBin, float frac)
const float theThresholdInE_BPix_L2
signal_map_type::const_iterator signal_map_const_iterator
const bool use_deadmodule_DB_
edm::ESHandle< SiPixelFedCablingMap > map_
const float electronsPerVCAL
const double theThresholdSmearing_FPix
SignalPoint(float x, float y, float sigma_x, float sigma_y, float t, const PSimHit &hit, float a=1.0)
const std::unique_ptr< SiPixelGainCalibrationOfflineSimService > theSiPixelGainCalibrationService_
std::vector< float > track
const std::vector< float > & individualampl() const
EnergyDepositUnit(float energy, Local3DPoint position)
const float tanLorentzAnglePerTesla_FPix
PixelEfficiencies pixelEfficiencies_
const std::unique_ptr< SiG4UniversalFluctuation > fluctuate
const int theAdcFullScale
SignalPoint(float x, float y, float sigma_x, float sigma_y, float t, float a=1.0)
const LocalPoint & position() const
const double theThresholdSmearing_BPix_L1
LocalVector drift(const StripGeomDetUnit *, const MagneticField &, const SiStripLorentzAngle &)
const float theThresholdInE_FPix
const double theThresholdSmearing_BPix
std::vector< std::vector< double > > thePUEfficiency
const float electronsPerVCAL_Offset
const bool addThresholdSmearing
Amplitude(float amp, float frac)
const bool fluctuateCharge
edm::ESHandle< TrackerGeometry > geom_
double calcQ(float x) const
const float GeVperElectron
const double theThresholdSmearing_BPix_L2
const bool use_ineff_from_db_
const SiPixel2DTemplateDBObject * dbobject_num
const bool use_LorentzAngle_DB_
const std::map< int, CalParameters, std::less< int > > calmap
double theInstLumiScaleFactor
const Parameters DeadModules
boost::multi_array< float, 2 > array_2d
std::map< uint32_t, std::vector< double > > PixelGeomFactorsROCBigPixels
const float theTofUpperCut
SiPixelTemplate2D templ2D
const bool use_module_killing_
std::vector< double > pu_scale
const bool PrintTemplates
const float theThresholdInE_BPix
const int NumberOfEndcapDisks
signal_map_type::iterator signal_map_iterator
SignalPoint & set_amplitude(float amp)
EnergyDepositUnit(float energy, float x, float y, float z)
const float theThresholdInE_BPix_L1
edm::ESHandle< SiPixelLorentzAngle > SiPixelLorentzAngle_
const float electronsPerVCAL_L1_Offset
boost::multi_array< float, 2 > array_2d
const bool UseReweighting
const bool doMissCalibrate
const bool AddPixelInefficiency
const bool addChargeVCALSmearing
const std::unique_ptr< GaussianTailNoiseGenerator > theNoiser
const float theNoiseInElectrons
std::map< int, Amplitude, std::less< int > > signal_map_type
const float theOffsetSmearing
std::vector< float > _frac
const float electronsPerVCAL_L1
edm::ESHandle< SiPixelQuality > SiPixelBadModule_
std::map< uint32_t, std::vector< double > > PixelGeomFactorsROCStdPixels
std::vector< bool > ydouble
GloballyPositioned< double > Frame
const float tanLorentzAnglePerTesla_BPix
std::vector< SiPixelTemplateStore2D > templateStores_
edm::ESHandle< SiPixelDynamicInefficiency > SiPixelDynamicInefficiency_
std::vector< bool > xdouble
std::vector< edm::ParameterSet > Parameters
const std::vector< SimHitInfoForLinks > & hitInfos() const
const float theGainSmearing
std::vector< SimHitInfoForLinks > _hitInfos
void operator+=(const float &)
static int position[264][3]
const float theReadoutNoise
std::map< uint32_t, double > ColGeomFactors
const float theTofLowerCut
std::map< uint32_t, double > PixelGeomFactors
std::map< uint32_t, double > ChipGeomFactors
const bool addNoisyPixels
const PixelAging pixelAging_
std::map< uint32_t, signal_map_type > signalMaps
std::map< uint32_t, size_t > iPU
const float theElectronPerADC
const bool makeDigiSimLinks_
const int NumberOfBarrelLayers