EventFileReader.hh

//
// This file is part of MARLEY (Model of Argon Reaction Low Energy Yields)
//
// MARLEY is free software: you can redistribute it and/or modify it under the
// terms of version 3 of the GNU General Public License as published by the
// Free Software Foundation.
//
// For the full text of the license please see COPYING or
// visit http://opensource.org/licenses/GPL-3.0
//
// Please respect the MCnet academic usage guidelines. See GUIDELINES
// or visit https://www.montecarlonet.org/GUIDELINES for details.
 
#pragma once
#include <deque>
#include <fstream>
#include <string>
 
#include "marley/JSON.hh"
#include "marley/OutputFile.hh"
 
namespace marley {
 
  // Forward-declare the Event class
  class Event;
 
  class EventFileReader {
 
    public:
 
      EventFileReader( const std::string& file_name );
 
      virtual ~EventFileReader() = default;
 
      virtual bool next_event( marley::Event& ev );
 
      double flux_averaged_xsec( bool natural_units = false );
 
      inline EventFileReader& operator>>( marley::Event& ev ) {
        next_event( ev );
        return *this;
      }
 
      virtual operator bool() const;
 
    protected:
 
      std::string file_name_;
 
      OutputFile::Format format_;
 
      std::ifstream in_;
 
      marley::JSON json_event_array_;
      marley::JSON::JSONWrapper< std::deque< marley::JSON > >
        json_event_array_wrapper_;
      std::deque<marley::JSON>::iterator json_event_iter_;
 
      double flux_avg_tot_xs_ = 0.;
 
      bool initialized_ = false;
 
      virtual bool deduce_file_format();
 
      virtual void initialize();
 
      void ensure_initialized();
  };
 
}