print "\n----- H -----\n"
tot_adjoint_elastic_cs = adjoint_data.getAdjointTotalElasticCrossSection()
reduced_cutoff_ratio = adjoint_data.getReducedCutoffCrossSectionRatios()
adjoint_cutoff_cs = adjoint_data.getAdjointCutoffElasticCrossSection()
adjoint_rutherford_cs = adjoint_data.getAdjointScreenedRutherfordElasticCrossSection(
)
adjoint_sr_index = adjoint_data.getAdjointScreenedRutherfordElasticCrossSectionThresholdEnergyIndex(
)
angular_energy_grid = adjoint_data.getAdjointElasticAngularEnergyGrid()
moment_cs_reduction = adjoint_data.getAdjointMomentPreservingCrossSectionReduction(
)
subshells = adjoint_data.getSubshells()
shell = subshells[0]
adjoint_ionization_cs = adjoint_data.getAdjointElectroionizationCrossSection(
shell)
adjoint_cutoff_dist = Collision.createLogLogLogCorrelatedCutoffElasticDistribution(
adjoint_data, 0.9, 1e-7)
adjoint_mp_reaction = Collision.createLogLogLogCorrelatedMomentPreservingElasticReaction(
adjoint_data, 0.9, 1e-7)
adjoint_hybrid_reaction = Collision.createLogLogLogCorrelatedHybridElasticReaction(
adjoint_data, 0.9, 1e-7)
adjoint_hybrid_dist = Collision.createLogLogLogCorrelatedHybridElasticDistribution(
adjoint_data, 0.9, 1e-7)
adjoint_cutoff_reaction = Collision.createLogLogLogCorrelatedCutoffElasticReaction(
adjoint_data, 0.9, 1e-7)
energies = [1e-5, 1e-3, 20.0]
#energies = [1.0, 10.0, 20.0]
for energy in energies:
index = 0
for i in range(0, adjoint_energy_grid.size):
if adjoint_energy_grid[i] <= energy:
screen_rutherford_cs = native_data.getScreenedRutherfordElasticCrossSection(
)
screen_rutherford_index = native_data.getScreenedRutherfordElasticCrossSectionThresholdEnergyIndex(
)
moment_cross_sections = Collision.createLogLogLogCorrelatedMomentPreservingElasticReaction(
native_data, 0.9, 1e-7)
# if z == 'Pb-Native':
# energies = [1e-5,4e-4,1e5]
# else:
# energies = [1e-5, 1e-3, 1e5 ]
for interp in interps:
print "\n--- ", interp, "Tests ---"
cutoff_dist = Collision.createLogLogLogCorrelatedCutoffElasticDistribution(
native_data, 0.9, 1e-7)
if interp == "LinLinLog":
cutoff_dist = Collision.createLinLinLogCorrelatedCutoffElasticDistribution(
native_data, 0.9, 1e-15)
elif interp == "LinLinLin":
cutoff_dist = Collision.createLinLinLinCorrelatedCutoffElasticDistribution(
native_data, 0.9, 1e-15)
###
### Moment Preserving Reaction Unit Test Check
###
for energy in energies:
index = 0
for i in range(0, energy_grid.size):
if energy_grid[i] <= energy: