Herwig7Hadronizer Class Reference

Inheritance diagram for Herwig7Hadronizer:

Public Member Functions
void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
const char *	classname () const
bool	decay ()
bool	declareSpecialSettings (const std::vector< std::string >)
bool	declareStableParticles (const std::vector< int > &pdgIds)
void	finalizeEvent ()
bool	generatePartonsAndHadronize ()
std::unique_ptr< GenLumiInfoHeader >	getGenLumiInfoHeader () const override
bool	hadronize ()
	Herwig7Hadronizer (const edm::ParameterSet &params)
bool	initializeForExternalPartons ()
bool	initializeForInternalPartons ()
void	randomizeIndex (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine)
bool	readSettings (int)
bool	residualDecay ()
void	statistics ()
	~Herwig7Hadronizer () override
Public Member Functions inherited from Herwig7Interface
	Herwig7Interface (const edm::ParameterSet &params)
void	setPEGRandomEngine (CLHEP::HepRandomEngine *)
	~Herwig7Interface () noexcept
Public Member Functions inherited from gen::BaseHadronizer
	BaseHadronizer (edm::ParameterSet const &ps)
void	cleanLHE ()
void	generateLHE (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine, unsigned int ncpu)
edm::Event &	getEDMEvent () const
std::unique_ptr< HepMC::GenEvent >	getGenEvent ()
std::unique_ptr< GenEventInfoProduct >	getGenEventInfo ()
GenRunInfoProduct &	getGenRunInfo ()
std::unique_ptr< lhef::LHEEvent >	getLHEEvent ()
const std::shared_ptr< lhef::LHERunInfo > &	getLHERunInfo () const
const std::string &	gridpackPath () const
int	randomIndex () const
const std::string &	randomInitConfigDescription () const
void	randomizeIndex (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine)
void	resetEvent (std::unique_ptr< HepMC::GenEvent > event)
void	resetEventInfo (std::unique_ptr< GenEventInfoProduct > eventInfo)
virtual bool	select (HepMC::GenEvent *) const
void	setEDMEvent (edm::Event &event)
void	setLHEEvent (std::unique_ptr< lhef::LHEEvent > event)
void	setLHERunInfo (std::unique_ptr< lhef::LHERunInfo > runInfo)
void	setRandomEngine (CLHEP::HepRandomEngine *v)
std::vector< std::string > const &	sharedResources () const
virtual	~BaseHadronizer () noexcept(false)

Private Member Functions
void	doSetRandomEngine (CLHEP::HepRandomEngine *v) override

Private Attributes
unsigned int	currentLumiBlock = 0
unsigned int	eventsToPrint
unsigned int	firstLumiBlock = 0
const std::string	handlerDirectory_
bool	haveEvt = false
edm::ParameterSet	paramSettings
std::shared_ptr< lhef::LHEProxy >	proxy_
const std::string	runFileName
ThePEG::EventPtr	thepegEvent

Additional Inherited Members
Public Attributes inherited from Herwig7Interface
ThePEG::EGPtr	eg_
Protected Member Functions inherited from Herwig7Interface
void	callHerwigGenerator ()
void	createInputFile (const edm::ParameterSet &params)
void	flushRandomNumberGenerator ()
bool	initGenerator ()
void	initRepository (const edm::ParameterSet &params)
Protected Member Functions inherited from gen::BaseHadronizer
std::unique_ptr< HepMC::GenEvent > &	event ()
std::unique_ptr< GenEventInfoProduct > &	eventInfo ()
lhef::LHEEvent *	lheEvent ()
lhef::LHERunInfo *	lheRunInfo ()
GenRunInfoProduct &	runInfo ()
Static Protected Member Functions inherited from Herwig7Interface
static std::auto_ptr< HepMC::GenEvent >	convert (const ThePEG::EventPtr &event)
static double	pthat (const ThePEG::EventPtr &event)
Protected Attributes inherited from Herwig7Interface
std::shared_ptr< Herwig::HerwigUIProvider >	HwUI_
std::auto_ptr< HepMC::IO_BaseClass >	iobc_
Protected Attributes inherited from gen::BaseHadronizer
std::string	lheFile_
int	randomIndex_

