def plt_A_B_bsgnedm(i, j): # Bsgamma-result and N-EDM in {A,B} plane
mass_axis1, mass_axis2 = i, j
print(ABarray4()[0], mass_axis1, len(ABarray4()))
print(ABarray4()[1], mass_axis2, len(ABarray4()))
resultb = []
resultn = []
resulte = []
#B>Xs+gamma SECTION
for n in np.arange(0, len(ABarray4())):
y3hdm= bsg.BR_B_Xs_gamma(mb,mw,mass_axis1,mass_axis2,\
exe.Y2(*ABarray4()[n] ), exe.complexyfunction(*ABarray4()[n] ),\
exe.Y3(*ABarray4()[n] ), exe.complexyfunction3(*ABarray4()[n] ))
resultb.append(y3hdm / (1e-4))
#Nedm SECTION
nedm3hdm = abs(dn(mass_axis1,mass_axis2, exe.complexyfunction(*ABarray4()[n]),\
exe.complexyfunction3(*ABarray4()[n]) ) / (5.06e13) )\
/ 1e-26
resultn.append(nedm3hdm)
#eedm SECTION
eedm3hdm = abs(de(mass_axis1,mass_axis2,exe.yconjz2(*ABarray4()[n]),\
exe.yconjz3(*ABarray4()[n]) ) /1e-29 )
resulte.append(eedm3hdm)
#########
ned = plt.contourf(exe.A, exe.B, \
np.resize(np.array(resultn).flatten() ,len(np.array(resultn).flatten() ) ).\
reshape(len(exe.B),len(exe.A)) ,\
levels = np.array([0.0,1.8]),colors = ['red'] )
#########
bsgamm = plt.contourf(exe.A, exe.B, \
np.resize(np.array(resultb).flatten() ,len(np.array(resultb).flatten() ) ).\
reshape(len(exe.B),len(exe.A)) ,\
levels = np.array([2.99,3.55]),colors = ['green'] )
#########
# eed = plt.contourf(exe.A, exe.B, \
# np.resize(np.array(resulte).flatten() ,len(np.array(resulte).flatten() ) ).\
# reshape(len(exe.B),len(exe.A)) ,\
# levels = np.array([0.0,1.1]),colors = ['blue'] )
plt.title('BR($\\bar{B} \\to X_{s} \gamma$) and NEDM in '\
+ str("%02d" % mass_axis1) +', ' + str("%02d"% mass_axis2) )
plt.xlabel(exe.readlist[int(exe.read1)])
plt.ylabel(exe.readlist[int(exe.read2)])
# plt.axis([0,60,-1.6,0]) #{tanbeta/tangamma,theta} plane
# plt.axis([0,2 * PI ,-1.6,0]) #{theta,delta} plane
plt.axis([0, 60, 0, 60]) # {tanbeta,tangamma} plane
plt.savefig(
str("%02d" % mass_axis1) + str("%02d" % mass_axis2) + 'bsg.png')
plt.show()
plt.close()
def NEUTRONEDMtext(): #save Neutron edm result in txt file
f = open('Neutron_EDM.txt', 'w')
f.write('%s %s %s %s %s %s\n' % (" N"," MH+1", " MH+2",\
" X_1Y_1*", " X_2Y_2*"," Result") )
for n in range(len(ABarray4())):
# print('dn',n, dn(80,170,exe.complexyfunction(*ABarray4()[n]).imag, \
# exe.complexyfunction3(*ABarray4()[n]).imag) )
f.write( "%5.0f %5.1f %5.1f %10.10e %10.10e %10.15e\n" % (n,80,170,\
- exe.complexyfunction(*ABarray4()[n]),\
- exe.complexyfunction3(*ABarray4()[n]), \
dn(80,170,- exe.complexyfunction(*ABarray4()[n]), \
- exe.complexyfunction3(*ABarray4()[n])) ) )
f.close()
return
def plt_A_B_cpsdiffer(i, j): #Delta-CPS-asymmetry in {A,B} plane
result_deltas = []
mass_axis1, mass_axis2 = i, j
for n in np.arange(0, len(ABarray4())):
cpsdif = bsg.newdifferacps(mb,mw,mass_axis1,mass_axis2,\
exe.Y2(*ABarray4()[n] ), exe.complexyfunction(*ABarray4()[n] ),\
# [0.0],[0.0])
exe.Y3(*ABarray4()[n] ), exe.complexyfunction3(*ABarray4()[n] ))
result_deltas.append(cpsdif)
result = plt.contourf(exe.A, exe.B, \
np.resize(np.array(result_deltas).flatten() ,\
len(np.array(result_deltas).flatten() ) ).\
reshape(len(exe.B),len(exe.A)), \
cmap = plt.cm.get_cmap('RdBu_r'))#levels = np.arange(-20,-8,2) )
plt.colorbar(result)
plt.title('$\\Delta_{X_s\gamma}$ with charged Higgs: '\
+ str("%02d" % mass_axis1) +', ' + str("%02d"% mass_axis2)+' GeV.' )
plt.xlabel(exe.readlist[int(exe.read1)])
plt.ylabel(exe.readlist[int(exe.read2)])
plt.axis([0, 6.5, -1.6, 0])
plt.savefig('cpsdiffer' + str("%02d" % mass_axis1) +
str("%02d" % mass_axis2) + '.png')
plt.show()
plt.close()
def plt_A_B_untag(i, j): #Untag-asymmetry in {A,B} plane
result_untag = []
mass_axis1, mass_axis2 = i, j
for n in np.arange(0, len(ABarray4())):
untagg = bsg.untag_cp(mb,mw,mass_axis1,mass_axis2,\
exe.Y2(*ABarray4()[n] ), exe.complexyfunction(*ABarray4()[n] ),\
# [0.0],[0.0])
exe.Y3(*ABarray4()[n] ), exe.complexyfunction3(*ABarray4()[n] ))
result_untag.append(untagg)
result = plt.contourf(exe.A, exe.B, \
np.resize(np.array(result_untag).flatten() ,\
len(np.array(result_untag).flatten() ) ).\
reshape(len(exe.B),len(exe.A)), \
cmap = plt.cm.get_cmap('RdBu_r'))#levels = np.arange(-20,-8,2) )
plt.colorbar(result)
plt.title('$A_{CP} (B \\to X_{s + d} \gamma )$ with charged Higgs: '\
+ str("%02d" % mass_axis1) +', ' + str("%02d"% mass_axis2)+' GeV.' )
plt.xlabel(exe.readlist[int(exe.read1)])
plt.ylabel(exe.readlist[int(exe.read2)])
plt.axis([0, 6.5, -1.6, 0])
plt.savefig('untag' + str("%02d" % mass_axis1) + str("%02d" % mass_axis2) +
'.png')
plt.show()
plt.close()
def plt_A_B_cps(i, j): #CP-asymmetry in {A,B} plane
result_cp = []
mass_axis1, mass_axis2 = i, j
for n in np.arange(0, len(ABarray4())):
cpasymetry = bsg.newa_cp(mb,mw,mass_axis1,mass_axis2,\
exe.Y2(*ABarray4()[n] ), exe.complexyfunction(*ABarray4()[n] ),\
# [0.0],[0.0])
exe.Y3(*ABarray4()[n] ), exe.complexyfunction3(*ABarray4()[n] ))
result_cp.append(cpasymetry)
cpresult = plt.contourf(exe.A, exe.B, \
np.resize(np.array(result_cp).flatten() ,\
len(np.array(result_cp).flatten() ) ).\
reshape(len(exe.B),len(exe.A)), \
cmap = plt.cm.get_cmap('RdBu_r') )# levels = np.array([-12,-10,-8,-6,-4,-2,0,2,4]) )
plt.colorbar(cpresult)
plt.title('$A_{CP}(B \\to X_{s}\gamma)$ with charged Higgs: '\
+ str("%02d" % mass_axis1) +', ' + str("%02d"% mass_axis2)+' GeV.' )
plt.xlabel(exe.readlist[int(exe.read1)])
plt.ylabel(exe.readlist[int(exe.read2)])
plt.axis([0, 6.5, -1.6, 0])
plt.savefig('cp' + str("%02d" % mass_axis1) + str("%02d" % mass_axis2) +
'.png')
plt.show()
plt.close()
def plt_A_B_bsg(i, j): # Bsgamma-result in {A,B} plane
mass_axis1, mass_axis2 = i, j
print(ABarray4()[0], mass_axis1, len(ABarray4()))
print(ABarray4()[1], mass_axis2, len(ABarray4()))
resultb = []
#B>Xs+gamma SECTION
for n in np.arange(0, len(ABarray4())):
y3hdm= bsg.BR_B_Xs_gamma(mb,mw,mass_axis1,mass_axis2,\
exe.Y2(*ABarray4()[n] ), exe.complexyfunction(*ABarray4()[n] ),\
exe.Y3(*ABarray4()[n] ), exe.complexyfunction3(*ABarray4()[n] ))
resultb.append(y3hdm / (1e-4))
#########
bsgamm = plt.contourf(exe.A, exe.B, \
np.resize(np.array(resultb).flatten() ,len(np.array(resultb).flatten() ) ).\
reshape(len(exe.B),len(exe.A)) ,\
levels = np.array([2.99,3.55]),colors = ['green'] )
plt.colorbar(bsgamm)
plt.title('BR($\\bar{B} \\to X_{s} \gamma$) in '\
+ str("%02d" % mass_axis1) +', ' + str("%02d"% mass_axis2) )
plt.xlabel(exe.readlist[int(exe.read1)])
plt.ylabel(exe.readlist[int(exe.read2)])
plt.axis([0, 6.5, -1.6, 0])
# plt.axis([1,60,-1.6,0])
# plt.axis([0,60,0,60])
plt.savefig(
str("%02d" % mass_axis1) + str("%02d" % mass_axis2) + 'bsg.png')
plt.show()
plt.close()
def numerical():
mass_axis = (80.0, 250.0)
result = []
for n in np.arange(0, len(ABarray4())):
y3hdm= bsg.BR_B_Xs_gamma(mb,mw,mass_axis[0],mass_axis[1],\
exe.Y2(*ABarray4()[n] ),- exe.complexyfunction(*ABarray4()[n] ),\
exe.Y3(*ABarray4()[n] ),- exe.complexyfunction3(*ABarray4()[n] ))
# print(y3hdm / (1e-4),n)
result.append(y3hdm / (1e-4))
return np.concatenate(result).ravel()
def plt_A_B_nedm(i, j): #[A,B] plane with MHP1 = i, MHp2 = j Nedm
resultn = []
#nedm result
for n in np.arange(0, len(ABarray4())):
nedm3hdm = abs(dn(i,j, exe.complexyfunction(*ABarray4()[n]),\
exe.complexyfunction3(*ABarray4()[n]) ) / (5.06e13) )
resultn.append(nedm3hdm)
#########
ned = plt.contourf(exe.A, exe.B, \
np.resize(np.array(resultn).flatten() ,len(np.array(resultn).flatten() ) ).\
reshape(len(exe.B),len(exe.A)) ,\
levels = np.array([0.0,1.8e-26]) )
# levels = np.arange(3.0,4.2,0.2),\
# colors = ['black','royalblue','purple','darkgreen',\
# 'brown','red','gray','orange','pink'])
plt.colorbar(ned)
plt.title('Neutron EDM with '\
+ str("%02d" % i) +', ' + str("%02d"% j) )
plt.xlabel(exe.readlist[int(exe.read1)])
plt.ylabel(exe.readlist[int(exe.read2)])
return