It is currently Sat Dec 04, 2021 6:42 am

All times are UTC - 8 hours [ DST ]




Post new topic Reply to topic  [ 2 posts ] 
Author Message
PostPosted: Thu Nov 17, 2011 1:25 am 
Offline

Joined: Tue Mar 01, 2011 6:58 pm
Posts: 9
Dear Prof. Strong,

Here ,I have some problems about energy loss rate of electron :

I have read your paper:
On possible interpretations of the high energy electron–positron spectrum
measured by the Fermi Large Area Telescope at Astroparticle Physics 32 (2009) 140–151

At the appendix A,I want to know if or not the parameter b(E) (energy loss reate ) is changed at the GALDEF file.I not find which is relative to b(E) at the GALDEF file.If or not I may modify the GALPROP code to do it.

If the equation (A-1) may be resolved by running the GALPROP or not ? The GALPROP 's propagation equation is different to A-1 and if I may modify the GALDEF parameters to fit the equation (A-1) .

I want to know how to change the b(E) by modifying some parameters or GALPROP code.


Top
 Profile  
 
PostPosted: Sun Nov 27, 2011 6:22 pm 
Offline

Joined: Tue Mar 01, 2011 6:58 pm
Posts: 9
I give the following code to resolve the problems.

if(oneParticle.n_spatial_dimensions==2)
{for(int ir=0; ir<oneParticle.n_rgrid; ir++)
{for(int iz=0; iz<oneParticle.n_zgrid; iz++)
{for(int ip=0; ip<oneParticle.n_pgrid; ip++)
{if(oneParticle.A==0)
{
oneParticle.dpdt.d2[ir] [iz].s[ip]=oneParticleEnergyLoss(oneParticle.Etot[ip]);
oneParticle.Dxx.d2[ir][iz].s[ip] = oneParticle.beta[ip]*oneParticleDXX(oneParticle.Etot[ip]);
}
} // p
} // z
} // r
}

oneParticleEnergyLoss and oneParticleDXX give the expression of some model.

for electron Etot=rigidity


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 2 posts ] 

All times are UTC - 8 hours [ DST ]


Who is online

Users browsing this forum: Google [Bot] and 1 guest


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
cron
Powered by phpBB® Forum Software © phpBB Group