No OneTemporary
Actions

Size

46 KB

Subscribers

None

View Options

	diff --git a/Models/Susy/SusyBase.cc b/Models/Susy/SusyBase.cc
	--- a/Models/Susy/SusyBase.cc
	+++ b/Models/Susy/SusyBase.cc
	@@ -1,763 +1,797 @@
	// -- C++ --
	//
	// SusyBase.cc is a part of Herwig++ - A multi-purpose Monte Carlo event generator
	// Copyright (C) 2002-2007 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 SusyBase class.
	//

	#include "SusyBase.h"
	#include "ThePEG/Interface/ClassDocumentation.h"
	#include "ThePEG/Persistency/PersistentOStream.h"
	#include "ThePEG/Persistency/PersistentIStream.h"
	#include "ThePEG/Interface/Switch.h"
	#include "ThePEG/Interface/Reference.h"
	#include "ThePEG/Interface/Parameter.h"
	#include "ThePEG/Repository/Repository.h"
	#include "ThePEG/Utilities/StringUtils.h"
	#include "ThePEG/Repository/CurrentGenerator.h"
	#include "ThePEG/PDT/EnumParticles.h"
	#include "ThePEG/PDT/MassGenerator.h"
	#include "ThePEG/PDT/WidthGenerator.h"
	#include "ThePEG/PDT/DecayMode.h"

	using namespace Herwig;

	SusyBase::SusyBase() : readFile_(false), topModesFromFile_(false),
	tolerance_(1e-6), MPlanck_(2.4e18*GeV),
	gravitino_(false),
	tanBeta_(0), mu_(ZERO),
	M1_(ZERO), M2_(ZERO), M3_(ZERO),
	mH12_(ZERO),mH22_(ZERO),
	meL_(ZERO),mmuL_(ZERO),mtauL_(ZERO),
	meR_(ZERO),mmuR_(ZERO),mtauR_(ZERO),
	mq1L_(ZERO),mq2L_(ZERO),mq3L_(ZERO),
	mdR_(ZERO),muR_(ZERO),msR_(ZERO),
	mcR_(ZERO),mbR_(ZERO),mtR_(ZERO),
	gluinoPhase_(1.)
	{}

	IBPtr SusyBase::clone() const {
	return new_ptr(*this);
	}

	IBPtr SusyBase::fullclone() const {
	return new_ptr(*this);
	}

	void SusyBase::doinit() {
	addVertex(WSFSFVertex_);
	addVertex(NFSFVertex_);
	addVertex(GFSFVertex_);
	addVertex(HSFSFVertex_);
	addVertex(CFSFVertex_);
	addVertex(GSFSFVertex_);
	addVertex(GGSQSQVertex_);
	addVertex(GSGSGVertex_);
	addVertex(NNZVertex_);
	if(NNPVertex_) addVertex(NNPVertex_);
	if(GNGVertex_) addVertex(GNGVertex_);
	addVertex(CCZVertex_);
	addVertex(CNWVertex_);
	addVertex(GOGOHVertex_);
	addVertex(WHHVertex_);
	addVertex(NCTVertex_);
	if(gravitino_) {
	if(GVNHVertex_) addVertex(GVNHVertex_);
	if(GVNVVertex_) addVertex(GVNVVertex_);
	if(GVFSVertex_) addVertex(GVFSVertex_);
	}
	StandardModel::doinit();
	}

	void SusyBase::persistentOutput(PersistentOStream & os) const {
	os << readFile_ << topModesFromFile_ << gravitino_
	<< NMix_ << UMix_ << VMix_ << WSFSFVertex_
	<< NFSFVertex_ << GFSFVertex_ << HSFSFVertex_ << CFSFVertex_
	<< GSFSFVertex_ << GGSQSQVertex_ << GSGSGVertex_
	<< NNZVertex_ << NNPVertex_ << CCZVertex_ << CNWVertex_
	<< GOGOHVertex_ << WHHVertex_ << GNGVertex_ << NCTVertex_
	<< GVNHVertex_ << GVNVVertex_ << GVFSVertex_
	<< tanBeta_ << ounit(mu_,GeV)
	<< ounit(M1_,GeV) << ounit(M2_,GeV) << ounit(M3_,GeV)
	<< ounit(mH12_,GeV2) << ounit(mH22_,GeV2)
	<< ounit(meL_,GeV) << ounit(mmuL_,GeV) << ounit(mtauL_,GeV)
	<< ounit(meR_,GeV) << ounit(mmuR_,GeV) << ounit(mtauR_,GeV)
	<< ounit(mq1L_,GeV) << ounit(mq2L_,GeV) << ounit(mq3L_,GeV)
	<< ounit(mdR_,GeV) << ounit(muR_,GeV) << ounit(msR_,GeV)
	<< ounit(mcR_,GeV) << ounit(mbR_,GeV) << ounit(mtR_,GeV)
	<< gluinoPhase_ << tolerance_ << ounit(MPlanck_,GeV);
	}

	void SusyBase::persistentInput(PersistentIStream & is, int) {
	is >> readFile_ >> topModesFromFile_ >> gravitino_
	>> NMix_ >> UMix_ >> VMix_ >> WSFSFVertex_
	>> NFSFVertex_ >> GFSFVertex_ >> HSFSFVertex_ >> CFSFVertex_
	>> GSFSFVertex_ >> GGSQSQVertex_ >> GSGSGVertex_
	>> NNZVertex_ >> NNPVertex_ >> CCZVertex_ >> CNWVertex_
	>> GOGOHVertex_ >> WHHVertex_ >> GNGVertex_ >> NCTVertex_
	>> GVNHVertex_ >> GVNVVertex_ >> GVFSVertex_
	>> tanBeta_ >> iunit(mu_,GeV)
	>> iunit(M1_,GeV) >> iunit(M2_,GeV) >> iunit(M3_,GeV)
	>> iunit(mH12_,GeV2) >> iunit(mH22_,GeV2)
	>> iunit(meL_,GeV) >> iunit(mmuL_,GeV) >> iunit(mtauL_,GeV)
	>> iunit(meR_,GeV) >> iunit(mmuR_,GeV) >> iunit(mtauR_,GeV)
	>> iunit(mq1L_,GeV) >> iunit(mq2L_,GeV) >> iunit(mq3L_,GeV)
	>> iunit(mdR_,GeV) >> iunit(muR_,GeV) >> iunit(msR_,GeV)
	>> iunit(mcR_,GeV) >> iunit(mbR_,GeV) >> iunit(mtR_,GeV)
	>> gluinoPhase_ >> tolerance_ >> iunit(MPlanck_,GeV);
	}

	ClassDescription<SusyBase> SusyBase::initSusyBase;
	// Definition of the static class description member.

	void SusyBase::Init() {

	static ClassDocumentation<SusyBase> documentation
	("This is the base class for any SUSY model.",
	"SUSY spectrum files follow the Les Houches accord"
	" \\cite{Skands:2003cj,Allanach:2008qq}.",
	" %\\cite{Skands:2003cj}\n"
	"\\bibitem{Skands:2003cj}\n"
	" P.~Skands {\\it et al.},\n"
	" ``SUSY Les Houches accord: Interfacing SUSY spectrum calculators, decay\n"
	" %packages, and event generators,''\n"
	" JHEP {\\bf 0407}, 036 (2004)\n"
	" [arXiv:hep-ph/0311123].\n"
	" %%CITATION = JHEPA,0407,036;%%\n"
	"%\\cite{Allanach:2008qq}\n"
	"\\bibitem{Allanach:2008qq}\n"
	" B.~Allanach {\\it et al.},\n"
	" %``SUSY Les Houches Accord 2,''\n"
	" Comput.\\ Phys.\\ Commun.\\ {\\bf 180}, 8 (2009)\n"
	" [arXiv:0801.0045 [hep-ph]].\n"
	" %%CITATION = CPHCB,180,8;%%\n"
	);

	static Switch<SusyBase,bool> interfaceTopModes
	("TopModes",
	"Whether ro use the Herwig++ SM top decays or those from the SLHA file",
	&SusyBase::topModesFromFile_, false, false, false);
	static SwitchOption interfaceTopModesFile
	(interfaceTopModes,
	"File",
	"Take the modes from the files",
	true);
	static SwitchOption interfaceTopModesHerwig
	(interfaceTopModes,
	"Herwig",
	"Use the SM ones",
	false);

	static Reference<SusyBase,Helicity::AbstractVSSVertex> interfaceVertexWSS
	("Vertex/WSFSF",
	"Reference to Susy W SF SF vertex",
	&SusyBase::WSFSFVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractFFSVertex> interfaceVertexNFSF
	("Vertex/NFSF",
	"Reference to the neutralino-fermion-sfermion vertex",
	&SusyBase::NFSFVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractFFSVertex> interfaceVertexGFSF
	("Vertex/GFSF",
	"Reference to the gluino-fermion-sfermion vertex",
	&SusyBase::GFSFVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractSSSVertex> interfaceVertexHSFSF
	("Vertex/HSFSF",
	"Reference to the Higgs-fermion-sfermion vertex",
	&SusyBase::HSFSFVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractFFSVertex> interfaceVertexCFSF
	("Vertex/CFSF",
	"Reference to the chargino-fermion-sfermion vertex",
	&SusyBase::CFSFVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractVSSVertex> interfaceVertexGSFSF
	("Vertex/GSFSF",
	"Reference to the gluon-sfermion-sfermion vertex",
	&SusyBase::GSFSFVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractVVSSVertex> interfaceVertexGGSS
	("Vertex/GGSQSQ",
	"Reference to the gluon-gluon-squark-squark vertex",
	&SusyBase::GGSQSQVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractFFVVertex> interfaceVertexGSGSG
	("Vertex/GSGSG",
	"Reference to the gluon-gluino-gluino vertex",
	&SusyBase::GSGSGVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractFFVVertex> interfaceVertexNNZ
	("Vertex/NNZ",
	"Reference to Z-~chi_i0-~chi_i0 vertex",
	&SusyBase::NNZVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractFFVVertex> interfaceVertexNNP
	("Vertex/NNP",
	"Reference to photon-~chi_i0-~chi_i0 vertex",
	&SusyBase::NNPVertex_, false, false, true, true, false);

	static Reference<SusyBase,Helicity::AbstractFFVVertex> interfaceVertexGNG
	("Vertex/GNG",
	"Reference to gluon-~chi_i0-gluino vertex",
	&SusyBase::GNGVertex_, false, false, true, true, false);

	static Reference<SusyBase,Helicity::AbstractFFVVertex> interfaceVertexCCZ
	("Vertex/CCZ",
	"Reference to ~chi_i+-~chi_i-Z vertex",
	&SusyBase::CCZVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractFFVVertex> interfaceVertexCNW
	("Vertex/CNW",
	"Reference to ~chi_i+-chi_i0-W vertex",
	&SusyBase::CNWVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractFFSVertex> interfaceVertexGOGOH
	("Vertex/GOGOH",
	"Reference to the gaugino-gaugino-higgs vertex",
	&SusyBase::GOGOHVertex_, false, false, true, false);

	static Reference<SusyBase,Helicity::AbstractVSSVertex> interfaceVertexWHH
	("Vertex/SSWHH",
	"Reference to Susy WHHVertex",
	&SusyBase::WHHVertex_, false, false, true, false);

	static Reference<SusyBase,AbstractFFSVertex> interfaceVertexNCT
	("Vertex/NCT",
	"Vertex for the flavour violating coupling of the top squark "
	"to the neutralino and charm quark.",
	&SusyBase::NCTVertex_, false, false, true, false, false);

	static Reference<SusyBase,AbstractRFSVertex> interfaceVertexGVNH
	("Vertex/GVNH",
	"Vertex for the interfaction of the gravitino-neutralino"
	" and Higgs bosons",
	&SusyBase::GVNHVertex_, false, false, true, true, false);

	static Reference<SusyBase,AbstractRFVVertex> interfaceVertexGVNV
	("Vertex/GVNV",
	"Vertex for the interfaction of the gravitino-neutralino"
	" and vector bosons",
	&SusyBase::GVNVVertex_, false, false, true, true, false);

	static Reference<SusyBase,AbstractRFSVertex> interfaceVertexGVFS
	("Vertex/GVFS",
	"Vertex for the interfaction of the gravitino-fermion"
	" and sfermion",
	&SusyBase::GVFSVertex_, false, false, true, true, false);

	static Parameter<SusyBase,double> interfaceBRTolerance
	("BRTolerance",
	"Tolerance for the sum of branching ratios to be difference from one.",
	&SusyBase::tolerance_, 1e-6, 1e-8, 0.01,
	false, false, Interface::limited);

	static Parameter<SusyBase,Energy> interfaceMPlanck
	("MPlanck",
	"The Planck mass for GMSB models",
	&SusyBase::MPlanck_, GeV, 2.4e18GeV, 1.e16GeV, 1.e20*GeV,
	false, false, Interface::limited);

	}

	void SusyBase::readSetup(istream & is) {
	string filename = dynamic_ptr_cast<istringstream*>(&is)->str();
	if(readFile_)
	throw SetupException()
	<< "A second SLHA file " << filename << " has been opened."
	<< "This is probably unintended and as it can cause crashes"
	<< " and other unpredictable behaviour it is not allowed."
	<< Exception::runerror;
	- ifstream file(filename.c_str());
	- if( !file ) throw SetupException()
	- << "SusyBase::readSetup - An error occurred in opening the "
	- << "spectrum file \"" << filename << "\". A SUSY model cannot be "
	- << "run without this."
	- << Exception::runerror;
	-
	+ CFileLineReader cfile;
	+ cfile.open(filename);
	+ if( !cfile ) throw SetupException()
	+ << "SusyBase::readSetup - An error occurred in opening the "
	+ << "spectrum file \"" << filename << "\". A SUSY model cannot be "
	+ << "run without this."
	+ << Exception::runerror;
	useMe();
	//Before reading the spectrum/decay files the SM higgs
	//decay modes, mass and width generators need to be turned off.
	PDPtr h0 = getParticleData(ParticleID::h0);
	h0->widthGenerator(WidthGeneratorPtr());
	h0->massGenerator(MassGenPtr());
	DecaySet::const_iterator dit = h0->decayModes().begin();
	DecaySet::const_iterator dend = h0->decayModes().end();
	for( ; dit != dend; ++dit ) {
	const InterfaceBase * ifb =
	BaseRepository::FindInterface(*dit, "BranchingRatio");
	ifb->exec(**dit, "set", "0.0");
	ifb = BaseRepository::FindInterface(*dit, "OnOff");
	ifb->exec(**dit, "set", "Off");
	}
	// if taking the top modes from the file
	// delete the SM stuff
	if(topModesFromFile_) {
	PDPtr top = getParticleData(ParticleID::t);
	top->widthGenerator(WidthGeneratorPtr());
	top->massGenerator(MassGenPtr());
	DecaySet::const_iterator dit = top->decayModes().begin();
	DecaySet::const_iterator dend = top->decayModes().end();
	for( ; dit != dend; ++dit ) {
	const InterfaceBase * ifb =
	BaseRepository::FindInterface(*dit, "BranchingRatio");
	ifb->exec(**dit, "set", "0.0");
	ifb = BaseRepository::FindInterface(*dit, "OnOff");
	ifb->exec(**dit, "set", "Off");
	- }
	-
	+ }
	}
	- string line;
	+ // read first line and chjeck if this is a Les Houches event file
	+ cfile.readline();
	+ bool lesHouches = cfile.find("<LesHouchesEvents");
	+ bool reading = !lesHouches;
	+ if(lesHouches) cfile.readline();
	//function pointer for putting all characters to lower case.
	int (*pf)(int) = tolower;
	- //Start to read in the file
	- while(getline(file, line)) {
	- // ignore comment lines
	- if(line[0] == '#') continue;
	- // make everything lower case
	- transform(line.begin(), line.end(), line.begin(), pf);
	- // start of a block
	- if(line.find("block") == 0) {
	- string name = StringUtils::car(StringUtils::cdr(line), " #");
	- // mixing matrix
	- if((name.find("mix") != string::npos &&
	- name.find("hmix") != 0)) {
	- unsigned int row(0),col(0);
	- MixingVector vals = readMatrix(file,row,col);
	- mixings_[name] = make_pair(make_pair(row,col),vals);
	- }
	- else if(name.find("au") == 0\|\|name.find("ad") == 0\|\|
	- name.find("ae") == 0 ) {
	- string test = StringUtils::car(line, "#");
	- while (test.find("=")!= string::npos) {
	- test = StringUtils::cdr(test, "=");
	- }
	- istringstream is(test);
	- double scale;
	- is >> scale;
	- if(scale>1e10) continue;
	- unsigned int row(0),col(0);
	- MixingVector vals = readMatrix(file,row,col);
	- mixings_[name] = make_pair(make_pair(row,col),vals);
	- }
	- else if( name.find("info") == string::npos)
	- readBlock(file,name,line);
	- }
	- // decays
	- else if( line.find("decay") == 0 )
	- readDecay(file, line);
	- }
	- // extract the relevant parameters
	- extractParameters();
	- // create the mixing matrices we need
	- createMixingMatrices();
	- // set the masses, this has to be done after the
	- // mixing matrices have been created
	- resetRepositoryMasses();
	- // have now read the file
	- readFile_=true;
	+ do {
	+ string line = cfile.getline();
	+ // check for start of slha block in SLHA files
	+ if(lesHouches && !reading) {
	+ if(line.find("<slha")==0) reading = true;
	+ if(!cfile.readline()) break;
	+ continue;
	+ }
	+ // ignore comment lines
	+ if(line[0] == '#') {
	+ if(!cfile.readline()) break;
	+ continue;
	+ }
	+ // make everything lower case
	+ transform(line.begin(), line.end(), line.begin(), pf);
	+ // start of a block
	+ if(line.find("block") == 0) {
	+ string name = StringUtils::car(StringUtils::cdr(line), " #");
	+ // mixing matrix
	+ if((name.find("mix") != string::npos &&
	+ name.find("hmix") != 0)) {
	+ unsigned int row(0),col(0);
	+ MixingVector vals = readMatrix(cfile,row,col);
	+ mixings_[name] = make_pair(make_pair(row,col),vals);
	+ }
	+ else if(name.find("au") == 0\|\|name.find("ad") == 0\|\|
	+ name.find("ae") == 0 ) {
	+ string test = StringUtils::car(line, "#");
	+ while (test.find("=")!= string::npos) {
	+ test = StringUtils::cdr(test, "=");
	+ }
	+ istringstream is(test);
	+ double scale;
	+ is >> scale;
	+ if(scale>1e10) continue;
	+ unsigned int row(0),col(0);
	+ MixingVector vals = readMatrix(cfile,row,col);
	+ mixings_[name] = make_pair(make_pair(row,col),vals);
	+ }
	+ else if( name.find("info") == string::npos) {
	+ readBlock(cfile,name,line);
	+ }
	+ else {
	+ if(!cfile.readline()) break;
	+ }
	+ continue;
	+ }
	+ // decays
	+ else if( line.find("decay") == 0 ) {
	+ readDecay(cfile, line);
	+ continue;
	+ }
	+ else if( lesHouches && line.find("</slha") == 0 ) {
	+ reading = false;
	+ break;
	+ }
	+ if(!cfile.readline()) break;
	+ }
	+ while(true);
	+ // extract the relevant parameters
	+ extractParameters();
	+ // create the mixing matrices we need
	+ createMixingMatrices();
	+ // set the masses, this has to be done after the
	+ // mixing matrices have been created
	+ resetRepositoryMasses();
	+ // have now read the file
	+ readFile_=true;
	}

	-void SusyBase::readBlock(ifstream & ifs,string name,string linein) {
	- if(!ifs)
	+void SusyBase::readBlock(CFileLineReader & cfile,string name,string linein) {
	+ if(!cfile)
	throw SetupException()
	<< "SusyBase::readBlock() - The input stream is in a bad state"
	<< Exception::runerror;
	- string line;
	+ // storage or the parameters
	+ string test = StringUtils::car(linein, "#");
	ParamMap store;
	bool set = true;
	- string test = StringUtils::car(linein, "#");
	// special for the alpha block
	if(name.find("alpha") == 0 ) {
	double alpha;
	- getline(ifs, line);
	+ cfile.readline();
	+ string line = cfile.getline();
	istringstream iss(line);
	iss >> alpha;
	store.insert(make_pair(1,alpha));
	}
	else {
	// extract the scale from the block if present
	if(test.find("=")!= string::npos) {
	while(test.find("=")!=string::npos)
	test= StringUtils::cdr(test,"=");
	istringstream is(test);
	double scale;
	is >> scale;
	// only store the lowest scale block
	if(parameters_.find(name)!=parameters_.end()) {
	set = scale < parameters_[name][-1];
	}
	else {
	store.insert(make_pair(-1,scale));
	}
	}
	- while(getline(ifs, line)) {
	- if(line[0] == '#') {
	- if( ifs.peek() == 'D' \|\| ifs.peek() == 'B' \|\|
	- ifs.peek() == 'd' \|\| ifs.peek() == 'b' ) break;
	- else continue;
	+ while(cfile.readline()) {
	+ string line = cfile.getline();
	+ // skip comments
	+ if(line[0] == '#') continue;
	+ // reached the end
	+ if( line[0] == 'B' \|\| line[0] == 'b' \|\|
	+ line[0] == 'D' \|\| line[0] == 'd' \|\|
	+ line[0] == '<' ) {
	+ cfile.resetline();
	+ break;
	}
	istringstream is(line);
	long index;
	double value;
	is >> index >> value;
	store.insert(make_pair(index, value));
	- if(ifs.peek() == 'B' \|\| ifs.peek() == 'D' \|\|
	- ifs.peek() == 'b' \|\| ifs.peek() == 'd' \|\|
	- ifs.peek() == '#') break;
	}
	}
	if(set) parameters_[name]=store;
	}

	-void SusyBase::readDecay(ifstream & ifs,
	+void SusyBase::readDecay(CFileLineReader & cfile,
	string decay) const{
	- if(!ifs)
	+ if(!cfile)
	throw SetupException()
	<<"SusyBase::readDecay - The input stream is in a bad state";
	- //if the particle is stable then next line will be the start of a
	- //new decaymode or a new block so just return if this is found
	- if( ifs.peek() == 'D' \|\| ifs.peek() == 'B' \|\|
	- ifs.peek() == 'd' \|\| ifs.peek() == 'b' ) return;
	+ // extract parent PDG code and width
	+ long parent(0);
	+ Energy width(ZERO);
	istringstream iss(decay);
	string dummy;
	- long parent(0);
	- Energy width(ZERO);
	iss >> dummy >> parent >> iunit(width, GeV);
	PDPtr inpart = getParticleData(parent);
	- if(!topModesFromFile_&&abs(parent)==ParticleID::t) return;
	- if(!inpart) {
	- throw SetupException()
	- << "SusyBase::readDecay() - A ParticleData object with the PDG code "
	- << parent << " does not exist. " << Exception::runerror;
	+ if(!topModesFromFile_&&abs(parent)==ParticleID::t) {
	+ cfile.readline();
	return;
	}
	+ if(!inpart) throw SetupException()
	+ << "SusyBase::readDecay() - "
	+ << "A ParticleData object with the PDG code "
	+ << parent << " does not exist. "
	+ << Exception::runerror;
	inpart->width(width);
	if( width > ZERO ) inpart->cTau(hbarc/width);
	inpart->widthCut(5.*width);
	- string prefix("decaymode " + inpart->name() + "->"), tag(""),line("");
	+ string prefix("decaymode " + inpart->name() + "->");
	double brsum(0.);
	unsigned int nmode = 0;
	- while(getline(ifs, line)) {
	- if(line[0] == '#') {
	- if( ifs.peek() == 'D' \|\| ifs.peek() == 'B' \|\|
	- ifs.peek() == 'd' \|\| ifs.peek() == 'b' ) break;
	- else continue;
	+ while(cfile.readline()) {
	+ string line = cfile.getline();
	+ // skip comments
	+ if(line[0] == '#') continue;
	+ // reached the end
	+ if( line[0] == 'B' \|\| line[0] == 'b' \|\|
	+ line[0] == 'D' \|\| line[0] == 'd' \|\|
	+ line[0] == '<' ) {
	+ cfile.resetline();
	+ break;
	}
	+ // read the mode
	+ string tag = prefix;
	istringstream is(line);
	double brat(0.);
	unsigned int nda(0),npr(0);
	is >> brat >> nda;
	brsum += brat;
	++nmode;
	while( true ) {
	long t;
	is >> t;
	if( is.fail() ) break;
	if( t == abs(parent) ) {
	throw SetupException()
	<< "An error occurred while read a decay of the "
	<< inpart->PDGName() << ". One of its products has the same PDG code "
	<< "as the parent particle. Please check the SLHA file.\n"
	<< Exception::runerror;
	}
	tcPDPtr p = getParticleData(t);
	if( !p ) {
	throw SetupException()
	<< "SusyBase::readDecay() - An unknown PDG code has been encounterd "
	<< "while reading a decay mode. ID: " << t
	<< Exception::runerror;
	}
	else {
	++npr;
	tag += p->name() + ",";
	}
	}
	if( npr != nda ) {
	throw SetupException()
	<< "SusyBase::readDecay - While reading a decay of the "
	<< inpart->PDGName() << " from an SLHA file, an inconsistency "
	<< "between the number of decay products and the value in "
	<< "the 'NDA' column was found. Please check if the spectrum "
	<< "file is correct.\n"
	<< Exception::warning;
	}
	if( npr > 1 ) {
	inpart->stable(false);
	tag.replace(tag.size() - 1, 1, ";");
	- createDecayMode(prefix + tag, brat);
	- tag.clear();
	+ createDecayMode(tag, brat);
	}
	- if( ifs.peek() == 'D' \|\| ifs.peek() == 'B' \|\|
	- ifs.peek() == 'd' \|\| ifs.peek() == 'b' ) break;
	}
	if( abs(brsum - 1.) > tolerance_ && nmode!=0 ) {
	cerr << "Warning: The total branching ratio for " << inpart->PDGName()
	<< " from the spectrum file does not sum to 1. The branching fractions"
	<< " will be rescaled.\n";
	cerr << setprecision(13) << abs(brsum - 1.) << "\n";
	}
	}

	const MixingVector
	-SusyBase::readMatrix(ifstream & ifs,
	+SusyBase::readMatrix(CFileLineReader & cfile,
	unsigned int & row, unsigned int & col) {
	- if(!ifs)
	+ if(!cfile)
	throw SetupException()
	<< "SusyBase::readMatrix() - The input stream is in a bad state."
	<< Exception::runerror;
	- string line;
	unsigned int rowmax(0), colmax(0);
	MixingVector values;
	- while(getline(ifs, line)) {
	- if(line[0] == '#') {
	- if( ifs.peek() == 'D' \|\| ifs.peek() == 'B' \|\|
	- ifs.peek() == 'd' \|\| ifs.peek() == 'b' ) break;
	- else continue;
	+ while(cfile.readline()) {
	+ string line = cfile.getline();
	+ // skip comments
	+ if(line[0] == '#') continue;
	+ // reached the end
	+ if( line[0] == 'B' \|\| line[0] == 'b' \|\|
	+ line[0] == 'D' \|\| line[0] == 'd' \|\|
	+ line[0] == '<' ) {
	+ cfile.resetline();
	+ break;
	}
	istringstream is(line);
	unsigned int index1, index2;
	double real(0.), imag(0.);
	is >> index1 >> index2 >> real >> imag;
	values.push_back(MixingElement(index1,index2,Complex(real, imag)));
	- if(index1 > rowmax) rowmax = index1;
	- if(index2 > colmax) colmax = index2;
	- if( ifs.peek() == 'B' \|\| ifs.peek() == 'D' \|\|
	- ifs.peek() == 'b' \|\| ifs.peek() == 'd' ) break;
	+ if(index1 > rowmax) rowmax = index1;
	+ if(index2 > colmax) colmax = index2;
	+
	}
	col=colmax;
	row=rowmax;
	return values;
	}

	void SusyBase::createDecayMode(string tag, double brat) const {
	ostringstream cmd;
	cmd << tag << string(" ")
	<< setprecision(13) << brat << string(" 1 /Herwig/Decays/Mambo");
	Repository::exec(cmd.str(), cerr);
	}

	void SusyBase::createMixingMatrix(MixingMatrixPtr & matrix,
	string name, const MixingVector & values,
	MatrixSize size) {
	matrix = new_ptr(MixingMatrix(size.first,size.second));
	for(unsigned int ix=0; ix < values.size(); ++ix)
	(*matrix)(values[ix].row-1,values[ix].col-1) = values[ix].value;

	if(name == "nmix") {
	vector<long> ids(4);
	ids[0] = 1000022; ids[1] = 1000023;
	ids[2] = 1000025; ids[3] = 1000035;
	matrix->setIds(ids);
	}
	else if(name == "nmnmix") {
	vector<long> ids(5);
	ids[0] = 1000022; ids[1] = 1000023;
	ids[2] = 1000025; ids[3] = 1000035;
	ids[4] = 1000045;
	matrix->setIds(ids);
	}
	else if(name == "umix") {
	vector<long> ids(2);
	ids[0] = 1000024; ids[1] = 1000037;
	matrix->setIds(ids);
	}
	else if(name == "vmix") {
	vector<long> ids(2);
	ids[0] = 1000024; ids[1] = 1000037;
	matrix->setIds(ids);
	}
	else if(name == "stopmix") {
	vector<long> ids(2);
	ids[0] = 1000006; ids[1] = 2000006;
	matrix->setIds(ids);
	}
	else if(name == "sbotmix") {
	vector<long> ids(2);
	ids[0] = 1000005; ids[1] = 2000005;
	matrix->setIds(ids);
	}
	else if(name == "staumix") {
	vector<long> ids(2);
	ids[0] = 1000015; ids[1] = 2000015;
	matrix->setIds(ids);
	}
	else if(name == "nmhmix") {
	vector<long> ids(3);
	ids[0] = 25; ids[1] = 35; ids[2] = 45;
	matrix->setIds(ids);
	}
	else if(name == "nmamix") {
	vector<long> ids(2);
	ids[0] = 36; ids[1] = 46;
	matrix->setIds(ids);
	}
	else
	throw SetupException() << "Cannot find correct title for mixing matrix "
	<< name << Exception::runerror;
	}

	void SusyBase::resetRepositoryMasses() {
	map<string,ParamMap>::const_iterator fit=parameters_.find("mass");
	if(fit==parameters_.end())
	throw Exception() << "BLOCK MASS not found in input file"
	<< " can't set masses of SUSY particles"
	<< Exception::runerror;
	ParamMap theMasses = fit->second;
	for(ParamMap::iterator it = theMasses.begin(); it != theMasses.end();
	++it) {
	long id = it->first;
	double mass = it->second;
	//a negative mass requires an adjustment to the
	//associated mixing matrix by a factor of i
	if(mass < 0.0) adjustMixingMatrix(id);
	PDPtr part = getParticleData(id);
	if(!part) throw SetupException()
	<< "SusyBase::resetRepositoryMasses() - Particle with PDG code " << id
	<< " not found." << Exception::warning;
	//Find interface nominal mass interface
	const InterfaceBase * ifb = BaseRepository::FindInterface(part, "NominalMass");
	ostringstream os;
	os << abs(it->second);
	ifb->exec(*part, "set", os.str());
	// switch on gravitino interactions?
	gravitino_ \|= id== ParticleID::SUSY_Gravitino;
	}
	theMasses.clear();
	}

	void SusyBase::adjustMixingMatrix(long id) {
	//get correct mixing matrix
	switch(id) {
	case 1000021 :
	gluinoPhase_ = Complex(0.,1.);
	break;
	case 1000022 :
	case 1000023 :
	case 1000025 :
	case 1000035 :
	case 1000045 :
	if(NMix_)
	NMix_->adjustPhase(id);
	else
	throw SetupException() << "SusyBase::adjustMixingMatrix - "
	<< "The neutralino mixing matrix pointer "
	<< "is null!" << Exception::runerror;
	break;
	case 1000024 :
	case 1000037 :
	if(UMix_)
	UMix_->adjustPhase(id);
	else
	throw SetupException() << "SusyBase::adjustMixingMatrix - "
	<< "The U-Type chargino mixing matrix pointer "
	<< "is null!" << Exception::runerror;
	if(VMix_)
	VMix_->adjustPhase(id);
	else
	throw SetupException() << "SusyBase::adjustMixingMatrix - "
	<< "The V-Type chargino mixing matrix pointer "
	<< "is null!" << Exception::runerror;
	break;
	default :
	throw SetupException()
	<< "SusyBase::adjustMixingMatrix - Trying to adjust mixing matrix "
	<< "phase for a particle that does not have a mixing matrix "
	<< "associated with it. " << id << " must have a negative mass in "
	<< "the spectrum file, this should only occur for particles that mix."
	<< Exception::runerror;
	}
	}

	void SusyBase::createMixingMatrices() {
	map<string,pair<MatrixSize, MixingVector > >::const_iterator it;
	for(it=mixings_.begin();it!=mixings_.end();++it) {
	string name=it->first;
	// create the gaugino mixing matrices
	if(name == "nmix" \|\| name == "nmnmix" ){
	createMixingMatrix(NMix_,name,it->second.second,it->second.first);
	}
	else if (name == "umix" ) {
	createMixingMatrix(UMix_,name,it->second.second,it->second.first);
	}
	else if (name == "vmix") {
	createMixingMatrix(VMix_,name,it->second.second,it->second.first);
	}
	}
	}

	void SusyBase::extractParameters(bool checkmodel) {
	map<string,ParamMap>::const_iterator pit;
	ParamMap::const_iterator it;
	// try and get tan beta from extpar first
	pit=parameters_.find("extpar");
	// extract tan beta
	tanBeta_ = -1.;
	if(pit!=parameters_.end()) {
	it = pit->second.find(25);
	if(it!=pit->second.end()) tanBeta_ = it->second;
	}
	// otherwise from minpar
	if(tanBeta_<0.) {
	pit=parameters_.find("minpar");
	if(pit!=parameters_.end()) {
	it = pit->second.find(3);
	if(it!=pit->second.end()) tanBeta_ = it->second;
	}
	}
	if(tanBeta_<0.)
	throw Exception() << "Can't find tan beta in BLOCK MINPAR"
	<< " or BLOCK EXTPAR " << Exception::runerror;
	// extract parameters from hmix
	pit=parameters_.find("hmix");
	if(pit==parameters_.end()) {
	cerr << "BLOCK HMIX not found setting mu to zero\n";
	mu_=ZERO;
	}
	else {
	it = pit->second.find(1);
	if(it==pit->second.end()) {
	cerr << "mu not found in BLOCK HMIX setting to zero\n";
	}
	else {
	mu_=it->second*GeV;
	}
	}
	pit = parameters_.find("msoft");
	if( pit == parameters_.end() )
	throw Exception() << "BLOCK MSOFT not found in "
	<< "SusyBase::extractParameters()"
	<< Exception::runerror;
	M1_ = pit->second.find(1 )->second*GeV;
	M2_ = pit->second.find(2 )->second*GeV;
	M3_ = pit->second.find(3 )->second*GeV;
	mH12_ = pit->second.find(21)->second*GeV2;
	mH22_ = pit->second.find(22)->second*GeV2;
	meL_ = pit->second.find(31)->second*GeV;
	mmuL_ = pit->second.find(32)->second*GeV;
	mtauL_ = pit->second.find(33)->second*GeV;
	meR_ = pit->second.find(34)->second*GeV;
	mmuR_ = pit->second.find(35)->second*GeV;
	mtauR_ = pit->second.find(36)->second*GeV;
	mq1L_ = pit->second.find(41)->second*GeV;
	mq2L_ = pit->second.find(42)->second*GeV;
	mq3L_ = pit->second.find(43)->second*GeV;
	muR_ = pit->second.find(44)->second*GeV;
	mcR_ = pit->second.find(45)->second*GeV;
	mtR_ = pit->second.find(46)->second*GeV;
	mdR_ = pit->second.find(47)->second*GeV;
	msR_ = pit->second.find(48)->second*GeV;
	mbR_ = pit->second.find(49)->second*GeV;
	if(checkmodel) {
	throw Exception() << "The SusyBase class should not be used as a "
	<< "Model class, use one of the models which inherit"
	<< " from it" << Exception::runerror;
	}
	}
	diff --git a/Models/Susy/SusyBase.h b/Models/Susy/SusyBase.h
	--- a/Models/Susy/SusyBase.h
	+++ b/Models/Susy/SusyBase.h
	@@ -1,712 +1,713 @@
	// -- C++ --
	//
	// SusyBase.h is a part of Herwig++ - A multi-purpose Monte Carlo event generator
	// Copyright (C) 2002-2007 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.
	//
	#ifndef HERWIG_SusyBase_H
	#define HERWIG_SusyBase_H
	//
	// This is the declaration of the SusyBase class.
	//

	#include "Herwig++/Models/StandardModel/StandardModel.h"
	#include "MixingMatrix.h"
	+#include "ThePEG/Utilities/CFileLineReader.h"
	#include "ThePEG/Helicity/Vertex/AbstractVSSVertex.h"
	#include "ThePEG/Helicity/Vertex/AbstractSSSVertex.h"
	#include "ThePEG/Helicity/Vertex/AbstractVVSSVertex.h"
	#include "ThePEG/Helicity/Vertex/AbstractRFSVertex.h"
	#include "ThePEG/Helicity/Vertex/AbstractRFVVertex.h"
	#include "SusyBase.fh"

	namespace Herwig {
	using namespace ThePEG;

	/@name Some convenient typedefs. /
	//@{
	/**
	* Map to hold key, parameter pairs.
	*/
	typedef map<long, double> ParamMap;
	//@}

	/** \ingroup Models
	* This class is designed to be a base class for SUSY models. There is
	* an interface to set the name of the spectrum file to read in
	* the necessary parameters for a SUSY model.
	*
	* @see \ref SusyBaseInterfaces "The interfaces"
	* defined for SusyBase.
	* @see StandardModel
	*/

	class SusyBase: public StandardModel {

	public:

	/**
	* The default constructor.
	*/
	SusyBase();

	public:

	/** @name Access functions. */
	//@{
	/**
	* Value of \f$\tan\beta\f$.
	*/
	double tanBeta() const { return tanBeta_; }

	/**
	* Value of \f$\mu\f$ parameter.
	*/
	Energy muParameter() const { return mu_; }

	/**
	* The neutralino mixing matrix
	*/
	const MixingMatrixPtr & neutralinoMix() const {
	return NMix_;
	}

	/**
	* The U-type chargino mixing matrix
	*/
	const MixingMatrixPtr & charginoUMix() const {
	return UMix_;
	}

	/**
	* The V-type chargino mixing matrix
	*/
	const MixingMatrixPtr & charginoVMix() const {
	return VMix_;
	}

	/**
	* The phase for gluino vertices
	*/
	const Complex & gluinoPhase() const {return gluinoPhase_;}
	//@}

	public:

	/** @name Functions used by the persistent I/O system. */
	//@{
	/**
	* Function used to write out object persistently.
	* @param os the persistent output stream written to.
	*/
	void persistentOutput(PersistentOStream & os) const;

	/**
	* Function used to read in object persistently.
	* @param is the persistent input stream read from.
	* @param version the version number of the object when written.
	*/
	void persistentInput(PersistentIStream & is, int version);
	//@}

	/**
	* The standard Init function used to initialize the interfaces.
	* Called exactly once for each class by the class description system
	* before the main function starts or
	* when this class is dynamically loaded.
	*/
	static void Init();

	public:

	/**
	* Soft breaking parameters
	*/
	//@{
	/**
	* The bilinear breaking mass term for the bino
	*/
	const Energy & M1() const {return M1_;}

	/**
	* The bilinear breaking mass term for the wino
	*/
	const Energy & M2() const {return M2_;}

	/**
	* The bilinear breaking mass term for the gluinos
	*/
	const Energy & M3() const {return M3_;}

	/**
	* The soft breaking mass squared for \f$H_1\f$
	*/
	const Energy2 & Mh12() const {return mH12_;}

	/**
	* The soft breaking mass squared for \f$H_2\f$
	*/
	const Energy2 & Mh22() const {return mH22_;}

	/**
	* Soft breaking mass for the first generation lepton doublet
	*/
	const Energy & MeL() const {return meL_;}

	/**
	* Soft breaking mass for the second generation lepton doublet
	*/
	const Energy & MmuL() const {return mmuL_;}

	/**
	* Soft breaking mass for the third generation lepton doublet
	*/
	const Energy & MtauL() const {return mtauL_;}

	/**
	* Soft breaking mass for the first generation lepton singlet
	*/
	const Energy & MeR() const {return meR_;}

	/**
	* Soft breaking mass for the second generation lepton singlet
	*/
	const Energy & MmuR() const {return mmuR_;}

	/**
	* Soft breaking mass for the third generation lepton singlet
	*/
	const Energy & MtauR() const {return mtauR_;}

	/**
	* Soft breaking mass for the first generation quark doublet
	*/
	const Energy & Mq1L() const {return mq1L_;}

	/**
	* Soft breaking mass for the second generation quark doublet
	*/
	const Energy & Mq2L() const {return mq2L_;}

	/**
	* Soft breaking mass for the third generation quark doublet
	*/
	const Energy & Mq3L() const {return mq3L_;}

	/**
	* Soft breaking mass for the down singlet
	*/
	const Energy & MdR() const {return mdR_;}

	/**
	* Soft breaking mass for the up singlet
	*/
	const Energy & MuR() const {return muR_;}

	/**
	* Soft breaking mass for the strange singlet
	*/
	const Energy & MsR() const {return msR_;}

	/**
	* Soft breaking mass for the charm singlet
	*/
	const Energy & McR() const {return mcR_;}

	/**
	* Soft breaking mass for the bottom singlet
	*/
	const Energy & MbR() const {return mbR_;}

	/**
	* Soft breaking mass for the top singlet
	*/
	const Energy & MtR() const {return mtR_;}
	//@}

	/**
	* Planck mass
	*/
	const Energy & MPlanck() const {return MPlanck_;}

	protected:

	/**
	* Function to read information from a setup file.
	* @param is istream object to read file.
	*/
	virtual void readSetup(istream & is);

	private:

	/*@name Functions to help file read-in. /
	//@{
	/**
	* Read block from LHA file
	* @param ifs input stream containg data
	* @param name The name of the block
	*/
	- void readBlock(ifstream & ifs,string name,string line);
	+ void readBlock(CFileLineReader & ifs,string name,string line);

	/**
	* Function to read mixing matrix from LHA file
	* @param ifs input stream containg data
	* @param row Number of rows
	* @param col Number of columns
	*/
	- const MixingVector readMatrix(ifstream & ifs, unsigned int & row,
	+ const MixingVector readMatrix(CFileLineReader & ifs, unsigned int & row,
	unsigned int & col);

	/**
	* Read decaymodes from LHA file
	* @param ifs input stream containg data
	* @param decay string containing name of parent and value of total width
	*/
	- void readDecay(ifstream & ifs, string decay) const;
	+ void readDecay(CFileLineReader & ifs, string decay) const;

	/**
	* Create a DecayMode object in the repository
	* @param tag The tag identifying the decay mode including the prefix
	* 'decaymode'
	* @param brat Branching ratio of this mode
	*/
	void createDecayMode(string tag, double brat) const;

	protected:

	/**
	* Create the mixing matrices for the model
	*/
	virtual void createMixingMatrices();

	/**
	* Extract the parameters from the input blocks
	*/
	virtual void extractParameters(bool checkmodel=true);

	/**
	* Create a object MixingMatrix in the repository
	* @param matrix Pointer to the mixing matrix
	* @param name Name of the mixing matrix, i.e. nmix, umix...
	* @param values Value of each entry in the matrix
	* @param size The size of the matrix
	*/
	void createMixingMatrix(MixingMatrixPtr & matrix, string name,
	const MixingVector & values,
	MatrixSize size);

	/**
	* Reset masses in the repository to values read from LHA file.
	*/
	void resetRepositoryMasses();

	/**
	* Adjust row of Mixing Matrix if a negative mass occurs in LHA file
	* @param id The PDG code of the particle with a negative mass
	*/
	virtual void adjustMixingMatrix(long id);
	//@}

	/**
	* Access to the mixings and parameters for the inheriting classes
	*/
	//@{
	/**
	* Parameter blocks
	*/
	const map<string,ParamMap> & parameters() const {
	return parameters_;
	}

	/**
	* Mixing blocks
	*/
	const map<string,pair<MatrixSize,MixingVector> > & mixings() const {
	return mixings_;
	}
	//@}

	/**
	* Reset neutralino mixing matrix
	*/
	void neutralinoMix(MixingMatrixPtr nm) { NMix_ = nm; }

	/**
	* Reset the U-type chargino mixing matrix
	*/
	void charginoUMix(MixingMatrixPtr um) { UMix_ = um; }

	/**
	* Reset the V-type chargino mixing matrix
	*/
	void charginoVMix(MixingMatrixPtr vm) { VMix_ = vm; }

	protected:

	/** @name Clone Methods. */
	//@{
	/**
	* Make a simple clone of this object.
	* @return a pointer to the new object.
	*/
	virtual IBPtr clone() const;

	/** Make a clone of this object, possibly modifying the cloned object
	* to make it sane.
	* @return a pointer to the new object.
	*/
	virtual IBPtr fullclone() const;
	//@}

	protected:

	/** @name Standard Interfaced functions. */
	//@{
	/**
	* Initialize this object after the setup phase before saving an
	* EventGenerator to disk.
	* @throws InitException if object could not be initialized properly.
	*/
	virtual void doinit();
	//@}

	private:

	/**
	* The static object used to initialize the description of this class.
	* Indicates that this is a concrete class with persistent data.
	*/
	static ClassDescription<SusyBase> initSusyBase;

	/**
	* The assignment operator is private and must never be called.
	* In fact, it should not even be implemented.
	*/
	SusyBase & operator=(const SusyBase &);

	private:

	/**
	* Whether or not the SLHA fiel has been read
	*/
	bool readFile_;

	/**
	* Whether or not to replace the top decay modes with those from
	* the SLHA files
	*/
	bool topModesFromFile_;

	/**
	* Tolerance for branching ratios
	*/
	double tolerance_;

	/**
	* Planck mass needed in GMSB models
	*/
	Energy MPlanck_;

	/**
	* Whether or not to include gravitino interactions
	*/
	bool gravitino_;

	/*
	* Storage of the parameters.
	*/
	//@{
	/**
	* Parameter blocks
	*/
	map<string,ParamMap> parameters_;

	/**
	* Mixing blocks
	*/
	map<string,pair<MatrixSize, MixingVector> > mixings_;

	/**
	* \f$\tan\beta\f$
	*/
	double tanBeta_;

	/**
	* \f$\mu\f$
	*/
	Energy mu_;
	//@}

	/**
	* Soft breaking parameters
	*/
	//@{
	/**
	* The bilinear breaking mass term for the bino
	*/
	Energy M1_;

	/**
	* The bilinear breaking mass term for the wino
	*/
	Energy M2_;

	/**
	* The bilinear breaking mass term for the gluinos
	*/
	Energy M3_;

	/**
	* The soft breaking mass squared for \f$H_1\f$
	*/
	Energy2 mH12_;

	/**
	* The soft breaking mass squared for \f$H_2\f$
	*/
	Energy2 mH22_;

	/**
	* Soft breaking mass for the first generation lepton doublet
	*/
	Energy meL_;

	/**
	* Soft breaking mass for the second generation lepton doublet
	*/
	Energy mmuL_;

	/**
	* Soft breaking mass for the third generation lepton doublet
	*/
	Energy mtauL_;

	/**
	* Soft breaking mass for the first generation lepton singlet
	*/
	Energy meR_;

	/**
	* Soft breaking mass for the second generation lepton singlet
	*/
	Energy mmuR_;

	/**
	* Soft breaking mass for the third generation lepton singlet
	*/
	Energy mtauR_;

	/**
	* Soft breaking mass for the first generation quark doublet
	*/
	Energy mq1L_;

	/**
	* Soft breaking mass for the second generation quark doublet
	*/
	Energy mq2L_;

	/**
	* Soft breaking mass for the third generation quark doublet
	*/
	Energy mq3L_;

	/**
	* Soft breaking mass for the down singlet
	*/
	Energy mdR_;

	/**
	* Soft breaking mass for the up singlet
	*/
	Energy muR_;

	/**
	* Soft breaking mass for the strange singlet
	*/
	Energy msR_;

	/**
	* Soft breaking mass for the charm singlet
	*/
	Energy mcR_;

	/**
	* Soft breaking mass for the bottom singlet
	*/
	Energy mbR_;

	/**
	* Soft breaking mass for the top singlet
	*/
	Energy mtR_;
	//@}

	/**
	* Phase for the gluino
	*/
	Complex gluinoPhase_;

	/**
	* Neutralino and Chargino mixing matrices
	*/
	//@{
	/**
	* Store pointers to the gaugino mixing matrices
	*/
	//@{
	/**
	* The neutralino mixing matrix
	*/
	MixingMatrixPtr NMix_;

	/**
	* The \f$U\f$ mixing matrix for the charginos
	*/
	MixingMatrixPtr UMix_;

	/**
	* The \f$V\f$ mixing matrix for the charginos
	*/
	MixingMatrixPtr VMix_;
	//@}

	/*@name Vertex pointers. /
	//@{
	/**
	* Pointer to the gauge boson sfermion-sfermion vertex
	*/
	AbstractVSSVertexPtr WSFSFVertex_;

	/**
	* Pointer to the neutralino-fermion-sfermion vertex
	*/
	AbstractFFSVertexPtr NFSFVertex_;

	/**
	* Pointer to the gluino-fermion-sfermion coupling
	*/
	AbstractFFSVertexPtr GFSFVertex_;

	/**
	* Pointer to the Higgs-sfermion-sfermion vertex
	*/
	AbstractSSSVertexPtr HSFSFVertex_;

	/**
	* Pointer to the \f$\tilde{\chi}^+\f$-fermion-sfermion vertex
	*/
	AbstractFFSVertexPtr CFSFVertex_;

	/**
	* Pointer to the gluon-sfermion-sfermion vertex
	*/
	AbstractVSSVertexPtr GSFSFVertex_;

	/**
	* Pointer to the gluon-gluon-squark-squark vertex;
	*/
	AbstractVVSSVertexPtr GGSQSQVertex_;

	/**
	* Pointer to the gluon-gluino-gluino vertex
	*/
	AbstractFFVVertexPtr GSGSGVertex_;

	/**
	* Pointer to the gluino-neutralino-gluon vertex
	*/
	AbstractFFVVertexPtr GNGVertex_;

	/**
	* Pointer to the neutralino-neutralino-Z vertex
	*/
	AbstractFFVVertexPtr NNZVertex_;

	/**
	* Pointer to the neutralino-neutralino-photon vertex
	*/
	AbstractFFVVertexPtr NNPVertex_;

	/**
	* Pointer to the vertex chargino-chargino-Z vertex
	*/
	AbstractFFVVertexPtr CCZVertex_;

	/**
	* Pointer to the vertex chargino-neutralino-Z vertex
	*/
	AbstractFFVVertexPtr CNWVertex_;

	/**
	* Pointer to the vertex gaugino-gaugino-higgs vertex
	*/
	AbstractFFSVertexPtr GOGOHVertex_;

	/**
	* Pointer to the vertex for a gauge boson and higgs
	*/
	AbstractVSSVertexPtr WHHVertex_;

	/**
	* Pointer to the vertex for flavour changing stop decay
	*/
	AbstractFFSVertexPtr NCTVertex_;

	/**
	* Pointer to the vertex for gravitino-neutralino vector boson
	*/
	AbstractRFVVertexPtr GVNVVertex_;

	/**
	* Pointer to the vertex for gravitino-neutralino Higgs boson
	*/
	AbstractRFSVertexPtr GVNHVertex_;

	/**
	* Pointer to the vertex for gravitino-fermion sfermion
	*/
	AbstractRFSVertexPtr GVFSVertex_;
	//@}
	};

	}

	#include "ThePEG/Utilities/ClassTraits.h"

	namespace ThePEG {

	/** @cond TRAITSPECIALIZATIONS */

	/** This template specialization informs ThePEG about the
	* base classes of SusyBase. */
	template <>
	struct BaseClassTrait<Herwig::SusyBase,1> {
	/** Typedef of the first base class of SusyBase. */
	typedef Herwig::StandardModel NthBase;
	};

	/** This template specialization informs ThePEG about the name of
	* the SusyBase class and the shared object where it is defined. */
	template <>
	struct ClassTraits<Herwig::SusyBase>
	: public ClassTraitsBase<Herwig::SusyBase> {
	/** Return a platform-independent class name */
	static string className() { return "Herwig::SusyBase"; }
	/**
	* The name of a file containing the dynamic library where the class
	* SusyBase is implemented. It may also include several, space-separated,
	* libraries if the class SusyBase depends on other classes (base classes
	* excepted). In this case the listed libraries will be dynamically
	* linked in the order they are specified.
	*/
	static string library() { return "HwSusy.so"; }
	};

	/** @endcond */

	}


	#endif /* HERWIG_SusyBase_H */

File Metadata

Mime Type: text/x-diff
Expires: Thu, Apr 24, 6:36 AM (1 d, 19 h)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 4887775
Default Alt Text: (46 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions