SimpleNoiseCalculator Class Reference

#include <SimpleNoiseCalculator.h>

Inheritance diagram for SimpleNoiseCalculator:

Public Member Functions
int	nevents () const
void	newEvent () override
ApvAnalysis::PedestalType	noise () const override
void	resetNoise () override
void	setStripNoise (ApvAnalysis::PedestalType &in) override
	SimpleNoiseCalculator ()
	SimpleNoiseCalculator (int evnt_ini, bool useDB)
ApvAnalysis::PedestalType	stripCMPSubtractedSignal () const
float	stripNoise (int in) const override
void	updateNoise (ApvAnalysis::PedestalType &in) override
void	updateStatus () override
	~SimpleNoiseCalculator () override
Public Member Functions inherited from TkNoiseCalculator
TkStateMachine *	status ()
virtual	~TkNoiseCalculator ()

Protected Member Functions
void	init ()

Protected Attributes
bool	alreadyUsedEvent
int	eventsRequiredToCalibrate_
int	numberOfEvents
ApvAnalysis::PedestalType	theCMPSubtractedSignal
std::vector< unsigned short >	theEventPerStrip
ApvAnalysis::PedestalType	theNoise
std::vector< double >	theNoiseSqSum
std::vector< double >	theNoiseSum
bool	useDB_
Protected Attributes inherited from TkNoiseCalculator
TkStateMachine	theStatus

Detailed Description

Concrete implementation of TkNoiseCalculator for Simple.

Definition at line 9 of file SimpleNoiseCalculator.h.

Constructor & Destructor Documentation

SimpleNoiseCalculator() [1/2]

SimpleNoiseCalculator::SimpleNoiseCalculator

(

)

Definition at line 10 of file SimpleNoiseCalculator.cc.

References gather_cfg::cout, and init().

                                             : numberOfEvents(0), alreadyUsedEvent(false) {
  if (false)
    cout << "Constructing SimpleNoiseCalculator " << endl;
  init();
}

SimpleNoiseCalculator() [2/2]

SimpleNoiseCalculator::SimpleNoiseCalculator	(	int	evnt_ini,
		bool	useDB
	)

Definition at line 16 of file SimpleNoiseCalculator.cc.

References gather_cfg::cout, eventsRequiredToCalibrate_, init(), and useDB_.

                                                                      : numberOfEvents(0), alreadyUsedEvent(false) {
  if (false)
    cout << "Constructing SimpleNoiseCalculator " << endl;
  useDB_ = use_DB;
  eventsRequiredToCalibrate_ = evnt_ini;
  //  eventsRequiredToUpdate_    = evnt_iter;
  //  cutToAvoidSignal_          = sig_cut;
  init();
}

~SimpleNoiseCalculator()

SimpleNoiseCalculator::~SimpleNoiseCalculator ( )

override

Definition at line 39 of file SimpleNoiseCalculator.cc.

References gather_cfg::cout.

                                              {
  if (false)
    cout << "Destructing SimpleNoiseCalculator " << endl;
}

Member Function Documentation

init()

void SimpleNoiseCalculator::init ( void  )

protected

Definition at line 28 of file SimpleNoiseCalculator.cc.

References theCMPSubtractedSignal, theEventPerStrip, theNoise, theNoiseSqSum, and theNoiseSum.

Referenced by SimpleNoiseCalculator().

                                 {
  theCMPSubtractedSignal.clear();
  theNoise.clear();
  theNoiseSum.clear();
  theNoiseSqSum.clear();
  theEventPerStrip.clear();
  // theStatus.setCalibrating();
}

nevents()

int SimpleNoiseCalculator::nevents ( ) const

inline

Definition at line 23 of file SimpleNoiseCalculator.h.

References numberOfEvents.

{ return numberOfEvents; }

newEvent()

void SimpleNoiseCalculator::newEvent ( )

overridevirtual

Tell noise calculator that a new event is available

Reimplemented from TkNoiseCalculator.

Definition at line 112 of file SimpleNoiseCalculator.cc.

References alreadyUsedEvent.

{ alreadyUsedEvent = false; }

noise()

ApvAnalysis::PedestalType SimpleNoiseCalculator::noise ( ) const

inlineoverridevirtual

Return reconstructed noise

Implements TkNoiseCalculator.

Definition at line 21 of file SimpleNoiseCalculator.h.

References theNoise.

{ return theNoise; }

resetNoise()

void SimpleNoiseCalculator::resetNoise ( )

inlineoverridevirtual

Implements TkNoiseCalculator.

Definition at line 26 of file SimpleNoiseCalculator.h.

References theNoise.

{ theNoise.clear(); }

setStripNoise()

void SimpleNoiseCalculator::setStripNoise ( ApvAnalysis::PedestalType &  in )

inlineoverridevirtual

Implements TkNoiseCalculator.

Definition at line 17 of file SimpleNoiseCalculator.h.

References recoMuon::in, and theNoise.

                                                           {
    theNoise.clear();
    theNoise = in;
  }

stripCMPSubtractedSignal()

ApvAnalysis::PedestalType SimpleNoiseCalculator::stripCMPSubtractedSignal ( ) const

inline

Definition at line 30 of file SimpleNoiseCalculator.h.

References theCMPSubtractedSignal.

{ return theCMPSubtractedSignal; }

stripNoise()

float SimpleNoiseCalculator::stripNoise ( int  in ) const

inlineoverridevirtual

Implements TkNoiseCalculator.

Definition at line 22 of file SimpleNoiseCalculator.h.

References recoMuon::in, and theNoise.

{ return theNoise[in]; }

updateNoise()

void SimpleNoiseCalculator::updateNoise ( ApvAnalysis::PedestalType &  )

overridevirtual

