#include <SimplePedestalCalculator.h>

Inheritance diagram for SimplePedestalCalculator:

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)

	SimplePedestalCalculator (int evnt_ini)

void	updatePedestal (ApvAnalysis::RawSignalType &in) override

void	updateStatus () override

	~SimplePedestalCalculator () override

Public Member Functions inherited from TkPedestalCalculator
virtual void	setNoise (ApvAnalysis::PedestalType &)

TkStateMachine *	status ()

virtual	~TkPedestalCalculator ()

Protected Attributes
bool	alreadyUsedEvent

int	eventsRequiredToCalibrate

int	numberOfEvents

std::vector< unsigned short >	theEventPerStrip

ApvAnalysis::PedestalType	thePedestal

std::vector< int >	thePedSqSum

std::vector< int >	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 Simple.

Definition at line 10 of file SimplePedestalCalculator.h.

Constructor & Destructor Documentation

◆ SimplePedestalCalculator()

SimplePedestalCalculator::SimplePedestalCalculator ( int evnt_ini )

Definition at line 8 of file SimplePedestalCalculator.cc.

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

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

◆ ~SimplePedestalCalculator()

SimplePedestalCalculator::~SimplePedestalCalculator ( )

override

Definition at line 30 of file SimplePedestalCalculator.cc.

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

References gather_cfg::cout.

Member Function Documentation

◆ init()

void SimplePedestalCalculator::init ( void )

private

Definition at line 19 of file SimplePedestalCalculator.cc.

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

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

Referenced by SimplePedestalCalculator().

◆ initializePedestal()

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

private

Definition at line 63 of file SimplePedestalCalculator.cc.

                                                                               {
   if (numberOfEvents == 1) {
     thePedSum.clear();
     thePedSqSum.clear();
     theEventPerStrip.clear();
  
     thePedSum.reserve(128);
     thePedSqSum.reserve(128);
     theEventPerStrip.reserve(128);
  
     thePedSum.resize(in.data.size(), 0);
     thePedSqSum.resize(in.data.size(), 0);
     theEventPerStrip.resize(in.data.size(), 0);
   }
  
   //eventsRequiredToCalibrate is considered the minimum number of events to be used
  
   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++) {
       // the pedestal is calculated as int, as required by FED.
       int avVal = (theEventPerStrip[ii]) ? thePedSum[ii] / theEventPerStrip[ii] : 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++;
     }
   }
 }

References eventsRequiredToCalibrate, mps_fire::i, cuy::ii, recoMuon::in, numberOfEvents, mathSSE::sqrt(), theEventPerStrip, thePedestal, thePedSqSum, thePedSum, and theRawNoise.

Referenced by updatePedestal().

◆ newEvent()

void SimplePedestalCalculator::newEvent ( )

overridevirtual

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

Reimplemented from TkPedestalCalculator.

Definition at line 148 of file SimplePedestalCalculator.cc.

148 { alreadyUsedEvent = false; }

References alreadyUsedEvent.

◆ pedestal()

ApvAnalysis::PedestalType SimplePedestalCalculator::pedestal ( ) const

inlineoverridevirtual

Return reconstructed pedestals

Implements TkPedestalCalculator.

Definition at line 27 of file SimplePedestalCalculator.h.

27 { return thePedestal; }

References thePedestal.

◆ rawNoise()

ApvAnalysis::PedestalType SimplePedestalCalculator::rawNoise ( ) const

inlineoverridevirtual

Implements TkPedestalCalculator.

Definition at line 26 of file SimplePedestalCalculator.h.

26 { return theRawNoise; }

References theRawNoise.

◆ refinePedestal()

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

private

Definition at line 113 of file SimplePedestalCalculator.cc.

                                                                           {
   // keep adding th adc count for any events
  
   unsigned int ii = 0;
   ApvAnalysis::RawSignalType::const_iterator i = in.data.begin();
   for (; i < in.data.end(); i++) {
     thePedSum[ii] += (*i).adc();
     thePedSqSum[ii] += ((*i).adc()) * ((*i).adc());
     theEventPerStrip[ii]++;
  
     ii++;
   }
  
   // calculate a new pedestal any events, so it will come for free when for the last event
  
   for (unsigned int iii = 0; iii < in.data.size(); iii++) {
     if (theEventPerStrip[iii] > 10) {
       int avVal = (theEventPerStrip[iii]) ? thePedSum[iii] / theEventPerStrip[iii] : 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);
       }
     }
   }
 }

References mps_fire::i, cuy::ii, recoMuon::in, mathSSE::sqrt(), theEventPerStrip, thePedestal, thePedSqSum, thePedSum, and theRawNoise.

Referenced by updatePedestal().

◆ resetPedestals()

void SimplePedestalCalculator::resetPedestals ( )

inlineoverridevirtual

Implements TkPedestalCalculator.

Definition at line 15 of file SimplePedestalCalculator.h.

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

References thePedestal, and theRawNoise.

◆ setPedestals()

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

inlineoverridevirtual

Implements TkPedestalCalculator.

Definition at line 19 of file SimplePedestalCalculator.h.

19 { thePedestal = in; }

References recoMuon::in, and thePedestal.

◆ setRawNoise()

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

inline

Definition at line 20 of file SimplePedestalCalculator.h.

20 { theRawNoise = in; }

References recoMuon::in, and theRawNoise.

◆ updatePedestal()

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

overridevirtual

Update pedestals with current event

Implements TkPedestalCalculator.

Definition at line 47 of file SimplePedestalCalculator.cc.

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

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

◆ updateStatus()

void SimplePedestalCalculator::updateStatus ( )

overridevirtual

Request that status flag be updated

Implements TkPedestalCalculator.

Definition at line 38 of file SimplePedestalCalculator.cc.

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