coupling_comment:='( (g3 " strong interaction ")  (g "

weak interaction ") )
 $
delta_qnumber_list:='nil
 $
pa_decay_list:='((pr3 (67 (pr2 pr1) 2) (15 (ef ef) 2)) (pr4
(64 (pr2 pr1) 2) (16 (ef ef) 2)) (pr5 (60 (pr2 pr1) 2)) (pr6
(55 (pr2 pr1) 2) (53 (pr3 pr1) 2)) (pr7 (48 (pr2 pr1) 2) (46
(pr3 pr1) 2) (17 (ef ef) 2)) (pr8 (40 (pr2 pr1) 2) (38 (pr3
pr1) 2) (36 (pr4 pr1) 2)) (pr9 (33 (pr2 pr1) 2) (31 (pr3 pr1
) 2) (29 (pr4 pr1) 2) (18 (ef ef) 2)) (pr10 (81 (pr2 pr1) 2)
(79 (pr3 pr1) 2) (77 (pr4 pr1) 2)) (pr11 (26 (pr2 pr1) 2) (
24 (pr3 pr1) 2) (22 (pr4 pr1) 2) (19 (ef ef) 2)) (pr12 (74 (
pr2 pr1) 2) (72 (pr3 pr1) 2) (70 (pr4 pr1) 2) (68 (pr5 pr1)
2)) (pr13 (80 (pr10 pr2) 2) (78 (pr10 pr3) 2) (76 (pr10 pr4)
2) (75 (pr10 pr5) 2) (73 (pr12 pr2) 2) (71 (pr12 pr3) 2) (69
(pr12 pr4) 2) (66 (pr3 pr2) 2) (65 (pr3 pr3) 2) (63 (pr4 pr2
) 2) (62 (pr4 pr3) 2) (61 (pr4 pr4) 2) (59 (pr5 pr2) 2) (58
(pr5 pr3) 2) (57 (pr5 pr4) 2) (56 (pr5 pr5) 2) (54 (pr6 pr2)
2) (52 (pr6 pr3) 2) (51 (pr6 pr4) 2) (50 (pr6 pr5) 2) (49 (
pr6 pr6) 2) (47 (pr7 pr2) 2) (45 (pr7 pr3) 2) (44 (pr7 pr4)
2) (43 (pr7 pr5) 2) (42 (pr7 pr6) 2) (41 (pr7 pr7) 2) (39 (
pr8 pr2) 2) (37 (pr8 pr3) 2) (35 (pr8 pr4) 2) (34 (pr8 pr5)
2) (32 (pr9 pr2) 2) (30 (pr9 pr3) 2) (28 (pr9 pr4) 2) (27 (
pr9 pr5) 2) (25 (pr11 pr2) 2) (23 (pr11 pr3) 2) (21 (pr11 
pr4) 2) (20 (pr11 pr5) 2)))
 $
two_pa_vertex_list:='nil
 $
stable_pal0:='(pr2 p nut numu nue ef qu qd qs mu pr1 pr3 qc
pr4 pr5 pr6 pr7 pr8 pr9 pr10 pr11 tau pr12 pr13 qb w qt)
 $
mass_order_list:='(h0 qt z w qb pr13 pr12 tau pr11 pr10 pr9
pr8 pr7 pr6 pr5 pr4 qc pr3 pr1 mu qs qd qu ef nue numu nut 
gsg p gs pr2)
 $
pal_order_list:='(pr13 pr1 pr3 pr2 pr4 pr5 pr7 w pr11 pr9 
pr6 h0 pr8 pr12 pr10 ef pa1 z pa5 pa4 p pa2 pa3)
 $
