CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes | Friends
HcalSimParameters Class Reference

#include <HcalSimParameters.h>

Inheritance diagram for HcalSimParameters:
CaloSimParameters

Public Member Functions

int delayQIE () const
 
bool doSiPMSmearing () const
 
bool doTimeSmear () const
 
double fCtoGeV (const DetId &detId) const
 
 HcalSimParameters (double simHitToPhotoelectrons, double samplingFactor, double timePhase, int readoutFrameSize, int binOfMaximum, bool doPhotostatistics, bool syncPhase, int firstRing, const std::vector< double > &samplingFactors, double sipmTau)
 
 HcalSimParameters (const edm::ParameterSet &p)
 
double photoelectronsToAnalog (const DetId &detId) const override
 
int pixels (const DetId &detId) const
 
virtual double samplingFactor (const DetId &detId) const
 
void setDbService (const HcalDbService *service)
 
unsigned int signalShape (const DetId &detId) const
 
double simHitToPhotoelectrons (const DetId &detId) const override
 
double sipmCrossTalk (const DetId &detId) const
 
double sipmDarkCurrentuA (const DetId &detId) const
 
std::vector< float > sipmNonlinearity (const DetId &detId) const
 
double sipmTau () const
 
double threshold_currentTDC () const
 
double timeSmearRMS (double ampl) const
 
 ~HcalSimParameters () override
 
- Public Member Functions inherited from CaloSimParameters
int binOfMaximum () const
 
 CaloSimParameters (double simHitToPhotoelectrons, double photoelectronsToAnalog, double samplingFactor, double timePhase, int readoutFrameSize, int binOfMaximum, bool doPhotostatistics, bool syncPhase=true)
 
 CaloSimParameters (const edm::ParameterSet &p, bool skipPe2Fc=false)
 
bool doPhotostatistics () const
 whether or not to apply Poisson statistics to photoelectrons More...
 
double photoelectronsToAnalog () const
 the factor which goes from photoelectrons to whatever gets read by ADCs More...
 
int readoutFrameSize () const
 for now, the LinearFrames and trhe digis will be one-to-one. More...
 
void setBinOfMaximum (int binOfMax)
 
void setReadoutFrameSize (int frameSize)
 some datamixing apps need this to be set dynamically More...
 
double simHitToPhotoelectrons () const
 
bool syncPhase () const
 
double timePhase () const
 the adjustment you need to apply to get the signal where you want it More...
 
virtual ~CaloSimParameters ()
 

Private Member Functions

void defaultTimeSmearing ()
 

Private Attributes

int delayQIE_
 
bool doTimeSmear_
 
const HcalDbServicetheDbService
 
int theFirstRing
 
std::vector< double > theSamplingFactors
 
const HcalSiPMCharacteristicstheSiPMcharacteristics
 
bool theSiPMSmearing
 
double theSiPMTau
 
HcalTimeSmearSettings theSmearSettings
 
double threshold_currentTDC_
 

Friends

class HcalSimParameterMap
 

Detailed Description

Definition at line 10 of file HcalSimParameters.h.

Constructor & Destructor Documentation

HcalSimParameters::HcalSimParameters ( double  simHitToPhotoelectrons,
double  samplingFactor,
double  timePhase,
int  readoutFrameSize,
int  binOfMaximum,
bool  doPhotostatistics,
bool  syncPhase,
int  firstRing,
const std::vector< double > &  samplingFactors,
double  sipmTau 
)

Definition at line 12 of file HcalSimParameters.cc.

References defaultTimeSmearing(), and theSiPMSmearing.

