CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
GEMSimpleModel Class Reference

#include <GEMSimpleModel.h>

Inheritance diagram for GEMSimpleModel:
GEMDigiModel

Public Member Functions

 GEMSimpleModel (const edm::ParameterSet &)
 
int getSimHitBx (const PSimHit *, CLHEP::HepRandomEngine *)
 
void setup () override
 
std::vector< std::pair< int, int > > simulateClustering (const GEMEtaPartition *, const PSimHit *, const int, CLHEP::HepRandomEngine *) override
 
void simulateNoise (const GEMEtaPartition *, CLHEP::HepRandomEngine *) override
 
void simulateSignal (const GEMEtaPartition *, const edm::PSimHitContainer &, CLHEP::HepRandomEngine *) override
 
 ~GEMSimpleModel () override
 
- Public Member Functions inherited from GEMDigiModel
void fillDigis (int rollDetId, GEMDigiCollection &)
 
const GEMDigiSimLinksgemDigiSimLinks () const
 
const GEMGeometrygetGeometry ()
 
void setGeometry (const GEMGeometry *geom)
 
const StripDigiSimLinksstripDigiSimLinks () const
 
virtual ~GEMDigiModel ()
 

Private Attributes

double averageEfficiency_
 
double averageNoiseRate_
 
double averageShapingTime_
 
int bxwidth_
 
bool digitizeOnlyMuons_
 
bool doBkgNoise_
 
bool doNoiseCLS_
 
bool fixedRollRadius_
 
double GE11ElecBkgParam0_
 
double GE11ElecBkgParam1_
 
double GE11ElecBkgParam2_
 
double GE11ModNeuBkgParam0_
 
double GE11ModNeuBkgParam1_
 
double GE11ModNeuBkgParam2_
 
double GE21ElecBkgParam0_
 
double GE21ElecBkgParam1_
 
double GE21ElecBkgParam2_
 
double GE21ModNeuBkgParam0_
 
double GE21ModNeuBkgParam1_
 
double GE21ModNeuBkgParam2_
 
double instLumi_
 
int maxBunch_
 
int minBunch_
 
double rateFact_
 
const double referenceInstLumi_
 
double resolutionX_
 
double signalPropagationSpeed_
 
bool simulateElectronBkg_
 
bool simulateIntrinsicNoise_
 
double timeJitter_
 
double timeResolution_
 

Additional Inherited Members

- Public Types inherited from GEMDigiModel
typedef edm::DetSet< GEMDigiSimLinkGEMDigiSimLinks
 
typedef edm::DetSet< StripDigiSimLinkStripDigiSimLinks
 
- Protected Types inherited from GEMDigiModel
typedef std::multimap< std::pair< unsigned int, int >, const PSimHit *, std::less< std::pair< unsigned int, int > > > DetectorHitMap
 
- Protected Member Functions inherited from GEMDigiModel
void addLinks (unsigned int strip, int bx)
 creates links from Digi to SimTrack More...
 
void addLinksWithPartId (unsigned int strip, int bx)
 
 GEMDigiModel (const edm::ParameterSet &)
 
- Protected Attributes inherited from GEMDigiModel
DetectorHitMap detectorHitMap_
 
const GEMGeometrygeometry_
 
StripDigiSimLinks stripDigiSimLinks_
 
std::set< std::pair< int, int > > strips_
 
GEMDigiSimLinks theGemDigiSimLinks_
 

Detailed Description

Class for the GEM strip response simulation based on a very simple model

Author
Sven Dildick by Roumyana Hadjiiska

Definition at line 22 of file GEMSimpleModel.h.

Constructor & Destructor Documentation

GEMSimpleModel::GEMSimpleModel ( const edm::ParameterSet config)

Definition at line 13 of file GEMSimpleModel.cc.

