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

Merge pull request #51 from hrosiak/kox

Browse Source 
kox formula CC change
This commit is contained in:
Andrej Prochazka 2019-01-17 18:22:33 +01:00 committed by GitHub
parent 273520a2e9 6c051b6f84
commit 84c2ef84e4
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
reactions.cpp
View File

@ -89,7 +89,7 @@ double SigmaR_Kox(int Ap, int Zp, double E, int At, int Zt){
double Rs = r0 * ((a*Ap13*At13)/(Ap13+At13)-c)+D; double Rs = r0 * ((a*Ap13*At13)/(Ap13+At13)-c)+D;
double Rv = r0 * (Ap13 + At13); double Rv = r0 * (Ap13 + At13);
double Ri = Rv + Rs; double Ri = Rv + Rs;
double Ecm = Ecm_from_T_relativistic(E,Ap,At); double Ecm = Ecm_from_T(E,Ap,At);
return 10.0*PI*Ri*Ri*(1-(Bc/Ecm)); return 10.0*PI*Ri*Ri*(1-(Bc/Ecm));
} }
#endif #endif

6
reactions.h
View File

@ -58,6 +58,7 @@ double SigmaR_Kox(int Ap, int Zp, double E, int At, int Zt);
inline double p_from_T(double T, double M=1.0){ inline double p_from_T(double T, double M=1.0){
return M*sqrt(T*T + 2*T*atomic_mass_unit); return M*sqrt(T*T + 2*T*atomic_mass_unit);
} }
/*
inline double Ecm_from_T_relativistic(double T, double Ap, double At){ inline double Ecm_from_T_relativistic(double T, double Ap, double At){
double mp = Ap*atomic_mass_unit; double mp = Ap*atomic_mass_unit;
double mt = At*atomic_mass_unit; double mt = At*atomic_mass_unit;
@ -67,6 +68,11 @@ inline double Ecm_from_T_relativistic(double T, double Ap, double At){
double pcm = plab * mt / ecm; double pcm = plab * mt / ecm;
return sqrt(pcm*pcm+mp*mp)-mp; return sqrt(pcm*pcm+mp*mp)-mp;
} }
*/
/// calculates Ecm fom T, return in MeV units
inline double Ecm_from_T(double T, double Ap, double At){
return T*Ap*At/(Ap+At);
}
#endif //NUREX #endif //NUREX
#endif #endif