def test_mg5_mssm_13(self):
"""Test a semi-complete list of mssm 1->3 processes"""
# Create a list of processes to check automatically
my_proc_list = me_comparator.create_proc_list_enhanced(
init_part_list=['t1', 't1~', 't2', 't2~', 'b2', 'b2~'],
final_part_list_1=['b1', 'b1~', 'b2', 'b2~', 'n1', 'n2'],
final_part_list_2=[
'u', 'u~', 'd', 'd~', 'c', 'c~', 'b', 'b~', 't', 't~'
],
initial=1,
final_1=1,
final_2=2,
charge_conservation=False)
my_proc_list += [
't2 > b1 d~ u', 't2 > b2 d~ u', 't1 > b1 d~ u', 't1 > b2 d~ u',
't2~ > b1~ d u~', 't2~ > b2~ d u~', 't1~ > b1~ d u~',
't1~ > b2~ d u~'
]
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list,
orders={
'QED': 4,
'QCD': 4
},
model="MSSM_SLHA2",
energy=2000,
filename="mssm_13_%s.log")
def test_mg5_mssm_22(self):
"""Test a semi-complete list of mssm 2->2 processes"""
# Create a list of processes to check automatically
sm_parts = [
'w+', 'w-', 'a', 'z', 'h1', 'h+', 'h-', 'g', 'u', 'u~', 'd', 'd~',
'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~'
]
mssm_parts = [
'dl', 'dl~', 'dr', 'dr~', 'ul', 'ul~', 'ur', 'ur~', 'b1', 'b1~',
'b2', 'b2~', 't1', 't1~', 'ta1-', 'ta1+', 'ta2-', 'ta2+', 'svt',
'svt~', 'x1-', 'x1+', 'x2-', 'x2+', 'go', 'n1'
]
# Generate 2 -> 2 processes, with MSSM particles in pairs in
# final state
my_proc_list = me_comparator.create_proc_list_enhanced(\
sm_parts, mssm_parts)
for i in range(len(my_proc_list) // 500):
print 'step %s/%s' % (i + 1, len(my_proc_list) // 500)
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list[500 * i:500 * (i + 1)],
orders={
'QED': 2,
'QCD': 2
},
model="MSSM_SLHA2",
energy=2000,
filename="mssm_22.log")
def test_mg5_heft_23(self):
"""Test a heft 2->3 processes"""
# Create a list of processes to check automatically
sm_parts = ['g', 'a', 'w+', 'w-']
heft_parts = ['h']
# Generate 2 -> 1+2 processes, with one Higgs in final state
# final state
my_proc_list = me_comparator.create_proc_list_enhanced(\
sm_parts, sm_parts, heft_parts)
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list,
orders = {'QED':2, 'QCD':0, 'HIG':1, 'HIW':1},
model = "heft",
energy = 500,
filename = "heft_23.log")
def test_mg5_sm_13(self):
"""Test a semi-complete list of sm 1->3 processes"""
# Create a list of processes to check automatically
my_proc_list = me_comparator.create_proc_list_enhanced(
init_part_list = ['t','t~','ta+','ta-'],
final_part_list_1 = ['b','b~','vt','vt~'],
final_part_list_2 = ['u', 'u~', 'd', 'd~', 'c','c~', 'e+', 've', 'e-', 've~','mu-','vm'],
initial=1, final_1=1, final_2=2, charge_conservation=True)
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list,
orders = {'QED':4, 'QCD':4},
model = "sm",
energy = 1000,
filename = "sm_13.log")
def test_mg5_heft_23(self):
"""Test a heft 2->3 processes"""
# Create a list of processes to check automatically
sm_parts = ['g', 'a', 'w+', 'w-']
heft_parts = ['h']
# Generate 2 -> 1+2 processes, with one Higgs in final state
# final state
my_proc_list = me_comparator.create_proc_list_enhanced(\
sm_parts, sm_parts, heft_parts)
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list,
orders = {'QED':2, 'QCD':0, 'HIG':1, 'HIW':1},
model = "heft",
energy = 500,
filename = "heft_23.log")
def test_mg5_sm_13(self):
"""Test a semi-complete list of sm 1->3 processes"""
# Create a list of processes to check automatically
my_proc_list = me_comparator.create_proc_list_enhanced(
init_part_list = ['t','t~','ta+','ta-'],
final_part_list_1 = ['b','b~','vt','vt~'],
final_part_list_2 = ['u', 'u~', 'd', 'd~', 'c','c~', 'e+', 've', 'e-', 've~','mu-','vm'],
initial=1, final_1=1, final_2=2, charge_conservation=True)
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list,
orders = {'QED':4, 'QCD':4},
model = "sm",
energy = 1000,
filename = "sm_13.log")
def test_mg5_mssm_13(self):
"""Test a semi-complete list of mssm 1->3 processes"""
# Create a list of processes to check automatically
my_proc_list = me_comparator.create_proc_list_enhanced(
init_part_list = ['t1','t1~','t2', 't2~','b2','b2~'],
final_part_list_1 = ['b1','b1~','b2','b2~','n1','n2'],
final_part_list_2 = ['u', 'u~', 'd', 'd~', 'c','c~', 'b', 'b~', 't','t~'],
initial=1, final_1=1, final_2=2, charge_conservation=False)
my_proc_list += ['t2 > b1 d~ u', 't2 > b2 d~ u', 't1 > b1 d~ u', 't1 > b2 d~ u',
't2~ > b1~ d u~', 't2~ > b2~ d u~', 't1~ > b1~ d u~', 't1~ > b2~ d u~']
print len(my_proc_list)
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list,
orders = {'QED':4, 'QCD':4},
model = "mssm",
energy = 2000,
filename = "mssm_13_%s.log")
def test_mg5_mssm_22(self):
"""Test a semi-complete list of mssm 2->2 processes"""
# Create a list of processes to check automatically
sm_parts = ['w+', 'w-', 'a', 'z', 'h1', 'h+', 'h-', 'g', 'u', 'u~',
'd', 'd~', 'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~']
mssm_parts = ['dl', 'dl~', 'dr', 'dr~', 'ul', 'ul~', 'ur', 'ur~', 'b1',
'b1~', 'b2', 'b2~', 't1', 't1~', 'ta1-', 'ta1+', 'ta2-',
'ta2+', 'svt', 'svt~', 'x1-', 'x1+', 'x2-', 'x2+',
'go', 'n1']
# Generate 2 -> 2 processes, with MSSM particles in pairs in
# final state
my_proc_list = me_comparator.create_proc_list_enhanced(\
sm_parts, mssm_parts)
for i in range(len(my_proc_list)//500):
print 'step %s/%s' %(i+1,len(my_proc_list)//500 )
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list[500*i:500*(i+1)],
orders = {'QED':2, 'QCD':2},
model = "mssm",
energy = 2000,
filename = "mssm_22.log")
