Exemple #1
0
 def Brems_spec(self,EdNdE_e,Ee,Eg):
     '''
     Calculates the bremsstrahlungs spectrum for given electron spectrum and ambient density 
     '''
     dlogEe=np.log(Ee[1]/Ee[0])*np.ones_like(Ee)
     dNe=EdNdE_e*dlogEe
     sigma=crs.sigma_brems(Eg,Ee,dNe,self.source.Zeff)
     dNdE_gamma=self.source.nHeff*const.C*const.m2cm*sigma
     return dNdE_gamma
Exemple #2
0
def plot_cross(typ=['ic','brems','pp']):
    Eg=np.logspace(6.,12.,300)
    Ee=np.logspace(7.,15.,100)
    Ne=spsh.plcut(Ee,[1.,9.,-3.,14.])
    Ne2=spsh.plcut(Ee,[1.,9.,-2.0,11.])
    ylist=[]
    xlist=[]
    taglist=[]
    elist=[]

    if 'brems' in typ:
        Ee/=10**9
        Eg/=10**9
        ylist.append(const.C*const.m2cm*crs.sigma_brems(Eg,Ee,Ne))
        print const.C*const.m2cm*crs.sigma_brems(Eg,Ee,Ne)
        #ylist.append(crs.sigma_brems(Eg,Ee,Ne2))
        ylist.append(const.C*const.m2cm*crs.sigma_brems_ee(Eg,Ee,Ne))
        #ylist.append(crs.sigma_brems_ee(Eg,Ee,Ne2))
        ylist.append(const.C*const.m2cm*gs.sigmaB_test(Ne,Eg,Ee))
        xlist.append(Eg)
        xlist.append(Eg)
        xlist.append(Eg)
        #xlist.append(Eg)
        #xlist.append(Eg)
        for t in['brems_N -1.5','brems_ee -1.5','brems dima']: 
            taglist.append(t)
            elist.append(0)
    if 'ic' in typ:
        Ee/=10**9
        Eg/=10**9
        #print Ee,Eg
        source=sed.source()
        source.set_radiation_densities(0.25,0.84,3.0e-3,0.9)#W51C densities Abdo et al 2009
        source.set_radiation_photons(-8.,4.,0.1)
        '''
        print 'CMB: ',source.phCMB
        print 'Dust : ',source.phDust
        print 'Star: ',source.phStar
        print 'tot: ',source.phTot
        '''
        for ph in [source.phTot,source.phCMB,source.phDust,source.phStar]:
            def sigma_fcn(x):
                sigma=crs.sigma_ic(x,source.Eph,Ee)
                #print 'sigma: ',np.size(sigma)
                return const.C*const.m2cm*np.dot(ph,np.dot(sigma,np.ones_like(Ee)))

            sigma=np.frompyfunc(sigma_fcn,1,1)
            ylist.append(sigma(Eg))
            xlist.append(Eg)
            elist.append(0)
        for t in['IC total','IC CMB','IC Dust','IC Star']: 
            taglist.append(t)
        plot([source.Eph,source.Eph,source.Eph,source.Eph],[source.phTot,source.phCMB,source.phDust,source.phStar],[],elist,['Tot','CMB','Dust','Star'],r'E [eV]',r'$N_{ph}$')
    
    #print ylist
    
    plot(xlist,ylist,[],elist,taglist,r'E [eV]',r'$N_{e}\times \sigma_{\mathrm{brems}}$')

    #plot([source.Eph,source.Eph,source.Eph,source.Eph],[source.phTot,source.phCMB,source.phDust,source.phStar],[],elist,['Tot','CMB','Dust','Star'],r'E [eV]',r'$N_{ph}$')
    
    plt.show()