gwm17/catima
1
0
Fork 0
mirror of https://github.com/gwm17/catima.git synced 2024-11-22 18:28:51 -05:00
Code Issues Projects Releases Wiki Activity

improved angular scattering

Browse Source
This commit is contained in:
hrocho 2020-12-16 23:00:02 +01:00
parent 73d86d925d
commit 732dddbab7
2 changed files with 3 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
calculations.h
View File

@ -100,7 +100,8 @@ namespace catima{
*/ */
double sezi_dedx_e(const Projectile &p, const Material &mat, const Config &c=default_config); double sezi_dedx_e(const Projectile &p, const Material &mat, const Config &c=default_config);
constexpr double Es2_FR =2*PI/fine_structure* electron_mass * electron_mass; //constexpr double Es2_FR =2*PI/fine_structure* electron_mass * electron_mass;
constexpr double Es2_FR = 198.81;
double angular_scattering_power(const Projectile &p, const Target &t, double Es2=Es2_FR); double angular_scattering_power(const Projectile &p, const Target &t, double Es2=Es2_FR);
double angular_scattering_power(const Projectile &p, const Material &material, double Es2=Es2_FR); double angular_scattering_power(const Projectile &p, const Material &material, double Es2=Es2_FR);

6
catima.cpp
View File

@ -57,11 +57,7 @@ double domega2dx(const Projectile &p, const Material &mat, const Config &c){
return sum; return sum;
} }
double da2dx(const Projectile &p, const Material &mat, const Config &c){ double da2dx(const Projectile &p, const Material &mat, const Config &c){
//double Es2 = Es2_FR;
//if(c.scattering == scattering_types::atima_scattering){
// Es2 = 198.81;
//}
const double Es2 = 198.81; const double Es2 = 198.81;
return angular_scattering_power(p,mat, Es2); return angular_scattering_power(p,mat, Es2);
} }