Changes between Version 12 and Version 13 of LoopInducedTimesTree

Timestamp:

Jun 16, 2017, 6:05:14 PM (7 years ago)

Author:

Valentin Hirschi

Comment:

--

LoopInducedTimesTree

@@ -208 +208 @@
 
 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.
+
+Notice that if the "background" tree-level structure is very different than the loop-induced ones you are interfering against, then it might be a good idea to add the corresponding tree-level topologies to the multi-channel.
+This can be achieved with the following small modifications. Around line 104 of the file '<MG_root_path>/madgraph/various/diagram_symmetry.py' , to the lines
+{{{
+        base_diagrams = base_objects.DiagramList(
+                   [(d.get_contracted_loop_diagram(base_model,FDStructRepo) if
+                   isinstance(d,loop_base_objects.LoopDiagram) else d) for d in
+               matrix_element.get('base_amplitude').get('loop_diagrams') \
+                                                            if d.get('type')>0])
+        diagrams = matrix_element.get_loop_diagrams()
+}}}
+Add the following ones right below:
+{{{
+        base_diagrams.extend([(d.get_contracted_loop_diagram(base_model,FDStructRepo) if
+                   isinstance(d,loop_base_objects.LoopUVCTDiagram) else d) for d in
+               matrix_element.get('base_amplitude').get('loop_UVCT_diagrams')])
+        diagrams.extend(matrix_element.get_loop_UVCT_diagrams())
+}}}
+And then around the line 317 of '<MG_root_path>/madgraph/iolibs/group_subprocs.py' , to the line
+{{{
+       diagrams = [(d.get_contracted_loop_diagram(model,FDStructRepo) if
+                   isinstance(d,loop_base_objects.LoopDiagram) else d) for d in
+                 me.get('base_amplitude').get('loop_diagrams') if d.get('type')>0]
+}}}
+Add the following one, right below:
+{{{
+      diagrams.extend([(d.get_contracted_loop_diagram(model,FDStructRepo) if
+                   isinstance(d,loop_base_objects.LoopUVCTDiagram) else d) for d in
+                  me.get('base_amplitude').get('loop_UVCT_diagrams')])
+}}}
+Finally, when using these modifications and in order to activate all these extra integration channels, you must set
+{{{
+    MG5_aMC>set group_subprocesses False
+}}}
+before the generation of your process (irrespectively of your MG5aMC distribution this time).