74 {
75
76
77 if (ps == 0.) return;
78
79
81
82
83 else if (id1Abs == 22) {
84 double chgI3 = (id2Abs%2 == 0) ? 0.5 : -0.5;
85 double chgY = (id2Abs < 9) ? 1. / 6. : -0.5;
87 widNow = preFac * alpEM *
pow2(chg) / 4.;
88 }
89
90 else if (id1Abs == 23) {
91 double chgI3 = (id2Abs%2 == 0) ? 0.5 : -0.5;
92 double chgY = (id2Abs < 9) ? 1. / 6. : -0.5;
96 }
97
98
99 else if (id1Abs == 24) widNow = preFac * (alpEM *
pow2(
m_coupF)
101
102
103
104 else {
105
106
107 if( id1Abs < 17 && id2Abs < 17 && id3Abs>0 && id3Abs < 17 ){
109 if( id3Abs < 10 ) widNow *= 3.;
110 if( id1Abs == id2Abs && id1Abs == id3Abs ){
111 if( std::abs(idRes)-4000000 < 10 ) widNow *= 4./3.;
112 else widNow *= 2.;
113 }
114 }
115 }
116
117 return;
118 }