2
INTEGER LMAXCONFIGSB_USED
3
PARAMETER (LMAXCONFIGSB_USED=1)
4
INTEGER MAX_BRANCHB_USED
5
PARAMETER (MAX_BRANCHB_USED=2)
6
INTEGER MAPCONFIG(0 : LMAXCONFIGSB_USED)
7
INTEGER IFOREST(2, -MAX_BRANCHB_USED:-1, LMAXCONFIGSB_USED)
8
INTEGER SPROP(-MAX_BRANCHB_USED:-1, LMAXCONFIGSB_USED)
9
INTEGER TPRID(-MAX_BRANCHB_USED:-1, LMAXCONFIGSB_USED)
10
LOGICAL GFORCEBW(-MAX_BRANCHB_USED : -1, LMAXCONFIGSB_USED)
12
C Here are the congifurations
13
C Process: t > b u d~ [ real = QED QCD ] QCD^2<=2 QED^2<=4
14
C Process: t > b c s~ [ real = QED QCD ] QCD^2<=2 QED^2<=4
1
15
C Diagram 1, Amplitude 1
2
16
DATA MAPCONFIG( 1)/ 1/
3
DATA (IFOREST(I, -1, 1),I=1,2)/ 4, 3/
17
DATA (IFOREST(IFR, -1, 1),IFR=1,2)/ 4, 3/
4
18
DATA SPROP( -1, 1)/ 24/
5
DATA (IFOREST(I, -2, 1),I=1,2)/ -1, 2/
19
DATA (IFOREST(IFR, -2, 1),IFR=1,2)/ -1, 2/
6
20
DATA SPROP( -2, 1)/ 6/
8
22
DATA MAPCONFIG(0)/ 1/
25
DATA GFORCEBW(-1,1)/.FALSE./
26
DATA GFORCEBW(-2,1)/.FALSE./