No OneTemporary
Actions

Size

37 KB

Subscribers

None

View Options

	diff --git a/Cards/param_card_TemplateBasis.dat b/Cards/param_card_TemplateBasis.dat
	index 9ec7f64..1d44392 100644
	--- a/Cards/param_card_TemplateBasis.dat
	+++ b/Cards/param_card_TemplateBasis.dat
	@@ -1,125 +1,130 @@
	######################################################################
	## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
	######################################################################
	## ##
	## Width set on Auto will be computed following the information ##
	## present in the decay.py files of the model. ##
	## See arXiv:1402.1178 for more details. ##
	## ##
	######################################################################

	###################################
	## INFORMATION FOR LOOP
	###################################
	Block loop
	1 9.118800e+01 # MU_R

	###################################
	## INFORMATION FOR MASS
	###################################
	Block mass
	5 4.700000e+00 # MB
	6 1.730000e+02 # MT
	15 1.730000e+02 # MTAU
	23 9.118800e+01 # MZ
	25 1.250000e+02 # MH
	## Dependent parameters, given by model restrictions.
	## Those values should be edited following the
	## analytical expression. MG5 ignores those values
	## but they are important for interfacing the output of MG5
	## to external program such as Pythia.
	1 0.000000 # d : 0.0
	2 0.000000 # u : 0.0
	3 0.000000 # s : 0.0
	4 0.000000 # c : 0.0
	11 0.000000 # e- : 0.0
	12 0.000000 # ve : 0.0
	13 0.000000 # mu- : 0.0
	14 0.000000 # vm : 0.0
	# 15 0.000000 # ta- : 0.0
	16 0.000000 # vt : 0.0
	21 0.000000 # g : 0.0
	22 0.000000 # a : 0.0
	###################################
	## INFORMATION FOR NEWCOUP
	###################################
	Block basis
	0 Template
	Block newcoup
	0 0.1 # a
	1 -0.2 # b
	2 0.3 # c

	###################################
	## INFORMATION FOR SMINPUTS
	###################################
	Block sminputs
	1 1.325070e+02 # aEWM1

	###################################
	## INFORMATION FOR YUKAWA
	###################################
	Block yukawa
	5 4.700000e+00 # ymb
	6 1.730000e+02 # ymt
	15 1.777000e+00 # ymtau

	###################################
	## INFORMATION FOR DECAY
	###################################
	DECAY 6 1.491500e+00 # WT
	DECAY 23 2.441404e+00 # WZ
	DECAY 24 2.047600e+00 # WW
	DECAY 25 6.382339e-03 # WH
	+# BR NDA ID1 ID2
	+ 0.e+01 2 13 13 # BR(H -> mu mu)
	+ 0.e+01 2 21 21 # BR(H -> a a)
	+ 0.e+01 2 5 5 # BR(H -> b b)
	+ 0.e+01 2 4 4 # BR(H -> c c)
	## Dependent parameters, given by model restrictions.
	## Those values should be edited following the
	## analytical expression. MG5 ignores those values
	## but they are important for interfacing the output of MG5
	## to external program such as Pythia.
	DECAY 1 0.000000 # d : 0.0
	DECAY 2 0.000000 # u : 0.0
	DECAY 3 0.000000 # s : 0.0
	DECAY 4 0.000000 # c : 0.0
	DECAY 5 0.000000 # b : 0.0
	DECAY 11 0.000000 # e- : 0.0
	DECAY 12 0.000000 # ve : 0.0
	DECAY 13 0.000000 # mu- : 0.0
	DECAY 14 0.000000 # vm : 0.0
	DECAY 15 0.000000 # ta- : 0.0
	DECAY 16 0.000000 # vt : 0.0
	DECAY 21 0.000000 # g : 0.0
	DECAY 22 0.000000 # a : 0.0
	DECAY 82 0.000000 # ghg : 0.0
	DECAY 9000001 0.000000 # gha : 0.0
	DECAY 9000002 2.495200 # ghz : WZ
	DECAY 9000003 2.085000 # ghwp : WW
	DECAY 9000004 2.085000 # ghwm : WW
	#===========================================================
	# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
	#===========================================================

	Block QNUMBERS 9000001 # gha
	1 0 # 3 times electric charge
	2 1 # number of spin states (2S+1)
	3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
	4 1 # Particle/Antiparticle distinction (0=own anti)
	Block QNUMBERS 9000002 # ghz
	1 0 # 3 times electric charge
	2 1 # number of spin states (2S+1)
	3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
	4 1 # Particle/Antiparticle distinction (0=own anti)
	Block QNUMBERS 9000003 # ghwp
	1 3 # 3 times electric charge
	2 1 # number of spin states (2S+1)
	3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
	4 1 # Particle/Antiparticle distinction (0=own anti)
	Block QNUMBERS 9000004 # ghwm
	1 -3 # 3 times electric charge
	2 1 # number of spin states (2S+1)
	3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
	4 1 # Particle/Antiparticle distinction (0=own anti)
	Block QNUMBERS 82 # ghg
	1 0 # 3 times electric charge
	2 1 # number of spin states (2S+1)
	3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
	4 1 # Particle/Antiparticle distinction (0=own anti)
	diff --git a/Rosetta/Basis.py b/Rosetta/Basis.py
	index 45420ec..73a3134 100644
	--- a/Rosetta/Basis.py
	+++ b/Rosetta/Basis.py
	@@ -1,527 +1,570 @@
	import StringIO
	import re, sys, os, math, datetime
	from collections import namedtuple,OrderedDict
	from query import query_yes_no as Y_or_N
	from __init__ import PID, default_inputs, default_masses, input_names, particle_names, input_to_PID
	from eHDECAY import eHDECAY
	####################################################################################################
	# Base Basis class
	class Basis(object):
	'''
	Base class from which to derive other Higgs EFT basis classes.

	Designed to be instantiated with an SLHA format parameter card using
	read_param_card(). The contents of the parameter card are stored in
	self.param_card, while relevant information from blocks "basis",
	"newcoup", "mass" and "sminput" are stored in self.name, self.par_dict,
	self.mass and self.input respectively.

	self.name - string storing the value of the 0th element of Block basis.
	self.par_dict - ``OrderedDict`` of (name,value) pairs with the name taken to be
	the first non-whitespace characters after a "#" character
	in a parameter definition within Block newcoup.
	self.mass - ``OrderedDict`` of (PID,value) pairs within Block mass
	self.input - ``OrderedDict`` as for self.par_dict but within Block sminput.

	self.independent and self.dependent should be defined in accordance
	with the contents of Block newcoup, as well as self.required_inputs and
	self.required_masses with Blocks sminput and mass respectively
	(the Z ahd Higgs masses are also stored in SM input). A number of checks
	related to these definitions are performed by check_param_data() on the
	data read in from the parameter card.

	Basis also stores the coefficients in self.coeffs as a ``namedtuple``
	which stores the coefficients as data members and from which an ``OrderedDict``
	can also be recovered with self.coeffs._asdict()

	type(self.par_dict)==type(self.coeffs._asdict()) # True

	self.par_dict['myvar'] = 10.
	self.coeffs._asdict()['myvar'] == 10. # True
	self.coeffs.myvar == 10. # True

	The user should define methods calculate_dependent() and translate() to
	firstly calculate any dependent parameters and then translate the coefficients
	into the mass basis and populate self.newpar, a template of which can be
	obtained by creating an instance of HEFTRosetta.MassBasis.MassBasis() without
	the optional param_card argument.

	from MassBasis import MassBasis
	new_instance = MassBasis() # Empty MassBasis instance with all coeffs set to 0
	new_dict = new_instance.coeffs._asdict()

	set_newcard() uses the contents of self.newpar to write a modified parameter card
	in the mass basis and write_param_card() saves it to file.

	- Users should call translate() and set_newcard() in the __init__ of their derived class.
	+ Users should call translate() and set_newcard() in the __init__ of their derived class
	+ should they choose to overload the constructor.

	Derived classes can be used by the command line script "translate"
	'''

	independent, dependent=[], []
	required_inputs, required_masses = set(),set()
	translate_to = {'mass'}

	def __init__(self, param_card=None, block_in='newcoup', block_out='newcoup', output_basis='mass', keep_old=True):
	# Check that target basis has a translation implemented
	if output_basis in self.translate_to:
	self.target_basis = output_basis
	else:
	raise ValueError('''
	{}.translate_to does not contain "{}".
	Rosetta doesn't know of an implemented translation in {}
	'''.format(self.__class__.__name__,output_basis, self.__class__))

	self.param_card = param_card
	self.block_in, self.block_out, self.keep_old = block_in, block_out, keep_old # other options
	self.card, self.newcard = StringIO.StringIO(), StringIO.StringIO()
	self.SLHA_sminputs = OrderedDict()
	self.input, self.mass, self.par_dict = OrderedDict(), OrderedDict(), OrderedDict()
	self.newinput, self.newmass, self.newpar = dict(), dict(), dict()
	self.newname = 'Basis'
	self.all_coeffs = self.independent + self.dependent

	if param_card is not None: # read param card (sets self.par_dict, self.input, self.mass, self.name, self.card)
	assert os.path.exists(self.param_card), '{} does not exist!'.format(self.param_card)

	self.read_param_card() # read
	self.check_param_data() # consistency check on inputs
	self.calculate_dependent() # calculate dependent parameters (User defined)
	self.check_calculated_data() # consistency check on dependent parameters
	coeffclass = namedtuple(self.name, self.all_coeffs) # declare coefficient namedtuple
	self.coeffs = coeffclass(**self.par_dict) # store coefficients
	self.translate() # translate to new basis (User defined)
	self.check_new() # consistency check between self.newinput and self.newmass
	self.set_new_masses() # set new masses in self.newinput

	try: # run eHDECAY if eHDECAY_input() is implemented
	self.BRs = eHDECAY(self)
	print self.BRs
	except NotImplementedError:
	pass

	self.set_newcard() # set new param_card

	else: # if param_card option not given, instatiate with class name and all coeffs set
	# to 0 (used for creating an empty basis instance for use in translate() method)
	self.name=self.__class__.__name__
	coeffclass = namedtuple(self.name, self.all_coeffs)
	self.coeffs = coeffclass(**{k:0. for k in self.all_coeffs})


	def read_param_card(self):
	'''
	Reads SLHA style param card and stores values of parameters.
	Looks for Blocks "basis", "mass", "newcoup" and "sminput" and
	updates self.name, self.mass (by PDG number), self.par_dict (name, value)
	and self.input (name, value)/self.SLHA_sminputs (SLHA ID, value) respectively.
	The whole param_card is also stored in self.card for later use.
	'''

	def read_pattern(plines, patt, pdict, ikey=2, ival=1, convert_key=str, convert_val=float):
	'''Takes a list of lines and looks for regexp pattern "patt" assuming the pattern captures 2 groups.
	"key" and "val" refer to the index of the resulting group (1 or 2).
	"convert_val" and "convert_key" specify a functions with which to treat the resulting string match e.g. convert to float.
	If pdict is None, function returns convert_val( match.group(val) ) (used to assign one variable).
	Otherwise, pdict is appended with a key value pair ( convert_key(match.group(key)), convert_val(match.group(val)) ).
	'''
	for pline in plines[:-1]:
	try:
	match = re.match(patt,pline)
	if pdict is not None:
	key, val = convert_key(match.group(ikey)), convert_val(match.group(ival))
	if key in pdict: # check if variable has already been assigned
	label = 'PDG mass' if type(key)==int else 'variable'
	print 'Warning: {} "{}" assigned more than once, kept value {}'.format( label,key,val )
	pdict[key]= val
	else:
	return convert_val(match.group(ival))
	except AttributeError:
	pass
	return None

	with open(self.param_card,'r') as card:
	relevant_blocks =['mass','basis',self.block_in,'sminputs']
	lines = iter(card)
	for line in lines:
	line, block = self.line_is_block(line)
	while block in relevant_blocks:
	if block=='mass': # Get masses
	print >> self.card, line.strip('\n')
	param_lines = self.read_until(lines,'Block','DECAY') # lines in block mass
	read_pattern(param_lines, r'\s(\d+)\s+(\S+)\s+.', self.mass, ikey=1, ival=2, convert_key=int) # by PDG number

	for pline in param_lines[:-1]: print >> self.card, pline.strip('\n')
	line, block = self.line_is_block(param_lines[-1]) # set last line read to current line
	relevant_blocks.remove('mass')

	elif block=='basis': # Get basis name
	print >> self.card, line.strip('\n')
	param_lines = self.read_until(lines,'Block','DECAY') # lines in block basis
	self.name = read_pattern(param_lines,r'\s0\s+(\S+).', None, convert_val=str).strip()

	for pline in param_lines[:-1]: print >> self.card, pline.strip('\n')
	line, block = self.line_is_block(param_lines[-1]) # set last line read to current line
	relevant_blocks.remove('basis')

	elif block==self.block_in: # Get basis coefficients
	print >> self.card, line.strip('\n')
	param_lines = self.read_until(lines,'Block', 'DECAY') # lines in block newcoup
	read_pattern(param_lines, r'\s*\d+\s+(\S+)\s+#+\s+(\S+)', self.par_dict)

	for pline in param_lines[:-1]: print >> self.card, pline.strip('\n')
	line, block = self.line_is_block(param_lines[-1]) # set last line read to current line
	relevant_blocks.remove(self.block_in)

	elif block=='sminputs': # Get SM inputs
	print >> self.card, line.strip('\n')
	param_lines = self.read_until(lines,'Block', 'DECAY') # lines in block sminputs
	read_pattern(param_lines, r'\s*(\d+)\s+(\S+)\s+#+\s+\S+', self.SLHA_sminputs, ikey=1, ival=2, convert_key = int)
	read_pattern(param_lines, r'\s*\d+\s+(\S+)\s+#+\s+(\S+)', self.input)

	for pline in param_lines[:-1]: print >> self.card, pline.strip('\n')
	line, block = self.line_is_block(param_lines[-1]) # set last line read to current line
	relevant_blocks.remove('sminputs')

	if not block or (block not in relevant_blocks): print >> self.card, line.strip('\n')

	if relevant_blocks: # checks if all relevant block were found
	if relevant_blocks==['basis']: # only optional basis block unread
	print 'Warning: block basis not found.'
	carry_on = Y_or_N('Continue assuming default value "{}"?'.format(self.__class__.__name__))
	if carry_on:
	self.name=self.__class__.__name__
	else:
	print 'Exit'
	sys.exit()
	else:
	raise IOError("Required block(s) ({}) not found in {}".format(', '.join(relevant_blocks), self.param_card))

	# Add MH,MW,MZ,mt,mb,mtau,aEWM1 to SLHA_sminputs
	for inID,PID in input_to_PID.iteritems():
	if inID in self.required_inputs:
	try:
	input_mass = self.mass[PID]
	self.input[input_names[inID]] = input_mass
	self.SLHA_sminputs[inID] =input_mass
	except KeyError:
	pass

	def check_param_data(self):
	'''
	Compares lists of coefficients declared in self.independent and self.dependent to those read in from self.param_card.
	1) Deals with unrecognised names by removing them from self.par_dict
	2) Prints a warning if coefficients declared as dependent are assigned values in param_card.
	If so, the user is asked whether or not they wish to continue.
	3) Checks if all coefficients declared as independent are assigned values.
	If not, the user is given the option to continue with them set to 0.
	4) Ensures all fermion masses declared in self.required_masses are defined.
	If not, the user is given the option to continue with them set to 0.
	5) Ensures all sm imputs declared in self.required_inputs are defined.
	If not, the user is given the option to continue with them set to default
	values defined in Rosetta.default_inputs.
	'''

	unknown_coeffs = set(self.par_dict.keys()).difference( self.all_coeffs )
	if unknown_coeffs: # Check for unrecognised coeff names
	print 'Warning: you have declared coefficients undefined in {}.'.format(self.__class__)
	print 'The following will be ignored: {}'.format(','.join(unknown_coeffs))
	for c in unknown_coeffs: del self.par_dict[c]

	defined_dependent = set(self.dependent).intersection( set(self.par_dict.keys()) )
	if defined_dependent: # Check if coefficients defined as dependent have been assigned values
	print 'Warning: you have assigned values to some coefficients defined as dependent in {}.'.format(self.__class__)
	print 'Coefficients: {}'.format(','.join(defined_dependent))
	print 'These may be overwritten by an implementation of {}.translate()'.format(self.__class__.__name__)
	carry_on = Y_or_N('Continue?')
	if carry_on:
	pass
	else:
	print 'Exit'
	sys.exit()

	missing_coeffs = set(self.independent).difference( set(self.par_dict.keys()) )
	if missing_coeffs: # Deal with unassigned independent coefficients
	print 'Warning: Set of independent coefficients read from {} does not match those specified in {}.'.format(self.param_card,self.__class__)
	print 'Undefined: {}'.format(', '.join(missing_coeffs))
	carry_on = Y_or_N('Continue assuming unspecified coefficients are Zero?')
	if carry_on:
	for m in missing_coeffs: self.par_dict[m]=0.
	else:
	print 'Exit'
	sys.exit()

	missing_masses = set(self.required_masses).difference(self.mass.keys())
	repr_default_masses = ['{}={:.5e}'.format(k,default_masses.get(k,0.)) for k in missing_masses]
	if missing_masses: # Deal with unassigned fermion masses
	print 'Warning: Not all required fermion masses are defined in {}.'.format(self.param_card)
	print 'Required PIDs: {}'.format(', '.join([str(x) for x in self.required_masses]))
	print 'Missing PIDs: {}'.format(', '.join([str(x) for x in missing_masses]))
	carry_on = Y_or_N('Continue assuming default values for unspecified masses? ({})'.format(', '.join(repr_default_masses)))
	if carry_on: # assigns a data member for unspecified inputs for writing to self.newcard
	self.missing_masses=dict()
	for m in missing_masses:
	self.mass[m]=default_masses[m]
	self.missing_masses[m]=default_masses[m]
	else:
	print 'Exit'
	sys.exit()

	missing_inputs = set(self.required_inputs).difference(self.SLHA_sminputs.keys())
	repr_default_inputs = ['{}={:.5e}'.format(input_names[k],default_inputs.get(k,0.)) for k in missing_inputs]
	if missing_inputs: # Deal with unassigned SM inputs
	print 'Warning: Not all required SM inputs are defined in {}.'.format(self.param_card)
	print 'Required inputs: {}'.format(', '.join(['{} ({})'.format(input_names[x],x) for x in self.required_inputs]))
	print 'Missing inputs: {}'.format(', '.join([str(x) for x in missing_inputs]))
	carry_on = Y_or_N('Continue with default values for unspecified inputs? ({})'.format(', '.join(repr_default_inputs)))
	if carry_on:# assigns a data member for unspecified inputs for writing to self.newcard
	self.missing_inputs=dict()
	for m in missing_inputs:
	self.SLHA_sminputs[m]=default_inputs[m]
	self.input[input_names[m]]=default_inputs[m]
	self.missing_inputs[m]=default_inputs[m]
	else:
	print 'Exit'
	sys.exit()
	print 'Param card data are OK.'

	def check_calculated_data(self):
	'''
	Compares self.dependent and self.par_dict keys to see if all dependent coefficients have been calculated.
	If not, asks whether the user wants to continue assuming they are zero.
	'''
	missing_dependents = set(self.dependent).difference(self.par_dict.keys())
	if missing_dependents and self.dependent:
	print 'Warning: Set of dependent coefficients calculated by {0}.calculate_dependent() does not match those specified in {0}.'.format(self.__class__)
	print 'Undefined: {}'.format(', '.join(missing_dependents))
	carry_on = Y_or_N('Continue assuming coefficients are Zero?')
	if carry_on:
	for m in missing_dependents: self.par_dict[m]=0.
	else:
	print 'Exit'
	sys.exit()
	print 'Calculated coefficients match those defined in {}.dependent.'.format(self.__class__.__name__)

	def check_new(self):
	'''
	Check consistency of modifications to mass parameters defined as possible SHLA inputs:
	MZ (PID = 23, SLHA key = 4)
	MW (PID = 24, SLHA key = 7)
	MH (PID = 25, SLHA key = 8)
	Compares self.newinput to self.newmass and makes sure user hasn't specified different new values
	'''
	for ID_input, ID_mass in input_to_PID.iteritems():
	try:
	ipar, mpar = self.newinput[ID_input],self.newmass[ID_mass]
	if not ipar==mpar:
	print 'Warning: modified mass {} in self.newinput ({}) does not match value in self.newmass ({})'.format(input_name[ID_input],ipar,mpar)
	carry_on = Y_or_N('Continue using values in self.newmass?')
	if carry_on:
	pass
	else:
	print 'Exit'
	sys.exit()
	except KeyError as e:
	pass

	def set_new_masses(self):
	'''
	If the user modifies a mass input parameter, appends the change to self.mass
	'''
	for inID, PID in input_to_PID.iteritems():
	if inID in self.newinput and PID not in self.newmass: self.newmass[PID] = self.newinput[inID]

	def set_newcard(self):
	'''
	Generates a new param card in self.newcard from self.card, adding the new set of coefficients in self.newpar after "Block newcoup".
	If self.keep_old is True, the original names and values of newcoup variables are included, commented out.
	'''
	print >> self.newcard, '######################################################################'
	print >> self.newcard, '############# COEFFICIENTS TRANSLATED BY ROSETTA MODULE #############'
	print >> self.newcard, '########### PARAM_CARD GENERATED {} ###########'.format(datetime.datetime.now().ctime().upper())
	print >> self.newcard, '######################################################################'

	- blocks_to_modify =('basis',self.block_in, 'mass', 'sminputs')
	+ blocks_to_modify =('basis', self.block_in, 'mass', 'sminputs')
	lines = iter(self.card.getvalue().splitlines()) # Lines of old param card
	for line in lines:
	+
	+

	line, block = self.line_is_block(line)
	if block==self.block_in: # Add line to self.newcard, modify block name
	print >> self.newcard, line.strip('\n').replace(self.block_in,self.block_out)
	+ elif self.line_is_decay(line)=='25': # special case for modified higgs width and BRs
	+ try:
	+ totalwidth = float(self.BRs['WTOT'])
	+ print >> self.newcard, 'DECAY 25 {} # New Higgs total width'.format(totalwidth)
	+ print >> self.newcard, '# BR NDA ID1 ID2'.format(totalwidth)
	+ particle_IDs = {v:k for k,v in particle_names.iteritems()}
	+ for channel,BR in self.BRs.iteritems():
	+ if channel!='WTOT':
	+ p1, p2 = channel[:len(channel)/2], channel[len(channel)/2:]
	+ id1, id2 = particle_IDs[p1], particle_IDs[p2]
	+ print >> self.newcard, ' {:.5e} 2 {: <2} {: <2} # BR(H -> {} {})'.format(BR, id1, id2, p1, p2)
	+ param_lines = self.read_until(lines, 'Block', 'DECAY')[:-1] # read old BR values
	+ for i,pline in enumerate(param_lines):
	+ if i==0:
	+ print >> self.newcard, '# Old BR info'
	+ print >> self.newcard, '#'+line.strip('\n')
	+ if self.keep_old:
	+ comment_out = (pline and not pline.startswith('#'))
	+ print >> self.newcard, '#' + pline.strip('\n') if comment_out else pline.strip('\n')
	+ else:
	+ if not pline.strip() or pline.startswith('#'): print >> self.newcard, pline.strip('\n')
	+ line, block = self.line_is_block(param_lines[-1]) # set last line read to current line
	+ except AttributeError:
	+ print >> self.newcard, line.strip('\n')
	+
	else: # Add line to self.newcard
	print >> self.newcard, line.strip('\n')

	+
	while block in blocks_to_modify: # Slightly different actions for each block

	if block=='basis': # When Block basis is reached
	print >> self.newcard, ' 0 {} # translated basis'.format(self.newname)
	param_lines = self.read_until(lines, 'Block', 'DECAY')
	for pline in param_lines[:-1]:
	if self.keep_old:
	comment_out = (pline and not pline.startswith('#'))
	print >> self.newcard, '#' + pline.strip('\n') if comment_out else pline.strip('\n')
	else:
	if not pline.strip() or pline.startswith('#'): print >> self.newcard, pline.strip('\n')
	line, block = self.line_is_block(param_lines[-1]) # set last line read to current line

	elif block==self.block_in: # When Block self.block_in is reached,
	for i,(par,val) in enumerate(self.newpar.items()): # write out new couplings
	print >> self.newcard, ' {} {:.5e} # {}'.format(i,val,par)
	param_lines = self.read_until(lines,'Block', 'DECAY')
	if self.keep_old: # write old parameters
	print >> self.newcard, ''
	print >> self.newcard, '###################################'
	print >> self.newcard, '## COEFFICIENTS IN {} BASIS'.format(self.name.upper())
	print >> self.newcard, '###################################'
	for pline in param_lines[:-1]:
	comment_out = (pline and not pline.startswith('#'))
	print >> self.newcard, '#' + pline.strip('\n') if comment_out else pline.strip('\n')
	else:
	for pline in param_lines[:-1]:
	if not pline.strip() or pline.startswith('#'): print >> self.newcard, pline.strip('\n')
	line, block = self.line_is_block(param_lines[-1]) # set last line read to current line

	elif block=='sminputs':# When block sminputs is reached

	try: # add missing inputs
	for i,(ID, inpt) in enumerate(self.missing_inputs.items()):
	if i==0: print >> self.newcard, '# missing inputs'
	if ID in self.newinput:
	print >> self.newcard, ' {} {:.5e} # {} '.format(ID, self.newinput[ID], input_names[ID])
	print >> self.newcard, '# {} {:.5e} # {} old value'.format(ID, inpt, input_names[ID])
	else:
	print >> self.newcard, ' {} {:.5e} # {} '.format(ID, inpt, input_names[ID])
	except AttributeError as e:
	pass

	# add additional new inputs
	input_names_reversed={v:k for k,v in input_names.iteritems()}
	for i,(ID, inpt) in enumerate(self.newinput.items()):
	if i==0: print >> self.newcard, '# additional inputs'
	if not ID in self.SLHA_sminputs:
	if type(ID)==int:
	print >> self.newcard, ' {} {:.5e} # {} '.format(ID, inpt, input_names.get(ID,''))
	elif type(ID)==str:
	print >> self.newcard, ' {} {:.5e} # {} '.format(input_names_reversed.get(ID,99), inpt,ID )

	print >> self.newcard, '# original inputs'

	param_lines = self.read_until(lines,'Block', 'DECAY')
	for pline in param_lines[:-1]:
	match = re.match(r'\s(\d+)\s+(\S+)\s+.',pline.strip()) # read old inputs
	try: # do if match
	ID, inpt = int(match.group(1)), match.group(2)
	try: # if a new value exists
	print >> self.newcard, pline.replace(inpt, '{:.5e}'.format(self.newinput[ID])).strip('\n')
	if self.keep_old: print >> self.newcard, '# '+pline.strip('\n')+' # old value'
	except KeyError as e:
	print >> self.newcard, pline.strip('\n')
	except AttributeError as e:
	print >> self.newcard, pline.strip('\n')

	line, block = self.line_is_block(param_lines[-1]) # set last line read to current line

	elif block=='mass': # When block mass is reached

	try: # add missing masses
	for i,(ID, mass) in enumerate(self.missing_masses.items()):
	if i==0: print >> self.newcard, '# missing masses'
	if ID in self.newmass:
	print >> self.newcard, ' {} {:.5e} # M{} '.format(ID, self.newmass[ID], particle_names[ID])
	print >> self.newcard, '# {} {:.5e} # M{} old value'.format(ID, mass, particle_names[ID])
	else:
	print >> self.newcard, ' {} {:.5e} # M{} '.format(ID, mass, particle_names[ID])
	except AttributeError as e: # self.missing_masses may not exist
	pass

	try: # add missing inputs that are masses
	for inID, PID in input_to_PID.iteritems():
	try:
	print >> self.newcard, ' {} {:.5e} # {} '.format(PID,self.missing_inputs[inID],input_names[inID])
	except KeyError:
	pass
	except AttributeError as e: # self.missing_inputs may not exist
	pass

	print >> self.newcard, '# original masses'

	param_lines = self.read_until(lines,'Block', 'DECAY')
	for pline in param_lines[:-1]:
	match = re.match(r'\s(\d+)\s+(\S+)\s+.',pline.strip()) # read old mass
	try: # do if match
	ID, mass = int(match.group(1)), match.group(2)
	try: # if a new value exists
	print >> self.newcard, pline.replace(mass, '{:.5e}'.format(self.newmass[ID])).strip('\n')
	if self.keep_old: print >> self.newcard, '# '+pline.strip('\n')+' # old value'
	except KeyError:
	print >> self.newcard, pline.strip('\n')
	except AttributeError:
	print >> self.newcard, pline.strip('\n')

	line, block = self.line_is_block(param_lines[-1]) # set last line read to current line

	if block==self.block_in:
	print >> self.newcard, line.strip('\n').replace(self.block_in,self.block_out) # Add line to self.newcard, modify block name
	else:
	- print >> self.newcard, line.strip('\n') # Add line to self.newcard
	+ print >> self.newcard, line.strip('\n') # Add line to self.newcard
	+
	return None

	def write_param_card(self,filename,overwrite=False):
	'''Write contents of self.newcard to filename'''
	contents = self.newcard.getvalue()
	if not contents:
	print '{}.newcard is empty, nothing done.'.format(self)
	return None
	if os.path.exists(filename) and not overwrite:
	print '{} already exists.'.format(filename)
	carry_on = Y_or_N('Overwrite?')
	else:
	carry_on=True
	if carry_on:
	with open(filename,'w') as param_card:
	param_card.write(contents)
	return True
	else:
	return False

	@staticmethod
	def read_until(lines, here, *args):
	'''Loops through an iterator of strings by calling next() until
	it reaches a line starting with a particular string.
	Case insensitive.
	Args:
	lines - iterator of strings
	here - string (plus any further argumnts). Reading will end if the line matches any of these.
	Return:
	lines_read - list of lines read
	line - last line that was read (containing string "here")
	'''
	end_strings = [here.lower()]+[a.lower() for a in args]
	lines_read = []
	line = ''
	while not any([line.strip().lower().startswith(x) for x in end_strings]):
	line = lines.next()
	lines_read.append(line.strip('\n'))
	return lines_read

	@staticmethod
	def line_is_block(line):
	'''
	Looks for "BLOCK XXXX" pattern in line.
	returns (line, XXXX) if line matches.
	returns (line, False) if line doesn't match
	'''
	try:
	return line, re.match(r'block\s+(\S+).*',line.strip().lower()).group(1)
	except AttributeError:
	return line, False
	+
	+ @staticmethod
	+ def line_is_decay(line):
	+ '''
	+ Looks for "DECAY XXXX YYYY" pattern in line.
	+ returns XXXX if line matches.
	+ returns None if line doesn't match
	+ '''
	+ match = re.match(r'decay\s+(\S+)\s+.*',line.strip().lower())
	+ try:
	+ return match.group(1)
	+ except AttributeError:
	+ return None

	def calculate_dependent(self):
	print 'Nothing done for {}.calculate_dependent()'.format(self.__class__.__name__)

	def translate(self): # default behaviour for translate()
	self.keep_old=False
	self.newpar = self.coeffs._asdict()
	def eHDECAY_inputs(self):
	raise NotImplementedError
	####################################################################################################
	\ No newline at end of file

File Metadata

Mime Type: text/x-diff
Expires: Tue, Nov 19, 6:26 PM (1 d, 21 h)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 3805585
Default Alt Text: (37 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions