def printHplusCrossSections(tanbetas=[10, 20, 30, 40], mu=defaultMu):
print "ttbar cross section %.1f pb" % ttCrossSection
print "mu %.1f" % mu
print
for tanbeta in [10, 20, 30, 40]:
print "=" * 98
print "tan(beta) = %d" % tanbeta
print "H+ M (GeV) | BR(t->bH+) | BR(H+->taunu) | sigma(tt->tbH+->tbtaunu) | sigma(tt->bbH+H-->bbtautau) |"
for mass in [90, 100, 120, 140, 150, 155, 160]:
br_tH = br.getBR_top2bHp(mass, tanbeta, mu)
br_Htaunu = br.getBR_Hp2tau(mass, tanbeta, mu)
print "%10d | %10f | %13f | %24f | %27f |" % (
mass, br_tH, br_Htaunu,
whTauNuCrossSectionMSSM(mass, tanbeta, mu),
hhTauNuCrossSectionMSSM(mass, tanbeta, mu))
def hhTauNuCrossSectionMSSM(mass, tanbeta, mu):
br_tH = br.getBR_top2bHp(mass, tanbeta, mu)
br_Htaunu = br.getBR_Hp2tau(mass, tanbeta, mu)
return hhTauNuCrossSection(br_tH, br_Htaunu)