100 G4double* pMicroRoughnessTransTable,
101 G4double* pmaxMicroRoughnessTransTable)
129 if (Nthetadim*NEdim > 0) {
190 std::ofstream dateir(
"MRrefl.dat",std::ios::out);
191 std::ofstream dateit(
"MRtrans.dat",std::ios::out);
200 IntIplus(E, fermipot, theta_i, AngNoTheta, AngNoPhi,
205 IntIminus(E, fermipot, theta_i, AngNoTheta, AngNoPhi,
223 std::ofstream dateic(
"MRcheck.dat",std::ios::out);
224 std::ofstream dateimr(
"MRmaxrefl.dat",std::ios::out);
225 std::ofstream dateimt(
"MRmaxtrans.dat",std::ios::out);
262 if (theta_i<theta_i_min || theta_i>
theta_i_max || Energy<Emin || Energy>
Emax)
288 if (theta_i<theta_i_min || theta_i>
theta_i_max || Energy<Emin || Energy>
Emax)
311 if (theta_i<theta_i_min || theta_i>
theta_i_max || Energy<Emin || Energy>
Emax)
332 Energy<Emin || Energy>
Emax) {
357 if (theta_i<theta_i_min || theta_i>
theta_i_max || Energy<Emin || Energy>
Emax)
378 Energy<Emin || Energy>
Emax) {
401 ProbIplus(Energy, fermipot, theta_i, theta_o, phi_o,
b,
w,
AngCut);
410 ProbIminus(Energy, fermipot,theta_i, theta_o, phi_o,
b,
w,
AngCut);
429 if (2*
b*k*std::cos(theta_i) < 1 && 2*
b*k_l < 1)
return true;
440 if (E*(std::cos(theta_i)*std::cos(theta_i)) < VFermi)
return false;
449 if (2*
b*std::sqrt(kS2)*std::cos(theta_i) < 1 && 2*
b*std::sqrt(k_l2) < 1)
return true;
G4double GetMRIntProbability(G4double, G4double)
void SetMicroRoughnessParameters(G4double, G4double, G4int, G4int, G4double, G4double, G4double, G4double, G4int, G4int, G4double)
G4double GetMRMaxTransProbability(G4double, G4double)
G4double GetConstProperty(const char *key) const
void InitMicroRoughnessTables()
static constexpr double nm
G4double * maxMicroRoughnessTable
void AddConstProperty(const char *key, G4double PropertyValue)
G4double GetMRIntTransProbability(G4double, G4double)
G4double GetMRProbability(G4double, G4double, G4double, G4double, G4double)
G4double GetMRTransProbability(G4double, G4double, G4double, G4double, G4double)
G4double * theMicroRoughnessTable
static constexpr double neutron_mass_c2
G4bool ConditionsValid(G4double E, G4double VFermi, G4double theta_i)
G4bool ConstPropertyExists(const char *key) const
static G4UCNMicroRoughnessHelper * GetInstance()
const XML_Char int const XML_Char * value
G4double * GetMicroRoughnessTransTable()
void SetMRMaxProbability(G4double, G4double, G4double)
G4UCNMaterialPropertiesTable()
static constexpr double eV
G4double GetMRMaxProbability(G4double, G4double)
G4double * theMicroRoughnessTransTable
G4double * GetMicroRoughnessTable()
void LoadMicroRoughnessTables(G4double *, G4double *, G4double *, G4double *)
virtual ~G4UCNMaterialPropertiesTable()
static constexpr double degree
G4GLOB_DLL std::ostream G4cout
void ComputeMicroRoughnessTables()
void RemoveConstProperty(const char *key)
G4double * maxMicroRoughnessTransTable
G4bool TransConditionsValid(G4double E, G4double VFermi, G4double theta_i)
static constexpr double hbarc_squared
void SetMRMaxTransProbability(G4double, G4double, G4double)