Detailed Description

Definition at line 36 of file Herwig7Hadronizer.cc.

Constructor & Destructor Documentation

Herwig7Hadronizer()

Herwig7Hadronizer::Herwig7Hadronizer

(

const edm::ParameterSet &

params

)

Definition at line 77 of file Herwig7Hadronizer.cc.

References Herwig7Interface::initRepository(), paramSettings, and muonDTDigis_cfi::pset.

    : Herwig7Interface(pset),
      BaseHadronizer(pset),
      eventsToPrint(pset.getUntrackedParameter<unsigned int>("eventsToPrint", 0)),
      handlerDirectory_(pset.getParameter<std::string>("eventHandlers")),
      runFileName(pset.getParameter<std::string>("run")) {
  initRepository(pset);
  paramSettings = pset;
}

~Herwig7Hadronizer()

Herwig7Hadronizer::~Herwig7Hadronizer ( )

override

Definition at line 87 of file Herwig7Hadronizer.cc.

{}

Member Function Documentation

beginLuminosityBlock()

void Herwig7Hadronizer::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

classname()

const char* Herwig7Hadronizer::classname ( ) const

inline

Definition at line 55 of file Herwig7Hadronizer.cc.

{ return "Herwig7Hadronizer"; }

decay()

bool Herwig7Hadronizer::decay

(

)

Definition at line 207 of file Herwig7Hadronizer.cc.

{ return true; }

declareSpecialSettings()

bool Herwig7Hadronizer::declareSpecialSettings ( const std::vector< std::string >  )

inline

Definition at line 45 of file Herwig7Hadronizer.cc.

{ return true; }

declareStableParticles()

bool Herwig7Hadronizer::declareStableParticles

(

const std::vector< int > &

pdgIds

)

Definition at line 119 of file Herwig7Hadronizer.cc.

{ return false; }

doSetRandomEngine()

void Herwig7Hadronizer::doSetRandomEngine ( CLHEP::HepRandomEngine *  v )

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 61 of file Herwig7Hadronizer.cc.

References Herwig7Interface::setPEGRandomEngine(), and findQualityFiles::v.

{ setPEGRandomEngine(v); }

finalizeEvent()

void Herwig7Hadronizer::finalizeEvent

(

)

Definition at line 192 of file Herwig7Hadronizer.cc.

References gen::BaseHadronizer::event(), gen::BaseHadronizer::eventInfo(), eventsToPrint, Herwig7Interface::iobc_, Herwig7Interface::pthat(), and thepegEvent.

                                      {
  eventInfo() = std::make_unique<GenEventInfoProduct>(event().get());
  eventInfo()->setBinningValues(std::vector<double>(1, pthat(thepegEvent)));

  if (eventsToPrint) {
    eventsToPrint--;
    event()->print();
  }

  if (iobc_.get())
    iobc_->write_event(event().get());

  edm::LogInfo("Generator|Herwig7Hadronizer") << "Event produced";
}

generatePartonsAndHadronize()

bool Herwig7Hadronizer::generatePartonsAndHadronize

(

)

Definition at line 128 of file Herwig7Hadronizer.cc.

References Herwig7Interface::convert(), Herwig7Interface::eg_, gen::BaseHadronizer::event(), cppFunctionSkipper::exception, and thepegEvent.

                                                    {
  edm::LogInfo("Generator|Herwig7Hadronizer") << "Start production";

  try {
    thepegEvent = eg_->shoot();
  } catch (std::exception& exc) {
    edm::LogWarning("Generator|Herwig7Hadronizer")
        << "EGPtr::shoot() thrown an exception, event skipped: " << exc.what();
    return false;
  }

  if (!thepegEvent) {
    edm::LogWarning("Generator|Herwig7Hadronizer") << "thepegEvent not initialized";
    return false;
  }

  event() = convert(thepegEvent);
  if (!event().get()) {
    edm::LogWarning("Generator|Herwig7Hadronizer") << "genEvent not initialized";
    return false;
  }

  return true;
}

