# if path=="ett" or path=="mtt":
# for i in range(4,7):
# output.SetBinContent(i,-100)
if scale is not None:
output.Scale(scale)
if title is not None:
output.SetTitle(title)
if style:
apply_style(output, style)
return output
if __name__ == "__main__":
# The input files
style1=GetStyleHtt()
style1.cd()
parser = argparse.ArgumentParser()
parser.add_argument('--prefit', action='store_true',
help="Don't use postfit")
parser.add_argument('--period', default="all",
choices=['7TeV', '8TeV', 'all'],
help="Which data taking period")
parser.add_argument('--MLfit', default="all",
choices=['channel', 'all'],
help="Which fit")
if path == "ett" or path == "mtt":
for i in range(4, 7):
output.SetBinContent(i, -100)
if scale is not None:
output.Scale(scale)
if title is not None:
output.SetTitle(title)
if style:
apply_style(output, style)
return output
if __name__ == "__main__":
# The input files
style1 = GetStyleHtt()
style1.cd()
parser = argparse.ArgumentParser()
parser.add_argument('--prefit',
action='store_true',
help="Don't use postfit")
parser.add_argument('--period',
default="all",
choices=['7TeV', '8TeV', 'all'],
help="Which data taking period")
parser.add_argument('--MLfit',
default="all",
if model!=1:
lumi.AddText("tan#beta = "+str(tanbeta))
return lumi
parser = argparse.ArgumentParser()
parser.add_argument('--model', type=int, default='1', help="Which type of 2HDM?")
parser.add_argument('--tanbeta', type=float, default='1', help="Which tan beta?")
parser.add_argument('--ma', type=int, default='1', help="Which pseudoscalar mass?")
args = parser.parse_args()
model=args.model
ma=args.ma
tanbeta=args.tanbeta
style=GetStyleHtt()
style.cd()
#model=1
#ma=40
#tanbeta=5
width=get_total_width(model,ma,tanbeta)
print "running alpha ",running_alpha(ma)
print "up ",gamma_quarks(tanbeta,ma,model,1)/width
print "down ",gamma_quarks(tanbeta,ma,model,2)/width
print "charm ",gamma_quarks(tanbeta,ma,model,3)/width
print "strange ",gamma_quarks(tanbeta,ma,model,4)/width
print "top ",gamma_quarks(tanbeta,ma,model,5)/width
print "bottom ",gamma_quarks(tanbeta,ma,model,6)/width
print "tau ",gamma_tau(tanbeta,ma,model)/width
print "muon ",gamma_mu(tanbeta,ma,model)/width
