Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-10-25 10:04:30

0001 #include "CLHEP/Units/GlobalSystemOfUnits.h"
0002 #include "DetectorDescription/Core/interface/DDSolid.h"
0003 #include "DetectorDescription/Core/interface/DDSolidShapes.h"
0004 #include "DetectorDescription/Core/interface/Sphere.h"
0005 #include "Utilities/Testing/interface/CppUnit_testdriver.icpp"
0006 #include "cppunit/extensions/HelperMacros.h"
0007 #include <G4Sphere.hh>
0008 #include <cmath>
0009 #include <limits>
0010 #include <string>
0011 
0012 using namespace std;
0013 
0014 class testSphere : public CppUnit::TestFixture {
0015   CPPUNIT_TEST_SUITE(testSphere);
0016   CPPUNIT_TEST(matched_g4_and_dd);
0017 
0018   CPPUNIT_TEST_SUITE_END();
0019 
0020 public:
0021   void matched_g4_and_dd(void);
0022 };
0023 
0024 void testSphere::matched_g4_and_dd(void) {
0025   double innerRadius = 10. * cm;
0026   double outerRadius = 15. * cm;
0027   double startPhi = 0. * deg;
0028   double deltaPhi = 90. * deg;
0029   double startTheta = 0. * deg;
0030   double deltaTheta = 180. * deg;
0031   string name("fred1");
0032 
0033   G4Sphere g4(name, innerRadius, outerRadius, startPhi, deltaPhi, startTheta, deltaTheta);
0034   DDI::Sphere dd(innerRadius, outerRadius, startPhi, deltaPhi, startTheta, deltaTheta);
0035   DDSphere dds = DDSolidFactory::sphere(name, innerRadius, outerRadius, startPhi, deltaPhi, startTheta, deltaTheta);
0036   cout << endl;
0037   dd.stream(cout);
0038   cout << endl;
0039 
0040   double g4v = g4.GetCubicVolume() / cm3;
0041   double ddv = dd.volume() / cm3;
0042   double ddsv = dds.volume() / cm3;
0043 
0044   cout << "\tg4 volume = " << g4v << " cm3" << endl;
0045   cout << "\tdd volume = " << ddv << " cm3" << endl;
0046   cout << "\tDD Information: " << dds << " vol=" << ddsv << " cm3" << endl;
0047 
0048   CPPUNIT_ASSERT(abs(g4v - ddv) < numeric_limits<float>::epsilon());
0049   CPPUNIT_ASSERT(abs(g4v - ddsv) < numeric_limits<float>::epsilon());
0050 }
0051 
0052 CPPUNIT_TEST_SUITE_REGISTRATION(testSphere);