IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 3c603763 authored by Felix Riehn's avatar Felix Riehn
Browse files

remove test for proton target, this seems to be broken in pythia/angantyr

parent db64dd18
No related branches found
No related tags found
1 merge request!549Draft: "Use Pythia for interactions"
...@@ -61,8 +61,8 @@ SECTION("pythia interaction: proton") { ...@@ -61,8 +61,8 @@ SECTION("pythia interaction: proton") {
SECTION("pythia interaction: PiPlus") { SECTION("pythia interaction: PiPlus") {
// test some combinations of valid target and projectile particles // test some combinations of valid target and projectile particles
// so far only hadron-hadron and hadron-Nucleus is allowed //
Code const target = GENERATE(Code::Proton, Code::Nitrogen, Code::Oxygen, Code::Argon); Code const target = GENERATE(Code::Nitrogen, Code::Oxygen, Code::Argon);
Code const projectile = Code::PiPlus; Code const projectile = Code::PiPlus;
CORSIKA_LOG_INFO("testing: {}-{}", projectile, target); CORSIKA_LOG_INFO("testing: {}-{}", projectile, target);
...@@ -82,7 +82,7 @@ SECTION("pythia interaction: PiPlus") { ...@@ -82,7 +82,7 @@ SECTION("pythia interaction: PiPlus") {
SECTION("pythia interaction: KPlus") { SECTION("pythia interaction: KPlus") {
// test some combinations of valid target and projectile particles // test some combinations of valid target and projectile particles
// so far only hadron-hadron and hadron-Nucleus is allowed // so far only hadron-hadron and hadron-Nucleus is allowed
Code const target = GENERATE(Code::Proton, Code::Nitrogen, Code::Oxygen, Code::Argon); Code const target = GENERATE(Code::Nitrogen, Code::Oxygen, Code::Argon);
Code const projectile = Code::KPlus; Code const projectile = Code::KPlus;
CORSIKA_LOG_INFO("testing: {}-{}", projectile, target); CORSIKA_LOG_INFO("testing: {}-{}", projectile, target);
...@@ -99,6 +99,26 @@ SECTION("pythia interaction: KPlus") { ...@@ -99,6 +99,26 @@ SECTION("pythia interaction: KPlus") {
REQUIRE(view.getSize() >= 2); REQUIRE(view.getSize() >= 2);
} }
SECTION("pythia interaction: angantyr fail") {
// test some combinations of valid target and projectile particles
// so far only hadron-hadron and hadron-Nucleus is allowed
Code const target = Code::Proton;
Code const projectile = GENERATE(Code::KPlus, Code::PiPlus);
CORSIKA_LOG_INFO("testing: {}-{}", projectile, target);
REQUIRE(
collision.getCrossSection(
projectile, target,
{calculate_total_energy(P0, get_mass(projectile)), {rootCS, {0_eV, 0_eV, P0}}},
{get_mass(target), {rootCS, {0_eV, 0_eV, 0_eV}}}) > 0_mb);
REQUIRE_THROWS(
collision.doInteraction(view, projectile, target,
{corsika::calculate_total_energy(P0, get_mass(projectile)),
{rootCS, {0_eV, 0_eV, P0}}},
{get_mass(target), {rootCS, {0_eV, 0_eV, 0_eV}}}));
}
SECTION("pythia too low energy") { SECTION("pythia too low energy") {
// this is a projectile neutron with very little energy // this is a projectile neutron with very little energy
...@@ -109,10 +129,11 @@ SECTION("pythia too low energy") { ...@@ -109,10 +129,11 @@ SECTION("pythia too low energy") {
corsika::pythia8::Interaction collision; corsika::pythia8::Interaction collision;
// 5 MeV lab is too low, 0 mb expected // 5 MeV lab is too low, 0 mb expected
REQUIRE(collision.getCrossSectionInelEla( REQUIRE_THROWS(
Code::Proton, Code::Proton, collision.getCrossSectionInelEla(
{calculate_total_energy(Proton::mass, 5_MeV), {rootCS, 0_eV, 0_eV, 5_MeV}}, Code::Proton, Code::Proton,
{Proton::mass, {rootCS, 0_eV, 0_eV, 0_eV}}) == std::tuple{0_mb, 0_mb}); {calculate_total_energy(Proton::mass, 5_MeV), {rootCS, 0_eV, 0_eV, 5_MeV}},
{Proton::mass, {rootCS, 0_eV, 0_eV, 0_eV}}) == std::tuple{0_mb, 0_mb});
REQUIRE_THROWS(collision.doInteraction( REQUIRE_THROWS(collision.doInteraction(
view, Code::Neutron, Code::Proton, view, Code::Neutron, Code::Proton,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment