← Back to branch summary

~maddevelopers/mg5amcnlo/2.7.1.3

« back to all changes in this revision

Viewing changes to Template/NLO/SubProcesses/BinothLHA_user.f

  • Committer: olivier Mattelaer
  • Date: 2016-05-12 11:00:18 UTC
  • mfrom: (262.1.150 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160512110018-sevb79f0wm4g8mpp
pass to 2.4.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/NLO/SubProcesses/BinothLHA_user.f
5
5
c also the born_wgt is passed to this subroutine.
6
6
c
7
7
C************************************************************************
8
 
c WARNING WARNING WARNING WARNING WARNING WARNING WARNING WARNING WARNING
9
 
C************************************************************************
10
 
c The Born in MadFKS -- and therefore also the virtual!-- should have a
11
 
c slightly adapted identical particle symmetry factor. The normal
12
 
c virtual weight as coming from the OLP should be divided by the number
13
 
c of gluons in the corresponding real-emission process (i.e.  the number
14
 
c of gluons in the Born plus one). This factor is passed to this
15
 
c subroutine in /numberofparticles/ common block, as "ngluons". So,
16
 
c divided virt_wgt by dble(ngluons) to get the correct virtual to be
17
 
c used in MadFKS. The born_wgt that is passed to this subroutine has
18
 
c already been divided by this factor.
19
 
C************************************************************************
20
8
c
21
9
      implicit none
22
10
      include "nexternal.inc"
27
15
      double precision virt_wgt,born_wgt,double,single,virt_wgts(3)
28
16
      double precision mu,ao2pi,conversion,alpha_S
29
17
      save conversion
30
 
      double precision fkssymmetryfactor,fkssymmetryfactorBorn,
31
 
     &     fkssymmetryfactorDeg
32
 
      integer ngluons,nquarks(-6:6)
33
 
      common/numberofparticles/fkssymmetryfactor,fkssymmetryfactorBorn,
34
 
     &                         fkssymmetryfactorDeg,ngluons,nquarks
35
18
      logical firsttime,firsttime_conversion
36
19
      data firsttime,firsttime_conversion /.true.,.true./
37
20
      integer           isum_hel
62
45
c======================================================================
63
46
c Replace the following line with the call to the one-loop code you 
64
47
c wish to use. virt_wgts contains finite part, single and double pole
65
 
c To understand why we need to divide them by ngluons, read the warning
66
 
c message at the beginning of this subroutine
67
48
c      
68
49
c      call sloopmatrix(p, virt_wgts)
69
 
c      virt_wgt= virt_wgts(1)/dble(ngluons)
70
 
c      single  = virt_wgts(2)/dble(ngluons)
71
 
c      double  = virt_wgts(3)/dble(ngluons)
 
50
c      virt_wgt= virt_wgts(1)
 
51
c      single  = virt_wgts(2)
 
52
c      double  = virt_wgts(3)
72
53
c
73
54
c======================================================================
74
55