def plot_cevns_rate_fixed_T():
fig2 = plt.figure()
Tmin = 0.001
Z = 32
N = 74.
e_arr = np.linspace(0., 1e7, 10000)
fig2.patch.set_facecolor('white')
plt.ylim((1e-44, 1e-42))
plt.semilogy(e_arr * 1e-6,
dsigmadT_cns(10., e_arr, Z, N),
'k:',
label='T=10 eV',
linewidth=2)
plt.semilogy(e_arr * 1e-6,
dsigmadT_cns(50., e_arr, Z, N),
'b-',
label='T=50 eV',
linewidth=2)
plt.semilogy(e_arr * 1e-6,
dsigmadT_cns(100., e_arr, Z, N),
'r--',
label='T=100 eV',
linewidth=2)
plt.semilogy(e_arr * 1e-6,
dsigmadT_cns(200., e_arr, Z, N),
'g-.',
label='T=200 eV',
linewidth=2)
plt.legend(prop={'size': 11})
plt.xlabel('Neutrino Energy (MeV)')
plt.ylabel('Differential XSec, cm^2/eV')
plt.title('Ge-74 Differential CEvNS XSec, Fixed T')
plt.savefig('plots/diff_xsec_fixed_T.png')
fig2.clf()
def plot_cevns_rate_fixed_Enu():
fig2 = plt.figure()
Tmin = 0.001
Z = 32
N = 74
t_arr = np.logspace(0, 4, 10000)
fig2.patch.set_facecolor('white')
plt.ylim((1e-44, 1e-42))
plt.loglog(t_arr,
dsigmadT_cns(t_arr, 1e6, Z, N),
'k:',
label='Enu = 1 MeV',
linewidth=2)
plt.loglog(t_arr,
dsigmadT_cns(t_arr, 2e6, Z, N),
'b-',
label='Enu = 2 MeV',
linewidth=2)
plt.loglog(t_arr,
dsigmadT_cns(t_arr, 4e6, Z, N),
'r--',
label='Enu = 4 MeV',
linewidth=2)
plt.loglog(t_arr,
dsigmadT_cns(t_arr, 6e6, Z, N),
'g-.',
label='Enu = 6 MeV',
linewidth=2)
plt.legend(prop={'size': 11})
plt.xlabel('Recoil Energy (eV)')
plt.ylabel('Differential XSec, cm^2/eV')
plt.title('Ge-74 Differential CEvNS XSec, Fixed Enu')
plt.savefig('plots/diff_xsec_fixed_Enu.png')
fig2.clf()
def plot_cevns_rate_vs_T_Enu(nu_spec):
fig = plt.figure()
fig.patch.set_facecolor('white')
e_arr = np.linspace(0., 1e7, 1000)
Tmin = 0.001
t_arr = np.logspace(0, 3, 1000)
Z = 32
N = 72.64 - 32.
T, E = np.meshgrid(t_arr, e_arr)
spec = dsigmadT_cns(T, E, Z, N)
smax = spec.max()
smin = smax * 1e-3
spec[spec < smin] = smin
im = plt.pcolor(T,
E * 1e-6,
spec,
norm=LogNorm(vmin=smin, vmax=smax),
cmap='PuBu_r')
fig.colorbar(im)
plt.xlabel("Recoil Energy T (eV)")
plt.ylabel("Neutrino Energy Enu (MeV)")
plt.title("Ge Differential XSec, cm^2/eV")
plt.savefig('plots/diff_xsec_vs_E_T.png')
fig = plt.figure()
fig.patch.set_facecolor('white')
e_arr = np.linspace(0., 1e7, 1000)
Tmin = 0.001
t_arr = np.logspace(0, 3, 1000)
Z = 32
N = 72.64 - 32.
T, E = np.meshgrid(t_arr, e_arr)
spec_flux = dsigmadT_cns(T, E, Z,
N) * nu_spec.d_phi_d_enu_ev(E) * s_per_day * 1e6
smax = spec_flux.max()
smin = smax * 1e-3
spec_flux[spec_flux < smin] = smin
im = plt.pcolor(T,
E * 1e-6,
spec_flux,
norm=LogNorm(vmin=smin, vmax=smax),
cmap=plt.get_cmap('PuBu_r'))
fig.colorbar(im)
plt.xlabel("Recoil Energy T (eV)")
plt.ylabel("Neutrino Energy Enu (MeV)")
plt.title("Ge Differential XSec * Reactor Flux, nu/(eV*MeV*day)")
plt.savefig('plots/diff_xsec_flux_vs_E_T.png')
def plot_cevns_rate_vs_T_Enu(nu_spec=None, nbins=1000):
fig = plt.figure()
fig.patch.set_facecolor('white')
e_arr = np.linspace(0., 1e7, nbins)
Tmin = 0.001
t_arr = np.logspace(0, 3, nbins)
Z = 32
N = 74
T, E = np.meshgrid(t_arr, e_arr)
spec = dsigmadT_cns(T, E, Z, N)
smax = spec.max()
smin = smax * 1e-3
spec[spec < smin] = smin
im = plt.pcolor(T,
E * 1e-6,
spec,
norm=LogNorm(vmin=smin, vmax=smax),
cmap='PuBu_r')
fig.colorbar(im)
plt.xlabel("Recoil Energy T (eV)")
plt.ylabel("Neutrino Energy Enu (MeV)")
plt.title("Ge-74 Differential XSec, cm^2/eV")
plt.savefig('plots/diff_xsec_vs_E_T.png')