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_Zphoton(
        Event const & event,
        double mur,
        Particle const & boson
    ) const;
    double virtual_corrections_Zphoton_single(
        Event const & event,
        double mur,
        Particle const & boson,
        bool emit_fwd
    ) const;
    std::vector<double> virtual_corrections_Zphoton_double(
        Event const & event,
        double mur,
        Particle const & boson
    ) const;
    std::vector<double> virtual_corrections_Zphoton_triple(
        Event const & event,
        double mur,
        Particle const & boson
    ) const;
 
    template<class InputIterator>
    double virtual_corrections_no_interference(
      InputIterator begin_parton, InputIterator end_parton,
      fastjet::PseudoJet & q,
      const double mur
    ) const;
 
    template<class InputIterator>
    std::vector<double> virtual_corrections_interference(
      InputIterator begin_parton, InputIterator end_parton,
      fastjet::PseudoJet & q_t,
      fastjet::PseudoJet & q_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;
    std::vector <double> tree_kin_photon(
        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