SiGaussianTailNoiseAdder Class Reference

#include <SiGaussianTailNoiseAdder.h>

Inheritance diagram for SiGaussianTailNoiseAdder:

Public Member Functions
void	addBaselineShift (std::vector< float > &, std::vector< bool > &) const override
void	addCMNoise (std::vector< float > &, float, std::vector< bool > &, CLHEP::HepRandomEngine *) const override
void	addNoise (std::vector< float > &, size_t &, size_t &, int, float, CLHEP::HepRandomEngine *) const override
void	addNoiseVR (std::vector< float > &, std::vector< float > &, CLHEP::HepRandomEngine *) const override
void	addPedestals (std::vector< float > &, std::vector< float > &) const override
	SiGaussianTailNoiseAdder (float)
	~SiGaussianTailNoiseAdder ()
Public Member Functions inherited from SiNoiseAdder
virtual	~SiNoiseAdder ()

Private Attributes
std::unique_ptr < GaussianTailNoiseGenerator >	genNoise
const float	threshold

Detailed Description

Definition at line 17 of file SiGaussianTailNoiseAdder.h.

Constructor & Destructor Documentation

SiGaussianTailNoiseAdder::SiGaussianTailNoiseAdder

(

float

th

)

Definition at line 4 of file SiGaussianTailNoiseAdder.cc.

                                                          :
  threshold(th),
  genNoise(new GaussianTailNoiseGenerator())
{
}

SiGaussianTailNoiseAdder::~SiGaussianTailNoiseAdder

(

)

Definition at line 10 of file SiGaussianTailNoiseAdder.cc.

                                                   {
}

Member Function Documentation

void SiGaussianTailNoiseAdder::addBaselineShift ( std::vector< float > &  in, std::vector< bool > &  badChannels  ) const

overridevirtual

Implements SiNoiseAdder.

Definition at line 60 of file SiGaussianTailNoiseAdder.cc.

Referenced by DigiSimLinkAlgorithm::run().

                                                                                                        {
  size_t nAPVs = in.size()/128;
  std::vector<float> vShift;
  double apvCharge, apvMult;
  
  size_t iChannel;
  for(size_t APVn =0; APVn < nAPVs; ++APVn){
    apvMult=0; apvCharge=0;
    for(iChannel=APVn*128; iChannel!=APVn*128+128; ++iChannel) {
       if(in[iChannel]>0){
         ++apvMult;
          apvCharge+= in[iChannel];
       }
       if(apvMult==0) vShift.push_back(0);
       else vShift.push_back(apvCharge/apvMult);
    }
  }
     
  for (iChannel=0; iChannel!=in.size(); ++iChannel) {
     if(!badChannels[iChannel]) in[iChannel] -= vShift[(int)(iChannel/128)];
  }
}

void SiGaussianTailNoiseAdder::addCMNoise ( std::vector< float > &  in, float  cmnRMS, std::vector< bool > &  badChannels, CLHEP::HepRandomEngine *  engine  ) const

overridevirtual

Implements SiNoiseAdder.

Definition at line 51 of file SiGaussianTailNoiseAdder.cc.

Referenced by DigiSimLinkAlgorithm::run().

                                                                                                                                              {
  int nAPVs = in.size()/128;
  std::vector<float> CMNv;
  for(int APVn =0; APVn < nAPVs; ++APVn) CMNv.push_back(CLHEP::RandGaussQ::shoot(engine, 0., cmnRMS));
  for (size_t iChannel=0; iChannel!=in.size(); iChannel++) {
     if(!badChannels[iChannel]) in[iChannel] += CMNv[(int)(iChannel/128)];
  }
}

void SiGaussianTailNoiseAdder::addNoise ( std::vector< float > &  in, size_t &  minChannel, size_t &  maxChannel, int  numStrips, float  noiseRMS, CLHEP::HepRandomEngine *  engine  ) const

overridevirtual

Implements SiNoiseAdder.

Definition at line 13 of file SiGaussianTailNoiseAdder.cc.

References genNoise, AlCaHLTBitMon_ParallelJobs::p, and threshold.

Referenced by DigiSimLinkAlgorithm::run().

                                                                            {
 
  std::vector<std::pair<int,float> > generatedNoise;
  genNoise->generate(numStrips,threshold,noiseRMS,generatedNoise, engine);
  
  // noise on strips with signal:
  for (size_t iChannel=minChannel; iChannel<maxChannel; iChannel++) {
    if(in[iChannel] != 0) {
      in[iChannel] += CLHEP::RandGaussQ::shoot(engine, 0., noiseRMS);
    }
  }

  // Noise on the other strips
  typedef std::vector<std::pair<int,float> >::const_iterator VI;  
  for(VI p = generatedNoise.begin(); p != generatedNoise.end(); p++){
    if(in[(*p).first] == 0) {
      in[(*p).first] += (*p).second;
    }
  }
}

void SiGaussianTailNoiseAdder::addNoiseVR ( std::vector< float > &  in, std::vector< float > &  noiseRMS, CLHEP::HepRandomEngine *  engine  ) const

overridevirtual

Implements SiNoiseAdder.

Definition at line 37 of file SiGaussianTailNoiseAdder.cc.

Referenced by DigiSimLinkAlgorithm::run().

                                                                                                                              {
  // Add noise
  // Full Gaussian noise is added everywhere
  for (size_t iChannel=0; iChannel!=in.size(); iChannel++) {
    if(noiseRMS[iChannel] > 0.) in[iChannel] += CLHEP::RandGaussQ::shoot(engine, 0., noiseRMS[iChannel]);
  }
}

void SiGaussianTailNoiseAdder::addPedestals ( std::vector< float > &  in, std::vector< float > &  ped  ) const

overridevirtual

Implements SiNoiseAdder.

Definition at line 45 of file SiGaussianTailNoiseAdder.cc.

Referenced by DigiSimLinkAlgorithm::run().

                                                                                             {
    for (size_t iChannel=0; iChannel!=in.size(); iChannel++) {
      if(ped[iChannel]>0.) in[iChannel] += ped[iChannel];      
    }
}

Member Data Documentation

std::unique_ptr<GaussianTailNoiseGenerator> SiGaussianTailNoiseAdder::genNoise

private

Definition at line 30 of file SiGaussianTailNoiseAdder.h.

Referenced by addNoise().

const float SiGaussianTailNoiseAdder::threshold

private

Definition at line 29 of file SiGaussianTailNoiseAdder.h.

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