v_pal_list:='((15 (pr3 -1 1) (ef 1 1)) (16 (pr4 -1 1) (ef 1
1)) (17 (pr7 -1 1) (ef 1 1)) (18 (pr9 -1 1) (ef 1 1)) (19 (
pr11 -1 1) (ef 1 1)) (14 (h0 -1 3)) (13 (h0 -1 4)) (12 (pa2
1 1) (h0 -1 1)) (11 (w 1 1) (pa1 -1 1)) (10 (pa5 1 0) (pa4 0
1) (w 0 1)) (9 (pa4 1 0) (pa5 0 1) (w 1 0)) (8 (pa3 1 1) (
pa1 -1 1)) (7 (w 2 2)) (6 (w 1 1) (p -1 1) (z -1 1)) (5 (w 1
1) (pa1 -1 2)) (4 (w 1 1) (h0 -1 1)) (3 (z -1 2) (h0 -1 1))
(2 (w 1 1) (h0 -1 2)) (1 (z -1 2) (h0 -1 2)) (20 (pr11 -1 1)
(pr5 -1 1) (pr13 -1 1)) (21 (pr11 -1 1) (pr4 -1 1) (pr13 -1
1)) (22 (pr11 -1 1) (pr4 -1 1) (pr1 -1 1)) (23 (pr11 -1 1) (
pr3 -1 1) (pr13 -1 1)) (24 (pr11 -1 1) (pr3 -1 1) (pr1 -1 1)
) (25 (pr11 -1 1) (pr2 -1 1) (pr13 -1 1)) (26 (pr11 -1 1) (
pr2 -1 1) (pr1 -1 1)) (27 (pr9 -1 1) (pr5 -1 1) (pr13 -1 1))
(28 (pr9 -1 1) (pr4 -1 1) (pr13 -1 1)) (29 (pr9 -1 1) (pr4 
-1 1) (pr1 -1 1)) (30 (pr9 -1 1) (pr3 -1 1) (pr13 -1 1)) (31
(pr9 -1 1) (pr3 -1 1) (pr1 -1 1)) (32 (pr9 -1 1) (pr2 -1 1)
(pr13 -1 1)) (33 (pr9 -1 1) (pr2 -1 1) (pr1 -1 1)) (34 (pr8
-1 1) (pr5 -1 1) (pr13 -1 1)) (35 (pr8 -1 1) (pr4 -1 1) (
pr13 -1 1)) (36 (pr8 -1 1) (pr4 -1 1) (pr1 -1 1)) (37 (pr8 
-1 1) (pr3 -1 1) (pr13 -1 1)) (38 (pr8 -1 1) (pr3 -1 1) (pr1
-1 1)) (39 (pr8 -1 1) (pr2 -1 1) (pr13 -1 1)) (40 (pr8 -1 1)
(pr2 -1 1) (pr1 -1 1)) (41 (pr7 -1 2) (pr13 -1 1)) (42 (pr7
-1 1) (pr6 -1 1) (pr13 -1 1)) (43 (pr7 -1 1) (pr5 -1 1) (
pr13 -1 1)) (44 (pr7 -1 1) (pr4 -1 1) (pr13 -1 1)) (45 (pr7
-1 1) (pr3 -1 1) (pr13 -1 1)) (46 (pr7 -1 1) (pr3 -1 1) (pr1
-1 1)) (47 (pr7 -1 1) (pr2 -1 1) (pr13 -1 1)) (48 (pr7 -1 1)
(pr2 -1 1) (pr1 -1 1)) (49 (pr6 -1 2) (pr13 -1 1)) (50 (pr6
-1 1) (pr5 -1 1) (pr13 -1 1)) (51 (pr6 -1 1) (pr4 -1 1) (
pr13 -1 1)) (52 (pr6 -1 1) (pr3 -1 1) (pr13 -1 1)) (53 (pr6
-1 1) (pr3 -1 1) (pr1 -1 1)) (54 (pr6 -1 1) (pr2 -1 1) (pr13
-1 1)) (55 (pr6 -1 1) (pr2 -1 1) (pr1 -1 1)) (56 (pr5 -1 2)
(pr13 -1 1)) (57 (pr5 -1 1) (pr4 -1 1) (pr13 -1 1)) (58 (pr5
-1 1) (pr3 -1 1) (pr13 -1 1)) (59 (pr5 -1 1) (pr2 -1 1) (
pr13 -1 1)) (60 (pr5 -1 1) (pr2 -1 1) (pr1 -1 1)) (61 (pr4 
-1 2) (pr13 -1 1)) (62 (pr4 -1 1) (pr3 -1 1) (pr13 -1 1)) (
63 (pr4 -1 1) (pr2 -1 1) (pr13 -1 1)) (64 (pr4 -1 1) (pr2 -1
1) (pr1 -1 1)) (65 (pr3 -1 2) (pr13 -1 1)) (66 (pr3 -1 1) (
pr2 -1 1) (pr13 -1 1)) (67 (pr3 -1 1) (pr2 -1 1) (pr1 -1 1))
(68 (pr12 -1 1) (pr5 -1 1) (pr1 -1 1)) (69 (pr12 -1 1) (pr4
-1 1) (pr13 -1 1)) (70 (pr12 -1 1) (pr4 -1 1) (pr1 -1 1)) (
71 (pr12 -1 1) (pr3 -1 1) (pr13 -1 1)) (72 (pr12 -1 1) (pr3
-1 1) (pr1 -1 1)) (73 (pr12 -1 1) (pr2 -1 1) (pr13 -1 1)) (
74 (pr12 -1 1) (pr2 -1 1) (pr1 -1 1)) (75 (pr10 -1 1) (pr5 
-1 1) (pr13 -1 1)) (76 (pr10 -1 1) (pr4 -1 1) (pr13 -1 1)) (
77 (pr10 -1 1) (pr4 -1 1) (pr1 -1 1)) (78 (pr10 -1 1) (pr3 
-1 1) (pr13 -1 1)) (79 (pr10 -1 1) (pr3 -1 1) (pr1 -1 1)) (
80 (pr10 -1 1) (pr2 -1 1) (pr13 -1 1)) (81 (pr10 -1 1) (pr2
-1 1) (pr1 -1 1)))
 $
