Weight_t.h

Go to the documentation of this file.

#ifndef _WEIGHT_T_H_
#define _WEIGHT_T_H_

#include "WeightCalc.h"

#include <numeric> // std::accumulate()

namespace evwgh {
  struct Weight_t {

    Weight_t()
      : fMinWeight(std::numeric_limits<double>::max())
      , fMaxWeight(std::numeric_limits<double>::min())
      , fAvgWeight(0)
      , fNcalls(0)
    {}

    std::vector<std::vector<double>> GetWeight(art::Event& e)
    {
      std::vector<std::vector<double>> wgh = fWeightCalc->GetWeight(e);
      for (unsigned int inu = 0; inu < wgh.size(); inu++) {
        double avgwgh = std::accumulate(wgh[inu].begin(), wgh[inu].end(), 0.0) / wgh[inu].size();
        fAvgWeight = (fAvgWeight * fNcalls + avgwgh) / float(fNcalls + 1);
        fMinWeight = std::min(fMinWeight, *std::min_element(wgh[inu].begin(), wgh[inu].end()));
        fMaxWeight = std::max(fMaxWeight, *std::max_element(wgh[inu].begin(), wgh[inu].end()));
        fNcalls++;
      }

      return wgh;
    }
    std::string fName;
    WeightCalc* fWeightCalc;
    std::string fWeightCalcType;
    double fMinWeight;
    double fMaxWeight;
    double fAvgWeight;
    long fNcalls;
    int fNmultisims;
  };

}

#endif // _WEIGHT_T_H_