IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 0af63322 authored by ralfulrich's avatar ralfulrich
Browse files

format

parent 5e465306
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Processes/QGSJetII/testQGSJetII.cc
+ 23
− 23
View file @ 0af63322
......@@ -45,37 +45,37 @@ auto sumMomentum(TStackView const& view, geometry::CoordinateSystem const& vCS)
TEST_CASE("QgsjetII", "[processes]") {
SECTION("QgsjetII -> Corsika") {
CHECK( particles::PiPlus::GetCode() == process::qgsjetII::ConvertFromQgsjetII(
process::qgsjetII::QgsjetIICode::PiPlus));
CHECK(particles::PiPlus::GetCode() == process::qgsjetII::ConvertFromQgsjetII(
process::qgsjetII::QgsjetIICode::PiPlus));
}
SECTION("Corsika -> QgsjetII") {
CHECK( process::qgsjetII::ConvertToQgsjetII(particles::PiMinus::GetCode()) ==
process::qgsjetII::QgsjetIICode::PiMinus);
CHECK( process::qgsjetII::ConvertToQgsjetIIRaw(particles::Proton::GetCode()) == 2);
CHECK(process::qgsjetII::ConvertToQgsjetII(particles::PiMinus::GetCode()) ==
process::qgsjetII::QgsjetIICode::PiMinus);
CHECK(process::qgsjetII::ConvertToQgsjetIIRaw(particles::Proton::GetCode()) == 2);
}
SECTION("canInteractInQgsjetII") {
CHECK( process::qgsjetII::CanInteract(particles::Proton::GetCode()));
CHECK( process::qgsjetII::CanInteract(particles::Code::KPlus));
CHECK( process::qgsjetII::CanInteract(particles::Nucleus::GetCode()));
CHECK(process::qgsjetII::CanInteract(particles::Proton::GetCode()));
CHECK(process::qgsjetII::CanInteract(particles::Code::KPlus));
CHECK(process::qgsjetII::CanInteract(particles::Nucleus::GetCode()));
// CHECK(process::qgsjetII::CanInteract(particles::Helium::GetCode()));
CHECK_FALSE( process::qgsjetII::CanInteract(particles::EtaC::GetCode()));
CHECK_FALSE( process::qgsjetII::CanInteract(particles::SigmaC0::GetCode()));
CHECK_FALSE(process::qgsjetII::CanInteract(particles::EtaC::GetCode()));
CHECK_FALSE(process::qgsjetII::CanInteract(particles::SigmaC0::GetCode()));
}
SECTION("cross-section type") {
CHECK( process::qgsjetII::GetQgsjetIIXSCode(particles::Code::Neutron) ==
process::qgsjetII::QgsjetIIXSClass::Baryons);
CHECK( process::qgsjetII::GetQgsjetIIXSCode(particles::Code::K0Long) ==
process::qgsjetII::QgsjetIIXSClass::Kaons);
CHECK( process::qgsjetII::GetQgsjetIIXSCode(particles::Code::Proton) ==
process::qgsjetII::QgsjetIIXSClass::Baryons);
CHECK( process::qgsjetII::GetQgsjetIIXSCode(particles::Code::PiMinus) ==
process::qgsjetII::QgsjetIIXSClass::LightMesons);
CHECK(process::qgsjetII::GetQgsjetIIXSCode(particles::Code::Neutron) ==
process::qgsjetII::QgsjetIIXSClass::Baryons);
CHECK(process::qgsjetII::GetQgsjetIIXSCode(particles::Code::K0Long) ==
process::qgsjetII::QgsjetIIXSClass::Kaons);
CHECK(process::qgsjetII::GetQgsjetIIXSCode(particles::Code::Proton) ==
process::qgsjetII::QgsjetIIXSClass::Baryons);
CHECK(process::qgsjetII::GetQgsjetIIXSCode(particles::Code::PiMinus) ==
process::qgsjetII::QgsjetIIXSClass::LightMesons);
}
}
......@@ -145,11 +145,11 @@ TEST_CASE("QgsjetIIInterface", "[processes]") {
[[maybe_unused]] const process::EProcessReturn ret = model.DoInteraction(projectile);
[[maybe_unused]] const GrammageType length = model.GetInteractionLength(particle);
CHECK( length / (1_g / square(1_cm)) == Approx(26568.64).margin(0.1) );
CHECK( view.GetSize() == 8 );
CHECK( sumCharge(view) == 2 );
CHECK(length / (1_g / square(1_cm)) == Approx(26568.64).margin(0.1));
CHECK(view.GetSize() == 8);
CHECK(sumCharge(view) == 2);
auto const secMomSum = sumMomentum(view, projectileMomentum.GetCoordinateSystem());
CHECK( (secMomSum - projectileMomentum).norm() / projectileMomentum.norm() ==
Approx(0).margin(1e-2) );
CHECK((secMomSum - projectileMomentum).norm() / projectileMomentum.norm() ==
Approx(0).margin(1e-2));
}
}
Processes/Sibyll/testSibyll.cc
+ 27
− 28
View file @ 0af63322
......@@ -30,41 +30,41 @@ using namespace corsika::units::si;
TEST_CASE("Sibyll", "[processes]") {
SECTION("Sibyll -> Corsika") {
CHECK( particles::Electron::GetCode() ==
process::sibyll::ConvertFromSibyll(process::sibyll::SibyllCode::Electron));
CHECK(particles::Electron::GetCode() ==
process::sibyll::ConvertFromSibyll(process::sibyll::SibyllCode::Electron));
}
SECTION("Corsika -> Sibyll") {
CHECK( process::sibyll::ConvertToSibyll(particles::Electron::GetCode()) ==
process::sibyll::SibyllCode::Electron);
CHECK( process::sibyll::ConvertToSibyllRaw(particles::Proton::GetCode()) == 13);
CHECK( process::sibyll::ConvertToSibyll(particles::XiStarC0::GetCode()) ==
process::sibyll::SibyllCode::XiStarC0);
CHECK(process::sibyll::ConvertToSibyll(particles::Electron::GetCode()) ==
process::sibyll::SibyllCode::Electron);
CHECK(process::sibyll::ConvertToSibyllRaw(particles::Proton::GetCode()) == 13);
CHECK(process::sibyll::ConvertToSibyll(particles::XiStarC0::GetCode()) ==
process::sibyll::SibyllCode::XiStarC0);
}
SECTION("canInteractInSibyll") {
CHECK( process::sibyll::CanInteract(particles::Proton::GetCode()));
CHECK( process::sibyll::CanInteract(particles::Code::XiCPlus));
CHECK(process::sibyll::CanInteract(particles::Proton::GetCode()));
CHECK(process::sibyll::CanInteract(particles::Code::XiCPlus));
CHECK_FALSE( process::sibyll::CanInteract(particles::Electron::GetCode()));
CHECK_FALSE( process::sibyll::CanInteract(particles::SigmaC0::GetCode()));
CHECK_FALSE(process::sibyll::CanInteract(particles::Electron::GetCode()));
CHECK_FALSE(process::sibyll::CanInteract(particles::SigmaC0::GetCode()));
CHECK_FALSE( process::sibyll::CanInteract(particles::Nucleus::GetCode()));
CHECK_FALSE( process::sibyll::CanInteract(particles::Helium::GetCode()));
CHECK_FALSE(process::sibyll::CanInteract(particles::Nucleus::GetCode()));
CHECK_FALSE(process::sibyll::CanInteract(particles::Helium::GetCode()));
}
SECTION("cross-section type") {
CHECK( process::sibyll::GetSibyllXSCode(particles::Code::Electron) == 0);
CHECK( process::sibyll::GetSibyllXSCode(particles::Code::K0Long) == 3);
CHECK( process::sibyll::GetSibyllXSCode(particles::Code::SigmaPlus) == 1);
CHECK( process::sibyll::GetSibyllXSCode(particles::Code::PiMinus) == 2);
CHECK(process::sibyll::GetSibyllXSCode(particles::Code::Electron) == 0);
CHECK(process::sibyll::GetSibyllXSCode(particles::Code::K0Long) == 3);
CHECK(process::sibyll::GetSibyllXSCode(particles::Code::SigmaPlus) == 1);
CHECK(process::sibyll::GetSibyllXSCode(particles::Code::PiMinus) == 2);
}
SECTION("sibyll mass") {
CHECK_FALSE( process::sibyll::GetSibyllMass(particles::Code::Electron) == 0_GeV);
CHECK_FALSE(process::sibyll::GetSibyllMass(particles::Code::Electron) == 0_GeV);
}
}
......@@ -197,15 +197,15 @@ TEST_CASE("SibyllInterface", "[processes]") {
// run checks
// lambda decays into proton and pi- or neutron and pi+
CHECK( stack.GetSize() == 3);
CHECK(stack.GetSize() == 3);
}
SECTION("DecayConfiguration") {
Decay model({particles::Code::PiPlus, particles::Code::PiMinus});
CHECK( model.IsDecayHandled(particles::Code::PiPlus));
CHECK( model.IsDecayHandled(particles::Code::PiMinus));
CHECK_FALSE( model.IsDecayHandled(particles::Code::KPlus));
CHECK(model.IsDecayHandled(particles::Code::PiPlus));
CHECK(model.IsDecayHandled(particles::Code::PiMinus));
CHECK_FALSE(model.IsDecayHandled(particles::Code::KPlus));
const std::vector<particles::Code> particleTestList = {
particles::Code::PiPlus, particles::Code::PiMinus, particles::Code::KPlus,
......@@ -213,16 +213,15 @@ TEST_CASE("SibyllInterface", "[processes]") {
// setup decays
model.SetHandleDecay(particleTestList);
for (auto& pCode : particleTestList)
CHECK( model.IsDecayHandled(pCode));
for (auto& pCode : particleTestList) CHECK(model.IsDecayHandled(pCode));
// individually
model.SetHandleDecay(particles::Code::KMinus);
// possible decays
CHECK_FALSE( model.CanHandleDecay(particles::Code::Proton));
CHECK_FALSE( model.CanHandleDecay(particles::Code::Electron));
CHECK( model.CanHandleDecay(particles::Code::PiPlus));
CHECK( model.CanHandleDecay(particles::Code::MuPlus));
CHECK_FALSE(model.CanHandleDecay(particles::Code::Proton));
CHECK_FALSE(model.CanHandleDecay(particles::Code::Electron));
CHECK(model.CanHandleDecay(particles::Code::PiPlus));
CHECK(model.CanHandleDecay(particles::Code::MuPlus));
}
}
Processes/UrQMD/testUrQMD.cc
+ 3
− 3
View file @ 0af63322
......@@ -159,7 +159,7 @@ TEST_CASE("UrQMD") {
SECTION("nucleus projectile") {
auto [env, csPtr, nodePtr] = setupEnvironment(particles::Code::Oxygen);
[[maybe_unused]] auto const& env_dummy = env; // against warnings
[[maybe_unused]] auto const& env_dummy = env; // against warnings
[[maybe_unused]] auto const& node_dummy = nodePtr; // against warnings
unsigned short constexpr A = 14, Z = 7;
......@@ -183,7 +183,7 @@ TEST_CASE("UrQMD") {
SECTION("\"special\" projectile") {
auto [env, csPtr, nodePtr] = setupEnvironment(particles::Code::Oxygen);
[[maybe_unused]] auto const& env_dummy = env; // against warnings
[[maybe_unused]] auto const& env_dummy = env; // against warnings
[[maybe_unused]] auto const& node_dummy = nodePtr; // against warnings
auto [stackPtr, secViewPtr] =
......@@ -209,7 +209,7 @@ TEST_CASE("UrQMD") {
SECTION("K0Long projectile") {
auto [env, csPtr, nodePtr] = setupEnvironment(particles::Code::Oxygen);
[[maybe_unused]] auto const& env_dummy = env; // against warnings
[[maybe_unused]] auto const& env_dummy = env; // against warnings
[[maybe_unused]] auto const& node_dummy = nodePtr; // against warnings
auto [stackPtr, secViewPtr] =
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

data privacy - impressum