HwDecayerBase.cc
No OneTemporary
Actions

Size

4 KB

Referenced Files

None

Subscribers

None

HwDecayerBase.cc
View Options

	// -- C++ --
	//
	// HwDecayerBase.cc is a part of Herwig++ - A multi-purpose Monte Carlo event generator
	// Copyright (C) 2002-2011 The Herwig Collaboration
	//
	// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
	// Please respect the MCnet academic guidelines, see GUIDELINES for details.
	//
	//
	// This is the implementation of the non-inlined, non-templated member
	// functions of the HwDecayerBase class.
	//

	#include "HwDecayerBase.h"
	#include "ThePEG/Interface/ClassDocumentation.h"
	#include "ThePEG/PDT/DecayMode.h"
	#include "ThePEG/Interface/Switch.h"
	#include "ThePEG/Persistency/PersistentOStream.h"
	#include "ThePEG/Persistency/PersistentIStream.h"
	#include "ThePEG/Repository/CurrentGenerator.h"
	#include "Herwig++/Shower/Base/HardTree.h"
	#include "Herwig++/Shower/Base/Branching.h"

	using namespace Herwig;

	bool HwDecayerBase::accept(const DecayMode & dm) const {
	// get the primary products
	tPDVector products=dm.orderedProducts();
	// add products for which the decay mode is all ready specified
	if(!dm.cascadeProducts().empty()) {
	for(ModeMSet::const_iterator mit=dm.cascadeProducts().begin();
	mit!=dm.cascadeProducts().end();++mit) {
	products.push_back(const_ptr_cast<PDPtr>((**mit).parent()));
	}
	}
	// can this mode be handled ?
	return accept(dm.parent(),products);
	}

	ParticleVector HwDecayerBase::decay(const DecayMode & dm,
	const Particle & p) const {
	// handling of the decay including the special features of the
	// DecayMode
	// get the primary products
	tPDVector products=dm.orderedProducts();
	// add products for which the decay mode is all ready specified
	if(!dm.cascadeProducts().empty()) {
	for(ModeMSet::const_iterator mit=dm.cascadeProducts().begin();
	mit!=dm.cascadeProducts().end();++mit) {
	products.push_back(const_ptr_cast<PDPtr>((**mit).parent()));
	}
	}
	// perform the primary decay
	ParticleVector output=decay(p,products);
	// perform the secondary decays
	if(!dm.cascadeProducts().empty()) {
	unsigned int iloc=dm.orderedProducts().size();
	for(ModeMSet::const_iterator mit=dm.cascadeProducts().begin();
	mit!=dm.cascadeProducts().end();++mit) {
	if(!(*mit)->decayer())
	throw Exception() << "Decay mode " << (**mit).tag()
	<< "does not have a decayer, can't perform"
	<< "decay in HwDecayerBase::decay()"
	<< Exception::eventerror;
	ParticleVector children=(mit)->decayer()->decay(mit,output[iloc]);
	for(unsigned int ix=0;ix<children.size();++ix) {
	output[iloc]->addChild(children[ix]);
	}
	++iloc;
	}
	}
	return output;
	}

	void HwDecayerBase::persistentOutput(PersistentOStream & os) const {
	os << _initialize << _dbOutput;
	}

	void HwDecayerBase::persistentInput(PersistentIStream & is, int) {
	is >> _initialize >> _dbOutput;
	}

	AbstractClassDescription<HwDecayerBase> HwDecayerBase::initHwDecayerBase;
	// Definition of the static class description member.

	void HwDecayerBase::Init() {

	static ClassDocumentation<HwDecayerBase> documentation
	("The HwDecayerBase class is the base class for Decayers in Hw++.");

	static Switch<HwDecayerBase,bool> interfaceInitialize
	("Initialize",
	"Initialization of the phase space calculation",
	&HwDecayerBase::_initialize, false, false, false);
	static SwitchOption interfaceInitializeon
	(interfaceInitialize,
	"Yes",
	"At initialisation find max weight and optimise the integration",
	true);
	static SwitchOption interfaceInitializeoff
	(interfaceInitialize,
	"No",
	"Use the maximum weight and channel weights supplied for the integration",
	false);

	static Switch<HwDecayerBase,bool> interfaceDatabaseOutput
	("DatabaseOutput",
	"Whether to print the database information",
	&HwDecayerBase::_dbOutput, false, false, false);
	static SwitchOption interfaceDatabaseOutputYes
	(interfaceDatabaseOutput,
	"Yes",
	"Output information on the decayer initialization",
	true);
	static SwitchOption interfaceDatabaseOutputNo
	(interfaceDatabaseOutput,
	"No",
	"Do not output information about the decayer initialization",
	false);
	}

	void HwDecayerBase::dofinish() {
	Decayer::dofinish();
	if(initialize() && databaseOutput()) {
	string fname = CurrentGenerator::current().filename() +
	string("-") + name() + string(".output");
	ofstream output(fname.c_str());
	dataBaseOutput(output,true);
	}
	}

	bool HwDecayerBase::softMatrixElementVeto(ShowerProgenitorPtr,
	ShowerParticlePtr,
	Branching) {
	return false;
	}

	HardTreePtr HwDecayerBase::generateHardest(ShowerTreePtr) {
	return HardTreePtr();
	}

File Metadata

Mime Type: text/x-c
Expires: Tue, Sep 30, 4:48 AM (10 h, 13 m)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 6480347
Default Alt Text: HwDecayerBase.cc (4 KB)

HwDecayerBase.ccNo OneTemporaryActions

HwDecayerBase.ccView Options

File Metadata

Event Timeline

HwDecayerBase.cc
No OneTemporary
Actions

HwDecayerBase.cc
View Options