← Back to branch summary

~maddevelopers/mg5amcnlo/2.7.1.3

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/long_ML_SMQCD_default/dux_mumvmxg/born_matrix.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:
tests/input_files/IOTestsComparison/long_ML_SMQCD_default/dux_mumvmxg/born_matrix.f
42
42
C     and helicities
43
43
C     for the point in phase space P(0:3,NEXTERNAL)
44
44
C     
45
 
C     Process: d u~ > m- vm~ g QED=2 QCD=1 [ virt = QCD ]
 
45
C     Process: d u~ > m- vm~ g QED<=2 QCD<=1 [ virt = QCD ]
46
46
C     
47
47
      IMPLICIT NONE
48
48
C     
119
119
      DO IHEL=1,NEXTERNAL
120
120
        JC(IHEL) = +1
121
121
      ENDDO
 
122
C     When spin-2 particles are involved, the Helicity filtering is
 
123
C      dangerous for the 2->1 topology.
 
124
C     This is because depending on the MC setup the initial PS points
 
125
C      have back-to-back initial states
 
126
C     for which some of the spin-2 helicity configurations are zero.
 
127
C      But they are no longer zero
 
128
C     if the point is boosted on the z-axis. Remember that HELAS
 
129
C      helicity amplitudes are no longer
 
130
C     lorentz invariant with expternal spin-2 particles (only the
 
131
C      helicity sum is).
 
132
C     For this reason, we simply remove the filterin when there is
 
133
C      only three external particles.
 
134
      IF (NEXTERNAL.LE.3) THEN
 
135
        DO IHEL=1,NCOMB
 
136
          GOODHEL(IHEL)=.TRUE.
 
137
        ENDDO
 
138
      ENDIF
122
139
      ANS = 0D0
123
140
      DO IHEL=1,NCOMB
124
141
        IF (USERHEL.EQ.-1.OR.USERHEL.EQ.IHEL) THEN
147
164
C     Returns amplitude squared summed/avg over colors
148
165
C     for the point with external lines W(0:6,NEXTERNAL)
149
166
C     
150
 
C     Process: d u~ > m- vm~ g QED=2 QCD=1 [ virt = QCD ]
 
167
C     Process: d u~ > m- vm~ g QED<=2 QCD<=1 [ virt = QCD ]
151
168
C     
152
169
      IMPLICIT NONE
153
170
C     
259
276
      RETURN
260
277
      END
261
278
 
 
279