(个人学习笔记,仅供参考)
t1 = np.arange(0.01,0.5,0.01) #从0.01开始步长为0.01,直到0.5
energy = 5 #给定能量值
DCS_pipi22 = [DCS_pipi(energy**2,i,lmd_pi_pf,lmd_pi_rf2) for i in t1]
DCS_pipi11 = [DCS_pipi(energy**2,i,lmd_pi_pf,lmd_pi_rf1) for i in t1]
DCS_pipi21 = [DCS_pipi(energy**2,i,lmd_pi_pf,math.sqrt(lmd_pi_rf2*lmd_pi_rf1)) for i in t1]
DCS_pip1 = [DCS_pip(energy**2,i,lmd_pi_pf,lmd_pi_rf1) for i in t1]
DCS_pip2 = [DCS_pip(energy**2,i,lmd_pi_pf,lmd_pi_rf2) for i in t1]
print(t1)
print(DCS_pipi22)
print(DCS_pipi11)
print(DCS_pipi21)
print(DCS_pip1)
print(DCS_pip2)
plt.plot(t1,DCS_pip1,c='red',label='$\pi^{+} p$',linestyle=(0,(3,3))) #(x轴上的数据, 曲线颜色, 标签, 设置线型"solid=虚线", 虚线的样式)
plt.plot(t1,DCS_pip2,c='blue',label='$\pi^{-} p$',linestyle=(0,(3,0)))
plt.plot(t1,DCS_pipi11,c='r',label='$\pi^{+}\pi^{+}$',linestyle=(0,(6,3,1,3)))
plt.plot(t1,DCS_pipi22,c='g',label='$\pi^{-}\pi^{-}$',linestyle=(0,(5,5)))
plt.plot(t1,DCS_pipi21,c='b',label='$\pi^{+}\pi^{-}$',linestyle=(0,(6,3,.5,3,.5,3)))
plt.yscale('log') #将y轴的刻度设置为对数刻度,使y轴上的数值按对数尺度展示
plt.legend(loc='best',numpoints=1) #设置图例的位置参数, 每个图例条目只显示一个标记点
plt.xlabel(r'$|t|\rm (GeV)$',fontsize=20) #x轴标签, 即字体大小
plt.ylabel(r'$\frac{d \sigma}{dt} \rm (mb/GeV^2)$',fontsize=20) #y轴标签, 即字体大小
plt.text(0.2,40,r'$\sqrt{s} =5 ~\rm GeV$',size='xx-large') #text函数在图形中添加文本,及文本位置
plt.savefig(r"E:\Code文件\Code note figure\predication_5GeV.pdf",bbox_inches='tight',dpi=600) #图片保存位置
plt.show()