GaussNoiseFP420.cc

Go to the documentation of this file.

// File: GaussNoiseFP420.cc
// Date: 12.2006
// Description: GaussNoiseFP420 for FP420
// Modifications: 
#include "SimRomanPot/SimFP420/interface/GaussNoiseFP420.h"
#include "SimRomanPot/SimFP420/interface/GaussNoiseProducerFP420.h"
#include "CLHEP/Random/RandGauss.h"

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


PileUpFP420::signal_map_type 
GaussNoiseFP420::addNoise(const PileUpFP420::signal_map_type& in){
  
  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;
  //
}
//