← Back to branch summary

~maddevelopers/mg5amcnlo/simple_unlops

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_group/matrix1.f

  • Committer: olivier-mattelaer
  • Date: 2021-11-12 09:29:31 UTC
  • mfrom: (967.1.15 3.3.0)
  • Revision ID: olivier-mattelaer-20211112092931-4ec9qfzgxkeovqog
version 3.3.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportV4IOTest/export_matrix_element_v4_madevent_group/matrix1.f
166
166
     $     .LE.MAXTRIES.OR.(ISUM_HEL.NE.0).OR.THIS_NTRY(IMIRROR).LE.10)
167
167
     $      THEN
168
168
            T=MATRIX1(P ,NHEL(1,I),JC(1),I)
 
169
 
169
170
            DO JJ=1,NINCOMING
170
171
              IF(POL(JJ).NE.1D0.AND.NHEL(JJ,I).EQ.INT(SIGN(1D0,POL(JJ))
171
172
     $         )) THEN
172
 
                T=T*ABS(POL(JJ))*NB_SPIN_STATE_IN(JJ)/2D0  ! NB_SPIN_STATE(JJ)/2d0 is added for polarised beam
 
173
                T=T*ABS(POL(JJ))
173
174
              ELSE IF(POL(JJ).NE.1D0)THEN
174
 
                T=T*(2D0-ABS(POL(JJ)))*NB_SPIN_STATE_IN(JJ)/2D0
 
175
                T=T*(2D0-ABS(POL(JJ)))
175
176
              ENDIF
176
177
            ENDDO
 
178
 
177
179
            IF (ISUM_HEL.NE.0.AND.DS_GET_DIM_STATUS('Helicity')
178
180
     $       .EQ.0.AND.ALLOW_HELICITY_GRID_ENTRIES) THEN
179
181
              CALL DS_ADD_ENTRY('Helicity',I,T)
229
231
 
230
232
        T=MATRIX1(P ,NHEL(1,I),JC(1),I)
231
233
 
 
234
 
232
235
        DO JJ=1,NINCOMING
233
236
          IF(POL(JJ).NE.1D0.AND.NHEL(JJ,I).EQ.INT(SIGN(1D0,POL(JJ))))
234
237
     $      THEN
237
240
            T=T*(2D0-ABS(POL(JJ)))
238
241
          ENDIF
239
242
        ENDDO
 
243
 
240
244
C       Always one helicity at a time
241
245
        ANS = T
242
246
C       Include the Jacobian from helicity sampling
343
347
      COMPLEX*16 DUM0,DUM1
344
348
      DATA DUM0, DUM1/(0D0, 0D0), (1D0, 0D0)/
345
349
 
 
350
      DOUBLE PRECISION FK_WZ
346
351
      DOUBLE PRECISION FK_ZERO
347
 
      DOUBLE PRECISION FK_WZ
 
352
      SAVE FK_WZ
348
353
      SAVE FK_ZERO
349
 
      SAVE FK_WZ
350
354
 
351
355
      LOGICAL FIRST
352
356
      DATA FIRST /.TRUE./