~maddevelopers/mg5amcnlo/2.9.4

« back to all changes in this revision

Viewing changes to Template/loop_material/StandAlone/Cards/MadLoopParams.dat

pass to v2.0.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
! This file is for the user to set the different parameters of MadLoop.
 
2
! The name of the variable to define must start with the '#' sign and then
 
3
! the value should be put immediately on the next line.
 
4
 
 
5
!     =================================================================================
 
6
!     The parameters below set the stability checks of MadLoop at run time
 
7
!     =================================================================================
 
8
 
 
9
!     Decide in which mode to set CutTools at runtime.
 
10
!     The possible modes are:
 
11
!
 
12
!     imode:|  actions performed by ctsxcut:                                     
 
13
!           |                                                                    
 
14
!       0   | (dp_dir,dp_inv)-> dp_Atest -> stable -> (only if stable=.false.) ->
 
15
!               | (mp_dir,mp_inv)-> mp_Atest -> stable                                
 
16
!       1   | (dp_dir)       -> dp_Ntest -> stable                               
 
17
!       2   | (dp_inv)       -> dp_Ntest -> stable                               
 
18
!       3   | (dp_dir,dp_inv)-> dp_Atest -> stable                               
 
19
!       4   | (mp_dir)       -> mp_Ntest -> stable                                 
 
20
!       5   | (mp_inv)       -> mp_Ntest -> stable                                
 
21
!       6   | (mp_dir,mp_inv)-> mp_Atest -> stable
 
22
!      -1   | Special mode, see below for the negative ones.
 
23
!     
 
24
!     Due to the architecture of the program, you are better off
 
25
!     rerunning the full PS point in quadruple precision than just a single loop 
 
26
!     because the two things would almost take the same time. So '-1' is always
 
27
!     very recommended.
 
28
#CTModeRun
 
29
-1
 
30
! Default :: -1
 
31
!     In the negative mode -1, MadLoop first evaluates each PS points in modes 1 and 2,
 
32
!     yielding results Res1 and Res2, and then check if:
 
33
!                    (Res1-Res2)/(2*(Res1+Res2)< MLStabThres
 
34
!     If it is not the case, MadLoop evaluates again the PS point in modes 4 and 5,
 
35
!     yielding results Res4 and Res5, and then check if:
 
36
!                    (Res4-Res5)/(2*(Res4+Res5)< MLStabThres
 
37
!     If it is the case then the unstable phase-space point could be cured. If it is
 
38
!     not the case, MadLoop outputs a warning. 
 
39
!     Notice that MLStabThres is used only when CTModeRun is negative.
 
40
#MLStabThres
 
41
1.0d-3
 
42
! Default :: 1.0d-3
 
43
!     You can add other evaluation method to check for the stability in DP and QP.
 
44
!     Below you can chose if you want to use zero, one or two rotations of the PS point 
 
45
!     in QP.
 
46
#NRotations_DP
 
47
1
 
48
! Default :: 1
 
49
#NRotations_QP
 
50
0
 
51
! Default :: 0
 
52
 
 
53
!     By default, MadLoop is allowed to slightly deform the Phase-Space point in input
 
54
!     so to insure perfect onshellness of the external particles and perfect energy-momentum
 
55
!     conservation. The deformation is minimal and such that it leaves the input PS point 
 
56
!     unchanged if it already satisfies the physical condiditions mentioned above.
 
57
!     This integer values select what is the method to be employed preferably to restore this
 
58
!     precision. It can take the following values:
 
59
!
 
60
!     -1 :: No method is used for double precision computations, and method 2 will be used
 
61
!           preferentially when quadruple precision (for which this precision improvement
 
62
!           is mandatory, otherwise quadruple precision is pointless)
 
63
!      1 :: This methods imitates what is done in PSMC, namely 
 
64
!           a) Set the space-like momentum of the last external particle to be the
 
65
!              opposite of the sum of the others (with a minus sign for the initial states).
 
66
!           b) Rescale all final state space-like momenta by a fixed value x computed such
 
67
!              that energy is conserved when particles are put exactly onshell. This value 
 
68
!              is determined numericaly via Ralph-Newton's method.
 
69
!           c) Set all energies to have particles exactly onshell.           
 
