-
Notifications
You must be signed in to change notification settings - Fork 6
/
Copy pathLesHouchesreg.f
41 lines (41 loc) · 1.38 KB
/
LesHouchesreg.f
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
subroutine gen_leshouches_reg
implicit none
include 'pwhg_math.h'
include 'nlegborn.h'
include 'pwhg_flst.h'
include 'pwhg_kn.h'
include 'pwhg_rad.h'
include 'LesHouches.h'
integer ireg
nup=nlegreal
scalup=sqrt(rad_pt2max)
do ireg=1,nup
c Remember: gluons are marked 0 here!
idup(ireg)=flst_regular(ireg,rad_realreg)
if(ireg.le.2) then
istup(ireg)=-1
mothup(1,ireg)=0
mothup(2,ireg)=0
else
istup(ireg)=1
mothup(1,ireg)=1
mothup(2,ireg)=2
endif
spinup(ireg)=9
vtimup(ireg)=0
enddo
c add resonances, perform decays, put particles on shell, etc.(or nothing!)
call finalize_lh
c no remnants for now!
c Don't forget to set scale for scalup equal to the pt of the
c radiation (whatever it is now!)
c set color connections for all particles
write(*,*) 'gen_leshouches_reg: dummy interface to regular'//
1 ' contributions'
write(*,*) ' Replace with your own process-dependent one, to be'
write(*,*) ' put in the process-specific directory (e.g. /W,'//'
# /Z, /VBF_H...)'
write(*,*) ' The Makefile will automatically compile the version'
write(*,*) ' in the process-specific directory'
call exit(1)
end