[pptt] in pp_tt.mad (and then gg_tt.mad) matrix1.f, add a sanity chec…

…k that the helicities from IMIROR=1,2 match (madgraph5#872 madgraph5#935 madgraph5#941)
valassi · Jul 31, 2024 · 681f6ac · 681f6ac
1 parent 1e43794
commit 681f6ac
Show file tree

Hide file tree

Showing 2 changed files with 36 additions and 0 deletions.
diff --git a/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx/matrix1.f b/epochX/cudacpp/gg_tt.mad/SubProcesses/P1_gg_ttx/matrix1.f
@@ -232,6 +232,24 @@ SUBROUTINE SMATRIX1(P, RHEL, RCOL, CHANNEL, IVEC, ANS, IHEL,
               WRITE (6,*) 'IMIRROR =', IMIRROR
               WRITE (6,*) 'NGOODHEL =', NGOODHEL(IMIRROR)
               WRITE (6,*) 'NCOMB =', NCOMB
+C           SANITY CHECK: check that the two lists of good helicities are identical
+C           (see madgraph4gpu #872, #935, #941)
+              IF (NGOODHEL(1).NE.-1 .AND. NGOODHEL(2).NE.-1) THEN
+                IF (NGOODHEL(1) .NE. NGOODHEL(2)) THEN
+                  WRITE (6,*) 'ERROR! IMIRROR=1,2 mismatch for NGOODHEL',
+     $              NGOODHEL(1), NGOODHEL(2)
+                  STOP
+                ENDIF
+                DO I=1,NCOMB
+                  IF (GOODHEL(I,1) .NEQV. GOODHEL(I,2)) THEN
+                    WRITE (6,*) 'ERROR! IMIRROR=1,2 mismatch for GOODHEL',
+     $                I, GOODHEL(I,1), GOODHEL(I,2)
+                    STOP
+                  ENDIF
+                END DO
+                WRITE (6,*) 'Helicity lists from IMIRROR=1,2 match',
+     $            NGOODHEL(1), NGOODHEL(2)
+              ENDIF
             ENDIF
           ENDIF
         ENDIF

diff --git a/epochX/cudacpp/pp_tt.mad/SubProcesses/P1_uux_ttx/matrix1.f b/epochX/cudacpp/pp_tt.mad/SubProcesses/P1_uux_ttx/matrix1.f
@@ -235,6 +235,24 @@ SUBROUTINE SMATRIX1(P, RHEL, RCOL, CHANNEL, IVEC, ANS, IHEL,
               WRITE (6,*) 'IMIRROR =', IMIRROR
               WRITE (6,*) 'NGOODHEL =', NGOODHEL(IMIRROR)
               WRITE (6,*) 'NCOMB =', NCOMB
+C           SANITY CHECK: check that the two lists of good helicities are identical
+C           (see madgraph4gpu #872, #935, #941)
+              IF (NGOODHEL(1).NE.-1 .AND. NGOODHEL(2).NE.-1) THEN
+                IF (NGOODHEL(1) .NE. NGOODHEL(2)) THEN
+                  WRITE (6,*) 'ERROR! IMIRROR=1,2 mismatch for NGOODHEL',
+     $              NGOODHEL(1), NGOODHEL(2)
+                  STOP
+                ENDIF
+                DO I=1,NCOMB
+                  IF (GOODHEL(I,1) .NEQV. GOODHEL(I,2)) THEN
+                    WRITE (6,*) 'ERROR! IMIRROR=1,2 mismatch for GOODHEL',
+     $                I, GOODHEL(I,1), GOODHEL(I,2)
+                    STOP
+                  ENDIF
+                END DO
+                WRITE (6,*) 'Helicity lists from IMIRROR=1,2 match',
+     $            NGOODHEL(1), NGOODHEL(2)
+              ENDIF
             ENDIF
           ENDIF
         ENDIF