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< HepMC3::GenEvent >	getGenEvent3 ()

std::unique_ptr< GenEventInfoProduct >	getGenEventInfo ()

std::unique_ptr< GenEventInfoProduct3 >	getGenEventInfo3 ()

GenRunInfoProduct &	getGenRunInfo ()

std::unique_ptr< lhef::LHEEvent >	getLHEEvent ()

const std::shared_ptr< lhef::LHERunInfo > &	getLHERunInfo () const

unsigned int	getVHepMC ()

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	resetEvent3 (std::unique_ptr< HepMC3::GenEvent > event3)

void	resetEventInfo (std::unique_ptr< GenEventInfoProduct > eventInfo)

void	resetEventInfo3 (std::unique_ptr< GenEventInfoProduct3 > 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< HepMC3::GenEvent > &	event3 ()

std::unique_ptr< GenEventInfoProduct > &	eventInfo ()

std::unique_ptr< GenEventInfoProduct3 > &	eventInfo3 ()

lhef::LHEEvent *	lheEvent ()

lhef::LHERunInfo *	lheRunInfo ()

GenRunInfoProduct &	runInfo ()

Static Protected Member Functions inherited from Herwig7Interface
static std::unique_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::unique_ptr< HepMC::IO_BaseClass >	iobc_

Protected Attributes inherited from gen::BaseHadronizer
unsigned int	ivhepmc = 2

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.

87 {}

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.

55 { return "Herwig7Hadronizer"; }

◆ decay()

bool Herwig7Hadronizer::decay ( )

Definition at line 207 of file Herwig7Hadronizer.cc.

207 { return true; }

◆ declareSpecialSettings()

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

inline

Definition at line 45 of file Herwig7Hadronizer.cc.

45 { return true; }

◆ declareStableParticles()

bool Herwig7Hadronizer::declareStableParticles ( const std::vector< int > & pdgIds )

Definition at line 119 of file Herwig7Hadronizer.cc.

119 { 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.

61 { setPEGRandomEngine(v); }

Herwig7Interface::setPEGRandomEngine

void setPEGRandomEngine(CLHEP::HepRandomEngine *)

Definition: Herwig7Interface.cc:73

findQualityFiles.v

v

Definition: findQualityFiles.py:179

◆ 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, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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.

41 { return true; }

◆ residualDecay()

bool Herwig7Hadronizer::residualDecay ( )

Definition at line 209 of file Herwig7Hadronizer.cc.

209 { 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().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Herwig7Hadronizer()

◆ ~Herwig7Hadronizer()

Member Function Documentation

◆ beginLuminosityBlock()

◆ classname()

◆ decay()

◆ declareSpecialSettings()

◆ declareStableParticles()

◆ doSetRandomEngine()

◆ finalizeEvent()

◆ generatePartonsAndHadronize()

◆ getGenLumiInfoHeader()

◆ hadronize()

◆ initializeForExternalPartons()

◆ initializeForInternalPartons()

◆ randomizeIndex()

◆ readSettings()

◆ residualDecay()

◆ statistics()

Member Data Documentation

◆ currentLumiBlock

◆ eventsToPrint

◆ firstLumiBlock

◆ handlerDirectory_

◆ haveEvt

◆ paramSettings

◆ proxy_

◆ runFileName

◆ thepegEvent