waveform modifier to add noise to waveform More...

#include <PndEmcShapingNoiseAdder.h>

Inheritance diagram for PndEmcShapingNoiseAdder:

Classes
class	MA

Public Member Functions
	PndEmcShapingNoiseAdder ()

	PndEmcShapingNoiseAdder (Double_t fShapingTime, Double_t sampleRate, Double_t noiseWidth, Double_t inputScale)

virtual void	Modify (PndEmcWaveform *wf)

virtual Double_t	GetScale ()

Protected Member Functions
std::vector< Double_t > &	GetWaveformReference (PndEmcWaveform *wf)

	ClassDef (PndEmcAbsWaveformModifier, 1)

Private Attributes
Double_t	fShapingTime

Double_t	fNoiseWidth

Double_t	fInputScale

Double_t	fSingleNoiseWidth

class PndEmcShapingNoiseAdder::MA	fMA_buffer

Detailed Description

waveform modifier to add noise to waveform

Author: Ph. Mahlberg mahlb.nosp@m.erg@.nosp@m.hiskp.nosp@m..uni.nosp@m.-bonn.nosp@m..de

Definition at line 15 of file PndEmcShapingNoiseAdder.h.

Constructor & Destructor Documentation

PndEmcShapingNoiseAdder::PndEmcShapingNoiseAdder ( )

Definition at line 8 of file PndEmcShapingNoiseAdder.cxx.

                                                  : 
         fShapingTime(0),
         fNoiseWidth(0.),
         fInputScale(0.0),
         fSingleNoiseWidth(0.0),
         fMA_buffer(0) {
 }

PndEmcShapingNoiseAdder::PndEmcShapingNoiseAdder	(	Double_t	fShapingTime,
		Double_t	sampleRate,
		Double_t	noiseWidth,
		Double_t	inputScale
	)

Definition at line 16 of file PndEmcShapingNoiseAdder.cxx.

References fInputScale, fMA_buffer, fNoiseWidth, fSingleNoiseWidth, PndEmcShapingNoiseAdder::MA::GetSize(), and CAMath::Sqrt().

                                                                                                                                     :
         fShapingTime(shapingTime),
         fNoiseWidth(noiseWidth),
         fInputScale(inputScale),
         //MA buffer size as long as no of samples forming rising edge of wf
 
         fMA_buffer((Int_t) (shapingTime*sampleRate)) {
                 Int_t bufferSize = fMA_buffer.GetSize();
 
                 if(!bufferSize) {
                         std::cout << "-E PndEmcShapedNoiseWaveformSimulator: pulse raises in less than one sample. Noise generation will fail" << std::endl;
                 }
 
                 //noise of individual sample is averaged over bufferSize samples --> fNoiseWidth² = bufferSize * fSingleNoiseWidth²
                 fSingleNoiseWidth = fNoiseWidth * TMath::Sqrt(bufferSize);                      //still in GeV
 
                 //Init buffer
                 while(bufferSize--) {                                                           //init buffer
                         fMA_buffer(gRandom->Gaus(0, fSingleNoiseWidth*fInputScale));
                 }
         }

Member Function Documentation

PndEmcAbsWaveformModifier::ClassDef	(	PndEmcAbsWaveformModifier	,
		1
	)

protectedinherited

virtual Double_t PndEmcAbsWaveformModifier::GetScale ( )

inlinevirtualinherited

Reimplemented in PndEmcWfRescaler, and PndEmcWaveformDigitizer.

Definition at line 19 of file PndEmcAbsWaveformModifier.h.

19 { return 1.; };

std::vector<Double_t>& PndEmcAbsWaveformModifier::GetWaveformReference ( PndEmcWaveform * wf )

inlineprotectedinherited

Definition at line 24 of file PndEmcAbsWaveformModifier.h.

References PndEmcWaveform::fSignal.

Referenced by PndEmcWaveformDigitizer::Modify(), Modify(), and PndEmcWfRescaler::Modify().

                                                                               {
                         if(PndEmcMultiWaveform* multiWf = dynamic_cast<PndEmcMultiWaveform*>(wf)) {
                                 Int_t activeWf = multiWf->GetActiveWaveform();
                                 return multiWf->fSignals.at(activeWf);
                         }
                         return wf->fSignal;
                 }

void PndEmcShapingNoiseAdder::Modify ( PndEmcWaveform * wf )

virtual

Implements PndEmcAbsWaveformModifier.

Definition at line 39 of file PndEmcShapingNoiseAdder.cxx.

References fInputScale, fMA_buffer, fSingleNoiseWidth, and PndEmcAbsWaveformModifier::GetWaveformReference().

                                                        {
 
         for(std::vector<Double_t>::iterator it=GetWaveformReference(wf).begin(); it!=GetWaveformReference(wf).end(); ++it) { 
         
                 (*it) += fMA_buffer(gRandom->Gaus(0, fSingleNoiseWidth*fInputScale)); //smear with shaping noise
         }
 }