22double H_cross_section(
double EgammaRyd ,
double EthRyd,
long n,
long l,
long nelem )
25 double rel_photon_energy;
30 rel_photon_energy = EgammaRyd / EthRyd;
31 rel_photon_energy =
MAX2( rel_photon_energy , 1. + FLT_EPSILON*2. );
33 cs =
H_photo_cs(rel_photon_energy , n, l, nelem + 1 );
35 ASSERT( cs > 0. && cs < 1.E-8 );