Python scipy - 指定自定义离散分布

我使用scipy.stats的各种连续分布(例如norm).所以如果我想找到P(Z <0.5),我会这样做:

from scipy.stats import norm
norm(0, 1).cdf(0.5)  # Z~N(0,1)

是否有工具(scipy.stats或statsmodels或其他)可用于描述离散分布,然后计算CDF/CMF等？我可以自己编写代码,但我想知道是否存在某些内容,例如:

pdf(x)= 1/3,x = 1,2,3; 否则0

然后我可以构造2个向量x = [1,2,3],p = [1/3,1/3,1/3]并将它们输入到一个库类中,然后提供.cdf()等？