TT6PedestalCalculator Class Reference

#include <TT6PedestalCalculator.h>

Inheritance diagram for TT6PedestalCalculator:

Public Member Functions
void	newEvent () override
ApvAnalysis::PedestalType	pedestal () const override
ApvAnalysis::PedestalType	rawNoise () const override
void	resetPedestals () override
void	setPedestals (ApvAnalysis::PedestalType &in) override
void	setRawNoise (ApvAnalysis::PedestalType &in)
	TT6PedestalCalculator (int evnt_ini, int evnt_iter, float sig_cut)
void	updatePedestal (ApvAnalysis::RawSignalType &in) override
void	updateStatus () override
	~TT6PedestalCalculator () override
Public Member Functions inherited from TkPedestalCalculator
virtual void	setNoise (ApvAnalysis::PedestalType &)
TkStateMachine *	status ()
virtual	~TkPedestalCalculator ()

Protected Attributes
bool	alreadyUsedEvent
float	cutToAvoidSignal
int	eventsRequiredToCalibrate
int	eventsRequiredToUpdate
int	numberOfEvents
std::vector< unsigned short >	theEventPerStrip
ApvAnalysis::PedestalType	thePedestal
std::vector< double >	thePedSqSum
std::vector< double >	thePedSum
ApvAnalysis::PedestalType	theRawNoise
Protected Attributes inherited from TkPedestalCalculator
TkStateMachine	theStatus

Private Member Functions
void	init ()
void	initializePedestal (ApvAnalysis::RawSignalType &in)
void	refinePedestal (ApvAnalysis::RawSignalType &in)

Detailed Description

Concrete implementation of TkPedestalCalculator for TT6.

Definition at line 10 of file TT6PedestalCalculator.h.

Constructor & Destructor Documentation

TT6PedestalCalculator::TT6PedestalCalculator	(	int	evnt_ini,
		int	evnt_iter,
		float	sig_cut
	)

Definition at line 7 of file TT6PedestalCalculator.cc.

References gather_cfg::cout, cutToAvoidSignal, eventsRequiredToCalibrate, eventsRequiredToUpdate, and init().

    : numberOfEvents(0), alreadyUsedEvent(false) {
  if (false)
    cout << "Constructing TT6PedestalCalculator " << endl;
  eventsRequiredToCalibrate = evnt_ini;
  eventsRequiredToUpdate = evnt_iter;
  cutToAvoidSignal = sig_cut;
  init();
}

TT6PedestalCalculator::~TT6PedestalCalculator ( )

override

Definition at line 30 of file TT6PedestalCalculator.cc.

References gather_cfg::cout.

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

Member Function Documentation

void TT6PedestalCalculator::init ( void  )

private

Definition at line 19 of file TT6PedestalCalculator.cc.

References TkStateMachine::setCalibrating(), theEventPerStrip, thePedestal, thePedSqSum, thePedSum, theRawNoise, and TkPedestalCalculator::theStatus.

Referenced by TT6PedestalCalculator().

                                 {
  theRawNoise.clear();
  thePedestal.clear();
  thePedSum.clear();
  thePedSqSum.clear();
  theEventPerStrip.clear();
  theStatus.setCalibrating();
}

void TT6PedestalCalculator::initializePedestal ( ApvAnalysis::RawSignalType &  in )

private

Definition at line 60 of file TT6PedestalCalculator.cc.

References edm::DetSet< T >::data, eventsRequiredToCalibrate, mps_fire::i, cuy::ii, numberOfEvents, mathSSE::sqrt(), theEventPerStrip, thePedestal, thePedSqSum, thePedSum, and theRawNoise.

