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 () override
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::SiGaussianTailNoiseAdder

(

float

th

)

Definition at line 4 of file SiGaussianTailNoiseAdder.cc.

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

~SiGaussianTailNoiseAdder()

SiGaussianTailNoiseAdder::~SiGaussianTailNoiseAdder ( )

override

Definition at line 7 of file SiGaussianTailNoiseAdder.cc.

{}

Member Function Documentation

addBaselineShift()

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

overridevirtual

Implements SiNoiseAdder.

Definition at line 67 of file SiGaussianTailNoiseAdder.cc.

References recoMuon::in, and createfilelist::int.

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)];
  }
}

addCMNoise()

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

overridevirtual

Implements SiNoiseAdder.

Definition at line 52 of file SiGaussianTailNoiseAdder.cc.

References recoMuon::in, and createfilelist::int.

Referenced by DigiSimLinkAlgorithm::run().

                                                                              {
  int nAPVs = in.size() / 128;
  std::vector<float> CMNv;
  CMNv.reserve(nAPVs);
  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)];
  }
}

addNoise()

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 9 of file SiGaussianTailNoiseAdder.cc.

References genNoise, recoMuon::in, 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;
    }
  }
}

addNoiseVR()

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

overridevirtual

Implements SiNoiseAdder.

Definition at line 34 of file SiGaussianTailNoiseAdder.cc.

References recoMuon::in.

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]);
  }
}

addPedestals()

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

overridevirtual

Implements SiNoiseAdder.

Definition at line 45 of file SiGaussianTailNoiseAdder.cc.

References recoMuon::in, and QIE10Task_cfi::ped.

Referenced by DigiSimLinkAlgorithm::run().

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

Member Data Documentation

genNoise

std::unique_ptr<GaussianTailNoiseGenerator> SiGaussianTailNoiseAdder::genNoise

private

Definition at line 30 of file SiGaussianTailNoiseAdder.h.

Referenced by addNoise().

threshold

const float SiGaussianTailNoiseAdder::threshold

private

Definition at line 29 of file SiGaussianTailNoiseAdder.h.

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