CDCFudgeFactorsForSigma.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: CDC group                                                *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <map>
#include <iostream>
#include <fstream>
#include <iomanip>
#include <TObject.h>
 
namespace Belle2 {
  class CDCFudgeFactorsForSigma: public TObject {
  public:
 
    CDCFudgeFactorsForSigma() {}
 
    void setGroupID(unsigned short mode)
    {
      m_groupID = mode;
    }
 
    void setFactors(unsigned short id, const std::vector<float>& factors)
    {
      m_ffs.insert(std::pair<unsigned short, std::vector<float>>(id, factors));
    }
 
    unsigned short getGroupID() const
    {
      return m_groupID;
    }
 
    unsigned short getEntries() const
    {
      return m_ffs.size();
    }
 
    std::map<unsigned short, std::vector<float>> getFactors() const
    {
      return m_ffs;
    }
 
    const std::vector<float>& getFactors(unsigned short id) const
    {
      std::map<unsigned short, std::vector<float>>::const_iterator it = m_ffs.find(id);
      if (it != m_ffs.end()) {
        return it->second;
      } else {
        B2FATAL("Specified id not found in getFactors !");
      }
    }
 
    void dump() const
    {
      std::cout << " " << std::endl;
      std::cout << "Fudge factor list" << std::endl;
      std::cout << "#entries= " << m_ffs.size() << std::endl;
      std::cout << "in order of id and factors" << std::endl;
 
      for (auto const& ent : m_ffs) {
        std::cout << ent.first;
        unsigned short np = (ent.second).size();
        for (unsigned short i = 0; i < np; ++i) {
          std::cout << "  " << (ent.second)[i];
        }
        std::cout << std::endl;
      }
    }
 
    /*
    void outputToFile(std::string fileName) const
    {
      std::ofstream fout(fileName);
 
      if (fout.bad()) {
        B2ERROR("Specified output file could not be opened!");
      } else {
        fout << m_twParamMode << "  " << m_nTwParams << std::endl;
        for (auto const& ent : m_tws) {
          fout << std::setw(3) << std::right << ent.first;
          for (unsigned short i = 0; i < m_nTwParams; ++i) {
            fout << "  " << std::setw(15) << std::scientific << std::setprecision(8) << ent.second[i];
          }
          fout << std::endl;
        }
        fout.close();
      }
    }
    */
 
  private:
    unsigned short m_groupID = 0;   
    std::map<unsigned short, std::vector<float>> m_ffs; 
    ClassDef(CDCFudgeFactorsForSigma, 1); 
  };
 
} // end namespace Belle2