pa_rv_list:='((pr10 (g (1 75 76 77 78 79 80 81))) (pr2 (g (1
25 26 32 33 39 40 47 48 54 55 59 60 63 64 66 67 73 74 80 81)
)) (pr1 (g (1 22 24 26 29 31 33 36 38 40 46 48 53 55 60 64 
67 68 70 72 74 77 79 81))) (pr13 (g (1 20 21 23 25 27 28 30
32 34 35 37 39 41 42 43 44 45 47 49 50 51 52 54 56 57 58 59
61 62 63 65 66 69 71 73 75 76 78 80))) (pr3 (g (1 15 23 24 
30 31 37 38 45 46 52 53 58 62 65 66 67 71 72 78 79))) (pr4 (
g (1 16 21 22 28 29 35 36 44 51 57 61 62 63 64 69 70 76 77))
) (pr5 (g (1 20 27 34 43 50 56 57 58 59 60 68 75))) (pr12 (g
(1 68 69 70 71 72 73 74))) (pr6 (g (1 42 49 50 51 52 53 54 
55))) (pr7 (g (1 17 41 42 43 44 45 46 47 48))) (pr8 (g (1 34
35 36 37 38 39 40))) (pr9 (g (1 18 27 28 29 30 31 32 33))) (
pr11 (g (1 19 20 21 22 23 24 25 26))) (z (g (2 6 1) (1 3)))
(h0 (g (2 13 2 1) (1 14 12 4 3))) (w (g (2 7 6 5 2) (1 11 10
9 4))) (pa3 (g (1 8))) (pa1 (g (2 5) (1 11 8))) (pa4 (g (1 
10 9))) (pa5 (g (1 10 9))) (pa2 (g (1 12))) (ef (g (1 15 16
17 18 19))) (p (g (2 6))))
 $
cvpd_list:='((15 0) (16 0) (17 0) (18 0) (19 0) (14 0) (13 0
) (11 1) (7 0) (6 0) (5 0) (4 0) (3 0) (2 0) (1 0) (20 0) (
21 0) (22 0) (23 0) (24 0) (25 0) (26 0) (27 0) (28 0) (29 0
) (30 0) (31 0) (32 0) (33 0) (34 0) (35 0) (36 0) (37 0) (
38 0) (39 2) (40 2) (41 0) (42 0) (43 0) (44 0) (45 0) (46 0
) (47 0) (48 0) (49 0) (50 0) (51 0) (52 0) (53 0) (54 2) (
55 2) (56 0) (57 0) (58 0) (59 2) (60 2) (61 0) (62 0) (63 0
) (64 0) (65 0) (66 0) (67 0) (68 2) (69 2) (70 2) (71 2) (
72 2) (73 2) (74 2) (75 2) (76 2) (77 2) (78 2) (79 2) (80 2
) (81 2))
 $
gbsymmetrylist:='(strange baryon len lmuon ltau)
 $