23  0.0,
25  timePhase,
29  syncPhase),
30  theDbService(nullptr),
31  theSiPMcharacteristics(nullptr),
32  theFirstRing(firstRing),
33  theSamplingFactors(samplingFactors),
34  theSiPMSmearing(false),
35  doTimeSmear_(true),
38 
39  edm::LogInfo("HcalSimParameters:") << " doSiPMsmearing = " << theSiPMSmearing;
40 }
const HcalSiPMCharacteristics * theSiPMcharacteristics
double sipmTau() const
CaloSimParameters(double simHitToPhotoelectrons, double photoelectronsToAnalog, double samplingFactor, double timePhase, int readoutFrameSize, int binOfMaximum, bool doPhotostatistics, bool syncPhase=true)
bool doPhotostatistics() const
whether or not to apply Poisson statistics to photoelectrons
double timePhase() const
the adjustment you need to apply to get the signal where you want it
double simHitToPhotoelectrons() const
std::vector< double > theSamplingFactors
virtual double samplingFactor(const DetId &detId) const
int readoutFrameSize() const
for now, the LinearFrames and trhe digis will be one-to-one.
Log< level::Info, false > LogInfo
const HcalDbService * theDbService
int binOfMaximum() const
bool syncPhase() const
HcalSimParameters::HcalSimParameters ( const edm::ParameterSet p)

Definition at line 42 of file HcalSimParameters.cc.

References defaultTimeSmearing(), and theSiPMSmearing.

43  : CaloSimParameters(p, true),
44  theDbService(nullptr),
45  theFirstRing(p.getParameter<int>("firstRing")),
46  theSamplingFactors(p.getParameter<std::vector<double> >("samplingFactors")),
47  theSiPMSmearing(p.getParameter<bool>("doSiPMSmearing")),
48  doTimeSmear_(p.getParameter<bool>("timeSmearing")),
49  theSiPMTau(p.getParameter<double>("sipmTau")),
50  threshold_currentTDC_(p.getParameter<double>("threshold_currentTDC")),
51  delayQIE_(p.getParameter<int>("delayQIE")) {
53 
54  edm::LogInfo("HcalSimParameters:") << " doSiPMsmearing = " << theSiPMSmearing;
55 }
CaloSimParameters(double simHitToPhotoelectrons, double photoelectronsToAnalog, double samplingFactor, double timePhase, int readoutFrameSize, int binOfMaximum, bool doPhotostatistics, bool syncPhase=true)
std::vector< double > theSamplingFactors
Log< level::Info, false > LogInfo
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
const HcalDbService * theDbService
HcalSimParameters::~HcalSimParameters ( )
inlineoverride

Definition at line 24 of file HcalSimParameters.h.

24 {}

Member Function Documentation

void HcalSimParameters::defaultTimeSmearing ( )
private

Definition at line 108 of file HcalSimParameters.cc.

References GeV2fC, and theSmearSettings.

Referenced by HcalSimParameters().

108  {
109  // GeV->ampl (fC), time (ns)
110  theSmearSettings.emplace_back(4.00 * GeV2fC, 4.050);
111  theSmearSettings.emplace_back(20.00 * GeV2fC, 3.300);
112  theSmearSettings.emplace_back(25.00 * GeV2fC, 2.925);
113  theSmearSettings.emplace_back(30.00 * GeV2fC, 2.714);
114  theSmearSettings.emplace_back(37.00 * GeV2fC, 2.496);
115  theSmearSettings.emplace_back(44.50 * GeV2fC, 2.278);
116  theSmearSettings.emplace_back(56.00 * GeV2fC, 2.138);
117  theSmearSettings.emplace_back(63.50 * GeV2fC, 2.022);
118  theSmearSettings.emplace_back(81.00 * GeV2fC, 1.788);
119  theSmearSettings.emplace_back(88.50 * GeV2fC, 1.695);
120  theSmearSettings.emplace_back(114.50 * GeV2fC, 1.716);
121  theSmearSettings.emplace_back(175.50 * GeV2fC, 1.070);
122  theSmearSettings.emplace_back(350.00 * GeV2fC, 1.564);
123  theSmearSettings.emplace_back(99999.00 * GeV2fC, 1.564);
124 }
static const double GeV2fC
HcalTimeSmearSettings theSmearSettings
int HcalSimParameters::delayQIE ( ) const
inline

Definition at line 43 of file HcalSimParameters.h.

