void
testPolyhedra::matched_g4_and_dd( void )
{
int sides = 3;
double phiStart = 45.*deg;
double phiTotal = 325.*deg;
double inner[] = { 0, 0, 0, 0, 0, 0, 0, 0, 0 };
vector<double> rInner( inner, inner + sizeof( inner ) / sizeof( double ));
double outer[] = { 0, 10, 10, 5, 5, 10, 10, 2, 2 };
vector<double> rOuter( outer, outer + sizeof( outer ) / sizeof( double ));
double pl[] = { 5, 7, 9, 11, 25, 27, 29, 31, 35 };
vector<double> z( pl, pl + sizeof( pl ) / sizeof( double ));
string name( "fred1" );
G4Polyhedra g4( name, phiStart, phiTotal, sides, z.size(), &z[0], &rInner[0], &rOuter[0] );
DDI::Polyhedra dd( sides, phiStart, phiTotal, z, rInner, rOuter );
DDPolyhedra dds = DDSolidFactory::polyhedra( name, sides, phiStart, phiTotal, z, rInner, rOuter );
cout << endl;
dd.stream( cout );
cout << endl;
double g4v = g4.GetCubicVolume()/cm3;
double ddv = dd.volume()/cm3;
double ddsv = dds.volume()/cm3;
cout << "\tg4 volume = " << g4v <<" cm3" << endl;
cout << "\tdd volume = " << ddv << " cm3" << endl;
cout << "\tDD Information: " << dds << " vol=" << ddsv << " cm3" << endl;
// FIXME: volumes do not match!
// CPPUNIT_ASSERT( abs( g4v - ddv ) < numeric_limits<float>::epsilon());
// CPPUNIT_ASSERT( abs( g4v - ddsv ) < numeric_limits<float>::epsilon());
}
void
doPolyhedra2( const std::string& name, int sides, double phiStart, double phiTotal,
const std::vector<double> & z,
const std::vector<double> & r )
{
G4Polyhedra g4( name, phiStart, phiTotal, sides, z.size(), &r[0], &z[0] );
DDI::Polyhedra dd( sides, phiStart, phiTotal, z, r );
DDPolyhedra dds = DDSolidFactory::polyhedra( name, sides, phiStart, phiTotal, z, r );
dd.stream(std::cout);
std::cout << std::endl;
std::cout << "\tg4 volume = " << g4.GetCubicVolume()/cm3 <<" cm3" << std::endl;
std::cout << "\tdd volume = " << dd.volume()/cm3 << " cm3"<< std::endl;
std::cout << "\tDD Information: " << dds << " vol= " << dds.volume() << std::endl;
}