jets.hh

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <cmath>
#include <complex>
#include <ostream>
#include <vector>
 
#include "CLHEP/Vector/LorentzVector.h"
 
namespace HEJ {
namespace currents {
  using COM = std::complex<double>;
  using current = std::array<COM, 4>;
  using HLV = CLHEP::HepLorentzVector;
 
 
  double ME_qQ(HLV const & p1out, HLV const & p1in,
               HLV const & p2out, HLV const & p2in);
 
 
  double ME_qQbar(HLV const & p1out, HLV const & p1in,
                  HLV const & p2out, HLV const & p2in);
 
 
  double ME_qbarQbar(HLV const & p1out, HLV const & p1in,
                     HLV const & p2out, HLV const & p2in);
 
 
  double ME_qg(HLV const & p1out, HLV const & p1in,
               HLV const & p2out, HLV const & p2in);
 
 
  double ME_qbarg(HLV const & p1out, HLV const & p1in,
                  HLV const & p2out, HLV const & p2in);
 
 
  double ME_gg(HLV const & p1out, HLV const & p1in,
               HLV const & p2out,  HLV const & p2in);
 
  // Unordered Backwards contributions:
 
  double ME_unob_qQ(HLV const & pg, HLV const & p1out,
                    HLV const & p1in, HLV const & p2out, HLV const & p2in);
 
 
  double ME_unob_qbarQ(HLV const & pg, HLV const & p1out,
                       HLV const & p1in, HLV const & p2out, HLV const & p2in);
 
 
  double ME_unob_qQbar(HLV const & pg, HLV const & p1out,
                       HLV const & p1in, HLV const & p2out, HLV const & p2in);
 
 
  double ME_unob_qbarQbar(HLV const & pg, HLV const & p1out,
                          HLV const & p1in, HLV const & p2out,
                          HLV const & p2in);
 
 
  double ME_unob_qg(HLV const & pg, HLV const & p1out,
                    HLV const & p1in, HLV const & p2out, HLV const & p2in);
 
 
  double ME_unob_qbarg(HLV const & pg, HLV const & p1out,
                       HLV const & p1in, HLV const & p2out, HLV const & p2in);
 
 
  double ME_Exqqbar_qbarqQ(HLV const & pgin, HLV const & pqout,
                           HLV const & pqbarout, HLV const & p2out,
                           HLV const & p2in);
 
 
  double ME_Exqqbar_qqbarQ(HLV const & pgin, HLV const & pqout,
                           HLV const & pqbarout, HLV const & p2out,
                           HLV const & p2in);
 
 
  double ME_Exqqbar_qbarqg(HLV const & pgin, HLV const & pqout,
                           HLV const & pqbarout, HLV const & p2out,
                           HLV const & p2in);
 
 
  double ME_Exqqbar_qqbarg(HLV const & pgin, HLV const & pqout,
                           HLV const & pqbarout, HLV const & p2out,
                           HLV const & p2in);
 
 
  double ME_Cenqqbar_qq(
    HLV const & pa, HLV const & pb,
    std::vector<HLV> const & partons,
    bool aqlinepa, bool aqlinepb,
    bool qqbarmarker, std::size_t nabove
  );
 
  double K_g(double p1minus, double paminus);
  double K_g(HLV const & pout, HLV const & pin);
} // namespace currents
} // namespace HEJ