def test_mg5_sm_22(self):
"""Test a semi-complete list of sm 2->2 processes"""
# Create a list of processes to check automatically
#my_proc_list = me_comparator.create_proc_list(\
# ['w+', 'w-', 'a', 'z', 'h', 'g', 'u', 'u~', 'd', 'd~',
# 'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~'],
# initial=2, final=2)
my_proc_list = me_comparator.create_proc_list(\
['w+', 'w-', 'a', 'z', 'h', 'g', 'u', 'u~', 'd', 'd~',
'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~'],
initial=2, final=2)
#my_proc_list = ['t t > t t']
# Store list of non-zero processes and results in file
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="sm",
energy=1000,
filename="sm_22.log")
self.assertTrue(len(my_proc_list) > 5000)
def test_gauge_2(self):
"""Test a semi-complete list of sm 2->2 processes"""
# Create a list of processes to check automatically
my_proc_list = me_comparator.create_proc_list(\
['a', 'g', 'u', 'u~', 'd', 'd~',
'b', 'b~', 'ta+', 'ta-', 'vt', 'vt~','ve','ve~'],
initial=2, final=2)
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list,
orders={
'QED': 99,
'QCD': 99
},
model="sm",
energy=90,
filename="sm_gauge_2_e90.log",
tolerance=1e-3)
# Do the test for high energy
self.compare_processes(my_proc_list,
orders={
'QED': 99,
'QCD': 99
},
model="sm",
energy=500,
filename="sm_gauge_2_e500.log",
tolerance=1e-3)
def test_mg5_minitest_sm(self):
"""Test a minimal list of sm 2->2 processes, mainly to test the test"""
# Create a list of processes to check automatically
my_proc_list = me_comparator.create_proc_list(\
['u'],
initial=2, final=2)
my_proc_list = ['e+ e- > a > e+ e-', 'h h > h h', 'e+ e+ > e- e-',
'u w+ > u w+']
# Store list of non-zero processes and results in file
#pickle_file = "mg4_sm_%sminitest.pkl" % self.suffix_name
self.compare_processes(my_proc_list, model='sm',
orders = {'QED':2, 'QCD':2},
filename = "sm_mini.log",
energy = 1000)
def test_mg5_minitest_sm(self):
"""Test a minimal list of sm 2->2 processes, mainly to test the test"""
# Create a list of processes to check automatically
my_proc_list = me_comparator.create_proc_list(\
['u'],
initial=2, final=2)
my_proc_list = ['e+ e- > a > e+ e-', 'h h > h h', 'e+ e+ > e- e-',
'u w+ > u w+']
# Store list of non-zero processes and results in file
#pickle_file = "mg4_sm_%sminitest.pkl" % self.suffix_name
self.compare_processes(my_proc_list, model='sm',
orders = {'QED':2, 'QCD':2},
filename = "sm_mini.log",
energy = 1000)
def test_gauge_2(self):
"""Test a semi-complete list of sm 2->2 processes"""
# Create a list of processes to check automatically
my_proc_list = me_comparator.create_proc_list(\
['a', 'g', 'u', 'u~', 'd', 'd~',
'b', 'b~', 'ta+', 'ta-', 'vt', 'vt~','ve','ve~'],
initial=2, final=2)
# Store list of non-zero processes and results in file
self.compare_processes(my_proc_list,
orders = {'QED':99, 'QCD':99},
model = "sm",
energy = 90,
filename = "sm_gauge_2_e90.log",
tolerance = 1e-3)
# Do the test for high energy
self.compare_processes(my_proc_list,
orders = {'QED':99, 'QCD':99},
model = "sm",
energy = 500,
filename = "sm_gauge_2_e500.log",
tolerance = 1e-3)
def test_mg5_sm_22(self):
"""Test a semi-complete list of sm 2->2 processes"""
# Create a list of processes to check automatically
#my_proc_list = me_comparator.create_proc_list(\
# ['w+', 'w-', 'a', 'z', 'h', 'g', 'u', 'u~', 'd', 'd~',
# 'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~'],
# initial=2, final=2)
my_proc_list = me_comparator.create_proc_list(\
['w+', 'w-', 'a', 'z', 'h', 'g', 'u', 'u~', 'd', 'd~',
'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~'],
initial=2, final=2)
#my_proc_list = ['e+ e- > e+ e-','e+ e- > e+ e- a']
# Store list of non-zero processes and results in file
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 = "sm",
energy = 1000,
filename = "sm_22.log")
from madgraph import MG4DIR
mg4_path = os.getcwd()
if '__main__' == __name__:
# Get full logging info
logging.config.fileConfig(os.path.join(mg5_path, 'tests', '.mg5_logging.conf'))
logging.root.setLevel(logging.INFO)
logging.getLogger('madgraph').setLevel(logging.INFO)
logging.getLogger('cmdprint').setLevel(logging.INFO)
logging.getLogger('tutorial').setLevel(logging.ERROR)
logging.basicConfig(level=logging.INFO)
#my_proc_list=['g g > h g', 'g g > h g g', 'g g > h g g g', 'g g > h g g g g']
my_proc_list = me_comparator.create_proc_list(['u', 'u~','t','t~','g','z','a', 'h'],
initial=2, final=2)
#my_proc_list = me_comparator.create_proc_list_enhanced(
# fermion, fermion, boson,
# initial=2, final_1=2, final_2 = 1)
# Set the model we are working with
model = 'sm'
# Create a MERunner object for MG4
my_mg4 = me_comparator.MG4Runner()
my_mg4.setup(mg4_path)
# Create a MERunner object for MG5
my_mg5 = me_comparator.MG5Runner()
my_mg5.setup(mg5_path, mg4_path)
import me_comparator
from madgraph import MG4DIR
mg4_path = os.getcwd()
if '__main__' == __name__:
# Get full logging info
logging.config.fileConfig(
os.path.join(mg5_path, 'tests', '.mg5_logging.conf'))
logging.root.setLevel(logging.INFO)
logging.getLogger('madgraph').setLevel(logging.INFO)
logging.getLogger('cmdprint').setLevel(logging.INFO)
logging.getLogger('tutorial').setLevel(logging.ERROR)
logging.basicConfig(level=logging.INFO)
#my_proc_list=['g g > h g', 'g g > h g g', 'g g > h g g g', 'g g > h g g g g']
my_proc_list = me_comparator.create_proc_list(
['u', 'u~', 't', 't~', 'g', 'z', 'a', 'h'], initial=2, final=2)
#my_proc_list = me_comparator.create_proc_list_enhanced(
# fermion, fermion, boson,
# initial=2, final_1=2, final_2 = 1)
# Set the model we are working with
model = 'sm'
# Create a MERunner object for MG4
my_mg4 = me_comparator.MG4Runner()
my_mg4.setup(mg4_path)
# Create a MERunner object for MG5
my_mg5 = me_comparator.MG5Runner()
my_mg5.setup(mg5_path, mg4_path)
from madgraph import MG4DIR
mg4_path = os.getcwd()
if "__main__" == __name__:
# Get full logging info
logging.config.fileConfig(os.path.join(mg5_path, "tests", ".mg5_logging.conf"))
logging.root.setLevel(logging.INFO)
logging.getLogger("madgraph").setLevel(logging.INFO)
logging.getLogger("cmdprint").setLevel(logging.INFO)
logging.getLogger("tutorial").setLevel(logging.ERROR)
logging.basicConfig(level=logging.INFO)
# my_proc_list=['g g > h g', 'g g > h g g', 'g g > h g g g', 'g g > h g g g g']
my_proc_list = me_comparator.create_proc_list(["u", "u~", "t", "t~", "g", "z", "a", "h"], initial=2, final=2)
# my_proc_list = me_comparator.create_proc_list_enhanced(
# fermion, fermion, boson,
# initial=2, final_1=2, final_2 = 1)
# Set the model we are working with
model = "sm"
# Create a MERunner object for MG4
my_mg4 = me_comparator.MG4Runner()
my_mg4.setup(mg4_path)
# Create a MERunner object for MG5
my_mg5 = me_comparator.MG5Runner()
my_mg5.setup(mg5_path, mg4_path)