Update noise with current event

Implements TkNoiseCalculator.

Definition at line 55 of file SimpleNoiseCalculator.cc.

References alreadyUsedEvent, gather_cfg::cout, eventsRequiredToCalibrate_, mps_fire::i, recoMuon::in, TkStateMachine::isCalibrating(), TkStateMachine::isUpdating(), numberOfEvents, mathSSE::sqrt(), theCMPSubtractedSignal, theEventPerStrip, theNoise, theNoiseSqSum, theNoiseSum, TkNoiseCalculator::theStatus, and updateStatus().

                                                                   {
  if (alreadyUsedEvent == false) {
    alreadyUsedEvent = true;
    numberOfEvents++;

    if (numberOfEvents == 1 && theNoise.size() != in.size()) {
      edm::LogError("SimpleNoiseCalculator:updateNoise")
          << " You did not initialize the Noise correctly prior to noise calibration.";
    }

    // Initialize sums used for estimating noise.
    if (numberOfEvents == 1) {
      theNoiseSum.clear();
      theNoiseSqSum.clear();
      theEventPerStrip.clear();

      theNoiseSum.resize(in.size(), 0.0);
      theNoiseSqSum.resize(in.size(), 0.0);
      theEventPerStrip.resize(in.size(), 0);
    }

    unsigned int i;

    // At every event Update sums used for estimating noise.
    for (i = 0; i < in.size(); i++) {
      theNoiseSum[i] += in[i];
      theNoiseSqSum[i] += in[i] * in[i];
      theEventPerStrip[i]++;
    }

    // Calculate noise.
    if ((theStatus.isCalibrating() && numberOfEvents == eventsRequiredToCalibrate_) || theStatus.isUpdating()) {
      theCMPSubtractedSignal.clear();
      theNoise.clear();

      for (i = 0; i < in.size(); i++) {
        double avVal = (theEventPerStrip[i]) ? theNoiseSum[i] / (theEventPerStrip[i]) : 0.0;
        double sqAvVal = (theEventPerStrip[i]) ? theNoiseSqSum[i] / (theEventPerStrip[i]) : 0.0;
        double corr_fac = (theEventPerStrip[i] > 1) ? (theEventPerStrip[i] / (theEventPerStrip[i] - 1)) : 1.0;
        double rmsVal = (sqAvVal - avVal * avVal > 0.0) ? sqrt(corr_fac * (sqAvVal - avVal * avVal)) : 0.0;

        theCMPSubtractedSignal.push_back(static_cast<float>(avVal));

        theNoise.push_back(static_cast<float>(rmsVal));

        if (false)
          cout << " SimpleNoiseCalculator::updateNoise " << theNoiseSum[i] << " " << theNoiseSqSum[i] << " "
               << theEventPerStrip[i] << " " << avVal << " " << sqAvVal << " " << (sqAvVal - avVal * avVal) << " "
               << rmsVal << endl;
      }
    }
    updateStatus();
  }
}

updateStatus()

void SimpleNoiseCalculator::updateStatus ( )

overridevirtual

Request that status flag be updated

Implements TkNoiseCalculator.

Definition at line 46 of file SimpleNoiseCalculator.cc.

References eventsRequiredToCalibrate_, TkStateMachine::isCalibrating(), numberOfEvents, TkStateMachine::setUpdating(), TkNoiseCalculator::theStatus, and useDB_.

Referenced by updateNoise().

                                         {
  if ((theStatus.isCalibrating() && numberOfEvents >= eventsRequiredToCalibrate_) ||
      (useDB_ == true && numberOfEvents == 1)) {
    theStatus.setUpdating();
  }
}

Member Data Documentation

alreadyUsedEvent

bool SimpleNoiseCalculator::alreadyUsedEvent

protected

Definition at line 41 of file SimpleNoiseCalculator.h.

Referenced by newEvent(), and updateNoise().

eventsRequiredToCalibrate_

int SimpleNoiseCalculator::eventsRequiredToCalibrate_

protected

Definition at line 44 of file SimpleNoiseCalculator.h.

Referenced by SimpleNoiseCalculator(), updateNoise(), and updateStatus().

numberOfEvents

int SimpleNoiseCalculator::numberOfEvents

protected

Definition at line 40 of file SimpleNoiseCalculator.h.

Referenced by nevents(), updateNoise(), and updateStatus().

theCMPSubtractedSignal

ApvAnalysis::PedestalType SimpleNoiseCalculator::theCMPSubtractedSignal

protected

Definition at line 37 of file SimpleNoiseCalculator.h.

Referenced by init(), stripCMPSubtractedSignal(), and updateNoise().

theEventPerStrip

std::vector<unsigned short> SimpleNoiseCalculator::theEventPerStrip

protected

Definition at line 39 of file SimpleNoiseCalculator.h.

Referenced by init(), and updateNoise().

theNoise

ApvAnalysis::PedestalType SimpleNoiseCalculator::theNoise

protected

Definition at line 36 of file SimpleNoiseCalculator.h.

Referenced by init(), noise(), resetNoise(), setStripNoise(), stripNoise(), and updateNoise().

theNoiseSqSum

std::vector<double> SimpleNoiseCalculator::theNoiseSqSum

protected

Definition at line 38 of file SimpleNoiseCalculator.h.

Referenced by init(), and updateNoise().

theNoiseSum

std::vector<double> SimpleNoiseCalculator::theNoiseSum

protected

Definition at line 38 of file SimpleNoiseCalculator.h.

Referenced by init(), and updateNoise().

useDB_

bool SimpleNoiseCalculator::useDB_

protected

Definition at line 42 of file SimpleNoiseCalculator.h.

Referenced by SimpleNoiseCalculator(), and updateStatus().