DDI::Sphere Class Reference

#include <Sphere.h>

Inheritance diagram for DDI::Sphere:

Public Member Functions
	Sphere (double innerRadius, double outerRadius, double startPhi, double deltaPhi, double startZ, double deltaZ)
void	stream (std::ostream &) const override
double	volume () const override
Public Member Functions inherited from DDI::Solid
const std::vector< double > &	parameters () const
void	setParameters (std::vector< double > const &p)
DDSolidShape	shape () const
	Solid ()
	Solid (DDSolidShape shape)
virtual	~Solid ()

Additional Inherited Members
Protected Attributes inherited from DDI::Solid
std::vector< double >	p_
DDSolidShape	shape_

Detailed Description

Definition at line 9 of file Sphere.h.

Constructor & Destructor Documentation

DDI::Sphere::Sphere	(	double	innerRadius,
		double	outerRadius,
		double	startPhi,
		double	deltaPhi,
		double	startZ,
		double	deltaZ
	)

Definition at line 14 of file Sphere.cc.

References ddsphere, and DDI::Solid::p_.

  : Solid(DDSolidShape::ddsphere)     
{
  p_.emplace_back(innerRadius);
  p_.emplace_back(outerRadius);
  p_.emplace_back(startPhi);
  p_.emplace_back(deltaPhi);
  p_.emplace_back(startTheta);
  p_.emplace_back(deltaTheta);
}    

Member Function Documentation

void DDI::Sphere::stream ( std::ostream &  os ) const

overridevirtual

Reimplemented from DDI::Solid.

Definition at line 30 of file Sphere.cc.

References geant_units::operators::convertMmToCm(), geant_units::operators::convertRadToDeg(), and DDI::Solid::p_.

{
  os << " innerRadius=" << convertMmToCm( p_[0] )
     << " outerRadius=" << convertMmToCm( p_[1] )
     << " startPhi=" << convertRadToDeg( p_[2] )
     << " deltaPhi=" << convertRadToDeg( p_[3] )
     << " startTheta=" << convertRadToDeg( p_[4] )
     << " deltaTheta=" << convertRadToDeg( p_[5] );
}

double DDI::Sphere::volume ( void  ) const

overridevirtual

Reimplemented from DDI::Solid.

Definition at line 40 of file Sphere.cc.

References funct::cos(), DDI::Solid::p_, and geant_units::piRadians().

{
  double volume(0.);
  if ( std::fabs(p_[3]) <= 2._pi && std::fabs(p_[5]) <= piRadians ) {
    volume = std::fabs((p_[1]*p_[1]*p_[1] - p_[0]*p_[0]*p_[0])/3. * (std::cos(p_[4]+p_[5]) - std::cos(p_[4]))*p_[3]);
  } else if (std::fabs(p_[3]) <= 2._pi && std::fabs(p_[5]) > piRadians ) {
    volume = std::fabs((p_[1]*p_[1]*p_[1] - p_[0]*p_[0]*p_[0])/3. * (std::cos(p_[4]+p_[5]-180._deg) - std::cos(p_[4]))*p_[3]);
  } else if (std::fabs(p_[3]) > 2._pi && std::fabs(p_[5]) <= piRadians ) {
    volume = std::fabs((p_[1]*p_[1]*p_[1] - p_[0]*p_[0]*p_[0])/3. * (std::cos(p_[4]+p_[5]) - std::cos(p_[4]))*(p_[3]-p_[2]));
  }
  return volume;
}