← Back to branch summary

~maddevelopers/mg5amcnlo/2.5.4_run_py8_at_evtgen

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/gg_wmtbx/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_optimized/gg_wmtbx/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: g g > w- t b~ QED=1 QCD=2 [ virt = QCD ]
 
45
C     Process: g g > w- t b~ QED<=1 QCD<=2 [ virt = QCD ]
46
46
C     
47
47
      IMPLICIT NONE
48
48
C     
135
135
      DO IHEL=1,NEXTERNAL
136
136
        JC(IHEL) = +1
137
137
      ENDDO
 
138
C     When spin-2 particles are involved, the Helicity filtering is
 
139
C      dangerous for the 2->1 topology.
 
140
C     This is because depending on the MC setup the initial PS points
 
141
C      have back-to-back initial states
 
142
C     for which some of the spin-2 helicity configurations are zero.
 
143
C      But they are no longer zero
 
144
C     if the point is boosted on the z-axis. Remember that HELAS
 
145
C      helicity amplitudes are no longer
 
146
C     lorentz invariant with expternal spin-2 particles (only the
 
147
C      helicity sum is).
 
148
C     For this reason, we simply remove the filterin when there is
 
149
C      only three external particles.
 
150
      IF (NEXTERNAL.LE.3) THEN
 
151
        DO IHEL=1,NCOMB
 
152
          GOODHEL(IHEL)=.TRUE.
 
153
        ENDDO
 
154
      ENDIF
138
155
      ANS = 0D0
139
156
      DO IHEL=1,NCOMB
140
157
        IF (USERHEL.EQ.-1.OR.USERHEL.EQ.IHEL) THEN
163
180
C     Returns amplitude squared summed/avg over colors
164
181
C     for the point with external lines W(0:6,NEXTERNAL)
165
182
C     
166
 
C     Process: g g > w- t b~ QED=1 QCD=2 [ virt = QCD ]
 
183
C     Process: g g > w- t b~ QED<=1 QCD<=2 [ virt = QCD ]
167
184
C     
168
185
      IMPLICIT NONE
169
186
C     
295
312
      RETURN
296
313
      END
297
314
 
 
315