GenEventInfoProduct.h

Go to the documentation of this file.

#ifndef SimDataFormats_GeneratorProducts_GenEventInfoProduct_h
#define SimDataFormats_GeneratorProducts_GenEventInfoProduct_h
 
#include <vector>
#include <memory>
 
#include "SimDataFormats/GeneratorProducts/interface/PdfInfo.h"
 
namespace HepMC {
  class GenEvent;
}  // namespace HepMC
 
class GenEventInfoProduct {
public:
  GenEventInfoProduct();
  GenEventInfoProduct(const HepMC::GenEvent *evt);
  GenEventInfoProduct(const GenEventInfoProduct &other);
  GenEventInfoProduct(GenEventInfoProduct &&other);
  virtual ~GenEventInfoProduct();
 
  GenEventInfoProduct &operator=(const GenEventInfoProduct &other);
  GenEventInfoProduct &operator=(GenEventInfoProduct &&other);
 
  typedef gen::PdfInfo PDF;
 
  // getters
 
  std::vector<double> &weights() { return weights_; }
  const std::vector<double> &weights() const { return weights_; }
 
  double weight() const { return weights_.empty() ? 1.0 : weights_[0]; }
 
  double weightProduct() const;
 
  unsigned int signalProcessID() const { return signalProcessID_; }
 
  double qScale() const { return qScale_; }
  double alphaQCD() const { return alphaQCD_; }
  double alphaQED() const { return alphaQED_; }
 
  const PDF *pdf() const { return pdf_.get(); }
  bool hasPDF() const { return pdf() != nullptr; }
 
  const std::vector<double> &binningValues() const { return binningValues_; }
  bool hasBinningValues() const { return !binningValues_.empty(); }
 
  const std::vector<float> &DJRValues() const { return DJRValues_; }
  bool hasDJRValues() const { return !DJRValues_.empty(); }
 
  int nMEPartons() const { return nMEPartons_; }
 
  int nMEPartonsFiltered() const { return nMEPartonsFiltered_; }
 
  // setters
 
  void setWeights(const std::vector<double> &weights) { weights_ = weights; }
 
  void setSignalProcessID(unsigned int procID) { signalProcessID_ = procID; }
 
  void setScales(double q = -1., double qcd = -1., double qed = -1.) { qScale_ = q, alphaQCD_ = qcd, alphaQED_ = qed; }
 
  void setPDF(const PDF *pdf) { pdf_.reset(pdf ? new PDF(*pdf) : nullptr); }
 
  void setBinningValues(const std::vector<double> &values) { binningValues_ = values; }
 
  void setDJR(const std::vector<float> &values) { DJRValues_ = values; }
 
  void setNMEPartons(int n) { nMEPartons_ = n; }
 
  void setNMEPartonsFiltered(int n) { nMEPartonsFiltered_ = n; }
 
private:
  // HepMC::GenEvent provides a list of weights
  std::vector<double> weights_;
 
  // generator-dependent process ID
  unsigned int signalProcessID_;
 
  // information about scales
  double qScale_;
  double alphaQCD_, alphaQED_;
 
  // optional PDF info
  std::unique_ptr<PDF> pdf_;
 
  // If event was produced in bis, this contains
  // the values that were used to define which
  // bin the event belongs in
  // This replaces the genEventScale, which only
  // corresponds to Pythia pthat.  The RunInfo
  // will contain the information what physical
  // quantity these values actually belong to
  std::vector<double> binningValues_;
  std::vector<float> DJRValues_;
  int nMEPartons_;
  int nMEPartonsFiltered_;
};
 
#endif  // SimDataFormats_GeneratorProducts_GenEventInfoProduct_h