limit_mmbb1=ymmbb1[i] mintanbeta=0.1 maxtanbeta=5 n=100 step=(maxtanbeta-mintanbeta)/n BRmm=1 BRtt=1 BRbb=1 a_ymmbb=[] a_x=[] a_ymmtt=[] for i in range(0,100): tanbeta=mintanbeta+step*i width=get_total_width(args.model,float(args.ma),tanbeta) BRmm=gamma_mu(tanbeta,float(args.ma),args.model)/width BRtt=gamma_tau(tanbeta,args.ma,args.model)/width BRbb=gamma_quarks(tanbeta,args.ma,args.model,6)/width a_ymmbb.append(limit_mmbb1*0.00017/(2*BRbb*BRmm)) a_ymmtt.append(limit_mmtt1/(BRtt*BRtt)) #print limit_mmbb1,limit_mmtt1 #print BRbb,BRmm,BRtt #print limit_mmbb1*0.00017/(2*BRbb*BRmm),limit_mmtt1/(BRtt*BRtt) a_x.append(tanbeta) x = array("d", a_x) ymmbb = array("d", a_ymmbb) ymmtt = array("d", a_ymmtt) gmmtt = ROOT.TGraph(len(x),x,ymmtt)
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?") args = parser.parse_args() style1=GetStyleHtt() style1.cd() #### h->aa->mmtt #### x_mmtt1, y_mmtt1 = np.loadtxt('mmtt.txt', unpack=True) x_mmtt=array("d",x_mmtt1) y_mmtt=array("d",y_mmtt1) for i in range(0,len(x_mmtt)): width=get_total_width(args.model,float(x_mmtt[i]),args.tanbeta) BRtt=gamma_tau(args.tanbeta,float(x_mmtt[i]),args.model)/width y_mmtt[i]=y_mmtt[i]/(BRtt*BRtt) gmmtt = ROOT.TGraph(len(x_mmtt), x_mmtt,y_mmtt) #### h->aa->mmbb #### x_mmbb, y_mmbb = np.loadtxt('mmbb.txt', unpack=True) for i in range(0,len(x_mmbb)): width=get_total_width(args.model,float(x_mmbb[i]),args.tanbeta) BRmm=gamma_mu(args.tanbeta,float(x_mmbb[i]),args.model)/width BRbb=gamma_quarks(args.tanbeta,float(x_mmbb[i]),args.model,6)/width y_mmbb[i]=y_mmbb[i]*0.00017 y_mmbb[i]=y_mmbb[i]/(2*BRmm*BRbb) gmmbb = ROOT.TGraph(len(x_mmbb), x_mmbb.flatten('C'),y_mmbb.flatten('C')) #### h->aa->tttt (HIG-14-019) ####
mintanbeta=0.1 maxtanbeta=5 n=100 step=(maxtanbeta-mintanbeta)/n BRmm=1 BRtt=1 BRbb=1 a_y1=[] a_y2=[] a_y3=[] a_y4=[] a_x=[] for i in range(0,100): tanbeta=mintanbeta+step*i width1=get_total_width(1,float(args.ma),tanbeta) width2=get_total_width(2,float(args.ma),tanbeta) width3=get_total_width(3,float(args.ma),tanbeta) width4=get_total_width(4,float(args.ma),tanbeta) BRmm1=gamma_mu(tanbeta,float(args.ma),1)/width1 BRtt1=gamma_tau(tanbeta,args.ma,1)/width1 BRbb1=gamma_quarks(tanbeta,args.ma,1,6)/width1 BRmm2=gamma_mu(tanbeta,float(args.ma),2)/width2 BRtt2=gamma_tau(tanbeta,args.ma,2)/width2 BRbb2=gamma_quarks(tanbeta,args.ma,2,6)/width2 BRmm3=gamma_mu(tanbeta,float(args.ma),3)/width3 BRtt3=gamma_tau(tanbeta,args.ma,3)/width3 BRbb3=gamma_quarks(tanbeta,args.ma,3,6)/width3 BRmm4=gamma_mu(tanbeta,float(args.ma),4)/width4 BRtt4=gamma_tau(tanbeta,args.ma,4)/width4 BRbb4=gamma_quarks(tanbeta,args.ma,4,6)/width4
minbeta=1.45 maxbeta=5 if (args.model==4): minbeta=0.28 maxbeta=0.9 x_mmtt1, y_mmtt1 = np.loadtxt('mmtt.txt', unpack=True) x_mmtt=array("d",x_mmtt1) y_mmtt=array("d",y_mmtt1) hist=ROOT.TH2F("hist","hist",len(x_mmtt)-1,x_mmtt[0],x_mmtt[len(x_mmtt)-1],binbeta,minbeta,maxbeta) for b in range(0,binbeta+1): tanbeta=0.001+minbeta+1.0*b*(maxbeta-minbeta)/(binbeta) for i in range(0,len(x_mmtt)): #for b in range(0,binbeta): #tanbeta=minbeta+1.0*b*(maxbeta-minbeta)/binbeta width=get_total_width(args.model,float(x_mmtt[i]),tanbeta) BRtt=gamma_tau(tanbeta,float(x_mmtt[i]),args.model)/width y=y_mmtt[i]/(BRtt*BRtt) print x_mmtt[i],tanbeta,y hist.Fill(x_mmtt[i],1.0*tanbeta,y) #hist.SetContour(500) canvas = MakeCanvas("asdf","asdf",800,800) canvas.SetRightMargin(0.2) canvas.SetLeftMargin(0.15) canvas.cd() canvas.SetLogz() hist.GetXaxis().SetTitle("m_{a} (GeV)") hist.GetYaxis().SetTitle("tan#beta") hist.GetZaxis().SetTitle("#frac{#sigma(h)}{#sigma_{SM}} #times B(h#rightarrowaa)") if (args.model==3):