70
!      2  :: This method applies a shift to the energy and the x and y components of the first 
 
71
!            initial state momentum in order to restore exact energy momentum conservation after
 
72
!            particles have been put exactly onshell via a shift of the z component of their 
 
73
!            momenta.
 
74
#ImprovePSPoint
 
75
2
 
76
! Default :: 2
 
77
!     =================================================================================
 
78
!     The parameters below set two CutTools internal parameters accessible to the user. 
 
79
!     =================================================================================
 
80
 
 
81
!     Choose here what library to chose for CutTools to compute the scalar loops of the
 
82
!     master integral basis. The choices are as follows:
 
83
!     2 | AVH
 
84
!     3 | QCDLoop
 
85
#CTLoopLibrary
 
86
2
 
87
! Default :: 2
 
88
 
 
89
!     Choose here the stability threshold used within CutTools to decide when to go to
 
90
!     higher precision.
 
91
#CTStabThres
 
92
1.0d-2
 
93
! Default :: 1.0d-2
 
94
 
 
95
!     =================================================================================
 
96
!     The parameters below set the general behavior of MadLoop for the initialization
 
97
!     =================================================================================
 
98
 
 
99
!     Decide in which mode to set CutTools when performing MadLoop's initialization of
 
100
!     the helicity (and possibly loop) filter. The possible modes are:
 
101
!
 
102
!     imode:|  actions performed by ctsxcut:                                     
 
103
!           |                                                                    
 
104
!       0   | (dp_dir,dp_inv)-> dp_Atest -> stable -> (only if stable=.false.) ->
 
105
!               | (mp_dir,mp_inv)-> mp_Atest -> stable                                
 
106
!       1   | (dp_dir)       -> dp_Ntest -> stable                               
 
107
!       2   | (dp_inv)       -> dp_Ntest -> stable                               
 
108
!       3   | (dp_dir,dp_inv)-> dp_Atest -> stable                               
 
109
!       4   | (mp_dir)       -> mp_Ntest -> stable                                 
 
110
!       5   | (mp_inv)       -> mp_Ntest -> stable                                
 
111
!       6   | (mp_dir,mp_inv)-> mp_Atest -> stable                               
 
112
!
 
113
#CTModeInit
 
114
1
 
115
! Default :: 0
 
116
 
 
117
!     CheckCycle sets on how many PS points trials the initialization filters must be 
 
118
!     obtained. As long as MadLoop does not find that many consecutive PS points for 
 
119
!     which the filters are the same, it will start over but only a maximum of
 
120
!     MaxAttempts times.
 
121
#CheckCycle
 
122
3
 
123
! Default :: 3
 
124
#MaxAttempts
 
125
10
 
126
! Default :: 10
 
127
 
 
128
!     Setting the threshold for deciding wether a numerical contribution is analytically 
 
129
!     zero or not.
 
130
#ZeroThres
 
131
1.0d-9
 
132
! Default :: 1.0d-9
 
133
 
 
134
!     The setting below is recommended to be on as it allows to systematically used the 
 
135
!     first PS point thrown at ML5 to be used for making sure that the helicity filter
 
136
!     read from HelFilter.dat is consistent as it might be no longer up to date with
 
137
!     certain changes of the paramaters by the user.
 
138
#DoubleCheckHelicityFilter
 
139
.TRUE.
 
140
! Default :: .TRUE.
 
141
 
 
142
!     Some loop contributions may be zero for some helicities which are however
 
143
!     contributing. In order to save their computing time, you can chose here to try
 
144
!     to filter them out. The gain is typically minimal, so it is turned off by default.
 
145
#UseLoopFilter
 
146
.FALSE.
 
147
! Default :: .FALSE.
 
148
 
 
149
!     This decides wether consecutive consistency for the loop filtering setup is also 
 
150
!     required.
 
151
#LoopInitStartOver
 
152
.FALSE.
 
153
! Default :: .FALSE.
 
154
 
 
155
!     This decides wether consecutive consistency for the helicity filtering setup is also 
 
156
!     required. Better to set it to false as it can cause problems for unstable processes.
 
157
#HelInitStartOver
 
158
.FALSE.
 
159
! Default :: .FALSE.
 
160
 
 
161
/* End of param file */