MatrixElement.hh

Go to the documentation of this file.

 
#pragma once
 
#include <functional>
#include <vector>
 
#include "fastjet/PseudoJet.hh"
 
#include "HEJ/Config.hh"
#include "HEJ/PDG_codes.hh"
#include "HEJ/Parameters.hh"
 
namespace CLHEP {
  class HepLorentzVector;
}
 
namespace HEJ {
  class Event;
  struct Particle;
 
  class MatrixElement{
  public:
    MatrixElement(
        std::function<double (double)> alpha_s,
        MatrixElementConfig conf
    );
 
    Weights operator()(Event const & event) const;
 
 
    Weights tree(Event const & event) const;
 
    std::vector<Weights> virtual_corrections(Event const & event) const;
 
    Weights tree_param(Event const & event) const;
 
    std::vector<double> tree_kin(Event const & event) const;
 
  private:
    double tree_param(
        Event const & event,
        double mur
    ) const;
 
    double virtual_corrections_W(
        Event const & event,
        double mur,
        Particle const & WBoson
    ) const;
    std::vector <double> virtual_corrections_WW(
        Event const & event,
        double mur
    ) const;
    std::vector <double> virtual_corrections_Z_qq(
        Event const & event,
        double mur,
        Particle const & ZBoson
    ) const;
    double virtual_corrections_Z_qg(
        Event const & event,
        double mur,
        Particle const & ZBoson,
        bool is_gq_event
    ) const;
 
    template<class InputIterator>
    std::vector <double> virtual_corrections_interference(
      InputIterator begin_parton, InputIterator end_parton,
      fastjet::PseudoJet const & q0_t,
      fastjet::PseudoJet const & q0_b,
      const double mur
    ) const;
 
    std::vector<double> virtual_corrections(
        Event const & event,
        double mur
    ) const;
 
    double omega0(
        double alpha_s, double mur,
        fastjet::PseudoJet const & q_j
    ) const;
 
    double tree_kin_jets(
        Event const & ev
    ) const;
    double tree_kin_W(
        Event const & ev
    ) const;
    std::vector <double> tree_kin_WW(
        Event const & ev
    ) const;
    std::vector <double> tree_kin_Z(
        Event const & ev
    ) const;
    double tree_kin_Higgs(
        Event const & ev
    ) const;
    double tree_kin_Higgs_first(
        Event const & ev
    ) const;
    double tree_kin_Higgs_last(
        Event const & ev
    ) const;
 
    double tree_kin_Higgs_between(
        Event const & ev
    ) const;
 
    double tree_param_partons(
        double alpha_s, double mur,
        std::vector<Particle> const & partons
    ) const;
 
    std::vector<int> in_extremal_jet_indices(
        std::vector<fastjet::PseudoJet> const & partons
    ) const;
 
    double MH2_backwardH(
      ParticleID type_forward,
      CLHEP::HepLorentzVector const & pa,
      CLHEP::HepLorentzVector const & pb,
      CLHEP::HepLorentzVector const & pH,
      CLHEP::HepLorentzVector const & pn
    ) const;
 
    double MH2_unob_forwardH(
      CLHEP::HepLorentzVector const & pa,
      CLHEP::HepLorentzVector const & pb,
      CLHEP::HepLorentzVector const & pg,
      CLHEP::HepLorentzVector const & p1,
      CLHEP::HepLorentzVector const & pH
    ) const;
 
    std::function<double (double)> alpha_s_;
 
    MatrixElementConfig param_;
  };
 
} // namespace HEJ