GenLumiInfoProduct.cc

Go to the documentation of this file.

#include <iostream>
#include <algorithm>
#include <map>
#include <utility>

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

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

using namespace edm;
using namespace std;

const bool operator<(const GenLumiInfoProduct::ProcessInfo& lhs, const GenLumiInfoProduct::ProcessInfo& rhs) {
  return (lhs.process() < rhs.process());
}

const bool operator!=(const GenLumiInfoProduct::ProcessInfo& lhs, const GenLumiInfoProduct::ProcessInfo& rhs) {
  bool condition = (lhs.process() != rhs.process()) || (lhs.lheXSec() != rhs.lheXSec());
  return condition;
}

const bool operator==(const GenLumiInfoProduct::ProcessInfo& lhs, const GenLumiInfoProduct::ProcessInfo& rhs) {
  bool condition =
      (lhs.process() == rhs.process() && lhs.lheXSec() == rhs.lheXSec() && lhs.nPassPos() == rhs.nPassPos() &&
       lhs.nPassNeg() == rhs.nPassNeg() && lhs.nTotalPos() == rhs.nTotalPos() && lhs.nTotalNeg() == rhs.nTotalNeg() &&
       lhs.tried() == rhs.tried() && lhs.selected() == rhs.selected() && lhs.killed() == rhs.killed() &&
       lhs.accepted() == rhs.accepted() && lhs.acceptedBr() == rhs.acceptedBr());
  return condition;
}

const bool operator!=(const GenLumiInfoProduct& lhs, const GenLumiInfoProduct& rhs) {
  std::vector<GenLumiInfoProduct::ProcessInfo> lhsVector = lhs.getProcessInfos();
  std::vector<GenLumiInfoProduct::ProcessInfo> rhsVector = rhs.getProcessInfos();
  std::sort(lhsVector.begin(), lhsVector.end());
  std::sort(rhsVector.begin(), rhsVector.end());
  unsigned int lhssize = lhsVector.size();
  unsigned int rhssize = rhsVector.size();
  bool condition = (lhs.getHEPIDWTUP() != rhs.getHEPIDWTUP()) || (lhssize != rhssize);
  bool fail = false;
  if (!condition) {
    for (unsigned int i = 0; i < lhssize; i++) {
      if (lhsVector[i] != rhsVector[i]) {
        fail = true;
        break;
      }
    }
  }
  return (condition || fail);
}

const bool operator==(const GenLumiInfoProduct& lhs, const GenLumiInfoProduct& rhs) {
  std::vector<GenLumiInfoProduct::ProcessInfo> lhsVector = lhs.getProcessInfos();
  std::vector<GenLumiInfoProduct::ProcessInfo> rhsVector = rhs.getProcessInfos();
  std::sort(lhsVector.begin(), lhsVector.end());
  std::sort(rhsVector.begin(), rhsVector.end());
  unsigned int lhssize = lhsVector.size();
  unsigned int rhssize = rhsVector.size();

  bool condition = (lhs.getHEPIDWTUP() == rhs.getHEPIDWTUP()) && (lhssize == rhssize);
  unsigned int passCounts = -999;
  if (condition) {
    for (unsigned int i = 0; i < lhssize; i++) {
      if (lhsVector[i] == rhsVector[i])
        passCounts++;
    }
  }
  return (condition && (passCounts == lhssize));
}

GenLumiInfoProduct::GenLumiInfoProduct() : hepidwtup_(-1) { internalProcesses_.clear(); }

GenLumiInfoProduct::GenLumiInfoProduct(const int id) : hepidwtup_(id) { internalProcesses_.clear(); }

GenLumiInfoProduct::GenLumiInfoProduct(GenLumiInfoProduct const& other)
    : hepidwtup_(other.hepidwtup_), internalProcesses_(other.internalProcesses_) {}

GenLumiInfoProduct::~GenLumiInfoProduct() {}

bool GenLumiInfoProduct::mergeProduct(GenLumiInfoProduct const& other) {
  std::map<int, ProcessInfo> processes;

  for (unsigned int i = 0; i < getProcessInfos().size(); i++) {
    int id = getProcessInfos()[i].process();
    ProcessInfo& x = processes[id];
    x = getProcessInfos()[i];
  }

  // the two GenLuminInfoProducts may not have the same number
  // of processes
  for (unsigned int i = 0; i < other.getProcessInfos().size(); i++) {
    int id = other.getProcessInfos()[i].process();
    ProcessInfo& x = processes[id];
    if (x.lheXSec().value() > 0)
      x.addOthers(other.getProcessInfos()[i]);
    else
      x = other.getProcessInfos()[i];
  }

  internalProcesses_.resize(processes.size());
  unsigned int i = 0;
  for (std::map<int, ProcessInfo>::const_iterator iter = processes.begin(); iter != processes.end(); ++iter, i++)
    internalProcesses_[i] = iter->second;
  return true;
}

void GenLumiInfoProduct::swap(GenLumiInfoProduct& other) {
  std::swap(hepidwtup_, other.hepidwtup_);
  internalProcesses_.swap(other.internalProcesses_);
}

bool GenLumiInfoProduct::isProductEqual(GenLumiInfoProduct const& other) const { return ((*this) == other); }