Referenced by updatePedestal().

                                                                           {
  if (numberOfEvents == 1) {
    thePedSum.clear();
    thePedSqSum.clear();
    theEventPerStrip.clear();

    thePedSum.reserve(128);
    thePedSqSum.reserve(128);
    theEventPerStrip.reserve(128);

    thePedSum.resize(in.data.size(), 0.0);
    thePedSqSum.resize(in.data.size(), 0.0);
    theEventPerStrip.resize(in.data.size(), 0);
  }
  if (numberOfEvents <= eventsRequiredToCalibrate) {
    edm::DetSet<SiStripRawDigi>::const_iterator i = in.data.begin();
    int ii = 0;
    for (; i != in.data.end(); i++) {
      thePedSum[ii] += (*i).adc();
      thePedSqSum[ii] += ((*i).adc()) * ((*i).adc());
      theEventPerStrip[ii]++;
      ii++;
    }
  }
  if (numberOfEvents == eventsRequiredToCalibrate) {
    thePedestal.clear();
    theRawNoise.clear();
    edm::DetSet<SiStripRawDigi>::const_iterator i = in.data.begin();
    int ii = 0;
    for (; i != in.data.end(); i++) {
      double avVal = (theEventPerStrip[ii]) ? thePedSum[ii] / theEventPerStrip[ii] : 0.0;
      double sqAvVal = (theEventPerStrip[ii]) ? thePedSqSum[ii] / theEventPerStrip[ii] : 0.0;
      double corr_fac = (theEventPerStrip[ii] > 1) ? (theEventPerStrip[ii] / (theEventPerStrip[ii] - 1)) : 1.0;
      double rmsVal = (sqAvVal - avVal * avVal > 0.0) ? sqrt(corr_fac * (sqAvVal - avVal * avVal)) : 0.0;

      thePedestal.push_back(static_cast<float>(avVal));
      theRawNoise.push_back(static_cast<float>(rmsVal));
      ii++;
    }
  }
}

void TT6PedestalCalculator::newEvent ( )

overridevirtual

Return raw noise, determined without CMN subtraction Tell pedestal calculator that a new event is available

Reimplemented from TkPedestalCalculator.

Definition at line 149 of file TT6PedestalCalculator.cc.

References alreadyUsedEvent.

{ alreadyUsedEvent = false; }

ApvAnalysis::PedestalType TT6PedestalCalculator::pedestal ( ) const

inlineoverridevirtual

Return reconstructed pedestals

Implements TkPedestalCalculator.

Definition at line 27 of file TT6PedestalCalculator.h.

References thePedestal.

{ return thePedestal; }

ApvAnalysis::PedestalType TT6PedestalCalculator::rawNoise ( ) const

inlineoverridevirtual

Implements TkPedestalCalculator.

Definition at line 26 of file TT6PedestalCalculator.h.

References theRawNoise.

{ return theRawNoise; }

void TT6PedestalCalculator::refinePedestal ( ApvAnalysis::RawSignalType &  in )

private

Definition at line 104 of file TT6PedestalCalculator.cc.

References cutToAvoidSignal, edm::DetSet< T >::data, eventsRequiredToCalibrate, eventsRequiredToUpdate, mps_fire::i, cuy::ii, numberOfEvents, mathSSE::sqrt(), theEventPerStrip, thePedestal, thePedSqSum, thePedSum, and theRawNoise.

Referenced by updatePedestal().

                                                                       {
  if (((numberOfEvents - eventsRequiredToCalibrate) % eventsRequiredToUpdate) == 1) {
    thePedSum.clear();
    thePedSqSum.clear();
    theEventPerStrip.clear();

    thePedSum.reserve(128);
    thePedSqSum.reserve(128);
    theEventPerStrip.reserve(128);

    thePedSum.resize(in.data.size(), 0.0);
    thePedSqSum.resize(in.data.size(), 0.0);
    theEventPerStrip.resize(in.data.size(), 0);
  }
  unsigned int ii = 0;
  ApvAnalysis::RawSignalType::const_iterator i = in.data.begin();
  for (; i < in.data.end(); i++) {
    if (fabs((*i).adc() - thePedestal[ii]) < cutToAvoidSignal * theRawNoise[ii]) {
      thePedSum[ii] += (*i).adc();
      thePedSqSum[ii] += ((*i).adc()) * ((*i).adc());
      theEventPerStrip[ii]++;
    }
    ii++;
  }
  if (((numberOfEvents - eventsRequiredToCalibrate) % eventsRequiredToUpdate) == 0) {
    for (unsigned int iii = 0; iii < in.data.size(); iii++) {
      if (theEventPerStrip[iii] > 10) {
        double avVal = (theEventPerStrip[iii]) ? thePedSum[iii] / theEventPerStrip[iii] : 0.0;
        double sqAvVal = (theEventPerStrip[iii]) ? thePedSqSum[iii] / theEventPerStrip[iii] : 0.0;
        double rmsVal = (sqAvVal - avVal * avVal > 0.0) ? sqrt(sqAvVal - avVal * avVal) : 0.0;

        if (avVal != 0) {
          thePedestal[iii] = static_cast<float>(avVal);
          theRawNoise[iii] = static_cast<float>(rmsVal);
        }
      }
    }
    thePedSum.clear();
    thePedSqSum.clear();
    theEventPerStrip.clear();
  }
}

