← Back to branch summary

~alifson/chiralityflow/trunk

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_wprod_fksew/%SubProcesses%P0_udx_veep%extra_cnt_wrapper.f

  • Committer: andrew.lifson at lu
  • Date: 2021-09-01 15:34:39 UTC
  • Revision ID: andrew.lifson@thep.lu.se-20210901153439-7fasjhav4cp4m88r
testing a new repository of a madgraph folder

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_wprod_fksew/%SubProcesses%P0_udx_veep%extra_cnt_wrapper.f
 
1
      SUBROUTINE EXTRA_CNT(P, ICNT, CNTS)
 
2
C     call the extra cnt corresponding to icnt
 
3
C     may be a dummy function, depending on the process
 
4
      IMPLICIT NONE
 
5
      INTEGER ICNT, I
 
6
      INCLUDE 'nexternal.inc'
 
7
      DOUBLE PRECISION P(0:3, NEXTERNAL-1)
 
8
      INCLUDE 'orders.inc'
 
9
      DOUBLE COMPLEX CNTS(2,NSPLITORDERS)
 
10
 
 
11
      DO I = 1, NSPLITORDERS
 
12
        CNTS(1,I) = DCMPLX(0D0,0D0)
 
13
        CNTS(2,I) = DCMPLX(0D0,0D0)
 
14
      ENDDO
 
15
 
 
16
      RETURN
 
17
      END
 
18
 
 
19
 
 
20
      INTEGER FUNCTION GET_EXTRA_CNT_COLOR(ICNT,IPART)
 
21
C     return the color of the ipart-th particle of 
 
22
C     counterterm icnt
 
23
      INTEGER ICNT IPART
 
24
      INCLUDE 'nexternal.inc'
 
25
      INTEGER NEXTERNALB
 
26
      PARAMETER(NEXTERNALB=NEXTERNAL-1)
 
27
      INTEGER NCNT,I
 
28
      PARAMETER (NCNT=1)
 
29
      INTEGER CNT_COLOR(NCNT,NEXTERNALB)
 
30
      DATA (CNT_COLOR(1,I), I=1,NEXTERNALB) / NEXTERNALB * 1 /
 
31
 
 
32
      IF (ICNT.GT.NCNT.OR.ICNT.LE.0) THEN
 
33
        WRITE(*,*) 'ERROR#1 in get_extra_cnt_color', ICNT
 
34
        STOP
 
35
      ENDIF
 
36
      IF (IPART.GE.NEXTERNAL.OR.IPART.LE.0) THEN
 
37
        WRITE(*,*) 'ERROR#2 in get_extra_cnt_color', IPART
 
38
        STOP
 
39
      ENDIF
 
40
 
 
41
      GET_EXTRA_CNT_COLOR=CNT_COLOR(ICNT,IPART)
 
42
 
 
43
      RETURN
 
44
      END
 
45
 
 
46
 
 
47
      INTEGER FUNCTION GET_EXTRA_CNT_PDG(ICNT,IPART)
 
48
C     return the pdg id of the ipart-th particle of 
 
49
C     counterterm icnt
 
50
      INTEGER ICNT IPART
 
51
      INCLUDE 'nexternal.inc'
 
52
      INTEGER NEXTERNALB
 
53
      PARAMETER(NEXTERNALB=NEXTERNAL-1)
 
54
      INTEGER NCNT,I
 
55
      PARAMETER (NCNT=1)
 
56
      INTEGER CNT_PDG(NCNT,NEXTERNALB)
 
57
      DATA (CNT_PDG(1,I), I=1,NEXTERNALB) / NEXTERNALB * 0 /
 
58
 
 
59
      IF (ICNT.GT.NCNT.OR.ICNT.LE.0) THEN
 
60
        WRITE(*,*) 'ERROR#1 in get_extra_cnt_pdg', ICNT
 
61
        STOP
 
62
      ENDIF
 
63
      IF (IPART.GE.NEXTERNAL.OR.IPART.LE.0) THEN
 
64
        WRITE(*,*) 'ERROR#2 in get_extra_cnt_pdg', IPART
 
65
        STOP
 
66
      ENDIF
 
67
 
 
68
      GET_EXTRA_CNT_PDG=CNT_PDG(ICNT,IPART)
 
69
 
 
70
      RETURN
 
71
      END
 
72
 
 
73
 
 
74
      DOUBLE PRECISION FUNCTION GET_EXTRA_CNT_CHARGE(ICNT,IPART)
 
75
C     return the charge id of the ipart-th particle of 
 
76
C     counterterm icnt
 
77
      INTEGER ICNT IPART
 
78
      INCLUDE 'nexternal.inc'
 
79
      INTEGER NEXTERNALB
 
80
      PARAMETER(NEXTERNALB=NEXTERNAL-1)
 
81
      INTEGER NCNT,I
 
82
      PARAMETER (NCNT=1)
 
83
      DOUBLE PRECISION CNT_CHARGE(NCNT,NEXTERNALB)
 
84
      DATA (CNT_CHARGE(1,I), I=1,NEXTERNALB) / NEXTERNALB * 0D0 /
 
85
 
 
86
      IF (ICNT.GT.NCNT.OR.ICNT.LE.0) THEN
 
87
        WRITE(*,*) 'ERROR#1 in get_extra_cnt_charge', ICNT
 
88
        STOP
 
89
      ENDIF
 
90
      IF (IPART.GE.NEXTERNAL.OR.IPART.LE.0) THEN
 
91
        WRITE(*,*) 'ERROR#2 in get_extra_cnt_charge', IPART
 
92
        STOP
 
93
      ENDIF
 
94
 
 
95
      GET_EXTRA_CNT_CHARGE=CNT_CHARGE(ICNT,IPART)
 
96
 
 
97
      RETURN
 
98
      END
 
99