Sample.cxx

Go to the documentation of this file.

// CAFAna includes
#include "CAFAna/Core/Sample.h"

#include "CAFAna/Core/Registry.h"

// c++ includes
#include <iostream>

using std::ostringstream;
using std::runtime_error;
using std::string;
using std::vector;

namespace ana {
  namespace covmx {

    //-------------------------------------------------------------------------
    Sample::Sample(Selection s, Polarity p, Detector d, Quantile q)
      : selection(s), polarity(p), detector(d), quantile(q), axis(nullptr),
        cut(nullptr), pot(-1), livetime(-1), prediction(nullptr),
        data(nullptr), cosmic(nullptr) {}

    //-------------------------------------------------------------------------
    Sample::Sample(Sample const& sample)
      : selection(sample.selection), polarity(sample.polarity),
        detector(sample.detector), quantile(sample.quantile),
        axis(sample.axis), cut(sample.cut), pot(sample.pot),
        livetime(sample.livetime), prediction(sample.prediction),
        data(sample.data), cosmic(sample.cosmic), systMap(sample.systMap)
    {}

    //-------------------------------------------------------------------------
    Sample::Sample(unsigned int id)
      : axis(nullptr), cut(nullptr), pot(-1), livetime(-1),\
        prediction(nullptr), data(nullptr), cosmic(nullptr) {
      selection = (covmx::Selection)((id & (3 << 7)) >> 7);
      polarity  = (covmx::Polarity) ((id & (3 << 5)) >> 5);
      detector  = (covmx::Detector) ((id & (3 << 3)) >> 3);
      quantile  = (covmx::Quantile) (id & 7);
    }

    //-------------------------------------------------------------------------
    Sample::~Sample() {}

    //-------------------------------------------------------------------------
    string Sample::GetSelection() const {
      vector<string> sel = { "numusel", "nuesel", "ncsel" };
      return sel[selection];
    }

    //-------------------------------------------------------------------------
    string Sample::GetPolarity() const {
      vector<string> pol = { "fhc", "rhc" };
      return pol[polarity];
    }

    //-------------------------------------------------------------------------
    string Sample::GetDetector() const {
      vector<string> pol = { "neardet", "fardet" };
      return pol[detector];
    }

    //-------------------------------------------------------------------------
    string Sample::GetQuantile() const {
      if (quantile == Quantile::kNoQ)
        throw runtime_error("This sample does not have quantiles.");
      vector<string> qnt = { "q0", "q1", "q2", "q3" };
      return qnt[quantile];
    }

    //-------------------------------------------------------------------------
    string Sample::GetName() const {
      ostringstream name;
      name << GetSelection() << " " << GetPolarity() << " " << GetDetector();
      if (quantile != Quantile::kNoQ)
        name << " " << GetQuantile();
      return name.str();
    }

    //-------------------------------------------------------------------------
    string Sample::GetLatexName() const {
      ostringstream name;
      vector<string> sel = { "CC $\\nu_{\\mu}$", "CC $\\nu_{e}$", "NC" };
      name << sel[selection];
      vector<string> pol = { "FHC", "RHC" };
      name << " " << pol[polarity];
      vector<string> det = { "ND", "FD" };
      name << " " << det[detector];
      if (quantile != kNoQ) name << " Quantile " << quantile+1;
      return name.str();
    }

    //-------------------------------------------------------------------------
    string Sample::GetTag() const {
      ostringstream tag;
      tag << GetSelection() << "_" << GetPolarity() << "_" << GetDetector();
      if (quantile != Quantile::kNoQ)
        tag << "_" << GetQuantile();
      return tag.str();
    }

    //-------------------------------------------------------------------------
    void Sample::SetPrediction(IPrediction* p) {
      prediction = p;
    }

    //-------------------------------------------------------------------------
    void Sample::SetCosmic(Spectrum* s) {
      cosmic = s;
    }

    //-------------------------------------------------------------------------
    void Sample::SetData(Spectrum* s) {
      data = s;
    }

    //-------------------------------------------------------------------------
    const HistAxis* Sample::GetAxis() const {
      if (!axis) throw runtime_error("Axis not set in sample!");
      return axis;
    }

    //-------------------------------------------------------------------------
    Binning Sample::GetBinning() const {
      if (!axis) throw runtime_error("Axis not set in sample!");
      return axis->GetBinnings()[0];
    }

    //-------------------------------------------------------------------------
    const Cut* Sample::GetCut() const {
      if (!cut) throw runtime_error("Cut not set in sample!");
      return cut;
    }

    //-------------------------------------------------------------------------
    double Sample::GetPOT() const {
      if (pot == -1) throw runtime_error("POT not set in sample!");
      return pot;
    }

    //-------------------------------------------------------------------------
    double Sample::GetLivetime() const {
      if (livetime == -1) throw runtime_error("Livetime not set in sample!");
      return livetime;
    }

    //-------------------------------------------------------------------------
    IPrediction* Sample::GetPrediction() const {
      // For what it's worth, I know having a const member function that
      // returns a bare pointer to a member is extremely not good. If you
      // feel strongly about this, please come yell at me and perhaps you'll
      // inspire me to be better
      if (!prediction) throw runtime_error("Prediction not set in sample!");
      return prediction;
    }

    //-------------------------------------------------------------------------
    Spectrum* Sample::GetCosmic() const {
      if (!cosmic) throw runtime_error("Cosmic spectrum not set in sample!");
      return cosmic;
    }

    //-------------------------------------------------------------------------
    Spectrum* Sample::GetData() const {
      if (!data) throw runtime_error("Data spectrum not set in sample!");
      return data;
    }

    //-------------------------------------------------------------------------
    SystShifts Sample::GetSystShifts(SystShifts shifts) const {
      SystShifts ret;
      for (const ISyst* syst : shifts.ActiveSysts()) {
        if (systMap.count(syst)) {
    if (systMap.at(syst)) { // If there's a nullptr here, skip it
      ret.SetShift(systMap.at(syst), shifts.GetShift(syst));
    }
        } else {
    ret.SetShift(syst, shifts.GetShift(syst));
        }
      }
      return ret;
    } // function Sample::GetSystShifts

    //-------------------------------------------------------------------------
    unsigned int Sample::GetID() const
    {
      unsigned int id = selection;
      id <<= 2;
      id += polarity;
      id <<= 2;
      id += detector;
      id <<= 3;
      id += quantile;
      return id;

    } // function Sample::GetID

    //-------------------------------------------------------------------------
  } // namespace covmx
} // namespace ana