Histogram Class Reference

Collaboration diagram for Histogram:

Public Member Functions
	Histogram ()
void	add (double relEnergy, double Expval, double stddev, Integer_t NdataInShell)
void	add (Histogram hist)
void	clear ()

Public Attributes
Integer_t	Ntot
std::vector< Integer_t >	Ndata
std::vector< Integer_t >	Nnan
std::vector< double >	ExpValue
std::vector< double >	Stddev
std::vector< double >	funcF
std::vector< double >	DOS

Private Attributes
Integer_t	m_Nbin

Constructor & Destructor Documentation

Histogram()

Histogram::Histogram ( )

inline

: m_Nbin(NBIN), Ntot(0), Ndata(NBIN,0), Nnan(NBIN,0), ExpValue(NBIN,0.0), Stddev(NBIN,0.0), funcF(NBIN,0.0), DOS(NBIN,0.0) {}

Member Function Documentation

add() [1/2]

void Histogram::add ( double  relEnergy, double  Expval, double  stddev, Integer_t  NdataInShell  )

inline

                                                                                     {
      Integer_t bin = (Integer_t)(m_Nbin*relEnergy); if(bin == m_Nbin) --bin;
      Ntot += 1;
      Ndata[bin]   += 1;
      ExpValue[bin] = ( (Ndata[bin]-1)*ExpValue[bin] + Expval )/(double)Ndata[bin];
      if( isinf(stddev) ) stddev = NAN;
      if( isnan(stddev) ) Nnan[bin] +=1;
      else {
        Integer_t Nstddev = Ndata[bin]-Nnan[bin];
        if(Nstddev != 0) {
          Stddev[bin] = ((Nstddev-1)*Stddev[bin]+stddev)/(double)Nstddev;
           funcF[bin] = ((Nstddev-1)* funcF[bin]+stddev*sqrt( (double)NdataInShell ))/(double)Nstddev;
        }
      }
    }

add() [2/2]

void Histogram::add ( Histogram  hist )

inline

                             {
      Ntot += 1;
      for(int bin = 0;bin < m_Nbin; ++bin) {
           Ndata[bin] += hist.Ndata[bin];
            Nnan[bin] +=  hist.Nnan[bin];
        ExpValue[bin] = ( (Ntot-1)*ExpValue[bin]+hist.ExpValue[bin] )/(double)Ntot;
          Stddev[bin] = ( (Ntot-1)*  Stddev[bin]+  hist.Stddev[bin] )/(double)Ntot;
           funcF[bin] = ( (Ntot-1)*   funcF[bin]+   hist.funcF[bin] )/(double)Ntot;
             DOS[bin] = ( (Ntot-1)*     DOS[bin]+     hist.DOS[bin] )/(double)Ntot;
        try{
          if( isinf(ExpValue[bin]) || isnan(ExpValue[bin]) ) throw std::runtime_error("ExpValue");
          if( isinf(  Stddev[bin]) || isnan(  Stddev[bin]) ) throw std::runtime_error("  Stddev");
          if( isinf(   funcF[bin]) || isnan(   funcF[bin]) ) throw std::runtime_error("   funcF");
          if( isinf(     DOS[bin]) || isnan(     DOS[bin]) ) throw std::runtime_error("     DOS");
        } catch(std::exception& e) {
          std::cerr << "Error:(bin=" << bin << ") " << e.what() << std::endl;
          std::cerr << "     :ExpValue=" << hist.ExpValue[bin] << ", "
                    << "Stddev=" << hist.Stddev[bin] << ", "
                    << "funcF="  << hist.funcF[bin]  << ", "
                    << "  DOS="  <<   hist.DOS[bin]  << std::endl;
          std::exit(EX_DATAERR);
        }
 
 
      }
    }

clear()

void Histogram::clear ( )

inline

                {
      Ntot = 0;
      std::fill(   Ndata.begin(),   Ndata.end(),0);
      std::fill(    Nnan.begin(),    Nnan.end(),0);
      std::fill(ExpValue.begin(),ExpValue.end(),0);
      std::fill(  Stddev.begin(),  Stddev.end(),0);
      std::fill(   funcF.begin(),   funcF.end(),0);
      std::fill(     DOS.begin(),     DOS.end(),0);
    }

Member Data Documentation

DOS

std::vector<double> Histogram::DOS

ExpValue

std::vector<double> Histogram::ExpValue

funcF

std::vector<double> Histogram::funcF

m_Nbin

Integer_t Histogram::m_Nbin

private

Ndata

std::vector<Integer_t> Histogram::Ndata

Nnan

std::vector<Integer_t> Histogram::Nnan

Ntot

Integer_t Histogram::Ntot

Stddev

std::vector<double> Histogram::Stddev

---

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

• /Users/shoki/GitHub/Locality/Processing/OpAverage.cpp