MatrixElement.hh

Go to the documentation of this file.

 
#pragma once
 
#include <functional>
#include <vector>
 
#include "fastjet/PseudoJet.hh"
 
#include "HEJ/PDG_codes.hh"
#include "HEJ/Weights.hh"
#include "HEJ/config.hh"
 
namespace CLHEP {
  class HepLorentzVector;
}
 
namespace HEJ{
  class Event;
  class 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;
 
    Weights virtual_corrections(Event const & event) const;
 
  private:
 
    Weights tree_param(
        Event const & event
    ) const;
 
    double tree_kin(Event const & event) const;
 
    double tree_param(
        Event const & event,
        double mur
    ) const;
 
    double virtual_corrections(
        Event const & event,
        double mur
    ) const;
 
    double omega0(
        double alpha_s, double mur,
        fastjet::PseudoJet const & q_j, double lambda
    ) const;
 
    double tree_kin_jets(
        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;
 
 
    std::vector<Particle> tag_extremal_jet_partons(
        Event const & ev
    ) const;
 
    double MH2_forwardH(
        CLHEP::HepLorentzVector p1out, CLHEP::HepLorentzVector p1in,
        pid::ParticleID type2,
        CLHEP::HepLorentzVector p2out, CLHEP::HepLorentzVector p2in,
        CLHEP::HepLorentzVector pH,
        double t1, double t2
    ) const;
 
    std::function<double (double)> alpha_s_;
 
    MatrixElementConfig param_;
  };
 
 
}