fit_para_list:='((81 f82) (80 f81) (79 f80) (78 f79) (77 f78
) (76 f77) (75 f76 f75) (74 f74) (73 f73) (72 f72) (71 f71)
(70 f70) (69 f69) (68 f68 f67) (67 f66) (66 f65) (65 f64) (
64 f63) (63 f62) (62 f61) (61 f60) (60 f59) (59 f58) (58 f57
f56) (57 f55 f54) (56 f53) (55 f52) (54 f51) (53 f50 f49) (
52 f48 f47) (51 f46 f45) (50 f44) (49 f43) (48 f42) (47 f41)
(46 f40) (45 f39) (44 f38) (43 f37 f36) (42 f35 f34) (41 f33
) (40 f32) (39 f31) (38 f30 f29) (37 f28 f27) (36 f26 f25) (
35 f24 f23) (34 f22) (33 f21) (32 f20) (31 f19) (30 f18) (29
f17) (28 f16) (27 f15 f14) (26 f13) (25 f12) (24 f11) (23 
f10) (22 f9) (21 f8) (20 f7 f6) (19 f5) (18 f4) (17 f3) (16
f2) (15 f1))
 $
palist:='(pr1 pr2 pr3 pr4 pr5 pr6 pr7 pr8 pr9 pr10 pr11 pr12
pr13 gs w z p gsg h0 qb qs qd qt qc qu tau mu ef nut numu 
nue)
 $
give_up_compact:='nil
 $
any_type_list:='((any_ext pr2 pr1) (any_nf pr1 pr2 pr3 pr4 
pr5 pr6 pr7 pr8 pr9 pr10 pr11 pr12 pr13 gs w z p h0) (any 
pr1 pr2 pr3 pr4 pr5 pr6 pr7 pr8 pr9 pr10 pr11 pr12 pr13 gs w
z p h0 qb qs qd qt qc qu tau mu ef nut numu nue) (any_f qb 
qs qd qt qc qu tau mu ef nut numu nue) (any_l nue numu nut 
ef mu tau) (any_q qu qc qt qd qs qb) (any_b h0 gs w z p) (
any_b1 gs w z p) (any_b0 h0))
 $
effective_numbers:='nil
 $
bound_particle_list:='nil
 $
bound_state_vertex_list:='nil
 $
charge_op_sign_list:='((12 1) (12 1) (19 -1) (18 -1) (17 -1)
(16 -1) (15 -1))
 $
color_vertices:='nil
 $
fcolor_vertices:='nil
 $
bound_pv_list:='nil
 $
fermionprog:='(sf)
 $
bosonprog:='(s11 s1 s0)
 $
propagatorl:='(s3 s11 s1 s0 sf)
 $
splist:='((s11 gs) (s1 pr2 pr3 pr4 pr5 pr6 pr7 pr8 pr9 pr11
w z p) (s0 pr1 pr13 gsg h0) (sf qt qc qu tau mu ef qb qs qd
nut numu nue) (s3 pr10 pr12))
 $
