File indexing completed on 2024-05-10 02:21:26
0001 #include <CLHEP/Units/SystemOfUnits.h>
0002 #include "DetectorDescription/Core/interface/DDSolid.h"
0003 #include "DetectorDescription/Core/interface/DDSolidShapes.h"
0004 #include "DetectorDescription/Core/interface/Torus.h"
0005 #include "Utilities/Testing/interface/CppUnit_testdriver.icpp"
0006 #include "cppunit/extensions/HelperMacros.h"
0007 #include <G4Torus.hh>
0008 #include <cmath>
0009 #include <limits>
0010 #include <string>
0011
0012 using namespace std;
0013 using CLHEP::cm;
0014 using CLHEP::cm3;
0015 using CLHEP::deg;
0016
0017 class testTorus : public CppUnit::TestFixture {
0018 CPPUNIT_TEST_SUITE(testTorus);
0019 CPPUNIT_TEST(matched_g4_and_dd);
0020
0021 CPPUNIT_TEST_SUITE_END();
0022
0023 public:
0024 void matched_g4_and_dd(void);
0025 };
0026
0027 void testTorus::matched_g4_and_dd(void) {
0028 double rMin = 10. * cm;
0029 double rMax = 15. * cm;
0030 double radius = 100. * cm;
0031 double sPhi = 0. * deg;
0032 double dPhi = 90. * deg;
0033 string name("fred1");
0034
0035 G4Torus g4(name, rMin, rMax, radius, sPhi, dPhi);
0036 DDI::Torus dd(rMin, rMax, radius, sPhi, dPhi);
0037 DDTorus dds = DDSolidFactory::torus(name, rMin, rMax, radius, sPhi, dPhi);
0038 cout << endl;
0039 dd.stream(cout);
0040 cout << endl;
0041
0042 double g4v = g4.GetCubicVolume() / cm3;
0043 double ddv = dd.volume() / cm3;
0044 double ddsv = dds.volume() / cm3;
0045
0046 cout << "\tg4 volume = " << g4v << " cm3" << endl;
0047 cout << "\tdd volume = " << ddv << " cm3" << endl;
0048 cout << "\tDD Information: " << dds << " vol= " << ddsv << " cm3" << endl;
0049
0050 CPPUNIT_ASSERT(abs(g4v - ddv) < numeric_limits<float>::epsilon());
0051 CPPUNIT_ASSERT(abs(g4v - ddsv) < numeric_limits<float>::epsilon());
0052 }
0053
0054 CPPUNIT_TEST_SUITE_REGISTRATION(testTorus);
