BIE Overview Quick Start Theoretical overview Available methods Parallel chains Computation engine Command Line Interface Graphical User Interface Assigning Output Visualization tool Data handling Software technology Parallel debugging Results to date Recent developments Links User guide Future features Project goals Project Team Copyright and license

/home/weinberg/src/BIE/include/Histogram1D.h // This is really -*- C++ -*- #ifndef Histogram1D_h #define Histogram1D_h #include "SmplHist.h" #include "Distribution.h" #include "Serializable.h" namespace BIE { class Bin: public Serializable { public: Bin() {}; double bot; double center; double top; double val; bool operator<(const double& x) const { return(x<top && x>bot); } // AUTO GENERATED BY ../persistence/autopersist.py private: friend class boost::serialization::access; template<class Archive> void serialize(Archive & ar, const unsigned int version) { this->pre_serialize(ar, version); try { ar & BOOST_SERIALIZATION_BASE_OBJECT_NVP(Serializable); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bot); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(center); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(top); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(val); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; class CumSampleHistogram : public SampleHistogram { private: bool accum; void compute_accum(void); int n; vector<Bin> cum; vector<double> bin; public: CumSampleHistogram(double low, double hi, double bucketWidth = -1.0) : SampleHistogram(low, hi, bucketWidth) { accum = false; } double quartile(double s); double fraction(double s); // AUTO GENERATED BY ../persistence/autopersist.py protected: CumSampleHistogram() {} private: friend class boost::serialization::access; template<class Archive> void serialize(Archive & ar, const unsigned int version) { this->pre_serialize(ar, version); try { ar & BOOST_SERIALIZATION_BASE_OBJECT_NVP(SampleHistogram); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(accum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(n); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(cum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; //+ CLICLASS Histogram1D SUPER BinnedDistribution class Histogram1D : public BinnedDistribution { private: CumSampleHistogram * ch; double low, high, width, toplimit; double val; bool AccumulateData(double v, vector<double> & x); bool AccumulateData(double v, RecordBuffer * datapoint); void initialize(double low, double high, double width); public: //+ CLICONSTR double double double RecordType* Histogram1D(double low, double hi, double w, RecordType * type); //+ CLICONSTR double double double string Histogram1D(double low, double hi, double w, string dataname); ~Histogram1D(); Histogram1D* New(); double PDF(vector<double>& x) {return ch->similarSamples(x[0])/width;} double logPDF(vector<double>& x) {return log(ch->similarSamples(x[0])/width);} double CDF(vector<double>& x) {return ch->fraction(x[0]);} vector<double> lower(void) {return vector<double>(1, ch->min());} vector<double> upper(void) {return vector<double>(1, ch->max());} vector<double> Mean(void) {return vector<double>(1, ch->mean());} vector<double> StdDev(void) {return vector<double>(1, ch->stdDev());} vector<double> Moments(int i) { if (i==0) return vector<double>(1, 1.0); if (i==2) return vector<double>(1, ch->stdDev()*ch->stdDev() + ch->mean()*ch->mean()); return vector<double>(1, ch->mean()); } vector<double> Sample(void) { cerr << "Histogram1D::Sample: Not implemented\n"; return vector<double>(1, 0.0); } void ComputeDistribution(); int numberData() { return ch->buckets()-2; } double getValue(int i) { if (i<0) i=0; if (i>=ch->buckets()-2) i=ch->buckets()-3; return (double)ch->inBucket(i+1); } vector<double> getLow(int i) { if (i<0) i=0; if (i>=ch->buckets()-2) i=ch->buckets()-3; return vector<double>(1, ch->bucketThreshold(i)); } vector<double> getHigh(int i) { if (i<0) i=0; if (i>=ch->buckets()-2) i=ch->buckets()-3; return vector<double>(1, ch->bucketThreshold(i+1)); } int getdim(int i) {return 1;} // AUTO GENERATED BY ../persistence/autopersist.py protected: Histogram1D() {} private: friend class boost::serialization::access; template<class Archive> void serialize(Archive & ar, const unsigned int version) { this->pre_serialize(ar, version); try { ar & BOOST_SERIALIZATION_BASE_OBJECT_NVP(BinnedDistribution); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ch); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(low); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(high); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(width); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(toplimit); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(val); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } // namespace BIE BIE_CLASS_TYPE_INFO(BIE::Bin) BIE_CLASS_TYPE_INFO(BIE::CumSampleHistogram) BIE_CLASS_TYPE_INFO(BIE::Histogram1D) BIE_CLASS_EXPORT_KEY(BIE::Bin) BIE_CLASS_EXPORT_KEY(BIE::CumSampleHistogram) BIE_CLASS_EXPORT_KEY(BIE::Histogram1D) #endif