getGenLumiInfoHeader()

std::unique_ptr< GenLumiInfoHeader > Herwig7Hadronizer::getGenLumiInfoHeader ( ) const

overridevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 211 of file Herwig7Hadronizer.cc.

References genWeightsTable_cfi::genLumiInfoHeader, and thepegEvent.

                                                                               {
  auto genLumiInfoHeader = BaseHadronizer::getGenLumiInfoHeader();

  if (thepegEvent) {
    int weights_number = thepegEvent->optionalWeights().size();

    if (weights_number > 1) {
      genLumiInfoHeader->weightNames().reserve(weights_number + 1);
      genLumiInfoHeader->weightNames().push_back("nominal");
      std::map<std::string, double> weights_map = thepegEvent->optionalWeights();
      for (std::map<std::string, double>::iterator it = weights_map.begin(); it != weights_map.end(); it++) {
        genLumiInfoHeader->weightNames().push_back(it->first);
      }
    }
  }

  return genLumiInfoHeader;
}

hadronize()

bool Herwig7Hadronizer::hadronize

(

)

Definition at line 153 of file Herwig7Hadronizer.cc.

References Herwig7Interface::convert(), Herwig7Interface::eg_, gen::BaseHadronizer::event(), lhef::LHEEvent::evtnum(), cppFunctionSkipper::exception, haveEvt, gen::BaseHadronizer::lheEvent(), and thepegEvent.

                                  {
  if (!haveEvt) {
    try {
      thepegEvent = eg_->shoot();
      haveEvt = true;
    } catch (std::exception& exc) {
      edm::LogWarning("Generator|Herwig7Hadronizer")
          << "EGPtr::shoot() thrown an exception, event skipped: " << exc.what();
      return false;
    }
  }
  int evtnum = lheEvent()->evtnum();
  if (evtnum == -1) {
    edm::LogError("Generator|Herwig7Hadronizer")
        << "Event number not set in lhe file, needed for correctly aligning Herwig and LHE events!";
    return false;
  }
  if (thepegEvent->number() < evtnum) {
    edm::LogError("Herwig7 interface") << "Herwig does not seem to be generating events in order, did you set "
                                          "/Herwig/EventHandlers/FxFxLHReader:AllowedToReOpen Yes?";
    return false;
  } else if (thepegEvent->number() == evtnum) {
    haveEvt = false;
    if (!thepegEvent) {
      edm::LogWarning("Generator|Herwig7Hadronizer") << "thepegEvent not initialized";
      return false;
    }

    event() = convert(thepegEvent);
    if (!event().get()) {
      edm::LogWarning("Generator|Herwig7Hadronizer") << "genEvent not initialized";
      return false;
    }
    return true;
  }
  edm::LogWarning("Generator|Herwig7Hadronizer") << "Event " << evtnum << " not generated (likely skipped in merging)";
  return false;
}

initializeForExternalPartons()

bool Herwig7Hadronizer::initializeForExternalPartons

(

)

Definition at line 104 of file Herwig7Hadronizer.cc.

References currentLumiBlock, beamvalidation::exit(), firstLumiBlock, Herwig7Interface::initGenerator(), Herwig7Interface::initRepository(), paramSettings, and runFileName.

                                                     {
  if (currentLumiBlock == firstLumiBlock) {
    std::ifstream runFile(runFileName + ".run");
    if (runFile.fail())  //required for showering of LHE files
    {
      initRepository(paramSettings);
    }
    if (!initGenerator()) {
      edm::LogInfo("Generator|Herwig7Hadronizer") << "No run step for Herwig chosen. Program will be aborted.";
      exit(0);
    }
  }
  return true;
}

initializeForInternalPartons()

bool Herwig7Hadronizer::initializeForInternalPartons

(

)

Definition at line 89 of file Herwig7Hadronizer.cc.

