/data/refman/pasoursint/CMSSW_5_3_6/src/SimDataFormats/GeneratorProducts/src/GenRunInfoProduct.cc

Go to the documentation of this file.

#include <iostream>
#include <algorithm> 

#include "FWCore/MessageLogger/interface/MessageLogger.h"

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

using namespace edm;
using namespace std;

GenRunInfoProduct::GenRunInfoProduct() :
        externalFilterEfficiency_(-1.)
{
}

GenRunInfoProduct::GenRunInfoProduct(GenRunInfoProduct const &other) :
        internalXSec_(other.internalXSec_),
        externalXSecLO_(other.externalXSecLO_),
        externalXSecNLO_(other.externalXSecNLO_),
        externalFilterEfficiency_(other.externalFilterEfficiency_)
{
}

GenRunInfoProduct::~GenRunInfoProduct()
{
}

bool GenRunInfoProduct::mergeProduct(GenRunInfoProduct const &other)
{
        // if external xsec and filter efficiency are identical we allow merging
        // unfortunately we cannot merge on-the-fly xsec, as we do not have
        // information about the corresponding statistics (i.e. uncertainty)
        // so we just ignore it.
        // So in the end this merger is a dummy (just a safety check)

        if (externalXSecLO_ == other.externalXSecLO_ &&
            externalXSecNLO_ == other.externalXSecNLO_ &&
            externalFilterEfficiency_ == other.externalFilterEfficiency_)
                return true;

        edm::LogWarning("GenRunInfoProduct|ProductsNotMergeable")
                << "You are merging runs with different cross-sections and/or "
                   "filter efficiencies (from GenRunInfoProduct)\n"
                   "The resulting cross-section will not be consistent." << std::endl;

        return true;
}

bool GenRunInfoProduct::isProductEqual(GenRunInfoProduct const &other) const
{
        return externalXSecLO_ == other.externalXSecLO_ &&
               externalXSecNLO_ == other.externalXSecNLO_ &&
               externalFilterEfficiency_ == other.externalFilterEfficiency_;
}