GaussNoiseFP420 Class Reference

#include <GaussNoiseFP420.h>

Inheritance diagram for GaussNoiseFP420:

Public Member Functions
PileUpFP420::signal_map_type	addNoise (const PileUpFP420::signal_map_type &) override
	GaussNoiseFP420 (int ns, float nrms, float th, bool aNpixel, int verbosity)
void	setNumPixels (int in)
void	setThreshold (float in)
Public Member Functions inherited from GNoiseFP420
virtual	~GNoiseFP420 ()

Private Attributes
bool	addNoisyPixels
float	noiseRMS
int	numPixels
float	threshold
int	verbosi

Detailed Description

Definition at line 7 of file GaussNoiseFP420.h.

Constructor & Destructor Documentation

GaussNoiseFP420::GaussNoiseFP420	(	int	ns,
		float	nrms,
		float	th,
		bool	aNpixel,
		int	verbosity
	)

Definition at line 11 of file GaussNoiseFP420.cc.

    : numPixels(ns), noiseRMS(nrms), threshold(th), addNoisyPixels(aNpixel), verbosi(verbosity) {}

Member Function Documentation

PileUpFP420::signal_map_type GaussNoiseFP420::addNoise ( const PileUpFP420::signal_map_type &  in )

overridevirtual

Implements GNoiseFP420.

Definition at line 14 of file GaussNoiseFP420.cc.

References addNoisyPixels, gather_cfg::cout, relval_steps::gen(), GaussNoiseProducerFP420::generate(), hgcalDigitizer_cfi::noise, noiseRMS, numPixels, AlCaHLTBitMon_ParallelJobs::p, threshold, and verbosi.

Referenced by FP420DigiMain::run().

                                                                                         {
  PileUpFP420::signal_map_type _signal;

  // Add noise on non-hit pixels
  std::map<int, float, std::less<int>> generatedNoise;

  //  int numberOfPixels = (numRows * numColumns);// numPixels=numberOfPixels

  GaussNoiseProducerFP420 gen;
  gen.generate(numPixels, threshold, noiseRMS,
               generatedNoise);  // threshold is thePixelThreshold

  // noise for channels with signal:
  // ----------------------------

  for (PileUpFP420::signal_map_type::const_iterator si = in.begin(); si != in.end(); si++) {
    if (verbosi > 0) {
      std::cout << " ***GaussNoiseFP420:  before noise:" << std::endl;
      std::cout << " for si->first=  " << si->first << "    _signal[si->first]=  " << _signal[si->first]
                << "        si->second=      " << si->second << std::endl;
    }
    // define Gaussian noise centered at 0. with sigma = noiseRMS:
    float noise(CLHEP::RandGauss::shoot(0., noiseRMS));
    //    float noise  = CLHEP::RandGaussQ::shoot(0.,theNoiseInElectrons) ;
    // add noise to each channel with signal:
    _signal[si->first] = si->second + noise;

    if (verbosi > 0) {
      std::cout << " ***GaussNoiseFP420: after noise added  = " << noise << std::endl;
      std::cout << "after noise added the _signal[si->first]=  " << _signal[si->first] << std::endl;
    }
  }

  //                                                                                                    //
  if (addNoisyPixels) {  // Option to skip noise in non-hit pixels
    // Noise on the other channels:
    typedef std::map<int, float, std::less<int>>::iterator MI;
    for (MI p = generatedNoise.begin(); p != generatedNoise.end(); p++) {
      if (_signal[(*p).first] == 0) {
        _signal[(*p).first] += (*p).second;
      }
    }  // for(MI
  }

  // or:
  //                                                                        //
  /*
    if(addNoisyPixels){  // Option to skip noise in non-hit pixels
    // Noise on the other channels:
    typedef std::map<int,float,std::less<int> >::iterator MI;
    for(MI p = generatedNoise.begin(); p != generatedNoise.end(); p++){
    int iy = ((*p).first) / numRows;
    int ix = ((*p).first) - (iy*numRows);
    // Keep for a while for testing.
    if( iy < 0 || iy > (numColumns-1) )
    LogWarning ("Pixel Geometry") << " error in iy " << iy ;
    if( ix < 0 || ix > (numRows-1) )
    LogWarning ("Pixel Geometry")  << " error in ix " << ix ;
    int chan = PixelDigi::pixelToChannel(ix, iy);
    LogDebug ("Pixel Digitizer")<<" Storing noise = " << (*mapI).first << " " <<
    (*mapI).second
    << " " << ix << " " << iy << " " << chan ;
    if(_signal[chan] == 0){
    _signal[(*p).first] += (*p).second;
    }//if
    }//for(MI
    }
  */
  //                                                                        //
  //                                                                                                    //

  //
  return _signal;
  //
}

void GaussNoiseFP420::setNumPixels ( int  in )

inline

Definition at line 12 of file GaussNoiseFP420.h.

References recoMuon::in, and numPixels.

{ numPixels = in; }

void GaussNoiseFP420::setThreshold ( float  in )

inline

Definition at line 13 of file GaussNoiseFP420.h.

References recoMuon::in, and threshold.

{ threshold = in; }

Member Data Documentation

bool GaussNoiseFP420::addNoisyPixels

private

Definition at line 19 of file GaussNoiseFP420.h.

Referenced by addNoise().

float GaussNoiseFP420::noiseRMS

private

Definition at line 17 of file GaussNoiseFP420.h.

Referenced by addNoise().

int GaussNoiseFP420::numPixels

private

Definition at line 16 of file GaussNoiseFP420.h.

Referenced by addNoise(), and setNumPixels().

float GaussNoiseFP420::threshold

private

Definition at line 18 of file GaussNoiseFP420.h.

Referenced by addNoise(), utils.StatisticalTest::get_status(), and setThreshold().

int GaussNoiseFP420::verbosi

private

Definition at line 20 of file GaussNoiseFP420.h.

Referenced by addNoise().