References delayQIE_.

Referenced by HcalAmplifier::amplify().

43 { return delayQIE_; }
bool HcalSimParameters::doSiPMSmearing ( ) const
inline

Definition at line 42 of file HcalSimParameters.h.

References theSiPMSmearing.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

42 { return theSiPMSmearing; }
bool HcalSimParameters::doTimeSmear ( ) const
inline

Definition at line 37 of file HcalSimParameters.h.

References doTimeSmear_.

Referenced by HcalTimeSlewSim::delay().

37 { return doTimeSmear_; }
double HcalSimParameters::fCtoGeV ( const DetId detId) const

{

Definition at line 75 of file HcalSimParameters.cc.

References cms::cuda::assert(), HcalDbService::getGain(), HcalDbService::getGainWidth(), HcalGain::getValue(), mps_fire::result, and theDbService.

Referenced by simHitToPhotoelectrons().

75  {
76  assert(theDbService != nullptr);
77  HcalGenericDetId hcalGenDetId(detId);
78  const HcalGain* gains = theDbService->getGain(hcalGenDetId);
79  const HcalGainWidth* gwidths = theDbService->getGainWidth(hcalGenDetId);
80  double result = 0.0;
81  if (!gains || !gwidths) {
82  edm::LogError("HcalAmplifier") << "Could not fetch HCAL conditions for channel " << hcalGenDetId;
83  } else {
84  // only one gain will be recorded per channel, so just use capID 0 for now
85  result = gains->getValue(0);
86  // if(doNoise_)
88  // result += CLHEP::RandGaussQ::shoot(0., gwidths->getValue(0));
89  // }
90  }
91  return result;
92 }
const HcalGainWidth * getGainWidth(const HcalGenericDetId &fId) const
Log< level::Error, false > LogError
float getValue(int fCapId) const
get value for capId = 0..3
Definition: HcalGain.h:21
assert(be >=bs)
tuple result
Definition: mps_fire.py:311
const HcalGain * getGain(const HcalGenericDetId &fId) const
const HcalDbService * theDbService
double HcalSimParameters::photoelectronsToAnalog ( const DetId detId) const
overridevirtual

Reimplemented from CaloSimParameters.

Definition at line 99 of file HcalSimParameters.cc.

References cms::cuda::assert(), HcalSiPMParameter::getFCByPE(), HcalDbService::getHcalSiPMParameter(), and theDbService.

Referenced by HcalSiPMHitResponse::addPEnoise().

99  {
100  //now always taken from database for HPDs or SiPMs (HB, HE, HO)
102  return theDbService->getHcalSiPMParameter(detId)->getFCByPE();
103 }
assert(be >=bs)
const HcalDbService * theDbService
const HcalSiPMParameter * getHcalSiPMParameter(const HcalGenericDetId &fId) const
float getFCByPE() const
get fcByPE
int HcalSimParameters::pixels ( const DetId detId) const

Definition at line 150 of file HcalSimParameters.cc.

References cms::cuda::assert(), HcalDbService::getHcalSiPMParameter(), HcalSiPMCharacteristics::getPixels(), HcalSiPMParameter::getType(), theDbService, and theSiPMcharacteristics.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

150  {
153  return theSiPMcharacteristics->getPixels(type);
154 }
const HcalSiPMCharacteristics * theSiPMcharacteristics
int getType() const
get SiPM type
assert(be >=bs)
const HcalDbService * theDbService
int getPixels(int type) const
get # of pixels
const HcalSiPMParameter * getHcalSiPMParameter(const HcalGenericDetId &fId) const
double HcalSimParameters::samplingFactor ( const DetId detId) const
virtual

the ratio of actual incident energy to deposited energy in the SimHit

Definition at line 94 of file HcalSimParameters.cc.

References HcalDetId::ietaAbs(), theFirstRing, and theSamplingFactors.

Referenced by simHitToPhotoelectrons().

94  {
95  HcalDetId hcalDetId(detId);
96  return theSamplingFactors.at(hcalDetId.ietaAbs() - theFirstRing);
97 }
std::vector< double > theSamplingFactors
void HcalSimParameters::setDbService ( const HcalDbService service)

Definition at line 57 of file HcalSimParameters.cc.

References cms::cuda::assert(), HcalDbService::getHcalSiPMCharacteristics(), theDbService, and theSiPMcharacteristics.

Referenced by HcalSimParameterMap::setDbService().

57  {
58  assert(service);
59  theDbService = service;
62 }
const HcalSiPMCharacteristics * theSiPMcharacteristics
assert(be >=bs)
const HcalSiPMCharacteristics * getHcalSiPMCharacteristics() const
const HcalDbService * theDbService
unsigned int HcalSimParameters::signalShape ( const DetId detId) const

Definition at line 172 of file HcalSimParameters.cc.

References cms::cuda::assert(), HcalDbService::getHcalMCParam(), HcalMCParam::signalShape(), and theDbService.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

172  {
174  return theDbService->getHcalMCParam(detId)->signalShape();
175 }
assert(be >=bs)
const HcalMCParam * getHcalMCParam(const HcalGenericDetId &fId) const
unsigned int signalShape() const
Definition: HcalMCParam.h:38
const HcalDbService * theDbService
double HcalSimParameters::simHitToPhotoelectrons ( const DetId detId) const
overridevirtual

Reimplemented from CaloSimParameters.

Definition at line 64 of file HcalSimParameters.cc.

References fCtoGeV(), HcalGenericDetId::HcalGenForward, HcalGenericDetId::HcalGenZDC, CaloSimParameters::photoelectronsToAnalog(), mps_fire::result, samplingFactor(), and CaloSimParameters::simHitToPhotoelectrons().

64  {
65  // the gain is in units of GeV/fC. We want a constant with pe/dGeV
66  // pe/dGeV = (GeV/dGeV) / (GeV/fC) / (fC/pe)
68  if (HcalGenericDetId(detId).genericSubdet() != HcalGenericDetId::HcalGenForward ||
69  HcalGenericDetId(detId).genericSubdet() != HcalGenericDetId::HcalGenZDC) {
70  result = samplingFactor(detId) / fCtoGeV(detId) / photoelectronsToAnalog(detId);
71  }
72  return result;
73 }
double fCtoGeV(const DetId &detId) const
tuple result
Definition: mps_fire.py:311
double simHitToPhotoelectrons() const
virtual double samplingFactor(const DetId &detId) const
double photoelectronsToAnalog() const
the factor which goes from photoelectrons to whatever gets read by ADCs
double HcalSimParameters::sipmCrossTalk ( const DetId detId) const

Definition at line 161 of file HcalSimParameters.cc.

References cms::cuda::assert(), HcalSiPMCharacteristics::getCrossTalk(), HcalDbService::getHcalSiPMParameter(), HcalSiPMParameter::getType(), theDbService, and theSiPMcharacteristics.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

161  {
164  return theSiPMcharacteristics->getCrossTalk(type);
165 }
const HcalSiPMCharacteristics * theSiPMcharacteristics
int getType() const
get SiPM type
assert(be >=bs)
float getCrossTalk(int type) const
get cross talk
const HcalDbService * theDbService
const HcalSiPMParameter * getHcalSiPMParameter(const HcalGenericDetId &fId) const
double HcalSimParameters::sipmDarkCurrentuA ( const DetId detId) const

Definition at line 156 of file HcalSimParameters.cc.

References cms::cuda::assert(), HcalSiPMParameter::getDarkCurrent(), HcalDbService::getHcalSiPMParameter(), and theDbService.

Referenced by HcalSiPMHitResponse::addPEnoise().

156  {
159 }
assert(be >=bs)
float getDarkCurrent() const
get dark current
const HcalDbService * theDbService
const HcalSiPMParameter * getHcalSiPMParameter(const HcalGenericDetId &fId) const
std::vector< float > HcalSimParameters::sipmNonlinearity ( const DetId detId) const

Definition at line 166 of file HcalSimParameters.cc.

References cms::cuda::assert(), HcalDbService::getHcalSiPMParameter(), HcalSiPMCharacteristics::getNonLinearities(), HcalSiPMParameter::getType(), theDbService, and theSiPMcharacteristics.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

166  {
170 }
const HcalSiPMCharacteristics * theSiPMcharacteristics
int getType() const
get SiPM type
assert(be >=bs)
std::vector< float > getNonLinearities(int type) const
get nonlinearity constants
const HcalDbService * theDbService
const HcalSiPMParameter * getHcalSiPMParameter(const HcalGenericDetId &fId) const
double HcalSimParameters::sipmTau ( ) const
inline

Definition at line 45 of file HcalSimParameters.h.

References theSiPMTau.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

45 { return theSiPMTau; }
double HcalSimParameters::threshold_currentTDC ( ) const
inline

Definition at line 44 of file HcalSimParameters.h.

References threshold_currentTDC_.

Referenced by HcalElectronicsSim::analogToDigital().

44 { return threshold_currentTDC_; }
double HcalSimParameters::timeSmearRMS ( double  ampl) const

Definition at line 126 of file HcalSimParameters.cc.

References first, mps_fire::i, and theSmearSettings.

Referenced by HcalTimeSlewSim::delay().

126  {
128  double smearsigma = 0;
129 
130  for (i = 0; i < theSmearSettings.size(); i++)
131  if (theSmearSettings[i].first > ampl)
132  break;
133 
134  // Smearing occurs only within the envelope definitions.
135  if (i != 0 && (i < theSmearSettings.size())) {
136  double energy1 = theSmearSettings[i - 1].first;
137  double sigma1 = theSmearSettings[i - 1].second;
138  double energy2 = theSmearSettings[i].first;
139  double sigma2 = theSmearSettings[i].second;
140 
141  if (energy2 != energy1)
142  smearsigma = sigma1 + ((sigma2 - sigma1) * (ampl - energy1) / (energy2 - energy1));
143  else
144  smearsigma = (sigma2 + sigma1) / 2.;
145  }
146 
147  return smearsigma;
148 }
uint16_t size_type
HcalTimeSmearSettings theSmearSettings

Friends And Related Function Documentation

friend class HcalSimParameterMap
friend

Definition at line 51 of file HcalSimParameters.h.

Member Data Documentation

int HcalSimParameters::delayQIE_
private

Definition at line 64 of file HcalSimParameters.h.

Referenced by delayQIE().

bool HcalSimParameters::doTimeSmear_
private

Definition at line 60 of file HcalSimParameters.h.

Referenced by doTimeSmear().

const HcalDbService* HcalSimParameters::theDbService
private
int HcalSimParameters::theFirstRing
private

Definition at line 57 of file HcalSimParameters.h.

Referenced by samplingFactor().

std::vector<double> HcalSimParameters::theSamplingFactors
private

Definition at line 58 of file HcalSimParameters.h.

Referenced by samplingFactor().

const HcalSiPMCharacteristics* HcalSimParameters::theSiPMcharacteristics
private

Definition at line 56 of file HcalSimParameters.h.

Referenced by pixels(), setDbService(), sipmCrossTalk(), and sipmNonlinearity().

bool HcalSimParameters::theSiPMSmearing
private

Definition at line 59 of file HcalSimParameters.h.

Referenced by doSiPMSmearing(), and HcalSimParameters().

double HcalSimParameters::theSiPMTau
private

Definition at line 62 of file HcalSimParameters.h.

Referenced by sipmTau().

HcalTimeSmearSettings HcalSimParameters::theSmearSettings
private

Definition at line 61 of file HcalSimParameters.h.

Referenced by defaultTimeSmearing(), and timeSmearRMS().

double HcalSimParameters::threshold_currentTDC_
private

Definition at line 63 of file HcalSimParameters.h.

Referenced by threshold_currentTDC().