/Users/barrand/private/dev/softinex/old/inexlib-1.2/inlib/inlib/sg/fit2plot

Go to the documentation of this file.

// Copyright (C) 2010, Guy Barrand. All rights reserved.
// See the file inlib.license for terms.

#ifndef inlib_sg_fit2plot
#define inlib_sg_fit2plot

// Inheritance :
#include "plottable"

#include "../smanip"

namespace inlib {
namespace sg {

class fit2plot : public virtual plottable {
public:
  static const std::string& s_class() {
    static const std::string s_v("inlib::sg::fit2plot");
    return s_v;
  }
public:
  virtual void* cast(const std::string& a_class) const {
    if(void* p = inlib::cmp_cast<fit2plot>(this,a_class)) {return p;}
    return plottable::cast(a_class);
  }
public: //plottable
  virtual plottable* copy() const {return new fit2plot(*this);}
  virtual bool is_valid() const {return true;}
  virtual std::string name(){return m_name;}
  virtual void set_name(const std::string& a_s) {m_name = a_s;}
  virtual std::string title(){return "";}
  virtual std::string legend(){return m_legend;}
  virtual void set_legend(const std::string& a_s) {m_legend = a_s;}

  virtual std::string infos(const std::string& a_opts){
    std::string f_lf("\n");
    std::string sinfos;
    std::vector<std::string> words;
    inlib::words(a_opts," ",false,words);
    std::vector<std::string>::const_iterator it;
  
    bool show_fit_ndf = false;
   {for(it=words.begin();it!=words.end();++it) {
      if((*it)=="fit_ndf") {show_fit_ndf = true;break;}
    }}
    bool show_fit_errors = false;
   {for(it=words.begin();it!=words.end();++it) {
      if((*it)=="fit_errors") {show_fit_errors = true;break;}
    }}

    for(it=words.begin();it!=words.end();++it) {
      if(((*it)=="name") && m_name.size()) {
        if(sinfos.size()) sinfos += f_lf;
        sinfos += "Name";
        sinfos += f_lf;
        sinfos += m_name;
  
      } else if((*it)=="fit_quality") {
        if(show_fit_ndf) {
          if(sinfos.size()) sinfos += f_lf;
          sinfos += "[h]^2! / ndf";
          sinfos += f_lf; 
          sinfos += to(m_output[0]);
          sinfos += " / "; 
          sinfos += to(m_output[1]);
        } else { //show chi2 only.
          if(sinfos.size()) sinfos += f_lf;
          sinfos += "[h]^2!";
          sinfos += f_lf;
          sinfos += to(m_output[0]);
        }

      } else if((*it)=="fit_parameters") {
        unsigned int nparam = m_names.size();
        for(unsigned int iparam=0;iparam<nparam;iparam++) {
          if(show_fit_errors) {
            if(sinfos.size()) sinfos += f_lf;
            sinfos += m_names[iparam];
            sinfos += f_lf;
            sinfos += to(m_output[2+4*iparam+0]); //value
            sinfos += " +&_ ";
            sinfos += to(m_output[2+4*iparam+1]); //error
          } else {
            if(sinfos.size()) sinfos += f_lf;
            sinfos += m_names[iparam];
            sinfos += f_lf;
            sinfos += to(m_output[2+4*iparam+0]);
          }
        }
      }
    }
  
    return sinfos;
  }
public:
  fit2plot(const std::vector<std::string>& a_names,const std::vector<double>& a_output)
  :m_names(a_names)
  ,m_output(a_output)
  {
    unsigned int nparam = (m_output.size()-2)/4;
    if(m_names.size()!=nparam) {
      //should issue a warning.
      m_names.clear();
      m_output.clear();
    }
  }
  virtual ~fit2plot(){}
protected:
  fit2plot(const fit2plot& a_from)
  : plottable(a_from)
  ,m_names(a_from.m_names)
  ,m_output(a_from.m_output)
  ,m_name(a_from.m_name)
  ,m_legend(a_from.m_legend)
  {}  
private:
  fit2plot& operator=(const fit2plot&){return *this;}
private:
  std::vector<std::string> m_names;
  std::vector<double> m_output;
  std::string m_name;
  std::string m_legend;
};

}}

#endif