G4TriangularPrism Class Reference

Triangular prism. More...

#include <G4TriangularPrism.h>

Public Member Functions
	G4TriangularPrism (const G4String &name, double r1, double phi1, double r2, double phi2, G4double halfZlen)
	Constructor. More...
	G4TriangularPrism (const G4TriangularPrism &)=delete
	Copy constructor (disabled).
G4TriangularPrism &	operator= (const G4TriangularPrism &)=delete
	Operator = (disabled).
	~G4TriangularPrism ()
	Destructor. More...
G4VSolid *	getSolid ()
	Get solid. More...

Private Attributes
G4Box *	m_box
	Box.
G4Tubs *	m_tube
	Tube.
G4IntersectionSolid *	m_is
	Intersection.

Detailed Description

Triangular prism.

Definition at line 35 of file G4TriangularPrism.h.

Constructor & Destructor Documentation

G4TriangularPrism()

G4TriangularPrism	(	const G4String &	name,
		double	r1,
		double	phi1,
		double	r2,
		double	phi2,
		G4double	halfZlen
	)

Constructor.

Parameters

[in]	name	Name of the volume.
[in]	r1	First point radius.
[in]	phi1	First point angle.
[in]	r2	Second point radius.
[in]	phi2	Second point angle.
[in]	halfZlen	Half of Z length.

Coordinates of the third point are (0, 0)

Definition at line 23 of file G4TriangularPrism.cc.

{
  double rl;
  double rs;
  double delta_phi;
  double tg_alpha;
  double alpha;
  double sin_alpha;
  double cos_alpha;
  G4Transform3D t;
  m_is = nullptr;
  try {
    m_tube = new G4Tubs("Tube_" + name, 0., std::max(r1, r2), halfZlen,
                        phi1, fabs(phi2 - phi1));
  } catch (std::bad_alloc& ba) {
    goto err_mem1;
  }
  if (r1 >= r2) {
    rl = r1;
    rs = r2;
  } else {
    rl = r2;
    rs = r1;
  }
  delta_phi = phi2 - phi1;
  tg_alpha = (rl / rs - cos(delta_phi)) / sin(delta_phi);
  alpha = atan(tg_alpha);
  sin_alpha = tg_alpha / sqrt(1.0 + tg_alpha * tg_alpha);
  cos_alpha = 1.0 / sqrt(1.0 + tg_alpha * tg_alpha);
  try {
    m_box = new G4Box("Box_" + name, rl * sin_alpha,
                      rl * cos_alpha, halfZlen);
  } catch (std::bad_alloc& ba) {
    goto err_mem2;
  }
  if (r1 >= r2)
    t = G4RotateZ3D((phi1 + alpha) * CLHEP::rad - 90.0 * CLHEP::deg);
  else
    t = G4RotateZ3D((phi2 - alpha) * CLHEP::rad + 90.0 * CLHEP::deg);
  try {
    m_is = new G4IntersectionSolid(name, m_tube, m_box, t);
  } catch (std::bad_alloc& ba) {
    goto err_mem3;
  }
  return;
err_mem3:
  delete m_box;
err_mem2:
  delete m_tube;
err_mem1:
  throw (std::bad_alloc());
}

~G4TriangularPrism()

~G4TriangularPrism

(

)

Destructor.

The resulting solid m_is (access by getSolid()) is not deleted here. It is meant to be deleted by Geant.

Definition at line 79 of file G4TriangularPrism.cc.

Member Function Documentation

getSolid()

G4VSolid* getSolid ( )

inline

Get solid.

Returns

Solid corresponding to the prism.

The resulting G4VSolid* is allocated via operator new and must be deleted when is is not necessary. If it is used in Geant geometry, do not delete it.

Definition at line 79 of file G4TriangularPrism.h.

---

The documentation for this class was generated from the following files:

• klm/eklm/geometry/include/G4TriangularPrism.h
• klm/eklm/geometry/src/G4TriangularPrism.cc