References currentLumiBlock, beamvalidation::exit(), firstLumiBlock, Herwig7Interface::initGenerator(), Herwig7Interface::initRepository(), paramSettings, and runFileName.

                                                     {
  if (currentLumiBlock == firstLumiBlock) {
    std::ifstream runFile(runFileName + ".run");
    if (runFile.fail())  //required for showering of LHE files
    {
      initRepository(paramSettings);
    }
    if (!initGenerator()) {
      edm::LogInfo("Generator|Herwig7Hadronizer") << "No run step for Herwig chosen. Program will be aborted.";
      exit(0);
    }
  }
  return true;
}

randomizeIndex()

void Herwig7Hadronizer::randomizeIndex	(	edm::LuminosityBlock const &	lumi,
		CLHEP::HepRandomEngine *	rengine
	)

Definition at line 230 of file Herwig7Hadronizer.cc.

References currentLumiBlock, and firstLumiBlock.

                                                                                                    {
  BaseHadronizer::randomizeIndex(lumi, rengine);

  if (firstLumiBlock == 0) {
    firstLumiBlock = lumi.id().luminosityBlock();
  }
  currentLumiBlock = lumi.id().luminosityBlock();
}

readSettings()

bool Herwig7Hadronizer::readSettings ( int  )

inline

Definition at line 41 of file Herwig7Hadronizer.cc.

{ return true; }

residualDecay()

bool Herwig7Hadronizer::residualDecay

(

)

Definition at line 209 of file Herwig7Hadronizer.cc.

{ return true; }

statistics()

void Herwig7Hadronizer::statistics

(

)

Definition at line 121 of file Herwig7Hadronizer.cc.

References Herwig7Interface::eg_, gen::BaseHadronizer::runInfo(), and GenRunInfoProduct::setInternalXSec().

                                   {
  if (eg_) {
    runInfo().setInternalXSec(
        GenRunInfoProduct::XSec(eg_->integratedXSec() / ThePEG::picobarn, eg_->integratedXSecErr() / ThePEG::picobarn));
  }
}

Member Data Documentation

currentLumiBlock

unsigned int Herwig7Hadronizer::currentLumiBlock = 0

private

Definition at line 74 of file Herwig7Hadronizer.cc.

Referenced by initializeForExternalPartons(), initializeForInternalPartons(), and randomizeIndex().

eventsToPrint

unsigned int Herwig7Hadronizer::eventsToPrint

private

Definition at line 63 of file Herwig7Hadronizer.cc.

Referenced by finalizeEvent().

firstLumiBlock

unsigned int Herwig7Hadronizer::firstLumiBlock = 0

private

Definition at line 73 of file Herwig7Hadronizer.cc.

Referenced by initializeForExternalPartons(), initializeForInternalPartons(), and randomizeIndex().

handlerDirectory_

const std::string Herwig7Hadronizer::handlerDirectory_

private

Definition at line 69 of file Herwig7Hadronizer.cc.

haveEvt

bool Herwig7Hadronizer::haveEvt = false

private

Definition at line 66 of file Herwig7Hadronizer.cc.

Referenced by hadronize().

paramSettings

edm::ParameterSet Herwig7Hadronizer::paramSettings

private

Definition at line 70 of file Herwig7Hadronizer.cc.

Referenced by Herwig7Hadronizer(), initializeForExternalPartons(), and initializeForInternalPartons().

proxy_

std::shared_ptr<lhef::LHEProxy> Herwig7Hadronizer::proxy_

private

Definition at line 68 of file Herwig7Hadronizer.cc.

runFileName

const std::string Herwig7Hadronizer::runFileName

private

Definition at line 71 of file Herwig7Hadronizer.cc.

Referenced by initializeForExternalPartons(), and initializeForInternalPartons().

thepegEvent

ThePEG::EventPtr Herwig7Hadronizer::thepegEvent

private

Definition at line 65 of file Herwig7Hadronizer.cc.

Referenced by finalizeEvent(), generatePartonsAndHadronize(), getGenLumiInfoHeader(), and hadronize().