File indexing completed on 2024-05-10 02:20:35
0001 #include <iomanip>
0002
0003 #include <cppunit/extensions/HelperMacros.h>
0004
0005 #include "DataFormats/Math/interface/GeantUnits.h"
0006 #include <CLHEP/Units/SystemOfUnits.h>
0007
0008 #include "cppunit/TestAssert.h"
0009 #include "cppunit/TestFixture.h"
0010
0011 using namespace geant_units;
0012 using namespace geant_units::operators;
0013 using namespace std;
0014
0015 class testDDUnits : public CppUnit::TestFixture {
0016 CPPUNIT_TEST_SUITE(testDDUnits);
0017 CPPUNIT_TEST(checkUnits);
0018 CPPUNIT_TEST_SUITE_END();
0019
0020 public:
0021 void setUp() override {}
0022 void tearDown() override {}
0023 void checkUnits();
0024 };
0025
0026 CPPUNIT_TEST_SUITE_REGISTRATION(testDDUnits);
0027
0028 void testDDUnits::checkUnits() {
0029 std::cout << "\nMy pi: " << std::setprecision(16) << piRadians << " == " << 1_pi << " == " << 1._pi << "\n";
0030
0031 CPPUNIT_ASSERT(M_PI == piRadians);
0032 CPPUNIT_ASSERT(M_PI == 1_pi);
0033 CPPUNIT_ASSERT(M_PI == 1._pi);
0034
0035 double twoPiAngle = 2_pi;
0036 std::cout << "My 2pi angle: " << twoPiAngle << " [rad] == " << convertRadToDeg(twoPiAngle) << " [deg]\n";
0037
0038 CPPUNIT_ASSERT(2 * M_PI == 2 * piRadians);
0039 CPPUNIT_ASSERT(2 * M_PI == twoPiAngle);
0040 CPPUNIT_ASSERT(2 * M_PI == 2_pi);
0041 CPPUNIT_ASSERT(2 * M_PI == 2._pi);
0042
0043 CPPUNIT_ASSERT(90_deg == 1_pi / 2);
0044 CPPUNIT_ASSERT(120_deg == 2_pi / 3);
0045 CPPUNIT_ASSERT(135_deg == 3_pi / 4);
0046 CPPUNIT_ASSERT(150_deg == 5_pi / 6);
0047
0048 double angle90 = convertRadToDeg(1_pi / 2);
0049 CPPUNIT_ASSERT(angle90 == 90);
0050
0051 double angle120 = convertRadToDeg(2_pi / 3);
0052 CPPUNIT_ASSERT(almostEqual(angle120, 120., 2));
0053
0054 double angle135 = convertRadToDeg(3_pi / 4);
0055 CPPUNIT_ASSERT(almostEqual(angle135, 135., 2));
0056
0057 double angle150 = convertRadToDeg(5_pi / 6);
0058 CPPUNIT_ASSERT(almostEqual(angle150, 150., 2));
0059
0060 cout << "Mass of 1 kg is " << 1._kg << " or " << 1 * CLHEP::kg << "\n";
0061 cout << "Mass of 1 g is " << 1._g << " or " << 1 * CLHEP::g << "\n";
0062 cout << "Ratio of 1._kg / CLHEP::kg is " << 1._kg / (1 * CLHEP::kg) << endl;
0063 cout << "Difference of 1._kg - CLHEP::kg is " << std::abs(1._kg - (1 * CLHEP::kg)) << endl;
0064 cout << "Ratio of 1._g / CLHEP::g is " << 1._g / (1 * CLHEP::g) << endl;
0065 cout << "Difference of 1._g - CLHEP::g is " << std::abs(1._g - (1 * CLHEP::g)) << endl;
0066 }