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:
index = i
for z in elements:
print "\n----------------------------"
print "-----", z, "Tests -----"
print "----------------------------"
data_list = cs_list.get(z)
file_name = datadir + data_list.get('electroatomic_file_path')
native_data = Native.ElectronPhotonRelaxationDataContainer(file_name)
energy_grid = native_data.getElectronEnergyGrid()
tot_elastic_cs = native_data.getTotalElasticCrossSection()
cutoff_cross_sections = native_data.getCutoffElasticCrossSection()
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":
