RPGaussianTailNoiseAdder Class Reference

#include <RPGaussianTailNoiseAdder.h>

Public Member Functions
simromanpot::strip_charge_map	addNoise (const simromanpot::strip_charge_map &theSignal)
	RPGaussianTailNoiseAdder (int numStrips, double theNoiseInElectrons, double theStripThresholdInE, CLHEP::HepRandomEngine &eng, int verbosity)

Private Attributes
int	numStrips_
CLHEP::HepRandomEngine &	rndEngine_
double	strips_above_threshold_prob_
double	theNoiseInElectrons
double	theStripThresholdInE
int	verbosity_

Detailed Description

Definition at line 7 of file RPGaussianTailNoiseAdder.h.

Constructor & Destructor Documentation

RPGaussianTailNoiseAdder::RPGaussianTailNoiseAdder	(	int	numStrips,
		double	theNoiseInElectrons,
		double	theStripThresholdInE,
		CLHEP::HepRandomEngine &	eng,
		int	verbosity
	)

Definition at line 12 of file RPGaussianTailNoiseAdder.cc.

References mathSSE::sqrt(), and strips_above_threshold_prob_.

    : numStrips_(numStrips),
      theNoiseInElectrons(theNoiseInElectrons),
      theStripThresholdInE(theStripThresholdInE),
      rndEngine_(eng),
      verbosity_(verbosity) {
  strips_above_threshold_prob_ = std::erfc(theStripThresholdInE / sqrt(2.0) / theNoiseInElectrons) / 2;
}

Member Function Documentation

simromanpot::strip_charge_map RPGaussianTailNoiseAdder::addNoise

(

const simromanpot::strip_charge_map &

theSignal

)

Definition at line 22 of file RPGaussianTailNoiseAdder.cc.

References mps_fire::i, dqmiolumiharvest::j, gpuVertexFinder::noise, numStrips_, rndEngine_, digitizers_cfi::strip, strips_above_threshold_prob_, theNoiseInElectrons, theStripThresholdInE, and verbosity_.

                                                                                                           {
  simromanpot::strip_charge_map the_strip_charge_map;

  // noise on strips with signal:
  for (simromanpot::strip_charge_map::const_iterator i = theSignal.begin(); i != theSignal.end(); ++i) {
    double noise = CLHEP::RandGauss::shoot(&(rndEngine_), 0.0, theNoiseInElectrons);
    the_strip_charge_map[i->first] = i->second + noise;
    if (verbosity_)
      edm::LogInfo("RPDigiProducer") << "noise added to signal strips: " << noise << "\n";
  }

  // noise on the other strips
  int strips_no_above_threshold =
      CLHEP::RandBinomial::shoot(&(rndEngine_), (long)numStrips_, strips_above_threshold_prob_);

  for (int j = 0; j < strips_no_above_threshold; j++) {
    int strip = CLHEP::RandFlat::shootInt(&(rndEngine_), numStrips_);
    if (the_strip_charge_map[strip] == 0) {
      the_strip_charge_map[strip] = 2 * theStripThresholdInE;
      //only binary decision later, no need to simulate the noise precisely,
      //enough to know that it exceeds the threshold
      if (verbosity_)
        edm::LogInfo("RPDigiProducer") << "nonsignal strips noise :" << strip << " " << the_strip_charge_map[strip]
                                       << "\n";
    }
  }

  return the_strip_charge_map;
}

Member Data Documentation

int RPGaussianTailNoiseAdder::numStrips_

private

Definition at line 17 of file RPGaussianTailNoiseAdder.h.

Referenced by addNoise().

CLHEP::HepRandomEngine& RPGaussianTailNoiseAdder::rndEngine_

private

Definition at line 21 of file RPGaussianTailNoiseAdder.h.

Referenced by addNoise().

double RPGaussianTailNoiseAdder::strips_above_threshold_prob_

private

Definition at line 20 of file RPGaussianTailNoiseAdder.h.

Referenced by addNoise(), and RPGaussianTailNoiseAdder().

double RPGaussianTailNoiseAdder::theNoiseInElectrons

private

Definition at line 18 of file RPGaussianTailNoiseAdder.h.

Referenced by addNoise().

double RPGaussianTailNoiseAdder::theStripThresholdInE

private

Definition at line 19 of file RPGaussianTailNoiseAdder.h.

Referenced by addNoise().

int RPGaussianTailNoiseAdder::verbosity_

private

Definition at line 22 of file RPGaussianTailNoiseAdder.h.

Referenced by addNoise().