/data/refman/pasoursint/CMSSW_5_3_9_patch3/src/SimDataFormats/GeneratorProducts/interface/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);
        }

        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