EventReweighter.hh

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <cstddef>
#include <memory>
#include <utility>
#include <vector>
 
#include "HEJ/Config.hh"
#include "HEJ/MatrixElement.hh"
#include "HEJ/PDF.hh"
#include "HEJ/PDG_codes.hh"
#include "HEJ/Parameters.hh"
#include "HEJ/ScaleFunction.hh"
#include "HEJ/StatusCode.hh"
#include "HEJ/event_types.hh"
 
namespace LHEF {
  class HEPRUP;
}
 
namespace HEJ {
  class Event;
  struct RNG;
 
 
  struct Beam{
    double E;                                
    std::array<ParticleID, 2> type;          
  };
 
  class EventReweighter{
    using EventType = event_type::EventType;
 
  public:
    EventReweighter(
        Beam const & beam,                    
        int pdf_id,                           
        ScaleGenerator scale_gen,             
        EventReweighterConfig conf,           
        std::shared_ptr<RNG> ran              
    );
 
    EventReweighter(
        LHEF::HEPRUP const & heprup,          
        ScaleGenerator scale_gen,             
        EventReweighterConfig conf,           
        std::shared_ptr<RNG> ran              
    );
 
    PDF const & pdf() const;
 
    EventTreatment treatment(EventType type) const;
 
 
    std::vector<Event> reweight(
        Event const & ev,
        std::size_t num_events
    );
 
 
    std::vector<StatusCode> const & status() const {
        return status_;
    }
 
  private:
 
    std::vector<Event> gen_res_events(
        Event const & ev, std::size_t phase_space_points
    );
    std::vector<Event> rescale(
        Event const & Born_ev, std::vector<Event> events
    ) const;
 
    bool jets_pass_resummation_cuts(Event const & ev) const;
 
    Weights pdf_factors(Event const & ev) const;
 
    Weights matrix_elements(Event const & ev) const;
 
    Weights fixed_order_scale_ME(Event const & ev) const;
 
    double tree_matrix_element(Event const & ev) const;
 
    EventReweighterConfig param_;
 
    double E_beam_;
 
    PDF pdf_;
 
    MatrixElement MEt2_;
    ScaleGenerator scale_gen_;
    std::shared_ptr<RNG> ran_;
    std::vector<StatusCode> status_;
  };
 
} // namespace HEJ