#include <Cons.h>

Inheritance diagram for DDI::Cons:

Public Member Functions
	Cons (double zhalf, double rInMinusZ, double rOutMinusZ, double rInPlusZ, double rOutPlusZ, double startPhi, double deltaPhi)

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 Cons.h.

Constructor & Destructor Documentation

◆ Cons()

DDI::Cons::Cons	(	double	zhalf,
		double	rInMinusZ,
		double	rOutMinusZ,
		double	rInPlusZ,
		double	rOutPlusZ,
		double	startPhi,
		double	deltaPhi
	)

Definition at line 14 of file Cons.cc.

     : Solid(DDSolidShape::ddcons) {
   p_.emplace_back(zhalf);
   p_.emplace_back(rInMinusZ);
   p_.emplace_back(rOutMinusZ);
   p_.emplace_back(rInPlusZ);
   p_.emplace_back(rOutPlusZ);
   p_.emplace_back(startPhi);
   p_.emplace_back(deltaPhi);
 }

References ddcons, SiPixelRawToDigiRegional_cfi::deltaPhi, and DDI::Solid::p_.

Member Function Documentation

◆ stream()

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

overridevirtual

Reimplemented from DDI::Solid.

Definition at line 31 of file Cons.cc.

                                          {
   os << " zhalf=" << convertMmToCm(p_[0]) << " rIn-Z=" << convertMmToCm(p_[1]) << " rOut-Z=" << convertMmToCm(p_[2])
      << " rIn+Z=" << convertMmToCm(p_[3]) << " rOut+Z=" << convertMmToCm(p_[4])
      << " startPhi=" << convertRadToDeg(p_[5]) << " deltaPhi=" << convertRadToDeg(p_[6]);
 }

References geant_units::operators::convertMmToCm(), and angle_units::operators::convertRadToDeg().

◆ volume()

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

overridevirtual

Reimplemented from DDI::Solid.

Definition at line 37 of file Cons.cc.

                              {
   /* zhalf is the half length of the cone,
      phiTo is always clockwise rotated from phiFrom 
      rInMinusZ is always smaller than rOutMinusZ (same for r*PlusZ)
      They are the distances relative to the rotation axes */
  
   /* calculation normalize from 0 to z */
  
   /* The function f=rInMinusZ+((rInPlusZ-rInMinusZ)/z)*x defines
      the radius of the the rotation from 0 to z. Raised to the power
      of 2 integrated on x from 0 to z. Multiplied by pi, gives the
      volume that needs to substracted from the other volume */
  
   /* f^2=rInMinusZ*rInMinusZ+2*rInMinusZ*((rInPlusZ-rInMinusZ)/z)*x+((rInPlusZ-rInMinusZ)*(rInPlusZ-rInMinusZ)*x*x)/(z*z) */
  
   /* primitive of f^2 is: rInMinusZ*rInMinusZ*x+rInMinusZ*((rInPlusZ-rInMinusZ)/z)*(x*x)+(rInPlusZ-rInMinusZ)*(rInPlusZ-rInMinusZ)*(x*x*x)/(3*z*z) */
  
   /*integration from 0 to z yields: pi*( rInMinusZ*rInMinusZ*z+rInMinusZ*(rInPlusZ-rInMinusZ)*z+((rInPlusZ-rInMinusZ)*(rInPlusZ-rInMinusZ)*z)/(3) ) */
  
   double zhalf = p_[0];
   double rInMinusZ = p_[1];
   double rOutMinusZ = p_[2];
   double rInPlusZ = p_[3];
   double rOutPlusZ = p_[4];
   //double phiFrom=p_[5]/rad;
   double deltaPhi = std::abs(p_[6]);
   double z = 2 * zhalf;
  
   double volume1 = 1_pi * (rInPlusZ * rInPlusZ + rInMinusZ * rInMinusZ + rInMinusZ * rInPlusZ) * z / 3;
  
   double volume2 = 1_pi * (rOutPlusZ * rOutPlusZ + rOutMinusZ * rOutMinusZ + rOutMinusZ * rOutPlusZ) * z / 3;
  
   double slice = deltaPhi / (2_pi);
   double volume = slice * (volume2 - volume1);
  
   return volume;
 }

References funct::abs(), SiPixelRawToDigiRegional_cfi::deltaPhi, and mergeAndRegister::slice.