From b4727f6ff21b2c320df1ed2bfbdeb69618f9b676 Mon Sep 17 00:00:00 2001 From: ralfulrich <ralf.ulrich@kit.edu> Date: Thu, 16 Jul 2020 22:00:01 +0200 Subject: [PATCH] format --- Environment/testShowerAxis.cc | 34 +++++++++++++++------------------- 1 file changed, 15 insertions(+), 19 deletions(-) diff --git a/Environment/testShowerAxis.cc b/Environment/testShowerAxis.cc index e86bb6754..b541e63f0 100644 --- a/Environment/testShowerAxis.cc +++ b/Environment/testShowerAxis.cc @@ -8,11 +8,11 @@ * the license. */ -#include <corsika/environment/ShowerAxis.h> #include <corsika/environment/DensityFunction.h> #include <corsika/environment/HomogeneousMedium.h> #include <corsika/environment/IMediumModel.h> #include <corsika/environment/NuclearComposition.h> +#include <corsika/environment/ShowerAxis.h> #include <corsika/environment/VolumeTreeNode.h> #include <corsika/geometry/Line.h> #include <corsika/geometry/RootCoordinateSystem.h> @@ -44,9 +44,8 @@ auto setupEnvironment(particles::Code vTargetCode) { using MyHomogeneousModel = environment::HomogeneousMedium<environment::IMediumModel>; theMedium->SetModelProperties<MyHomogeneousModel>( - density, - environment::NuclearComposition(std::vector<particles::Code>{vTargetCode}, - std::vector<float>{1.})); + density, environment::NuclearComposition(std::vector<particles::Code>{vTargetCode}, + std::vector<float>{1.})); auto const* nodePtr = theMedium.get(); universe.AddChild(std::move(theMedium)); @@ -54,7 +53,6 @@ auto setupEnvironment(particles::Code vTargetCode) { return std::make_tuple(std::move(env), &cs, nodePtr); } - TEST_CASE("Homogeneous Density") { auto [env, csPtr, nodePtr] = setupEnvironment(particles::Code::Nitrogen); auto const& cs = *csPtr; @@ -65,26 +63,24 @@ TEST_CASE("Homogeneous Density") { auto const injectionHeight = 10_km; auto const t = -observationHeight + injectionHeight; Point const showerCore{cs, 0_m, 0_m, observationHeight}; - Point const injectionPos = - showerCore + - Vector<dimensionless_d>{cs, {0, 0, 1}} * t; - - environment::ShowerAxis const showerAxis{injectionPos, - (showerCore - injectionPos), *env, 10000}; + Point const injectionPos = showerCore + Vector<dimensionless_d>{cs, {0, 0, 1}} * t; + + environment::ShowerAxis const showerAxis{injectionPos, (showerCore - injectionPos), + *env, 10000}; + + CHECK(showerAxis.steplength() == 1_m); - CHECK( showerAxis.steplength() == 1_m ); - - CHECK( showerAxis.maximumX()/(10_km * density) == Approx(1).epsilon(1e-8) ); + CHECK(showerAxis.maximumX() / (10_km * density) == Approx(1).epsilon(1e-8)); - CHECK( showerAxis.minimumX() == 0_g / square(1_cm) ); + CHECK(showerAxis.minimumX() == 0_g / square(1_cm)); const Point p{cs, 10_km, 20_km, 9_km}; - CHECK( showerAxis.projectedX(p)/(1_km * density) == Approx(1).epsilon(1e-8) ); + CHECK(showerAxis.projectedX(p) / (1_km * density) == Approx(1).epsilon(1e-8)); const units::si::LengthType d = 1_km; - CHECK( showerAxis.X(d)/(d * density) == Approx(1).epsilon(1e-8) ); + CHECK(showerAxis.X(d) / (d * density) == Approx(1).epsilon(1e-8)); const Vector<dimensionless_d> dir{cs, {0, 0, -1}}; - CHECK( showerAxis.GetDirection().GetComponents(cs) == dir.GetComponents(cs) ); - CHECK( showerAxis.GetStart().GetCoordinates() == injectionPos.GetCoordinates() ); + CHECK(showerAxis.GetDirection().GetComponents(cs) == dir.GetComponents(cs)); + CHECK(showerAxis.GetStart().GetCoordinates() == injectionPos.GetCoordinates()); } -- GitLab