13  :
14 GEMDigiModel(config)
15 , averageEfficiency_(config.getParameter<double> ("averageEfficiency"))
16 , averageShapingTime_(config.getParameter<double> ("averageShapingTime"))
17 , timeResolution_(config.getParameter<double> ("timeResolution"))
18 , timeJitter_(config.getParameter<double> ("timeJitter"))
19 , averageNoiseRate_(config.getParameter<double> ("averageNoiseRate"))
20 , signalPropagationSpeed_(config.getParameter<double> ("signalPropagationSpeed"))
21 , bxwidth_(config.getParameter<int> ("bxwidth"))
22 , minBunch_(config.getParameter<int> ("minBunch"))
23 , maxBunch_(config.getParameter<int> ("maxBunch"))
24 , digitizeOnlyMuons_(config.getParameter<bool> ("digitizeOnlyMuons"))
25 , doBkgNoise_(config.getParameter<bool> ("doBkgNoise"))
26 , doNoiseCLS_(config.getParameter<bool> ("doNoiseCLS"))
27 , fixedRollRadius_(config.getParameter<bool> ("fixedRollRadius"))
28 , simulateElectronBkg_(config.getParameter<bool> ("simulateElectronBkg"))
29 , instLumi_(config.getParameter<double>("instLumi"))
30 , rateFact_(config.getParameter<double>("rateFact"))
31 , referenceInstLumi_(config.getParameter<double>("referenceInstLumi"))
32 , resolutionX_(config.getParameter<double>("resolutionX"))
33 , GE11ElecBkgParam0_(config.getParameter<double>("GE11ElecBkgParam0"))
34 , GE11ElecBkgParam1_(config.getParameter<double>("GE11ElecBkgParam1"))
35 , GE11ElecBkgParam2_(config.getParameter<double>("GE11ElecBkgParam2"))
36 , GE21ElecBkgParam0_(config.getParameter<double>("GE21ElecBkgParam0"))
37 , GE21ElecBkgParam1_(config.getParameter<double>("GE21ElecBkgParam1"))
38 , GE21ElecBkgParam2_(config.getParameter<double>("GE21ElecBkgParam2"))
39 , GE11ModNeuBkgParam0_(config.getParameter<double>("GE11ModNeuBkgParam0"))
40 , GE11ModNeuBkgParam1_(config.getParameter<double>("GE11ModNeuBkgParam1"))
41 , GE11ModNeuBkgParam2_(config.getParameter<double>("GE11ModNeuBkgParam2"))
42 , GE21ModNeuBkgParam0_(config.getParameter<double>("GE11ModNeuBkgParam0"))
43 , GE21ModNeuBkgParam1_(config.getParameter<double>("GE11ModNeuBkgParam1"))
44 , GE21ModNeuBkgParam2_(config.getParameter<double>("GE11ModNeuBkgParam2"))
45 {
46 }
T getParameter(std::string const &) const
double signalPropagationSpeed_
double GE21ElecBkgParam1_
double GE11ModNeuBkgParam2_
double GE11ElecBkgParam2_
double averageNoiseRate_
double GE11ElecBkgParam0_
double GE21ElecBkgParam2_
double GE21ModNeuBkgParam1_
double averageEfficiency_
double GE21ModNeuBkgParam2_
double GE11ModNeuBkgParam1_
double GE21ModNeuBkgParam0_
double GE11ModNeuBkgParam0_
double GE21ElecBkgParam0_
const double referenceInstLumi_
GEMDigiModel(const edm::ParameterSet &)
Definition: GEMDigiModel.h:61
bool simulateElectronBkg_
double GE11ElecBkgParam1_
double averageShapingTime_
double timeResolution_
GEMSimpleModel::~GEMSimpleModel ( )
override

Definition at line 48 of file GEMSimpleModel.cc.

49 {
50 }

Member Function Documentation

int GEMSimpleModel::getSimHitBx ( const PSimHit simhit,
CLHEP::HepRandomEngine *  engine 
)

Definition at line 100 of file GEMSimpleModel.cc.

References averageShapingTime_, bxwidth_, gather_cfg::cout, debug, PSimHit::detUnitId(), MillePedeFileConverter_cfg::e, GEMGeometry::etaPartition(), Exception, GEMDigiModel::geometry_, triggerObjects_cff::id, PSimHit::localPosition(), signalPropagationSpeed_, mathSSE::sqrt(), timeJitter_, PSimHit::timeOfFlight(), timeResolution_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by simulateSignal().

101 {
102  int bx = -999;
103  const LocalPoint simHitPos(simhit->localPosition());
104  const float tof(simhit->timeOfFlight());
105  // random Gaussian time correction due to electronics jitter
106  float randomJitterTime = CLHEP::RandGaussQ::shoot(engine, 0., timeJitter_);
107  const GEMDetId id(simhit->detUnitId());
108  const GEMEtaPartition* roll(geometry_->etaPartition(id));
109  if (!roll)
110  {
111  throw cms::Exception("Geometry")<< "GEMSimpleModel::getSimHitBx() - GEM simhit id does not match any GEM roll id: " << id << "\n";
112  return 999;
113  }
114  if (roll->id().region() == 0)
115  {
116  throw cms::Exception("Geometry") << "GEMSimpleModel::getSimHitBx() - this GEM id is from barrel, which cannot happen: " << roll->id() << "\n";
117  }
118  const double cspeed = 299792458; // signal propagation speed in vacuum in [m/s]
119  const int nstrips = roll->nstrips();
120  float middleStrip = nstrips/2.;
121  const LocalPoint& middleOfRoll = roll->centreOfStrip(middleStrip);
122  const GlobalPoint& globMiddleRol = roll->toGlobal(middleOfRoll);
123  double muRadius = sqrt(globMiddleRol.x()*globMiddleRol.x() + globMiddleRol.y()*globMiddleRol.y() +globMiddleRol.z()*globMiddleRol.z());
124  double timeCalibrationOffset_ = (muRadius *1e+9)/(cspeed*1e+2); //[ns]
125 
126  const TrapezoidalStripTopology* top(dynamic_cast<const TrapezoidalStripTopology*> (&(roll->topology())));
127  const float halfStripLength(0.5 * top->stripLength());
128  const float distanceFromEdge(halfStripLength - simHitPos.y());
129 
130  // signal propagation speed in material in [cm/ns]
131  double signalPropagationSpeedTrue = signalPropagationSpeed_ * cspeed * 1e-7; // 1e+2 * 1e-9;
132 
133  // average time for the signal to propagate from the SimHit to the top of a strip
134  const float averagePropagationTime(distanceFromEdge / signalPropagationSpeedTrue);
135  // random Gaussian time correction due to the finite timing resolution of the detector
136  float randomResolutionTime = CLHEP::RandGaussQ::shoot(engine, 0., timeResolution_);
137  const float simhitTime(tof + averageShapingTime_ + randomResolutionTime + averagePropagationTime + randomJitterTime);
138  float referenceTime = 0.;
139  referenceTime = timeCalibrationOffset_ + halfStripLength / signalPropagationSpeedTrue + averageShapingTime_;
140  const float timeDifference(simhitTime - referenceTime);
141  // assign the bunch crossing
142  bx = static_cast<int> (std::round((timeDifference) / bxwidth_));
143 
144  // check time
145  const bool debug(false);
146  if (debug)
147  {
148  std::cout << "checktime " << "bx = " << bx << "\tdeltaT = " << timeDifference << "\tsimT = " << simhitTime
149  << "\trefT = " << referenceTime << "\ttof = " << tof << "\tavePropT = " << averagePropagationTime
150  << "\taveRefPropT = " << halfStripLength / signalPropagationSpeedTrue << std::endl;
151  }
152  return bx;
153 }
double signalPropagationSpeed_
T y() const
Definition: PV3DBase.h:63
const GEMEtaPartition * etaPartition(GEMDetId id) const
Return a GEMEtaPartition given its id.
Definition: GEMGeometry.cc:99
float timeOfFlight() const
Definition: PSimHit.h:69
Local3DPoint localPosition() const
Definition: PSimHit.h:44
T sqrt(T t)
Definition: SSEVec.h:18
T z() const
Definition: PV3DBase.h:64
#define debug
Definition: HDRShower.cc:19
const GEMGeometry * geometry_
Definition: GEMDigiModel.h:63
double averageShapingTime_
double timeResolution_
T x() const
Definition: PV3DBase.h:62
unsigned int detUnitId() const
Definition: PSimHit.h:93
void GEMSimpleModel::setup ( )
overridevirtual

Implements GEMDigiModel.

Definition at line 52 of file GEMSimpleModel.cc.

53 {
54  return;
55 }
std::vector< std::pair< int, int > > GEMSimpleModel::simulateClustering ( const GEMEtaPartition roll,
const PSimHit simHit,
const int  bx,
CLHEP::HepRandomEngine *  engine 
)
overridevirtual

Implements GEMDigiModel.

Definition at line 270 of file GEMSimpleModel.cc.

References funct::abs(), PSimHit::entryPoint(), PSimHit::exitPoint(), mps_fire::i, resolutionX_, GEMEtaPartition::strip(), and PV3DBase< T, PVType, FrameType >::x().

Referenced by simulateSignal().

273 {
274  const LocalPoint& hit_entry(simHit->entryPoint());
275  const LocalPoint& hit_exit(simHit->exitPoint());
276 
277  float hit_entry_smeardX;
278  float hit_exit_smeardX;
279  if (hit_entry.x()>hit_exit.x()) {
280  hit_entry_smeardX = hit_entry.x()+std::abs(CLHEP::RandGaussQ::shoot(engine, 0, resolutionX_));
281  hit_exit_smeardX = hit_exit.x()-std::abs(CLHEP::RandGaussQ::shoot(engine, 0, resolutionX_));
282  }
283  else {
284  hit_entry_smeardX = hit_entry.x()-std::abs(CLHEP::RandGaussQ::shoot(engine, 0, resolutionX_));
285  hit_exit_smeardX = hit_exit.x()+std::abs(CLHEP::RandGaussQ::shoot(engine, 0, resolutionX_));
286  }
287 
288  LocalPoint inPoint(hit_entry_smeardX, hit_entry.y(), hit_entry.z());
289  LocalPoint outPoint(hit_exit_smeardX, hit_exit.y(), hit_exit.z());
290 
291  int clusterStart = roll->strip(inPoint);
292  int clusterEnd = roll->strip(outPoint);
293 
294  std::vector < std::pair<int, int> > cluster_;
295  cluster_.clear();
296  for (int i = clusterStart; i<= clusterEnd ; i++) {
297  cluster_.emplace_back(i, bx);
298  }
299 
300  return cluster_;
301 }
Local3DPoint exitPoint() const
Exit point in the local Det frame.
Definition: PSimHit.h:38
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
float strip(const LocalPoint &lp) const
T x() const
Definition: PV3DBase.h:62
Local3DPoint entryPoint() const
Entry point in the local Det frame.
Definition: PSimHit.h:35
void GEMSimpleModel::simulateNoise ( const GEMEtaPartition roll,
CLHEP::HepRandomEngine *  engine 
)
overridevirtual

Implements GEMDigiModel.

Definition at line 155 of file GEMSimpleModel.cc.

References averageEfficiency_, averageNoiseRate_, bxwidth_, doBkgNoise_, doNoiseCLS_, MillePedeFileConverter_cfg::e, Exception, fixedRollRadius_, GE11ElecBkgParam0_, GE11ElecBkgParam1_, GE11ElecBkgParam2_, GE11ModNeuBkgParam0_, GE11ModNeuBkgParam1_, GE11ModNeuBkgParam2_, GE21ElecBkgParam0_, GE21ElecBkgParam1_, GE21ElecBkgParam2_, GE21ModNeuBkgParam0_, GE21ModNeuBkgParam1_, GE21ModNeuBkgParam2_, mps_fire::i, GEMEtaPartition::id(), instLumi_, gen::k, maxBunch_, minBunch_, GEMEtaPartition::nstrips(), GEMEtaPartition::pitch(), rateFact_, referenceInstLumi_, simulateElectronBkg_, simulateIntrinsicNoise_, GEMDigiModel::strips_, and GEMEtaPartition::topology().

156 {
157  if (!doBkgNoise_)
158  return;
159  const GEMDetId& gemId(roll->id());
160  const int nstrips(roll->nstrips());
161  double trArea(0.0);
162  double trStripArea(0.0);
163  if (gemId.region() == 0)
164  {
165  throw cms::Exception("Geometry")
166  << "GEMSynchronizer::simulateNoise() - this GEM id is from barrel, which cannot happen.";
167  }
168  const TrapezoidalStripTopology* top_(dynamic_cast<const TrapezoidalStripTopology*>(&(roll->topology())));
169  const float striplength(top_->stripLength());
170  trStripArea = (roll->pitch()) * striplength;
171  trArea = trStripArea * nstrips;
172  const int nBxing(maxBunch_ - minBunch_ + 1);
173  const float rollRadius(fixedRollRadius_ ? top_->radius() :
174  top_->radius() + CLHEP::RandFlat::shoot(engine, -1.*top_->stripLength()/2., top_->stripLength()/2.));
175 
176  //calculate noise from model
177  double averageNeutralNoiseRatePerRoll = 0.;
178  double averageNoiseElectronRatePerRoll = 0.;
179  double averageNoiseRatePerRoll = 0.;
180  if (gemId.station() == 1)
181  {
182  //simulate neutral background for GE1/1
183  averageNeutralNoiseRatePerRoll = (GE11ModNeuBkgParam0_
184  + GE11ModNeuBkgParam1_ * rollRadius
185  + GE11ModNeuBkgParam2_ * rollRadius * rollRadius); //simulate electron background for GE1/1
187  averageNoiseElectronRatePerRoll = (GE11ElecBkgParam0_
188  + GE11ElecBkgParam1_ * rollRadius
189  + GE11ElecBkgParam2_ * rollRadius * rollRadius);
190 
191  // Scale up/down for desired instantaneous lumi (reference is 5E34, double from config is in units of 1E34)
192  averageNoiseRatePerRoll = averageNeutralNoiseRatePerRoll + averageNoiseElectronRatePerRoll;
193  averageNoiseRatePerRoll *= instLumi_*rateFact_*1.0/referenceInstLumi_;
194  }
195  if (gemId.station() == 2)
196  {
197  //simulate neutral background for GE2/1
198  averageNeutralNoiseRatePerRoll = (GE21ModNeuBkgParam0_
199  + GE21ModNeuBkgParam1_ * rollRadius
200  + GE21ModNeuBkgParam2_ * rollRadius * rollRadius);
201  //simulate electron background for GE2/1
203  averageNoiseElectronRatePerRoll = (GE21ElecBkgParam0_
204  + GE21ElecBkgParam1_ * rollRadius
205  + GE21ElecBkgParam2_ * rollRadius * rollRadius);
206 
207  // Scale up/down for desired instantaneous lumi (reference is 5E34, double from config is in units of 1E34)
208  averageNoiseRatePerRoll = averageNeutralNoiseRatePerRoll + averageNoiseElectronRatePerRoll;
209  averageNoiseRatePerRoll *= instLumi_*rateFact_*1.0/referenceInstLumi_;
210  }
211 
212  //simulate intrinsic noise
214  {
215  const double aveIntrinsicNoisePerStrip(averageNoiseRate_ * nBxing * bxwidth_ * trStripArea * 1.0e-9);
216  for(int j = 0; j < nstrips; ++j)
217  {
218  CLHEP::RandPoissonQ randPoissonQ(*engine, aveIntrinsicNoisePerStrip);
219  const int n_intrHits(randPoissonQ.fire());
220 
221  for (int k = 0; k < n_intrHits; k++ )
222  {
223  const int time_hit(static_cast<int>(CLHEP::RandFlat::shoot(engine, nBxing)) + minBunch_);
224  std::pair<int, int> digi(k+1,time_hit);
225  strips_.emplace(digi);
226  }
227  }
228  }//end simulate intrinsic noise
229 
230  //simulate bkg contribution
231  const double averageNoise(averageNoiseRatePerRoll * nBxing * bxwidth_ * trArea * 1.0e-9);
232  CLHEP::RandPoissonQ randPoissonQ(*engine, averageNoise);
233  const int n_hits(randPoissonQ.fire());
234  for (int i = 0; i < n_hits; ++i)
235  {
236  const int centralStrip(static_cast<int> (CLHEP::RandFlat::shoot(engine, 1, nstrips)));
237  const int time_hit(static_cast<int>(CLHEP::RandFlat::shoot(engine, nBxing)) + minBunch_);
238  if (doNoiseCLS_)
239  {
240  std::vector < std::pair<int, int> > cluster_;
241  cluster_.clear();
242  cluster_.emplace_back(centralStrip, time_hit);
243  int clusterSize((CLHEP::RandFlat::shoot(engine)) <= 0.53 ? 1 : 2);
244  if (clusterSize == 2)
245  {
246  if(CLHEP::RandFlat::shoot(engine) < 0.5)
247  {
248  if (CLHEP::RandFlat::shoot(engine) < averageEfficiency_ && (centralStrip - 1 > 0))
249  cluster_.emplace_back(centralStrip - 1, time_hit);
250  }
251  else
252  {
253  if (CLHEP::RandFlat::shoot(engine) < averageEfficiency_ && (centralStrip + 1 <= nstrips))
254  cluster_.emplace_back(centralStrip + 1, time_hit);
255  }
256  }
257  for (const auto& digi : cluster_)
258  {
259  strips_.emplace(digi);
260  }
261  } //end doNoiseCLS_
262  else
263  {
264  strips_.emplace(centralStrip, time_hit);
265  }
266  }
267  return;
268 }
double GE21ElecBkgParam1_
double GE11ModNeuBkgParam2_
double GE11ElecBkgParam2_
double averageNoiseRate_
double GE11ElecBkgParam0_
const Topology & topology() const override
GEMDetId id() const
double GE21ElecBkgParam2_
int nstrips() const
number of readout strips in partition
bool simulateIntrinsicNoise_
double GE21ModNeuBkgParam1_
double averageEfficiency_
double GE21ModNeuBkgParam2_
double GE11ModNeuBkgParam1_
double GE21ModNeuBkgParam0_
double GE11ModNeuBkgParam0_
float pitch() const
double GE21ElecBkgParam0_
int k[5][pyjets_maxn]
const double referenceInstLumi_
bool simulateElectronBkg_
double GE11ElecBkgParam1_
std::set< std::pair< int, int > > strips_
Definition: GEMDigiModel.h:65
void GEMSimpleModel::simulateSignal ( const GEMEtaPartition roll,
const edm::PSimHitContainer simHits,
CLHEP::HepRandomEngine *  engine 
)
overridevirtual

