LHECommonBlocks.h

Go to the documentation of this file.

#ifndef SimDataFormats_GeneratorProducts_LHECommonBlocks_h
#define SimDataFormats_GeneratorProducts_LHECommonBlocks_h

#include "SimDataFormats/GeneratorProducts/interface/LesHouches.h"

extern "C" {
extern struct HEPRUP_ {
  int idbmup[2];
  double ebmup[2];
  int pdfgup[2];
  int pdfsup[2];
  int idwtup;
  int nprup;
  double xsecup[100];
  double xerrup[100];
  double xmaxup[100];
  int lprup[100];
} heprup_;

extern struct HEPEUP_ {
  int nup;
  int idprup;
  double xwgtup;
  double scalup;
  double aqedup;
  double aqcdup;
  int idup[500];
  int istup[500];
  int mothup[500][2];
  int icolup[500][2];
  double pup[500][5];
  double vtimup[500];
  double spinup[500];
} hepeup_;
}  // extern "C"

namespace lhef {

  class CommonBlocks {
  public:
    static void fillHEPRUP(const HEPRUP *heprup) {
      heprup_.idbmup[0] = heprup->IDBMUP.first;
      heprup_.idbmup[1] = heprup->IDBMUP.second;
      heprup_.ebmup[0] = heprup->EBMUP.first;
      heprup_.ebmup[1] = heprup->EBMUP.second;
      heprup_.pdfgup[0] = heprup->PDFGUP.first;
      heprup_.pdfgup[1] = heprup->PDFGUP.second;
      heprup_.pdfsup[0] = heprup->PDFSUP.first;
      heprup_.pdfsup[1] = heprup->PDFSUP.second;
      heprup_.idwtup = heprup->IDWTUP;
      heprup_.nprup = heprup->NPRUP;
      for (int i = 0; i < heprup->NPRUP; i++) {
        heprup_.xsecup[i] = heprup->XSECUP[i];
        heprup_.xerrup[i] = heprup->XERRUP[i];
        heprup_.xmaxup[i] = heprup->XMAXUP[i];
        heprup_.lprup[i] = heprup->LPRUP[i];
      }
    }

    static void fillHEPEUP(const HEPEUP *hepeup) {
      hepeup_.nup = hepeup->NUP;
      hepeup_.idprup = hepeup->IDPRUP;
      hepeup_.xwgtup = hepeup->XWGTUP;
      hepeup_.scalup = hepeup->SCALUP;
      hepeup_.aqedup = hepeup->AQEDUP;
      hepeup_.aqcdup = hepeup->AQCDUP;
      for (int i = 0; i < hepeup->NUP; i++) {
        hepeup_.idup[i] = hepeup->IDUP[i];
        hepeup_.istup[i] = hepeup->ISTUP[i];
        hepeup_.mothup[i][0] = hepeup->MOTHUP[i].first;
        hepeup_.mothup[i][1] = hepeup->MOTHUP[i].second;
        hepeup_.icolup[i][0] = hepeup->ICOLUP[i].first;
        hepeup_.icolup[i][1] = hepeup->ICOLUP[i].second;
        for (unsigned int j = 0; j < 5; j++)
          hepeup_.pup[i][j] = hepeup->PUP[i][j];
        hepeup_.vtimup[i] = hepeup->VTIMUP[i];
        hepeup_.spinup[i] = hepeup->SPINUP[i];
      }
    }

    static void readHEPRUP(HEPRUP *heprup) {
      heprup->resize(heprup_.nprup);
      heprup->IDBMUP.first = heprup_.idbmup[0];
      heprup->IDBMUP.second = heprup_.idbmup[1];
      heprup->EBMUP.first = heprup_.ebmup[0];
      heprup->EBMUP.second = heprup_.ebmup[1];
      heprup->PDFGUP.first = heprup_.pdfgup[0];
      heprup->PDFGUP.second = heprup_.pdfgup[1];
      heprup->PDFSUP.first = heprup_.pdfsup[0];
      heprup->PDFSUP.second = heprup_.pdfsup[1];
      heprup->IDWTUP = heprup_.idwtup;
      for (int i = 0; i < heprup->NPRUP; i++) {
        heprup->XSECUP[i] = heprup_.xsecup[i];
        heprup->XERRUP[i] = heprup_.xerrup[i];
        heprup->XMAXUP[i] = heprup_.xmaxup[i];
        heprup->LPRUP[i] = heprup_.lprup[i];
      }
    }

    static void readHEPEUP(HEPEUP *hepeup) {
      hepeup->resize(hepeup_.nup);
      hepeup->IDPRUP = hepeup_.idprup;
      hepeup->XWGTUP = hepeup_.xwgtup;
      hepeup->SCALUP = hepeup_.scalup;
      hepeup->AQEDUP = hepeup_.aqedup;
      hepeup->AQCDUP = hepeup_.aqcdup;
      for (int i = 0; i < hepeup->NUP; i++) {
        hepeup->IDUP[i] = hepeup_.idup[i];
        hepeup->ISTUP[i] = hepeup_.istup[i];
        hepeup->MOTHUP[i].first = hepeup_.mothup[i][0];
        hepeup->MOTHUP[i].second = hepeup_.mothup[i][1];
        hepeup->ICOLUP[i].first = hepeup_.icolup[i][0];
        hepeup->ICOLUP[i].second = hepeup_.icolup[i][1];
        for (unsigned int j = 0; j < 5; j++)
          hepeup->PUP[i][j] = hepeup_.pup[i][j];
        hepeup->VTIMUP[i] = hepeup_.vtimup[i];
        hepeup->SPINUP[i] = hepeup_.spinup[i];
      }
    }

  private:
    CommonBlocks() = delete;
    ~CommonBlocks() = delete;
  };

}  // namespace lhef

#endif  // SimDataFormats_GeneratorProducts_LHECommonBlocks_h