G4StatDouble Class Reference

#include <G4StatDouble.hh>

Public Member Functions
	G4StatDouble ()
	G4StatDouble (G4double)
virtual	~G4StatDouble ()=default
	G4StatDouble (const G4StatDouble &)=default
G4StatDouble &	operator= (const G4double &rhs)
G4StatDouble &	operator= (const G4StatDouble &rhs)=default
G4StatDouble &	operator+= (const G4double &rhs)
G4StatDouble &	operator+= (const G4StatDouble &rhs)
void	reset ()
void	fill (G4double x, G4double weight=1.)
void	scale (G4double)
G4double	mean () const
G4double	rms ()
G4double	mean (G4double ext_sum_w) const
G4double	rms (G4double ext_sum_w, G4int ext_n)
void	add (const G4StatDouble *)
G4int	n () const
G4double	sum_w () const
G4double	sum_w2 () const
G4double	sum_wx () const
G4double	sum_wx2 () const

Protected Member Functions
G4double	rms (G4double sum_wx, G4double sum_wx2, G4double sum_w, G4int n)

Protected Attributes
G4double	m_sum_wx = 0.0
G4double	m_sum_wx2 = 0.0
G4int	m_n = 0
G4double	m_sum_w = 0.0
G4double	m_sum_w2 = 0.0
G4double	m_scale = 0.0

Detailed Description

Definition at line 40 of file G4StatDouble.hh.

Constructor & Destructor Documentation

G4StatDouble() [1/3]

G4StatDouble::G4StatDouble

(

)

Definition at line 33 of file G4StatDouble.cc.

{ reset(); }

G4StatDouble() [2/3]

G4StatDouble::G4StatDouble

(

G4double

x

)

Definition at line 35 of file G4StatDouble.cc.

{ fill(x); }

~G4StatDouble()

virtual G4StatDouble::~G4StatDouble ( )

virtualdefault

G4StatDouble() [3/3]

G4StatDouble::G4StatDouble ( const G4StatDouble & )

default

Member Function Documentation

add()

void G4StatDouble::add

(

const G4StatDouble *

ptr

)

Definition at line 149 of file G4StatDouble.cc.

{
  m_n += ptr->n();
  m_sum_w += ptr->sum_w();
  m_sum_w2 += ptr->sum_w2();
  m_sum_wx += ptr->sum_wx();
  m_sum_wx2 += ptr->sum_wx2();
}

Referenced by operator+=().

fill()

void G4StatDouble::fill	(	G4double	x,
		G4double	weight = 1. )

Definition at line 47 of file G4StatDouble.cc.

{
  m_sum_wx += value * weight;
  m_sum_wx2 += value * value * weight;
  if(m_n < INT_MAX)
  {
    ++m_n;
  }
  m_sum_w += weight;
  m_sum_w2 += weight * weight;
 
  if(weight <= 0.)
  {
    G4cout << "[G4StatDouble::fill] WARNING: weight<=0. " << weight << G4endl;
  }
}

Referenced by G4StatDouble(), operator+=(), and operator=().

mean() [1/2]

G4double G4StatDouble::mean

(

)

const

Definition at line 66 of file G4StatDouble.cc.

{
  G4double mean_val = 0.;
  if(m_sum_w > 0.)
  {
    mean_val = m_sum_wx / m_sum_w;
  }
  return m_scale * mean_val;
}

Referenced by mean().

mean() [2/2]

G4double G4StatDouble::mean

(

G4double

ext_sum_w

)

const

Definition at line 76 of file G4StatDouble.cc.

{
  G4double factor = 0.;
  // factor to rescale the Mean for the requested number
  // of events (or sum of weights) ext_sum_w
 
  if(ext_sum_w > 0)
  {
    factor = m_sum_w;
    factor /= ext_sum_w;
  }
  return mean() * factor;
}

n()

G4int G4StatDouble::n ( ) const

inline

Definition at line 85 of file G4StatDouble.hh.

{ return m_n; }

Referenced by add().

operator+=() [1/2]

G4StatDouble & G4StatDouble::operator+= ( const G4double & rhs )

inline

Definition at line 56 of file G4StatDouble.hh.

  {
    fill(rhs);
    return *this;
  }

operator+=() [2/2]

G4StatDouble & G4StatDouble::operator+= ( const G4StatDouble & rhs )

inline

Definition at line 61 of file G4StatDouble.hh.

  {
    add(&rhs);
    return *this;
  }

operator=() [1/2]

G4StatDouble & G4StatDouble::operator= ( const G4double & rhs )

inline

Definition at line 49 of file G4StatDouble.hh.

  {
    reset();
    fill(rhs);
    return *this;
  }

operator=() [2/2]

G4StatDouble & G4StatDouble::operator= ( const G4StatDouble & rhs )

default

reset()

void G4StatDouble::reset

(

)

Definition at line 37 of file G4StatDouble.cc.

{
  m_sum_wx  = 0.;
  m_sum_wx2 = 0.;
  m_n       = 0;
  m_sum_w   = 0.;
  m_sum_w2  = 0.;
  m_scale   = 1.;
}

Referenced by G4StatDouble(), and operator=().

rms() [1/3]

G4double G4StatDouble::rms

(

)

Definition at line 128 of file G4StatDouble.cc.

{
  // this method computes the RMS with "all internal" parameters:
  // all the sums are the internal ones: m_sum_wx, m_sum_wx2, m_sum_w, m_n
 
  return rms(m_sum_wx, m_sum_wx2, m_sum_w, m_n);
}

Referenced by rms(), and rms().

rms() [2/3]

G4double G4StatDouble::rms	(	G4double	ext_sum_w,
		G4int	ext_n )

Definition at line 136 of file G4StatDouble.cc.

{
  // this method computes the RMS with sum_w and n coming from outside:
  // ext_sum_w and ext_n:
  // this means that the result is normalised to the external events
  // it is useful when, given a number ext_n of events with sum of the weights
  // ext_sum_w, only m_n (with sum of weights m_sum_w) are actually accumulated
  // in the internal summation (e.g. for a dose variable in a volume, because
  // only a few particles reach that volume)
 
  return rms(m_sum_wx, m_sum_wx2, ext_sum_w, ext_n);
}

rms() [3/3]

G4double G4StatDouble::rms ( G4double sum_wx, G4double sum_wx2, G4double sum_w, G4int n )

protected

Definition at line 90 of file G4StatDouble.cc.

{
  G4double vrms = 0.0;
  if(nn > 1)
  {
    G4double vmean = ssum_wx / ssum_w;
    G4double xn    = nn;
    G4double tmp =
      // from GNU Scientific Library. This part is equivalent to N/(N-1)
      // when w_i = w
      // ((m_sum_w * m_sum_w) / (m_sum_w * m_sum_w - m_sum_w2))
 
      // from NIST "DATAPLOT Reference manual", Page 2-66
      // http://www.itl.nist.gov/div898/software/dataplot/refman2/ch2/weightsd.pdf
      // rewritten based on: SUM[w(x-m)^2]/SUM[w] = SUM[wx^2]/SUM[w] - m^2
      // and dividing it by sqrt[n] to go from rms of distribution to the
      // rms of the mean value
 
      (xn / (xn - 1)) * ((ssum_wx2 / ssum_w) - (vmean * vmean));
 
    tmp  = std::max(tmp, 0.0);  // this avoids observed computation problem
    vrms = std::sqrt(tmp);
    //  G4cout << "[G4StatDoubleElement::rms] m_sum_wx: " << m_sum_wx
    //         << "  m_sum_wx2: " << m_sum_wx2 << "  m_sum_w: " << m_sum_w
    //         << "  m_n: " << m_n << "  tmp: " << tmp<< "  rms: " << rms
    //         << G4endl;
    //  G4cout << "[G4StatDoubleElement::rms] (m_n / (m_n - 1)): " << (xn/(xn -
    //  1))
    //     << "  (m_sum_wx2 / m_sum_w): " << (m_sum_wx2 / m_sum_w)
    //     << "  (mean * mean): " << (mean * mean)
    //     << "  ((m_sum_wx2 / m_sum_w) - (mean * mean)): "
    //         << ((m_sum_wx2 / m_sum_w) - (mean * mean))
    //     << G4endl;
  }
  return vrms * m_scale;
}

scale()

void G4StatDouble::scale

(

G4double

value

)

Definition at line 64 of file G4StatDouble.cc.

{ m_scale = m_scale * value; }

sum_w()

G4double G4StatDouble::sum_w ( ) const

inline

Definition at line 86 of file G4StatDouble.hh.

{ return m_sum_w; }

Referenced by add().

sum_w2()

G4double G4StatDouble::sum_w2 ( ) const

inline

Definition at line 87 of file G4StatDouble.hh.

{ return m_sum_w2; }

Referenced by add().

sum_wx()

G4double G4StatDouble::sum_wx ( ) const

inline

Definition at line 88 of file G4StatDouble.hh.

{ return m_sum_wx; }

Referenced by add().

sum_wx2()

G4double G4StatDouble::sum_wx2 ( ) const

inline

Definition at line 89 of file G4StatDouble.hh.

{ return m_sum_wx2; }

Referenced by add().

Member Data Documentation

m_n

G4int G4StatDouble::m_n = 0

protected

Definition at line 97 of file G4StatDouble.hh.

Referenced by add(), fill(), n(), reset(), and rms().

m_scale

G4double G4StatDouble::m_scale = 0.0

protected

Definition at line 100 of file G4StatDouble.hh.

Referenced by mean(), reset(), rms(), and scale().

m_sum_w

G4double G4StatDouble::m_sum_w = 0.0

protected

Definition at line 98 of file G4StatDouble.hh.

Referenced by add(), fill(), mean(), mean(), reset(), rms(), and sum_w().

m_sum_w2

G4double G4StatDouble::m_sum_w2 = 0.0

protected

Definition at line 99 of file G4StatDouble.hh.

Referenced by add(), fill(), reset(), and sum_w2().

m_sum_wx

G4double G4StatDouble::m_sum_wx = 0.0

protected

Definition at line 95 of file G4StatDouble.hh.

Referenced by add(), fill(), mean(), reset(), rms(), rms(), and sum_wx().

m_sum_wx2

G4double G4StatDouble::m_sum_wx2 = 0.0

protected

Definition at line 96 of file G4StatDouble.hh.

Referenced by add(), fill(), reset(), rms(), rms(), and sum_wx2().

---

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

• G4StatDouble.hh
• G4StatDouble.cc