Weighter.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
 
// Standard library includes
#include <memory>
#include <string>
#include <vector>
 
// MARLEY includes
#include "marley/EventProcessor.hh"
 
namespace marley {
 
  // Forward-declare the WeightCalculator class
  class WeightCalculator;

  class Weighter : public EventProcessor {
 
    public:
 
      Weighter( const marley::JSON& config,
        marley::Generator& gen );
 
      inline virtual ~Weighter() = default;
 
      virtual void process_event( HepMC3::GenEvent& event,
        marley::Generator& gen ) override;
 
      std::vector< double > compute_weights( HepMC3::GenEvent& event,
        marley::Generator& gen ) const;

      std::vector< std::string > get_weight_names() const;

      void set_use_cv_weight( bool use_it );

      inline std::vector< std::shared_ptr< WeightCalculator > >&
        get_weight_calculators() { return calc_vec_; }
 
    protected:

      static const std::string CV_WEIGHT_NAME;

      std::vector< std::shared_ptr< WeightCalculator > > calc_vec_;

      std::shared_ptr< marley::WeightCalculator > make_cv_weight_calc();
  };
 
}