| GALPROP https://galprop.stanford.edu/forum/ |
|
| some problems about energy loss rate of electron https://galprop.stanford.edu/forum/viewtopic.php?f=20&t=135 |
Page 1 of 1 |
| Author: | jhbcn [ Thu Nov 17, 2011 1:25 am ] |
| Post subject: | some problems about energy loss rate of electron |
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. |
|
| Author: | jhbcn [ Sun Nov 27, 2011 6:22 pm ] |
| Post subject: | Re: some problems about energy loss rate of electron |
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 |
|
| Page 1 of 1 | All times are UTC - 8 hours [ DST ] |
| Powered by phpBB® Forum Software © phpBB Group https://www.phpbb.com/ |
|