Package madgraph :: Package various :: Module process_checks :: Class LoopMatrixElementEvaluator
[hide private]
[frames] | no frames]

Class LoopMatrixElementEvaluator

source code

            object --+    
                     |    
MatrixElementEvaluator --+
                         |
                        LoopMatrixElementEvaluator
Known Subclasses:

Class taking care of matrix element evaluation for loop processes.

Instance Methods [hide private]
 
__init__(self, cuttools_dir=None, output_path=None, tir_dir={}, cmd=FakeInterface(), *args, **kwargs)
Allow for initializing the MG5 root where the temporary fortran output for checks is placed.
source code
 
evaluate_matrix_element(self, matrix_element, p=None, options=None, gauge_check=False, auth_skipping=None, output='m2', PS_name=None, MLOptions={})
Calculate the matrix element and evaluate it for a phase space point Output can only be 'm2.
source code
 
fix_MadLoopParamCard(self, dir_name, mp=False, loop_filter=False, DoubleCheckHelicityFilter=False, MLOptions={})
Set parameters in MadLoopParams.dat suited for these checks.MP stands for multiple precision and can either be a bool or an integer to specify the mode.
source code
 
get_me_value(self, proc, proc_id, working_dir, PSpoint=[], PS_name=None, verbose=True)
Compile and run ./check, then parse the output and return the result for process with id = proc_id and PSpoint if specified.
source code
 
setup_ward_check(self, working_dir, file_names, mp=False)
Modify loop_matrix.f so to have one external massless gauge boson polarization vector turned into its momentum.
source code

Inherited from MatrixElementEvaluator: get_momenta

Inherited from object: __delattr__, __format__, __getattribute__, __hash__, __new__, __reduce__, __reduce_ex__, __repr__, __setattr__, __sizeof__, __str__, __subclasshook__

Class Methods [hide private]
 
fix_PSPoint_in_check(cls, dir_path, read_ps=True, npoints=1, hel_config=-1, mu_r=0.0, split_orders=-1)
Set check_sa.f to be reading PS.input assuming a working dir dir_name.
source code
 
parse_check_output(cls, output, format='tuple')
Parse the output string and return a pair where first four values are the finite, born, single and double pole of the ME and the fourth is the GeV exponent and the second value is a list of 4 momenta for all particles involved.
source code
Properties [hide private]

Inherited from object: __class__

Method Details [hide private]

__init__(self, cuttools_dir=None, output_path=None, tir_dir={}, cmd=FakeInterface(), *args, **kwargs)
(Constructor)

source code 

Allow for initializing the MG5 root where the temporary fortran output for checks is placed.

Overrides: object.__init__

evaluate_matrix_element(self, matrix_element, p=None, options=None, gauge_check=False, auth_skipping=None, output='m2', PS_name=None, MLOptions={})

source code 

Calculate the matrix element and evaluate it for a phase space point Output can only be 'm2. The 'jamp' and 'amp' returned values are just empty lists at this point. If PS_name is not none the written out PS.input will be saved in the file PS.input_<PS_name> as well.

Overrides: MatrixElementEvaluator.evaluate_matrix_element

fix_PSPoint_in_check(cls, dir_path, read_ps=True, npoints=1, hel_config=-1, mu_r=0.0, split_orders=-1)
Class Method

source code 

Set check_sa.f to be reading PS.input assuming a working dir dir_name. if hel_config is different than -1 then check_sa.f is configured so to evaluate only the specified helicity. If mu_r > 0.0, then the renormalization constant value will be hardcoded directly in check_sa.f, if is is 0 it will be set to Sqrt(s) and if it is < 0.0 the value in the param_card.dat is used. If the split_orders target (i.e. the target squared coupling orders for the computation) is != -1, it will be changed in check_sa.f via the subroutine CALL SET_COUPLINGORDERS_TARGET(split_orders).

get_me_value(self, proc, proc_id, working_dir, PSpoint=[], PS_name=None, verbose=True)

source code 

Compile and run ./check, then parse the output and return the result for process with id = proc_id and PSpoint if specified. If PS_name is not none the written out PS.input will be saved in the file PS.input_<PS_name> as well

parse_check_output(cls, output, format='tuple')
Class Method

source code 

Parse the output string and return a pair where first four values are the finite, born, single and double pole of the ME and the fourth is the GeV exponent and the second value is a list of 4 momenta for all particles involved. Return the answer in two possible formats, 'tuple' or 'dict'.

setup_ward_check(self, working_dir, file_names, mp=False)

source code 

Modify loop_matrix.f so to have one external massless gauge boson polarization vector turned into its momentum. It is not a pretty and flexible solution but it works for this particular case.