GFException.h

Go to the documentation of this file.

/* Copyright 2008-2010, Technische Universitaet Muenchen,
   Authors: Christian Hoeppner & Sebastian Neubert

   This file is part of GENFIT.

   GENFIT is free software: you can redistribute it and/or modify
   it under the terms of the GNU Lesser General Public License as published
   by the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.

   GENFIT 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 Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public License
   along with GENFIT.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef GFEXCEPTION_H
#define GFEXCEPTION_H

#include <exception>
#include <iomanip> // std::setw()
#include <ios>     // std::ios::fmtflags
#include <sstream>
#include <string>
#include <vector>

#include "RtypesCore.h"
#include "TMatrixT.h"
#include "TVector3.h"

class GFException : public std::exception {
private:
  static bool fQuiet;

  std::string fExcString;
  int fLine;
  std::string fFile;

  std::string fNumbersLabel;
  std::string fMatricesLabel;
  std::vector<double> fNumbers;
  std::vector<TMatrixT<Double_t>> fMatrices;

  bool fFatal;

public:
  GFException(std::string, int, std::string);
  virtual ~GFException() throw();

  GFException& setFatal(bool b = true)
  {
    fFatal = b;
    return *this;
  }
  bool isFatal() { return fFatal; }
  GFException& setNumbers(std::string, const std::vector<double>&);
  GFException& setMatrices(std::string, const std::vector<TMatrixT<Double_t>>&);

  void info();

  const char* what() const noexcept override;

  std::string getExcString() { return fExcString; }

  static void quiet(bool b = true) { fQuiet = b; }
};

namespace genf {
  //------------------------------------------------------------------------------
  template <class ROOTOBJ>
  void PrintROOTobject(std::ostream&, const ROOTOBJ&);

  template <>
  void PrintROOTobject(std::ostream&, const TVector3& v);

  template <typename T>
  void PrintROOTmatrix(std::ostream& out, const TMatrixT<T>& m);

  template <class ROOTOBJ>
  std::string ROOTobjectToString(const ROOTOBJ& obj)
  {
    std::ostringstream sstr;
    PrintROOTobject(sstr, obj);
    return sstr.str();
  }

} // namespace genf

//------------------------------------------------------------------------------
// template definitions
//
template <class ROOTOBJ>
void genf::PrintROOTobject(std::ostream&, const ROOTOBJ& obj)
{
  obj.Print();
}

template <typename T>
void genf::PrintROOTmatrix(std::ostream& out, const TMatrixT<T>& m)
{

  constexpr std::streamsize fw = 11;
  constexpr std::streamsize ifw = 4 + (fw & 1);
  const Int_t rb = m.GetRowLwb(), cb = m.GetColLwb();

  const Int_t R = m.GetNrows(), C = m.GetNcols();
  out << R << "x" << C << " matrix is as follows";

  std::streamsize swidth = out.width(4);
  std::ios::fmtflags sflags = out.flags();
  out.unsetf(std::ios_base::floatfield); // out << std::defaultfloat;

  // header: column number
  std::string index_pad((fw - ifw) / 2, ' ');
  out << "\n" << std::string(ifw, ' ') << " |";
  for (Int_t c = 0; c < C; ++c)
    out << index_pad << std::setw(ifw) << (cb + c) << index_pad << "|";

  // dashed line
  out << "\n" << std::string((C + 1) * (fw + 1), '-');

  // content, row by row
  for (Int_t r = 0; r < R; ++r) {
    // header: row number
    out << "\n" << std::setw(ifw) << (rb + r) << " |";
    for (Int_t c = 0; c < C; ++c)
      out << std::setw(fw) << m(rb + r, cb + c) << " ";
  } // for r
  out << "\n\n";

  // restore the stream features
  out.flags(sflags);
  out.width(swidth);
} // genf::PrintROOTmatrix<TMatrixT<T>>()

#endif