ParameterStorage_ROOT.h

Go to the documentation of this file.

#ifndef _ParameterStorage_ROOT_h_
#define _ParameterStorage_ROOT_h_

#ifndef __CINT__
#  include <utl/ParameterStorage.h>
#endif
#include <io/EventIO.h>
#include <Rtypes.h>
#include <map>


namespace io {

  class ParameterStorage_ROOT {

  public:
    ParameterStorage_ROOT() { }
#ifndef __CINT__
    template<class ParameterStorage>
    ParameterStorage_ROOT(const ParameterStorage& h)
    {
      const auto pars = h.GetEnumVector();
      for (const auto p : pars)
        fParameterMap[p] = std::make_pair(h.GetParameter(p), h.GetParameterLockStatus(p));

      const auto ppars = h.GetCovarianceEnumVector();
      for (const auto& pp : ppars) {
        const auto& i = pp.first;
        const auto& j = pp.second;
        fParameterCovarianceMap[pp] =
          std::make_pair(h.GetParameterCovariance(i, j), h.GetParameterCovarianceLockStatus(i, j));
      }
    }

    template<class ParameterStorage>
    void
    operator>>(ParameterStorage& h)
      const
    {
      for (const auto& ivl : fParameterMap) {
        const auto i = static_cast<typename ParameterStorage::Parameter>(ivl.first);
        const auto& vl = ivl.second;
        const auto& v = vl.first;
        const auto& l = vl.second;
        h.SetParameter(i, v, l);
      }

      for (const auto& ijvl : fParameterCovarianceMap) {
        const auto& ij = ijvl.first;
        const auto i = static_cast<typename ParameterStorage::Parameter>(ij.first);
        const auto j = static_cast<typename ParameterStorage::Parameter>(ij.second);
        const auto& vl = ijvl.second;
        const auto& v = vl.first;
        const auto& l = vl.second;
        h.SetParameterCovariance(i, j, v, l);
      }
    }
#endif

    std::map<int, std::pair<double, bool>> fParameterMap;
    std::map<std::pair<int, int>, std::pair<double, bool>> fParameterCovarianceMap;
    std::string fParameterStorageName;

    ClassDefNV(ParameterStorage_ROOT, DONT_CHANGE_THIS_NUMBER);

  };

}


#endif