opdet::FlashHypothesis Class Reference

#include <FlashHypothesis.h>

Public Member Functions
	FlashHypothesis ()
	FlashHypothesis (size_t s)
	FlashHypothesis (std::vector< float > const &vector, std::vector< float > const &vector_error=std::vector< float >())
std::vector< float > const &	GetHypothesisVector () const
std::vector< float > const &	GetHypothesisErrorVector () const
void	SetHypothesisVector (std::vector< float > v)
void	SetHypothesisErrorVector (std::vector< float > v)
void	SetHypothesisVectorAndErrorVector (std::vector< float > v, std::vector< float > err=std::vector< float >(0))
float const &	GetHypothesis (size_t i_opdet) const
float const &	GetHypothesisError (size_t i_opdet) const
void	SetHypothesis (size_t i_opdet, float pe)
void	SetHypothesisError (size_t i_opdet, float err)
void	SetHypothesisAndError (size_t i_opdet, float pe, float err=-999)
float	GetTotalPEs () const
float	GetTotalPEsError () const
size_t	GetVectorSize () const
void	Normalize (float const &totalPE_target)
void	Print ()
FlashHypothesis	operator+ (const FlashHypothesis &fh)

Private Attributes
std::vector< float >	_NPEs_Vector
std::vector< float >	_NPEs_ErrorVector

Detailed Description

Definition at line 18 of file FlashHypothesis.h.

Constructor & Destructor Documentation

opdet::FlashHypothesis::FlashHypothesis ( )

inline

Definition at line 21 of file FlashHypothesis.h.

{}

opdet::FlashHypothesis::FlashHypothesis ( size_t  s )

inline

Definition at line 22 of file FlashHypothesis.h.

      { _NPEs_Vector = std::vector<float>(s,0.0); _NPEs_ErrorVector = std::vector<float>(s,0.0); }

opdet::FlashHypothesis::FlashHypothesis ( std::vector< float > const &  vector, std::vector< float > const &  vector_error = std::vector<float>()  )

inline

Definition at line 24 of file FlashHypothesis.h.

      { SetHypothesisVectorAndErrorVector(vector,vector_error); }

Member Function Documentation

float const& opdet::FlashHypothesis::GetHypothesis ( size_t  i_opdet ) const

inline

Definition at line 34 of file FlashHypothesis.h.

{ return _NPEs_Vector.at(i_opdet); }

float const& opdet::FlashHypothesis::GetHypothesisError ( size_t  i_opdet ) const

inline

Definition at line 35 of file FlashHypothesis.h.

{ return _NPEs_ErrorVector.at(i_opdet); }

std::vector<float> const& opdet::FlashHypothesis::GetHypothesisErrorVector ( ) const

inline

Definition at line 29 of file FlashHypothesis.h.

{ return _NPEs_ErrorVector; }

std::vector<float> const& opdet::FlashHypothesis::GetHypothesisVector ( ) const

inline

Definition at line 28 of file FlashHypothesis.h.

{ return _NPEs_Vector; }

float opdet::FlashHypothesis::GetTotalPEs ( ) const

inline

Definition at line 41 of file FlashHypothesis.h.

    { return std::accumulate(_NPEs_Vector.begin(),_NPEs_Vector.end(),0.0); }

float opdet::FlashHypothesis::GetTotalPEsError ( ) const

inline

Definition at line 43 of file FlashHypothesis.h.

    { return std::sqrt( std::inner_product(_NPEs_ErrorVector.begin(),_NPEs_ErrorVector.end(),_NPEs_ErrorVector.begin(),0.0) ); }

size_t opdet::FlashHypothesis::GetVectorSize ( ) const

inline

Definition at line 46 of file FlashHypothesis.h.

{ return _NPEs_Vector.size(); }

void opdet::FlashHypothesis::Normalize

(

float const &

totalPE_target

)

Definition at line 30 of file FlashHypothesis.cxx.

                                                               {

  if( GetTotalPEs() < std::numeric_limits<float>::epsilon() ) return;

  const float PE_ratio = totalPE_target/GetTotalPEs();
  for(size_t i_opdet=0; i_opdet<_NPEs_Vector.size(); i_opdet++){
    _NPEs_Vector[i_opdet] *= PE_ratio;
    _NPEs_ErrorVector[i_opdet] = std::sqrt(_NPEs_Vector[i_opdet]);
  }

}

FlashHypothesis opdet::FlashHypothesis::operator+ ( const FlashHypothesis &  fh )

inline

Definition at line 52 of file FlashHypothesis.h.

                                                        {

      if( _NPEs_Vector.size() != fh.GetVectorSize() )
        throw std::runtime_error("ERROR in FlashHypothesisAddition: Cannot add hypothesis of different size");

      FlashHypothesis flashhyp(_NPEs_Vector.size());
      for(size_t i=0; i<_NPEs_Vector.size(); i++){
        flashhyp._NPEs_Vector[i] = _NPEs_Vector[i] + fh._NPEs_Vector[i];
        flashhyp._NPEs_ErrorVector[i] =
          std::sqrt(this->_NPEs_ErrorVector[i]*this->_NPEs_ErrorVector[i] +
                    fh._NPEs_ErrorVector[i]*fh._NPEs_ErrorVector[i]);
      }
      return flashhyp;
    }

void opdet::FlashHypothesis::Print

(

void

)

Definition at line 42 of file FlashHypothesis.cxx.

{
  std::cout << "TotalPEs: " << GetTotalPEs() << " +/- " << GetTotalPEsError() << std::endl;
  std::cout << "Vector size: " << GetVectorSize() << std::endl;
  for(size_t i=0; i<GetVectorSize(); i++)
    std::cout << "\t" << i << ": " << GetHypothesis(i) << " +/- " << GetHypothesisError(i) << std::endl;
}

void opdet::FlashHypothesis::SetHypothesis ( size_t  i_opdet, float  pe  )

inline

Definition at line 36 of file FlashHypothesis.h.

{ _NPEs_Vector.at(i_opdet)=pe; }

void opdet::FlashHypothesis::SetHypothesisAndError	(	size_t	i_opdet,
		float	pe,
		float	err = -999
	)

Definition at line 23 of file FlashHypothesis.cxx.

{
  SetHypothesis(i_opdet,pe);
  if(err>0) SetHypothesisError(i_opdet,err);
  else SetHypothesisError(i_opdet,std::sqrt(pe));
}

void opdet::FlashHypothesis::SetHypothesisError ( size_t  i_opdet, float  err  )

inline

Definition at line 37 of file FlashHypothesis.h.

{ _NPEs_ErrorVector.at(i_opdet) = err; }

void opdet::FlashHypothesis::SetHypothesisErrorVector ( std::vector< float >  v )

inline

Definition at line 31 of file FlashHypothesis.h.

{ _NPEs_ErrorVector = v; _NPEs_Vector.resize(v.size()); }

void opdet::FlashHypothesis::SetHypothesisVector ( std::vector< float >  v )

inline

Definition at line 30 of file FlashHypothesis.h.

{ _NPEs_Vector=v; _NPEs_ErrorVector.resize(v.size()); }

void opdet::FlashHypothesis::SetHypothesisVectorAndErrorVector	(	std::vector< float >	v,
		std::vector< float >	err = std::vector<float>(0)
	)

Definition at line 6 of file FlashHypothesis.cxx.

{
  if(err.size()!=0 && err.size()!=v.size())
    throw std::runtime_error("ERROR in FlashHypothesisVectorSetter: Vector sizes not equal");

  _NPEs_Vector = v;
  if(err.size()==0){
    _NPEs_Vector = v;
    _NPEs_ErrorVector.resize(v.size());
    for(size_t i=0; i<_NPEs_Vector.size(); i++)
      _NPEs_ErrorVector[i] = std::sqrt(_NPEs_Vector[i]);
  }
  else
    _NPEs_ErrorVector = err;

}

Member Data Documentation

std::vector<float> opdet::FlashHypothesis::_NPEs_ErrorVector

private

Definition at line 69 of file FlashHypothesis.h.

std::vector<float> opdet::FlashHypothesis::_NPEs_Vector

private

Definition at line 68 of file FlashHypothesis.h.

---

The documentation for this class was generated from the following files:

• larana/larana/OpticalDetector/FlashHypothesis.h
• larana/larana/OpticalDetector/FlashHypothesis.cxx