Changes between Version 10 and Version 11 of LoopInducedTimesTree

Timestamp:

Mar 25, 2017, 9:37:22 PM (7 years ago)

Author:

Valentin Hirschi

Comment:

--

LoopInducedTimesTree

@@ -146 +146 @@
 }}}
 >d.7) MadGraph5_aMC@NLO will then ask you questions related to the details of either the MadLoop standalone run you are doing or the event generation run. 
+------------
+= Optional
+------------
+
+You can decide to use a normal model instead of the special one suggested here. Then, one can use an ad-hoc new command in MG5aMC to create a copy of all base interaction and orders so as to generate their 'BKG_' equivalent.
+This can be done by adding the following in the class 'MadGraphCmd' in '<MG_root_path>/madgraph/loop/loop_diagram_generation.py':
+
+{{{
+    def do_prepare_model_for_loopInducedXTrees(self, line):
+        """Commands for adding the necessary TREE UV interactions to the current model."""
+        
+        # Short-hand to the current model that will be modified
+        model = self._curr_model
+
+        # First add the necessary BKG coupling_orders to the available coupling_orders 
+        # in the model
+        if model['order_hierarchy']:
+            for coupling_order in model['order_hierarchy'].keys():
+                model['order_hierarchy']['BKG_%s'%coupling_order]=model['order_hierarchy'][coupling_order]
+                model['perturbation_couplings'].append('BKG_%s'%coupling_order)
+
+        new_order_hierarchy = dict(model['order_hierarchy'])
+        new_perturbation_couplings = list(model['perturbation_couplings'])
+
+        # Dubplicate base interactions to make UV ones
+        new_interactions=[]
+        id_offset = 1000000
+        for inter in model['interactions'].get_type('base'):
+            new_interactions.append(copy.deepcopy(inter))
+            # Deepcopy intorrectly treats colorstring object. Need to do it by hand.
+            new_interactions[-1].set('color',[cs.create_copy() for cs in inter.get('color')])
+            id_offset += 1
+            new_interactions[-1].set('id',id_offset)            
+            new_interactions[-1].set('type','UVtree')
+            new_interactions[-1].set('loop_particles',[[]])
+            new_orders = {}
+            for order in new_interactions[-1]['orders']:
+                new_orders['BKG_%s'%order]=new_interactions[-1]['orders'][order]
+            new_interactions[-1].set('orders',new_orders)
+
+
+        model.set('interactions',base_objects.InteractionList(model.get('interactions')+new_interactions))
+        # Reset the model dictionaries to sync changes
+        model.reset_dictionaries()
+        # Refresh dictionaries manually (would be done automatically anyway)
+        model.actualize_dictionaries()
+       
+        model.set('order_hierarchy', new_order_hierarchy)
+        model.set('perturbation_couplings', new_perturbation_couplings)
+
+        logger.info('Model %s successfully patched for (finite) Loops x Trees interference computation.'%model.get('name'))
+        return
+}}}
+
+And then add the following in the class 'Switcher' of '<MG_root_path>/madgraph/interface/master_interface.py':
+
+{{{
+    def do_prepare_model_for_loopInducedXTrees(self, *args, **opts):
+        return self.cmd.do_prepare_model_for_loopInducedXTrees(self, *args, **opts)
+}}}
+
+You can then simply run the command 'prepare_model_for_loopInducedXTrees' to preprocess your model before running loop-induced x tree computation as instructed above.