FakeUniverse.cxx

Go to the documentation of this file.

//
// Created by Brian Rebel on 11/15/19.
//
#include "messagefacility/MessageLogger/MessageLogger.h"

#include "CovarianceMatrixFit/dataProducts/FakeUniverse.h"

namespace cmf {
  //------------------------------------------------------------
  FakeUniverse::FakeUniverse()
  {
  }

  //------------------------------------------------------------
  FakeUniverse::FakeUniverse(cmf::Location const& loc,
                             cmf::Spectrum      & asimov,
                             cmf::Spectrum      & spect)
  {
    fAsimovSpectrum .swap(asimov);
    fPoissonSpectrum.swap(spect);

    // fill the oscillation parameters for this universe, just use
    // the FD values and ignore parameters we can't/don't measure, ie
    // rho, L, Dmsq21, Th12

    // We want to order the oscillation parameters for future use according to
    // which analysis the prediction is for:
    // For three flavor analyses, order by Th23, Dmsq32, Th13, dCP
    // For four flavor analyses, order by Th24, Dmsq41, Th34, Th23, Dmsq32

    std::vector<std::string> parOrder;
    if(loc.FDLocation().find("Th24") != loc.FDLocation().end() ){
      parOrder = std::vector<std::string>({"Th24", "Dmsq41", "Th34", "Th23", "Dmsq32"});
    }
      // look for an NSI parameter - this one needs to be fixed up
    else if(loc.FDLocation().find("Eps_mumu") != loc.FDLocation().end()){
      parOrder = std::vector<std::string>({"Eps_mumu", "Eps_ee", "Th23", "Dmsq32"});
    }
    else{
      parOrder = std::vector<std::string>({"Th23", "Dmsq32", "Th13", "dCP"});
    }

    for(auto const& itr : parOrder){
      if(loc.FDLocation().count(itr) > 0)
        fOscParVec.emplace_back(cmf::OscParm_t(loc.FDLocation().find(itr)->second.Key()),
                                float(loc.FDLocation().find(itr)->second.Value()));
    }

  }

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

  //------------------------------------------------------------
  std::ostream& operator << (std::ostream & o, cmf::FakeUniverse const& fu)
  {
    o << "universe location ";

    for(auto const& itr : fu.OscParMap())
      o << cmf::cOscParams_Strings[itr.first]
        << " "
        << itr.second
        << " ";

    o << std::endl;

//    for(size_t b = 0; b < fu.AsimovSpectrum().size(); ++b){
//      o << "fake universe spectrum bin "
//        << b
//        << " poisson: "
//        << fu.PoissonSpectrum()[b]
//        << " asimov: "
//        << fu.AsimovSpectrum()[b]
//        << "\n";
//    }

    return o;
  }

  //------------------------------------------------------------
  std::vector<float> const& FakeUniverse::AsimovSpectrum() const
  {
    return fAsimovSpectrum;
  }

  //------------------------------------------------------------
  std::vector<float> const& FakeUniverse::PoissonSpectrum() const
  {
    return fPoissonSpectrum;
  }

  //------------------------------------------------------------
  cmf::OscillationParameterMap const FakeUniverse::OscParMap() const
  {
    //loop over the far detector point and create an oscillation parameter map
    cmf::OscillationParameterMap oscPars;
    for(auto const& itr : fOscParVec){
      oscPars.emplace(itr.Param(), itr.Value());
    }

    return std::move(oscPars);
  }

} // end namespace