hvvvl:='((81 ((vvv x 81 r1 r2 vc v vb va) (b3 pr10 x 0 r1 v
va vb) (b1 pr2 x 0 r2 vc) (b0 pr1 x 0 r3))) (80 ((vvv x 80 
r1 r2 vc v vb va) (b3 pr10 x 0 r1 v va vb) (b1 pr2 x 0 r2 vc
) (b0 pr13 x 0 r3))) (79 ((vvv x 79 r1 r2 vc v vb va) (b3 
pr10 x 0 r1 v va vb) (b1 pr3 x 0 r2 vc) (b0 pr1 x 0 r3))) (
78 ((vvv x 78 r1 r2 vc v vb va) (b3 pr10 x 0 r1 v va vb) (b1
pr3 x 0 r2 vc) (b0 pr13 x 0 r3))) (77 ((vvv x 77 r1 r2 vc v
vb va) (b3 pr10 x 0 r1 v va vb) (b1 pr4 x 0 r2 vc) (b0 pr1 x
0 r3))) (76 ((vvv x 76 r1 r2 vc v vb va) (b3 pr10 x 0 r1 v 
va vb) (b1 pr4 x 0 r2 vc) (b0 pr13 x 0 r3))) (75 ((vvv x 75
r2 r1 vb va v vc) (b3 pr10 x 0 r1 v va vb) (b1 pr5 x 0 r2 vc
) (b0 pr13 x 0 r3))) (74 ((vvv x 74 r1 r2 vc v vb va) (b3 
pr12 x 0 r1 v va vb) (b1 pr2 x 0 r2 vc) (b0 pr1 x 0 r3))) (
73 ((vvv x 73 r1 r2 vc v vb va) (b3 pr12 x 0 r1 v va vb) (b1
pr2 x 0 r2 vc) (b0 pr13 x 0 r3))) (72 ((vvv x 72 r1 r2 vc v
vb va) (b3 pr12 x 0 r1 v va vb) (b1 pr3 x 0 r2 vc) (b0 pr1 x
0 r3))) (71 ((vvv x 71 r1 r2 vc v vb va) (b3 pr12 x 0 r1 v 
va vb) (b1 pr3 x 0 r2 vc) (b0 pr13 x 0 r3))) (70 ((vvv x 70
r1 r2 vc v vb va) (b3 pr12 x 0 r1 v va vb) (b1 pr4 x 0 r2 vc
) (b0 pr1 x 0 r3))) (69 ((vvv x 69 r1 r2 vc v vb va) (b3 
pr12 x 0 r1 v va vb) (b1 pr4 x 0 r2 vc) (b0 pr13 x 0 r3))) (
68 ((vvv x 68 r2 r1 vb va v vc) (b3 pr12 x 0 r1 v va vb) (b1
pr5 x 0 r2 vc) (b0 pr1 x 0 r3))) (67 ((vvv x 67 r2 r1 va v)
(b1 pr3 x 0 r1 v) (b1 pr2 x 0 r2 va) (b0 pr1 x 0 r3))) (66 (
(vvv x 66 r2 r1 va v) (b1 pr3 x 0 r1 v) (b1 pr2 x 0 r2 va) (
b0 pr13 x 0 r3))) (65 ((vvv x 65 r2 r1 va v) (b1 pr3 x 0 r1
v) (b1 pr3 x 0 r2 va) (b0 pr13 x 0 r3))) (64 ((vvv x 64 r2 
r1 va v) (b1 pr4 x 0 r1 v) (b1 pr2 x 0 r2 va) (b0 pr1 x 0 r3
))) (63 ((vvv x 63 r2 r1 va v) (b1 pr4 x 0 r1 v) (b1 pr2 x 0
r2 va) (b0 pr13 x 0 r3))) (62 ((vvv x 62 r2 r1 va v) (b1 pr4
x 0 r1 v) (b1 pr3 x 0 r2 va) (b0 pr13 x 0 r3))) (61 ((vvv x
61 r2 r1 va v) (b1 pr4 x 0 r1 v) (b1 pr4 x 0 r2 va) (b0 pr13
x 0 r3))) (60 ((vvv x 60 r3 r2 r1 va v) (b1 pr5 x 0 r1 v) (
b1 pr2 x 0 r2 va) (b0 pr1 x 0 r3))) (59 ((vvv x 59 r3 r2 r1
va v) (b1 pr5 x 0 r1 v) (b1 pr2 x 0 r2 va) (b0 pr13 x 0 r3))
) (58 ((vvv x 58 r2 r1 v va) (b1 pr5 x 0 r1 v) (b1 pr3 x 0 
r2 va) (b0 pr13 x 0 r3))) (57 ((vvv x 57 r2 r1 v va) (b1 pr5
x 0 r1 v) (b1 pr4 x 0 r2 va) (b0 pr13 x 0 r3))) (56 ((vvv x
56 r2 r1 va v) (b1 pr5 x 0 r1 v) (b1 pr5 x 0 r2 va) (b0 pr13
x 0 r3))) (55 ((vvv x 55 r3 r2 r1 va v) (b1 pr6 x 0 r1 v) (
b1 pr2 x 0 r2 va) (b0 pr1 x 0 r3))) (54 ((vvv x 54 r3 r2 r1
va v) (b1 pr6 x 0 r1 v) (b1 pr2 x 0 r2 va) (b0 pr13 x 0 r3))
) (53 ((vvv x 53 r2 r1 v va) (b1 pr6 x 0 r1 v) (b1 pr3 x 0 
r2 va) (b0 pr1 x 0 r3))) (52 ((vvv x 52 r2 r1 v va) (b1 pr6
x 0 r1 v) (b1 pr3 x 0 r2 va) (b0 pr13 x 0 r3))) (51 ((vvv x
51 r2 r1 v va) (b1 pr6 x 0 r1 v) (b1 pr4 x 0 r2 va) (b0 pr13
x 0 r3))) (50 ((vvv x 50 r2 r1 va v) (b1 pr6 x 0 r1 v) (b1 
pr5 x 0 r2 va) (b0 pr13 x 0 r3))) (49 ((vvv x 49 r2 r1 va v)
(b1 pr6 x 0 r1 v) (b1 pr6 x 0 r2 va) (b0 pr13 x 0 r3))) (48
((vvv x 48 r2 r1 va v) (b1 pr7 x 0 r1 v) (b1 pr2 x 0 r2 va)
(b0 pr1 x 0 r3))) (47 ((vvv x 47 r2 r1 va v) (b1 pr7 x 0 r1
v) (b1 pr2 x 0 r2 va) (b0 pr13 x 0 r3))) (46 ((vvv x 46 r2 
r1 va v) (b1 pr7 x 0 r1 v) (b1 pr3 x 0 r2 va) (b0 pr1 x 0 r3
))) (45 ((vvv x 45 r2 r1 va v) (b1 pr7 x 0 r1 v) (b1 pr3 x 0
r2 va) (b0 pr13 x 0 r3))) (44 ((vvv x 44 r2 r1 va v) (b1 pr7
x 0 r1 v) (b1 pr4 x 0 r2 va) (b0 pr13 x 0 r3))) (43 ((vvv x
43 r2 r1 v va) (b1 pr7 x 0 r1 v) (b1 pr5 x 0 r2 va) (b0 pr13
x 0 r3))) (42 ((vvv x 42 r2 r1 v va) (b1 pr7 x 0 r1 v) (b1 
pr6 x 0 r2 va) (b0 pr13 x 0 r3))) (41 ((vvv x 41 r2 r1 va v)
(b1 pr7 x 0 r1 v) (b1 pr7 x 0 r2 va) (b0 pr13 x 0 r3))) (40
((vvv x 40 r3 r2 r1 va v) (b1 pr8 x 0 r1 v) (b1 pr2 x 0 r2 
va) (b0 pr1 x 0 r3))) (39 ((vvv x 39 r3 r2 r1 va v) (b1 pr8
x 0 r1 v) (b1 pr2 x 0 r2 va) (b0 pr13 x 0 r3))) (38 ((vvv x
38 r2 r1 v va) (b1 pr8 x 0 r1 v) (b1 pr3 x 0 r2 va) (b0 pr1
x 0 r3))) (37 ((vvv x 37 r2 r1 v va) (b1 pr8 x 0 r1 v) (b1 
pr3 x 0 r2 va) (b0 pr13 x 0 r3))) (36 ((vvv x 36 r2 r1 v va)
(b1 pr8 x 0 r1 v) (b1 pr4 x 0 r2 va) (b0 pr1 x 0 r3))) (35 (
(vvv x 35 r2 r1 v va) (b1 pr8 x 0 r1 v) (b1 pr4 x 0 r2 va) (
b0 pr13 x 0 r3))) (34 ((vvv x 34 r2 r1 va v) (b1 pr8 x 0 r1
v) (b1 pr5 x 0 r2 va) (b0 pr13 x 0 r3))) (33 ((vvv x 33 r2 
r1 va v) (b1 pr9 x 0 r1 v) (b1 pr2 x 0 r2 va) (b0 pr1 x 0 r3
))) (32 ((vvv x 32 r2 r1 va v) (b1 pr9 x 0 r1 v) (b1 pr2 x 0
r2 va) (b0 pr13 x 0 r3))) (31 ((vvv x 31 r2 r1 va v) (b1 pr9
x 0 r1 v) (b1 pr3 x 0 r2 va) (b0 pr1 x 0 r3))) (30 ((vvv x 
30 r2 r1 va v) (b1 pr9 x 0 r1 v) (b1 pr3 x 0 r2 va) (b0 pr13
x 0 r3))) (29 ((vvv x 29 r2 r1 va v) (b1 pr9 x 0 r1 v) (b1 
pr4 x 0 r2 va) (b0 pr1 x 0 r3))) (28 ((vvv x 28 r2 r1 va v)
(b1 pr9 x 0 r1 v) (b1 pr4 x 0 r2 va) (b0 pr13 x 0 r3))) (27
((vvv x 27 r2 r1 v va) (b1 pr9 x 0 r1 v) (b1 pr5 x 0 r2 va)
(b0 pr13 x 0 r3))) (26 ((vvv x 26 r2 r1 va v) (b1 pr11 x 0 
r1 v) (b1 pr2 x 0 r2 va) (b0 pr1 x 0 r3))) (25 ((vvv x 25 r2
r1 va v) (b1 pr11 x 0 r1 v) (b1 pr2 x 0 r2 va) (b0 pr13 x 0
r3))) (24 ((vvv x 24 r2 r1 va v) (b1 pr11 x 0 r1 v) (b1 pr3
x 0 r2 va) (b0 pr1 x 0 r3))) (23 ((vvv x 23 r2 r1 va v) (b1
pr11 x 0 r1 v) (b1 pr3 x 0 r2 va) (b0 pr13 x 0 r3))) (22 ((
vvv x 22 r2 r1 va v) (b1 pr11 x 0 r1 v) (b1 pr4 x 0 r2 va) (
b0 pr1 x 0 r3))) (21 ((vvv x 21 r2 r1 va v) (b1 pr11 x 0 r1
v) (b1 pr4 x 0 r2 va) (b0 pr13 x 0 r3))) (20 ((vvv x 20 r2 
r1 v va) (b1 pr11 x 0 r1 v) (b1 pr5 x 0 r2 va) (b0 pr13 x 0
r3))) (1 ((vvv x 1 v va) (b1 z x 0 r4 va) (b1 z x 0 r3 v) (
b0 h0 x 0 r1) (b0 h0 x 0 r2))) (2 ((vvv x 2 v va) (b1 w x -1
r4 va) (b1 w x 1 r3 v) (b0 h0 x 0 r2) (b0 h0 x 0 r1))) (3 ((
vvv x 3 v va) (b1 z x 0 r3 va) (b1 z x 0 r2 v) (b0 h0 x 0 r1
))) (4 ((vvv x 4 v va) (b1 w x -1 r3 va) (b1 w x 1 r2 v) (b0
h0 x 0 r1))) (5 ((vvv x 5 pa1 v va vb vc) (b1 w x 1 r1 v) (
b1 w x -1 r2 va) (b1 pa1 x 0 r4 vc) (b1 pa1 x 0 r3 vb))) (6
((vvv x 6 v va vb vc) (b1 w x 1 r1 v) (b1 w x -1 r2 va) (b1
p x 0 r4 vc) (b1 z x 0 r3 vb))) (7 ((vvv x 7 v va vb vc) (b1
w x -1 r4 vc) (b1 w x -1 r3 vb) (b1 w x 1 r2 va) (b1 w x 1 
r1 v))) (8 ((vvv x 8 l pa1 pa3 v) (psi pa3 x -1 r1 id1) (psi
pa3 x 1 r2 id1) (b1 pa1 x 0 r3 v))) (9 ((vvv x 9 l pa5 pa4 v
) (psi pa4 x -1 r1 id1) (psi pa5 x 1 r2 id1) (b1 w x -1 r3 v
))) (10 ((vvv x 10 l pa4 pa5 v) (psi pa5 x -1 r1 id1) (psi 
pa4 x 1 r2 id1) (b1 w x 1 r3 v))) (11 ((vvv x 11 pa1 v va vb
r1 r2 r3) (b1 w x 1 r1 va) (b1 w x -1 r2 vb) (b1 pa1 x 0 r3
v))) (12 ((vvv x 12 l pa2) (psi pa2 x -1 r1 id1) (psi pa2 x
1 r2 id1) (b0 h0 x 0 r3))) (13 ((vvv x 13) (b0 h0 x 0 r4) (
b0 h0 x 0 r3) (b0 h0 x 0 r2) (b0 h0 x 0 r1))) (14 ((vvv x 14
) (b0 h0 x 0 r3) (b0 h0 x 0 r2) (b0 h0 x 0 r1))) (19 ((vvv x
19 l v) (b1 pr11 x 0 r1 v) (psi ef x -1 r2) (psi ef x 1 r3))
) (18 ((vvv x 18 l v) (b1 pr9 x 0 r1 v) (psi ef x -1 r2) (
psi ef x 1 r3))) (17 ((vvv x 17 l v) (b1 pr7 x 0 r1 v) (psi
ef x -1 r2) (psi ef x 1 r3))) (16 ((vvv x 16 l v) (b1 pr4 x
0 r1 v) (psi ef x -1 r2) (psi ef x 1 r3))) (15 ((vvv x 15 l
v) (b1 pr3 x 0 r1 v) (psi ef x -1 r2) (psi ef x 1 r3))))
 $