Implements GEMDigiModel.

Definition at line 57 of file GEMSimpleModel.cc.

References funct::abs(), averageEfficiency_, edm::DetSet< T >::clear(), GEMDigiModel::detectorHitMap_, digitizeOnlyMuons_, MillePedeFileConverter_cfg::e, JetChargeProducer_cfi::exp, getSimHitBx(), GEMEtaPartition::id(), cmsBatch::log, DetId::rawId(), simulateClustering(), GEMDigiModel::stripDigiSimLinks_, GEMDigiModel::strips_, and GEMDigiModel::theGemDigiSimLinks_.

58 {
60  detectorHitMap_.clear();
64  bool digiMuon = false;
65  bool digiElec = false;
66  for (const auto& hit : simHits)
67  {
68  if (std::abs(hit.particleType()) != 13 && digitizeOnlyMuons_)
69  continue;
70  double elecEff = 0.;
71  double partMom = hit.pabs();
72  double checkMuonEff = CLHEP::RandFlat::shoot(engine, 0., 1.);
73  double checkElecEff = CLHEP::RandFlat::shoot(engine, 0., 1.);
74  if (std::abs(hit.particleType()) == 13 && checkMuonEff < averageEfficiency_)
75  digiMuon = true;
76  if (std::abs(hit.particleType()) != 13) //consider all non muon particles with gem efficiency to electrons
77  {
78  if (partMom <= 1.95e-03)
79  elecEff = 1.7e-05 * std::exp(2.1 * partMom * 1000.);
80  if (partMom > 1.95e-03 && partMom < 10.e-03)
81  elecEff = 1.34 * log(7.96e-01 * partMom * 1000. - 5.75e-01)
82  / (1.34 + log(7.96e-01 * partMom * 1000. - 5.75e-01));
83  if (partMom > 10.e-03)
84  elecEff = 1.;
85  if (checkElecEff < elecEff)
86  digiElec = true;
87  }
88  if (!(digiMuon || digiElec))
89  continue;
90  const int bx(getSimHitBx(&hit, engine));
91  const std::vector<std::pair<int, int> >& cluster(simulateClustering(roll, &hit, bx, engine));
92  for (const auto & digi : cluster)
93  {
94  detectorHitMap_.emplace(digi,&hit);
95  strips_.emplace(digi);
96  }
97  }
98 }
edm::DetSet< GEMDigiSimLink > GEMDigiSimLinks
Definition: GEMDigiModel.h:37
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:47
int getSimHitBx(const PSimHit *, CLHEP::HepRandomEngine *)
GEMDigiSimLinks theGemDigiSimLinks_
Definition: GEMDigiModel.h:80
GEMDetId id() const
double averageEfficiency_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::vector< std::pair< int, int > > simulateClustering(const GEMEtaPartition *, const PSimHit *, const int, CLHEP::HepRandomEngine *) override
std::set< std::pair< int, int > > strips_
Definition: GEMDigiModel.h:65
void clear()
Definition: DetSet.h:69
DetectorHitMap detectorHitMap_
Definition: GEMDigiModel.h:78
StripDigiSimLinks stripDigiSimLinks_
Definition: GEMDigiModel.h:79
edm::DetSet< StripDigiSimLink > StripDigiSimLinks
Definition: GEMDigiModel.h:36

