rivetoutput.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2021 Christian Holm Christensen
 *    Copyright (c) 2021-2022
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_SMASH_RIVETOUTPUT_H_
#define SRC_INCLUDE_SMASH_RIVETOUTPUT_H_
 
#include <filesystem>
#include <memory>
#include <string>
#include <utility>
 
#include "Rivet/AnalysisHandler.hh"
 
#include "smash/fpenvironment.h"
#include "smash/hepmcinterface.h"
 
namespace smash {
 
class RivetOutput : public HepMcInterface {
 public:
  RivetOutput(const std::filesystem::path& path, std::string name,
              const bool full_event, const RivetOutputParameters& rivet_par);
  ~RivetOutput();
  void at_eventend(const Particles& particles, const int32_t event_number,
                   const EventInfo& event) override;
  void add_analysis(const std::string& name);
  void add_path(const std::string& path);
  void add_preload(const std::string& file);
  void set_ignore_beams(bool ignore = true);
  void set_log_level(const std::string& name, const std::string& level);
  void set_cross_section(double xs, double xserr);
 
 private:
  void setup(const RivetOutputParameters& params);
 
  struct Proxy {
    using Handler = std::shared_ptr<Rivet::AnalysisHandler>;
    explicit Proxy(Handler h) : h_(h), g_(0) {}
    Proxy(Proxy&& p) : h_(p.h_), g_(std::move(p.g_)) {}
    Handler operator->() { return h_; }
    Handler h_;
    DisableFloatTraps g_;
  };
  Proxy analysis_handler_proxy() { return Proxy(handler_); }
 
  std::shared_ptr<Rivet::AnalysisHandler> handler_;
  std::filesystem::path filename_;
  bool need_init_;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_RIVETOUTPUT_H_