UniverseOracle.h

Go to the documentation of this file.

#pragma once

#include "sbnana/CAFAna/Core/SystShifts.h"

#include <map>
#include <string>
#include <vector>

namespace ana
{
  /// Desired match type in UniverseOracle::ClosestShiftIndex
  enum class ESide{
    kAbove, kBelow, kEither
  };

  class UniverseOracle
  {
  public:
    static UniverseOracle& Instance();

    bool SystExists(const std::string& name) const;
    /// List of all known syst names
    std::vector<std::string> Systs() const;

    /// List of shifts for this syst in each universe
    const std::vector<float>& ShiftsForSyst(const std::string& name) const;

    /// Which index in the weights array corresponds to this parameter set?
    unsigned int ParameterSetIndex(const std::string& name) const;

    /// Which index in the weights array corresponds to the shifting of just this syst (in any parameter set)?
    unsigned int SystIndex(const std::string& name) const;

    /// Within that entry, which index corresponds most closely to 'shift'?
    unsigned int ClosestShiftIndex(const std::string& name,
                                   double shift,
                                   ESide side = ESide::kEither,
                                   double* trueShift = 0) const;
  protected:
    UniverseOracle();

    std::map<std::string, unsigned int> fPSetIdxs;
    std::map<std::string, unsigned int> fSystIdxs;
    std::map<std::string, std::vector<float>> fShiftVals;
  };
}