PedestalSub Class Reference

#include <PedestalSub.h>

Public Member Functions
void	calculate (const std::vector< double > &inputCharge, const std::vector< double > &inputPedestal, const std::vector< double > &inputNoise, std::vector< double > &corrCharge, int soi, int nSample) const
double	getCorrection (const std::vector< double > &inputCharge, const std::vector< double > &inputPedestal, const std::vector< double > &inputNoise, int soi, int nSample) const
	PedestalSub ()
	~PedestalSub ()

Detailed Description

Definition at line 7 of file PedestalSub.h.

Constructor & Destructor Documentation

PedestalSub::PedestalSub

(

)

Definition at line 8 of file PedestalSub.cc.

                         {
}

PedestalSub::~PedestalSub

(

)

Definition at line 11 of file PedestalSub.cc.

                          { 
}

Member Function Documentation

void PedestalSub::calculate	(	const std::vector< double > &	inputCharge,
		const std::vector< double > &	inputPedestal,
		const std::vector< double > &	inputNoise,
		std::vector< double > &	corrCharge,
		int	soi,
		int	nSample
	)		const

Definition at line 15 of file PedestalSub.cc.

References getCorrection(), and mps_fire::i.

                                                                                                                                                                                                        {

  double bseCorr=PedestalSub::getCorrection(inputCharge, inputPedestal, inputNoise, soi, nTS);
  for (auto i=0; i<nTS; i++) {
    corrCharge.push_back(inputCharge[i]-inputPedestal[i]-bseCorr);
  }
}

double PedestalSub::getCorrection	(	const std::vector< double > &	inputCharge,
		const std::vector< double > &	inputPedestal,
		const std::vector< double > &	inputNoise,
		int	soi,
		int	nSample
	)		const

Definition at line 23 of file PedestalSub.cc.

References mps_fire::i.

Referenced by calculate().

                                                                                                                                                                              {

  double baseline=0;

    for (auto i=0; i<nTS; i++) {
      if (i==soi || i==(soi+1)) continue;
      if ( (inputCharge[i]-inputPedestal[i])<3*inputNoise[i]) {
    baseline+=(inputCharge[i]-inputPedestal[i]);
      }
      else {
    baseline+=3*inputNoise[i];
      }
    }
    baseline/=(nTS-2);
  return baseline;
  
}