IR singularities with squared amplitudes produced by MadLoop
Hello,
I’m trying to use MadLoop to generate a standalone library for some processes with EW loop corrections, and have a few questions:
1. I’m confused by the IR poles in the results. For instance, for the simple QED process e- e+ -> mu- mu+, I did
define exc = g0 g+ g- h gha ghz ghwp ghwm z u c d e- b s mu- ta+ ta- w+
generate e- e+ > mu- mu+ QCD=0 [virt=QED] /exc
which produced a loop graph with a top quark loop on the photon propagator. When I did
launch -f
I got the result
|| Total(*) Born contribution (GeV^0):
| Born = 1.1012911820047
|| Total(*) virtual contribution normalized with born*alpha_
| Finite = -2.184221175345
| Single pole = -1.168534349724
| Double pole = 0.0000000000000
As far as I know, this loop graph should not have an IR singularity. Do you know the reason for the non-zero single pole?
2. Can I use the subtraction method in MadFKS to compute photon radiative corrections? It seems that only QCD NLO radiative corrections are available.
3. When I run MadLoop in standalone mode, how can I set a dynamical renormalization scale? The process I’m considering is again EW loop graphs. I don’t find the setscal.f file in my output.
Thanks,
Xiang Lv
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- Hua-Sheng Shao Edit question
- Solved by:
- Hua-Sheng Shao
- Solved:
- Last query:
- Last reply: