def mygbwkftau(mass, fullgam, rmsbeam, yint, phi, w): "yint = 0.20, 0.37, 0.27; phi = 0" if w > mass + 200.: return 0.076 / (w - mass) return gbwkftau.gbwkf(mass, fullgam, rmsbeam, yint, phi, w - mass)
def y3sfunc(area, mass, wide, back, en): had = area * gbwkf.gbwkf(mass, 0.0263, wide, 0.018, 0., en*2000. - mass) tau = area * (0.57 * 0.0239) * gbwkftau.gbwkf(mass, 0.053, wide, 0.27, 0., en*2000. - mass) cont = back * 4.717**2 / en**2 return had + tau + cont