// ReactorNuOscProc.cc
// Contact person: Stefan-Alexandru Nae <stefan@lip.pt>
// See ReactorNuOscProc.hh for more details.
//---------------------------------------------//

#include <string>
#include <sstream>                          // istringstream
#include <cmath>
#include <climits>                          // INT_MAX
#include <locale>                           // std::locale, std::tolower

#include <CLHEP/Random/Randomize.h>         // HepUniformRand()

#include <RAT/ReactorNuOscProc.hh>
#include <RAT/DS/Entry.hh>
#include <RAT/DS/MCParticle.hh>             // GetMetaInfo(), ...
#include <RAT/Log.hh>

namespace RAT
{

  ReactorNuOscProc::ReactorNuOscProc() : Processor("ReactorNuOscProc") { }

  ReactorNuOscProc::~ReactorNuOscProc() { }

  void ReactorNuOscProc::BeginOfRun(DS::Run &) {
    fReactorNuOsc = DU::Utility::Get()->GetReactorNuOsc();
  }

  Processor::Result ReactorNuOscProc::DSEvent(DS::Run &, DS::Entry &ds) {
    // Get antineutrino info for generated event
    const RAT::DS::MC &mc = ds.GetMC();
    const RAT::DS::MCParticle &mcParticle = mc.GetMCParent(0);
    const double nuKE = mcParticle.GetKineticEnergy();
    std::string coreData(mcParticle.GetMetaInfo());

    bool sample = fReactorNuOsc.ReactorOsc(nuKE, coreData);

    if(sample == true)
      return Processor::OKTRUE;
    else
      return Processor::OKFALSE;
  }

} // namespace RAT