htionlist:='((g (1 81 80 79 78 77 76 75 74 73 72 71 70 69 68
67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48
47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28
27 26 25 24 23 22 21 20 3 4 8 9 10 11 12 14 19 18 17 16 15)
(2 1 2 5 6 7 13)))
 $
couplist:='(g3 g)
 $
maslist2:='((gsg mgsg 0) (z zm (quotient wm (cos theta))) (
gs mgs 0))
 $
paramlist:='(((c 0 1) (quotient (expt hm 2) 2)) ((c 0 5) (
quotient (minus (times (expt g 2) (expt hm 2))) (times 8 (
expt wm 2)))) (wangm3 (quotient (times (sqrt 2) fmb g) (
times 2 wm))) (wangm2 (quotient (times (sqrt 2) fms g) (
times 2 wm))) (wangm1 (quotient (times (sqrt 2) fmd g) (
times 2 wm))) (wangm6 (quotient (times (sqrt 2) fmtau g) (
times 2 wm))) (wangm5 (quotient (times (sqrt 2) fmmu g) (
times 2 wm))) (wangm4 (quotient (times (sqrt 2) fme g) (
times 2 wm))) (wangm9 (quotient (times (sqrt 2) fmt g) (
times 2 wm))) (wangm8 (quotient (times (sqrt 2) fmc g) (
times 2 wm))) (wangm7 (quotient (times (sqrt 2) fmu g) (
times 2 wm))) (v0 (quotient (times 2 wm) g)))
 $