void TT6PedestalCalculator::resetPedestals ( )

inlineoverridevirtual

Implements TkPedestalCalculator.

Definition at line 15 of file TT6PedestalCalculator.h.

References thePedestal, and theRawNoise.

                                 {
    thePedestal.clear();
    theRawNoise.clear();
  }

void TT6PedestalCalculator::setPedestals ( ApvAnalysis::PedestalType &  in )

inlineoverridevirtual

Implements TkPedestalCalculator.

Definition at line 19 of file TT6PedestalCalculator.h.

References recoMuon::in, and thePedestal.

{ thePedestal = in; }

void TT6PedestalCalculator::setRawNoise ( ApvAnalysis::PedestalType &  in )

inline

Definition at line 20 of file TT6PedestalCalculator.h.

References recoMuon::in, and theRawNoise.

{ theRawNoise = in; }

void TT6PedestalCalculator::updatePedestal ( ApvAnalysis::RawSignalType &  in )

overridevirtual

Update pedestals with current event

Implements TkPedestalCalculator.

Definition at line 45 of file TT6PedestalCalculator.cc.

References alreadyUsedEvent, initializePedestal(), TkStateMachine::isCalibrating(), TkStateMachine::isUpdating(), numberOfEvents, refinePedestal(), TkPedestalCalculator::theStatus, and updateStatus().

                                                                       {
  if (alreadyUsedEvent == false) {
    alreadyUsedEvent = true;
    numberOfEvents++;
    if (theStatus.isCalibrating()) {
      initializePedestal(in);
    } else if (theStatus.isUpdating()) {
      refinePedestal(in);
    }
    updateStatus();
  }
}

void TT6PedestalCalculator::updateStatus ( )

overridevirtual

Request that status flag be updated

Implements TkPedestalCalculator.

Definition at line 37 of file TT6PedestalCalculator.cc.

References eventsRequiredToCalibrate, TkStateMachine::isCalibrating(), numberOfEvents, TkStateMachine::setUpdating(), and TkPedestalCalculator::theStatus.

Referenced by updatePedestal().

                                         {
  if (theStatus.isCalibrating() && numberOfEvents >= eventsRequiredToCalibrate) {
    theStatus.setUpdating();
  }
}

Member Data Documentation

bool TT6PedestalCalculator::alreadyUsedEvent

protected

Definition at line 45 of file TT6PedestalCalculator.h.

Referenced by newEvent(), and updatePedestal().

float TT6PedestalCalculator::cutToAvoidSignal

protected

Definition at line 44 of file TT6PedestalCalculator.h.

Referenced by refinePedestal(), and TT6PedestalCalculator().

int TT6PedestalCalculator::eventsRequiredToCalibrate

protected

Definition at line 42 of file TT6PedestalCalculator.h.

Referenced by initializePedestal(), refinePedestal(), TT6PedestalCalculator(), and updateStatus().

int TT6PedestalCalculator::eventsRequiredToUpdate

protected

Definition at line 43 of file TT6PedestalCalculator.h.

Referenced by refinePedestal(), and TT6PedestalCalculator().

int TT6PedestalCalculator::numberOfEvents

protected

Definition at line 41 of file TT6PedestalCalculator.h.

Referenced by initializePedestal(), refinePedestal(), updatePedestal(), and updateStatus().

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

protected

Definition at line 40 of file TT6PedestalCalculator.h.

Referenced by init(), initializePedestal(), and refinePedestal().

ApvAnalysis::PedestalType TT6PedestalCalculator::thePedestal

protected

Definition at line 37 of file TT6PedestalCalculator.h.

Referenced by init(), initializePedestal(), pedestal(), refinePedestal(), resetPedestals(), and setPedestals().

std::vector<double> TT6PedestalCalculator::thePedSqSum

protected

Definition at line 39 of file TT6PedestalCalculator.h.

Referenced by init(), initializePedestal(), and refinePedestal().

std::vector<double> TT6PedestalCalculator::thePedSum

protected

Definition at line 39 of file TT6PedestalCalculator.h.

Referenced by init(), initializePedestal(), and refinePedestal().

ApvAnalysis::PedestalType TT6PedestalCalculator::theRawNoise

protected

Definition at line 38 of file TT6PedestalCalculator.h.

Referenced by init(), initializePedestal(), rawNoise(), refinePedestal(), resetPedestals(), and setRawNoise().