LesHouches.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Copyright (C) 1999-2007 Leif Lonnblad
//
// Modified by C.Saout
//
#ifndef SimDataFormats_GeneratorProducts_LesHouches_h
#define SimDataFormats_GeneratorProducts_LesHouches_h
 
#include <utility>
#include <vector>
 
namespace lhef {
 
  class HEPRUP {
  public:
    HEPRUP() : IDWTUP(0), NPRUP(0) {}
 
  public:
    bool operator==(const HEPRUP &other) const {
      return IDBMUP == other.IDBMUP && EBMUP == other.EBMUP && PDFGUP == other.PDFGUP && PDFSUP == other.PDFSUP &&
             IDWTUP == other.IDWTUP && NPRUP == other.NPRUP && XSECUP == other.XSECUP && XERRUP == other.XERRUP &&
             XMAXUP == other.XMAXUP && LPRUP == other.LPRUP;
    }
 
    void resize(int nrup) {
      NPRUP = nrup;
      resize();
    }
 
    void resize() {
      XSECUP.resize(NPRUP);
      XERRUP.resize(NPRUP);
      XMAXUP.resize(NPRUP);
      LPRUP.resize(NPRUP);
    }
 
    void swap(HEPRUP &other) {
      IDBMUP.swap(other.IDBMUP);
      EBMUP.swap(other.EBMUP);
      PDFGUP.swap(other.PDFGUP);
      PDFSUP.swap(other.PDFSUP);
      std::swap(IDWTUP, other.IDWTUP);
      std::swap(NPRUP, other.NPRUP);
      XSECUP.swap(other.XSECUP);
      XERRUP.swap(other.XERRUP);
      XMAXUP.swap(other.XMAXUP);
      LPRUP.swap(other.LPRUP);
    }
 
    std::pair<int, int> IDBMUP;
 
    std::pair<double, double> EBMUP;
 
    std::pair<int, int> PDFGUP;
 
    std::pair<int, int> PDFSUP;
 
    int IDWTUP;
 
    int NPRUP;
 
    std::vector<double> XSECUP;
 
    std::vector<double> XERRUP;
 
    std::vector<double> XMAXUP;
 
    std::vector<int> LPRUP;
  };
 
  class HEPEUP {
  public:
    HEPEUP() : NUP(0), IDPRUP(0), XWGTUP(0.0), XPDWUP(0.0, 0.0), SCALUP(0.0), AQEDUP(0.0), AQCDUP(0.0) {}
 
  public:
    struct FiveVector {
      double operator[](unsigned int i) const { return x[i]; }
      double &operator[](unsigned int i) { return x[i]; }
 
      double x[5];
    };
 
    void resize(int nup) {
      NUP = nup;
      resize();
    }
 
    void resize() {
      IDUP.resize(NUP);
      ISTUP.resize(NUP);
      MOTHUP.resize(NUP);
      ICOLUP.resize(NUP);
      PUP.resize(NUP);
      VTIMUP.resize(NUP);
      SPINUP.resize(NUP);
    }
 
    int NUP;
 
    int IDPRUP;
 
    double XWGTUP;
 
    std::pair<double, double> XPDWUP;
 
    double SCALUP;
 
    double AQEDUP;
 
    double AQCDUP;
 
    std::vector<int> IDUP;
 
    std::vector<int> ISTUP;
 
    std::vector<std::pair<int, int> > MOTHUP;
 
    std::vector<std::pair<int, int> > ICOLUP;
 
    std::vector<FiveVector> PUP;
 
    std::vector<double> VTIMUP;
 
    std::vector<double> SPINUP;
  };
 
}  // namespace lhef
 
#endif  // SimDataFormats_GeneratorProducts_LesHouches_h