Member Data Documentation

double GEMSimpleModel::averageEfficiency_
private

Definition at line 43 of file GEMSimpleModel.h.

Referenced by simulateNoise(), and simulateSignal().

double GEMSimpleModel::averageNoiseRate_
private

Definition at line 47 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::averageShapingTime_
private

Definition at line 44 of file GEMSimpleModel.h.

Referenced by getSimHitBx().

int GEMSimpleModel::bxwidth_
private

Definition at line 49 of file GEMSimpleModel.h.

Referenced by getSimHitBx(), and simulateNoise().

bool GEMSimpleModel::digitizeOnlyMuons_
private

Definition at line 52 of file GEMSimpleModel.h.

Referenced by simulateSignal().

bool GEMSimpleModel::doBkgNoise_
private

Definition at line 53 of file GEMSimpleModel.h.

Referenced by simulateNoise().

bool GEMSimpleModel::doNoiseCLS_
private

Definition at line 54 of file GEMSimpleModel.h.

Referenced by simulateNoise().

bool GEMSimpleModel::fixedRollRadius_
private

Definition at line 55 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE11ElecBkgParam0_
private

Definition at line 64 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE11ElecBkgParam1_
private

Definition at line 65 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE11ElecBkgParam2_
private

Definition at line 66 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE11ModNeuBkgParam0_
private

Definition at line 71 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE11ModNeuBkgParam1_
private

Definition at line 72 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE11ModNeuBkgParam2_
private

Definition at line 73 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE21ElecBkgParam0_
private

Definition at line 67 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE21ElecBkgParam1_
private

Definition at line 68 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE21ElecBkgParam2_
private

Definition at line 69 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE21ModNeuBkgParam0_
private

Definition at line 74 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE21ModNeuBkgParam1_
private

Definition at line 75 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::GE21ModNeuBkgParam2_
private

Definition at line 76 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::instLumi_
private

Definition at line 58 of file GEMSimpleModel.h.

Referenced by simulateNoise().

int GEMSimpleModel::maxBunch_
private

Definition at line 51 of file GEMSimpleModel.h.

Referenced by simulateNoise().

int GEMSimpleModel::minBunch_
private

Definition at line 50 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::rateFact_
private

Definition at line 59 of file GEMSimpleModel.h.

Referenced by simulateNoise().

const double GEMSimpleModel::referenceInstLumi_
private

Definition at line 60 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::resolutionX_
private

Definition at line 61 of file GEMSimpleModel.h.

Referenced by simulateClustering().

double GEMSimpleModel::signalPropagationSpeed_
private

Definition at line 48 of file GEMSimpleModel.h.

Referenced by getSimHitBx().

bool GEMSimpleModel::simulateElectronBkg_
private

Definition at line 57 of file GEMSimpleModel.h.

Referenced by simulateNoise().

bool GEMSimpleModel::simulateIntrinsicNoise_
private

Definition at line 56 of file GEMSimpleModel.h.

Referenced by simulateNoise().

double GEMSimpleModel::timeJitter_
private

Definition at line 46 of file GEMSimpleModel.h.

Referenced by getSimHitBx().

double GEMSimpleModel::timeResolution_
private

Definition at line 45 of file GEMSimpleModel.h.

Referenced by getSimHitBx().