/data/refman/pasoursint/CMSSW_5_3_10/src/SimDataFormats/GeneratorProducts/interface/GenRunInfoProduct.h

Go to the documentation of this file.

#ifndef SimDataFormats_GeneratorProducts_GenRunInfoProduct_h
#define SimDataFormats_GeneratorProducts_GenRunInfoProduct_h

class GenRunInfoProduct {
    public:
        // a few forward declarations
        struct XSec;

        // constructors, destructors
        GenRunInfoProduct();
        GenRunInfoProduct(const GenRunInfoProduct &other);
        virtual ~GenRunInfoProduct();

        // getters

        const XSec &internalXSec() const { return internalXSec_; }
        const XSec &externalXSecLO() const { return externalXSecLO_; }
        const XSec &externalXSecNLO() const { return externalXSecNLO_; }
        double filterEfficiency() const { return externalFilterEfficiency_; }

        // setters

        void setInternalXSec(const XSec &xsec) { internalXSec_ = xsec; }
        void setExternalXSecLO(const XSec &xsec) { externalXSecLO_ = xsec; }
        void setExternalXSecNLO(const XSec &xsec) { externalXSecNLO_ = xsec; }
        void setFilterEfficiency(double effic) { externalFilterEfficiency_ = effic; }

        // struct definitions
        struct XSec {
            public:
                XSec() : value_(-1.), error_(-1.) {}
                XSec(double value, double error = -1.) :
                        value_(value), error_(error) {}
                XSec(const XSec &other) :
                        value_(other.value_), error_(other.error_) {}

                double value() const { return value_; }
                double error() const { return error_; }

                bool isSet() const { return value_ >= 0.; }
                bool hasError() const { return error_ >= 0.; }

                operator double() const { return value_; }
                operator bool() const { return isSet(); }

                bool operator == (const XSec &other) const
                { return value_ == other.value_ && error_ == other.error_; }
                bool operator != (const XSec &other) const { return !(*this == other); }

            private:
                double value_, error_;
        };

        // convenience (return the value, prefer externally specified over internal one)
        double crossSection() const
        { return externalXSecLO_ ? externalXSecLO_.value() : internalXSec_.value(); }

        // methods used by EDM
        virtual bool mergeProduct(const GenRunInfoProduct &other);
        virtual bool isProductEqual(const GenRunInfoProduct &other) const;

    private:
        // cross sections
        XSec    internalXSec_;  // the one computed during cmsRun
        XSec    externalXSecLO_, externalXSecNLO_;      // from config file
        double  externalFilterEfficiency_;              // from config file
};

#endif // SimDataFormats_GeneratorProducts_GenRunInfoProduct_h