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: