<SYSTEM_TASK:>
Given the following text description, write Python code to implement the functionality described below step by step
<END_TASK>
<USER_TASK:>
Description:
 Setup a simulated data to see if we can make it look right
Step1:  This will then serve as the background
Can the simplest model get started?
<ASSISTANT_TASK:>
Python Code:
np.random.seed(8675309)
sim_pa = np.arange(20,175)
sim_c = 890*np.sin(np.deg2rad(sim_pa))**0.8
# at each point draw a poisson variable with that mean
sim_c_n = np.asarray([np.random.poisson(v) for v in sim_c ])
prob=0.1
sim_c_n2 = np.asarray([np.random.negative_binomial((v*prob)/(1-prob), prob) for v in sim_c ])
# Two subplots, the axes array is 1-d
f, axarr = plt.subplots(3, sharex=True, sharey=True, figsize=(7,9))
axarr[0].plot(sim_pa, sim_c, lw=2)
axarr[0].set_ylabel('Truth')
axarr[0].set_xlim((0,180))
axarr[0].set_yscale('log')
axarr[1].plot(sim_pa, sim_c_n, lw=2)
axarr[1].set_ylabel('Measured\nPoisson')
axarr[1].set_xlim((0,180))
axarr[1].set_yscale('log')
axarr[2].plot(sim_pa, sim_c_n2, lw=2)
axarr[2].set_ylabel('Measured\nNegBin')
axarr[2].set_xlim((0,180))
axarr[2].set_yscale('log')
axarr[2].set_ylim(axarr[1].get_ylim())
Explanation: Setup a simulated data to see if we can make it look right
End of explanation
# generate some data
with pm.Model() as model:
    bkg = pm.NegativeBinomial('bkg', mu=pm.Uniform('m_bkg', 0, 1e5,shape=len(sim_pa), testval=1e3), 
                              alpha=0.1, 
                              observed=sim_c_n2, shape=len(sim_pa))
    
    
    
#     truth_mc = pm.Uniform('truth', 0, 100, shape=dat_len)
#     noisemean_mc = pm.Uniform('noisemean', 0, 100)
#     noise_mc = pm.Poisson('noise', noisemean_mc, observed=obs[1:20])
#     real_n_mc = pm.Poisson('real_n', truth_mc+noisemean_mc, shape=dat_len)
#     psf = pm.Uniform('psf', 0, 5, observed=det)
#     obs_mc = pm.Normal('obs', (truth_mc+noisemean_mc)*psf.max(), 1/5**2, observed=obs, shape=dat_len)
    
    trace = pm.sample(50000)
pm.traceplot(trace)
pm.summary(trace)
Explanation: This will then serve as the background
Can the simplest model get started?
End of explanation
<END_TASK>