In [83]:

    

%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
from matplotlib import rc, rcParams
rc('text',usetex=True)
rc('font',**{'family':'serif','serif':['Computer Modern'], 'size' : 18})
cols = ['#29A2C6','#FF6D31','#FFCB18','#73B66B','#EF597B', '#333333']


    

In [112]:

    

import neutrinocommon.astro.body as bd
import neutrinocommon.physconst.physicsconstants as PC
param = PC.PhysicsConstants()

sun = bd.Sun()
Rf = sun.Radius*param.km

ne = []
for i,r in enumerate(edens[:,0]):
    track = sun.track(Rf*r,Rf)
    ye = sun.ye(track)
    ne.append(sun.density(track)*ye/edens[i][1])
    
ne = np.array(ne)


    

In [113]:

    

edens = np.genfromtxt("../output/edens.dat")


    

In [114]:

    

plt.figure(figsize = (8,6))
plt.plot(edens[:,0],edens[:,2], lw = 2, color = "k")
plt.plot(edens[:,0],ne, lw = 2, color = "red")
plt.ylabel(r"$ne$")
plt.xlabel(r"$r")
plt.grid()
#plt.semilogx()


    

In [94]:

    

osc_prob = np.genfromtxt("../output/prob.dat")


    

In [95]:

    

plt.figure(figsize = (8,6))
plt.plot(osc_prob[:,0],osc_prob[:,1], lw = 2, color = "k")
plt.ylabel(r"$P_\odot(\nu_e\to\nu_e)$")
plt.xlabel(r"$E_\nu$ [MeV]")
plt.grid()
plt.semilogx()


    

    
Out[95]:





[]






    

In [87]:

    

stuff_dm = np.genfromtxt("../output/out_dm.dat")
stuff_std = np.genfromtxt("../output/out_std.dat")


    

In [89]:

    

plt.figure(figsize = (8,6))
plt.plot(stuff_dm[:,0],stuff_dm[:,2], lw = 2, color = "k")
plt.plot(stuff_std[:,0],stuff_std[:,2], lw = 2, color = "red")
plt.ylabel(r"$P_\odot(\nu_e\to\nu_e)$")
plt.xlabel(r"$E_\nu$ [MeV]")
plt.grid()
plt.xlim(0.2,14)
plt.semilogx()


    

    
Out[89]:





[]






    

In [91]:

    

plt.figure(figsize = (8,6))
plt.plot(stuff_dm[:,0],stuff_dm[:,1], lw = 2, color = "k")
plt.plot(stuff_std[:,0],stuff_std[:,1], lw = 2, color = "red")
plt.ylabel(r"$Fluxes$")
plt.xlabel(r"$E_\nu$ [MeV]")
plt.ylim(1.0e20,1.0e27)
plt.grid()
plt.loglog()


    

    
Out[91]:





[]






    

In [98]:

    

plt.figure(figsize = (8,6))
ratio = stuff_dm[:,1]/stuff_std[:,1]
plt.plot(stuff_dm[:,0],ratio, lw = 2, color = "k")
plt.ylabel(r"$Ratio$")
plt.xlabel(r"$E_\nu$ [MeV]")
#plt.ylim(1.0e20,1.0e27)
plt.grid()
#plt.loglog()


    

In [39]:

    

b8 = np.genfromtxt("../output/b8.dat")
plt.figure(figsize = (8,6))
plt.plot(b8[:,0],b8[:,1], lw = 2, color = "k")
plt.ylabel(r"$Fluxes$")
plt.xlabel(r"$E_\nu$ [MeV]")
plt.ylim(1.0e20,1.0e27)
plt.grid()
plt.loglog()


    

    
Out[39]:





[]






    

In [ ]: