ResultHistogram.hh

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway, 
University of London 2001 - 2023.
 
This file is part of BDSIM.
 
BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.
 
BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef COMPRESULTHISTOGRAM_H
#define COMPRESULTHISTOGRAM_H
 
#include "Result.hh"
 
#include <iomanip>
#include <sstream>
#include <string>
 
class ResultHistogram: public Result
{ 
public:
  double h1Entries;
  double h2Entries; 
  int    h1NXBins;
  int    h2NXBins;
  double h1XMean;
  double h2XMean;
  double h1XRms;
  double h2XRms;
  double h1Integral;
  double h2Integral;
  double chi2;
  double tolerance;
 
  virtual std::string print() const
  {
    const int nw = 15;
    std::stringstream ss;
    ss << Result::print();
    ss << "Parameter" << std::setw(nw) << std::right << "file1"    << " / "
       << std::left   << std::setw(nw) << "file2"    << "\n";
    ss << "Entries  " << std::setw(nw) << std::right << h1Entries  << " / "
       << std::left   << std::setw(nw) << h2Entries  << "\n";
    ss << "Bins     " << std::setw(nw) << std::right << h1NXBins   << " / "
       << std::left   << std::setw(nw) << h2NXBins   << "\n";
    ss << "X Mean   " << std::setw(nw) << std::right << h1XMean    << " / "
       << std::left   << std::setw(nw) << h2XMean    << "\n";
    ss << "X RMS    " << std::setw(nw) << std::right << h1XRms     << " / "
       << std::left   << std::setw(nw) << h2XRms     << "\n";
    ss << "Integral " << std::setw(nw) << std::right << h1Integral << " / "
       << std::left   << std::setw(nw) << h2Integral << "\n";
    ss << "Chi^2    " << chi2 << " vs tolerance of " << tolerance;
    return ss.str();
  }
 
  virtual ~ResultHistogram() {}
};
 
#endif