phyinput:='(p1m w1d p3m w3d p4m w4d p5m w5d p6m w6d p7m w7d
p8m w8d p9m w9d p10m w10d p11m w11d p12m w12d p13m w13d g3 g
theta wm hm fme fmmu fmtau fmd fmu fms fmc fmb fmt)
 $
colorlist:='(pa6 pa4 pa5 pa2 pa3 qd qs qb qu qc qt gs)
 $
hnsuml:='(pr13 pr12 pr11 pr10 pr9 pr8 pr7 pr6 pr5 pr4 pr3 
pr2 pr1 h0 gsg w gs)
 $
hsuml:='(pa3 pa2 pa5 pa4 pa6 pa1)
 $
hparameter:='((spacetime x) (sum pa3 pa2 pa5 pa4 pa6 pa1) (
nsum pr13 pr12 pr11 pr10 pr9 pr8 pr7 pr6 pr5 pr4 pr3 pr2 pr1
h0 gsg w gs) (single p z qt qc qu qb qs qd nut numu nue tau
mu ef) (index3 ic1 ic2 ic3 ic4 ic5) (index2 id1 id2) (index1
v va vb vc vd vf va1 va2))
 $
pnametype:='(sum nsum single)
 $
hway:='(number dotnumber wayorder nameorder dotnameorder)
 $
