Geometry/test/Cons_.cpp

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/Cons.h"
0005 #include "Utilities/Testing/interface/CppUnit_testdriver.icpp"
0006 #include "cppunit/extensions/HelperMacros.h"
0007 #include <G4Cons.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 testCons : public CppUnit::TestFixture {
0018   CPPUNIT_TEST_SUITE(testCons);
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 testCons::matched_g4_and_dd(void) {
0028   double rIn1 = 10. * cm;
0029   double rOut1 = 15. * cm;
0030   double rIn2 = 20. * cm;
0031   double rOut2 = 25. * cm;
0032   double zhalf = 20. * cm;
0033   double startPhi = 0. * deg;
0034   double deltaPhi = 90. * deg;
0035   string name("fred1");
0036
0037   G4Cons g4(name, rIn1, rOut1, rIn2, rOut2, zhalf, startPhi, deltaPhi);
0038   DDI::Cons dd(zhalf, rIn1, rOut1, rIn2, rOut2, startPhi, deltaPhi);
0039   DDCons dds = DDSolidFactory::cons(name, zhalf, rIn1, rOut1, rIn2, rOut2, startPhi, deltaPhi);
0040   cout << endl;
0041   dd.stream(cout);
0042   cout << endl;
0043
0044   double g4v = g4.GetCubicVolume() / cm3;
0045   double ddv = dd.volume() / cm3;
0046   double ddsv = dds.volume() / cm3;
0047
0048   cout << "\tg4 volume = " << g4v << " cm3" << endl;
0049   cout << "\tdd volume = " << ddv << " cm3" << endl;
0050   cout << "\tDD Information: " << dds << " vol=" << ddsv << " cm3" << endl;
0051
0052   CPPUNIT_ASSERT(abs(g4v - ddv) < numeric_limits<float>::epsilon());
0053   CPPUNIT_ASSERT(abs(g4v - ddsv) < numeric_limits<float>::epsilon());
0054 }
0055
0056 CPPUNIT_TEST_SUITE_REGISTRATION(testCons);