Esempio n. 1
0
 def test_mg5_mssm_23_p1(self):
     """Test a semi-complete list of mssm 2->3 processes"""
     # Create a list of processes to check automatically
     
     init_part_list1 = ['w+', 'a', 'z', 'h1', 'h+', 'g', 'u~',
         'd~', 'b~', 't~', 'ta+', 'vt~']
     init_part_list2 = ['w-', 'a', 'z', 'h1', 'h-', 'g', 'u',
         'd', 'b', 't', 'ta+', 'vt']        
     final_part_list1 = ['w-', 'w+', 'a', 'z', 'h1', 'h-', 'h+', 'g'] 
     final_part_list2 = ['dl~', 'dr~', 'ul~', 'ur~']
     final_part_list3 = ['dl', 'dr', 'ul', 'ur']
     
     my_proc_list = me_comparator.create_proc_list_2_3(
                 init_part_list1, init_part_list2, final_part_list1,
                   final_part_list2,final_part_list3)                                   
                                                               
     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+1 processes, with MSSM particles in pairs in
     # final state
     #my_proc_list = me_comparator.create_proc_list_enhanced(\
     #    sm_parts, mssm_parts, sm_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':3, 'QCD':3},
                          model = "mssm",
                          energy = 2000,
                          filename = "mssm_23_p1.log")    
 def test_mg5_mssm_23_p1(self):
     """Test a semi-complete list of mssm 2->3 processes"""
     # Create a list of processes to check automatically
     
     init_part_list1 = ['w+', 'a', 'z', 'h1', 'h+', 'g', 'u~',
         'd~', 'b~', 't~', 'ta+', 'vt~']
     init_part_list2 = ['w-', 'a', 'z', 'h1', 'h-', 'g', 'u',
         'd', 'b', 't', 'ta+', 'vt']        
     final_part_list1 = ['w-', 'w+', 'a', 'z', 'h1', 'h-', 'h+', 'g'] 
     final_part_list2 = ['dl~', 'dr~', 'ul~', 'ur~']
     final_part_list3 = ['dl', 'dr', 'ul', 'ur']
     
     my_proc_list = me_comparator.create_proc_list_2_3(
                 init_part_list1, init_part_list2, final_part_list1,
                   final_part_list2,final_part_list3)                                   
                                                               
     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+1 processes, with MSSM particles in pairs in
     # final state
     #my_proc_list = me_comparator.create_proc_list_enhanced(\
     #    sm_parts, mssm_parts, sm_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':3, 'QCD':3},
                          model = "mssm",
                          energy = 2000,
                          filename = "mssm_23_p1.log")    
Esempio n. 3
0
    def test_mg5_mssm_23_p2(self):
        """Test a semi-complete list of mssm 2->3 processes"""
        # Create a list of processes to check automatically

        init_part_list1 = [
            'w+', 'a', 'z', 'h1', 'h+', 'g', 'u~', 'd~', 'b~', 't~', 'ta+',
            'vt~'
        ]
        init_part_list2 = [
            'w-', 'a', 'z', 'h1', 'h-', 'g', 'u', 'd', 'b', 't', 'ta+', 'vt'
        ]
        final_part_list1 = ['w-', 'w+', 'a', 'z', 'h1', 'h-', 'h+', 'g']
        final_part_list2 = ['x1-', 'x1+', 'x2-', 'x2+', 'go', 'n1']
        final_part_list3 = ['x1-', 'x1+', 'x2-', 'x2+', 'go', 'n1']

        my_proc_list = me_comparator.create_proc_list_2_3(
            init_part_list1, init_part_list2, final_part_list1,
            final_part_list2, final_part_list3)

        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': 3,
                                       'QCD': 3
                                   },
                                   model="MSSM_SLHA2",
                                   energy=2000,
                                   filename="mssm_23_p2.log")
    def test_mg5_mssm_23_p2(self):
        """Test a semi-complete list of mssm 2->3 processes"""
        # Create a list of processes to check automatically
        
        init_part_list1 = ['w+', 'a', 'z', 'h1', 'h+', 'g', 'u~',
            'd~', 'b~', 't~', 'ta+', 'vt~']
        init_part_list2 = ['w-', 'a', 'z', 'h1', 'h-', 'g', 'u',
            'd', 'b', 't', 'ta+', 'vt']        
        final_part_list1 = ['w-', 'w+', 'a', 'z', 'h1', 'h-', 'h+', 'g'] 
        final_part_list2 = ['x1-', 'x1+', 'x2-', 'x2+', 'go', 'n1']
        final_part_list3 = ['x1-', 'x1+', 'x2-', 'x2+', 'go', 'n1']
        
        my_proc_list = me_comparator.create_proc_list_2_3(
                    init_part_list1, init_part_list2, final_part_list1,
                      final_part_list2,final_part_list3)

        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':3, 'QCD':3},
                             model = "mssm",
                             energy = 2000,
                             filename = "mssm_23_p2.log")