hparaway:='((index3 wayorder) (index2 wayorder) (index1 
wayorder) (spacetime number) (sum dotnumber) (nsum 
dotwayorder) (single mdotnumber))
 $
hpsumrange:='((pa3 ef mu tau nue numu nut qd qs qb qu qc qt)
(pa2 qt qc qu qb qs qd tau mu ef) (pa5 nut numu nue qt qc qu
) (pa4 tau mu ef qb qs qd) (pa6 qu qc qt) (pa1 z p))
 $
hpindex:='(spacetime indices sum)
 $
singlist:='(pr1 pr2 pr3 pr4 pr5 pr6 pr7 pr8 pr9 pr10 pr11 
pr12 pr13 gs w z p gsg h0 qb qs qd qt qc qu tau mu ef nut 
numu nue)
 $
rangetable:='((pa5 pa6 1) (pa2 pa6 1) (pa2 pa4 1) (pa5 pa2 
pa6) (pa3 pa6 1) (pa3 pa4 1) (pa3 pa5 1) (pa3 pa2 1))
 $
hconstr:='((((0 pa3 pa1) (1 z nut numu nue) (2 pa3 z) (2 pa2
p) (1 pa1 pa2 qt qc qu qb qs qd tau mu ef pa4 pa6)) 8) (((0
pa5 pa4) (0 tau nut) (0 mu numu) (0 ef nue) (0 qb qt) (0 qs
qc) (0 qd qu)) 9 10))
 $
hconstr1:='(((pa3 pa1) 8) ((pa5 pa4) 9 10))
 $
patable1:='((b1 s11) (b1 s1) (b0 s0) (psi sf) (b3 s3))
 $
end;
 $
