### Forward Elastic Unit Test Check
### -------------------------------------------------------------------------- ##
name = 'Pb-Native'
data_list = cs_list.get(name)
native_file_name = datadir + data_list.get('electroatomic_file_path')
native_data = Native.ElectronPhotonRelaxationDataContainer(native_file_name)
interps = ["LinLinLin", "LinLinLog"]
energies = [1e-4]
for interp in interps:
print "\n----------------------------"
print "--- ", interp, name, "Tests ---"
print "----------------------------"
hybrid_dist = Collision.createHybridElasticDistribution(
native_data, 0.9, interp, True, 1e-14)
###
### Sample
###
print '\n--- Sampling ---'
for energy in energies:
print "\n\t Energy = ", energy
random_numbers = [0.5, 0.9, 0.95, 1.0 - 1e-15]
Prng.RandomNumberGenerator.setFakeStream(random_numbers)
for rnd in random_numbers:
energy, angle = hybrid_dist.sample(energy)
print "\tsample[", energy, ",", '%.16e' % rnd, "] = ", '%.16e' % angle
### -------------------------------------------------------------------------- ##
# -------------------------------------------------------------------------- ##
data_list = cs_list.get('H-Native')
file_name = datadir + data_list.get('electroatomic_file_path')
native_data = Native.ElectronPhotonRelaxationDataContainer(file_name)
energy_grid = native_data.getElectronEnergyGrid()
elastic_energy_grid = native_data.getElasticAngularEnergyGrid()
cs_reduction = native_data.getMomentPreservingCrossSectionReduction()
hybrid_reaction = Collision.createHybridElasticReaction(
native_data, 0.9, False, True, 1e-15)
cutoff_reaction = Collision.createCutoffElasticReaction(
native_data, 1.0, False, True, 1e-15)
analog_reaction = Collision.createAnalogElasticReaction(
native_data, False, True, 1e-15)
hybrid_dist = Collision.createHybridElasticDistribution(
native_data, 0.9, False, True, 1e-15)
cutoff_dist = Collision.createCutoffElasticDistribution(
native_data, 1.0, False, True, 1e-15)
analog_dist = Collision.createAnalogElasticDistribution(
native_data, False, True, 1e-15)
bin_index = 4
energy = elastic_energy_grid[bin_index]
mp_cs_reduction = cs_reduction[bin_index]
cutoff_cs = native_data.getCutoffElasticCrossSection()
moment_cs = native_data.getMomentPreservingCrossSection()
moment_index = native_data.getMomentPreservingCrossSectionThresholdEnergyIndex(
)
index = 0