JGEOMETRY3D::JOmega3D Class Reference

Direction set covering (part of) solid angle. More...

#include <JOmega3D.hh>

Inheritance diagram for JGEOMETRY3D::JOmega3D:

Public Types
typedef std::pair< double, double >	range_type

Public Member Functions
	JOmega3D ()
	Default constructor. More...
	JOmega3D (const JAngle3D &dir)
	Constructor. More...
	JOmega3D (const double grid)
	Constructor. More...
	JOmega3D (const JAngle3D &dir, const range_type &theta, const double grid)
	Constructor. More...
void	configure (const JAngle3D &dir, const range_type &theta, const double grid)
	Configure direction set. More...
int	find (const JAngle3D &dir) const
	Get index of direction closest to given direction. More...

Detailed Description

Direction set covering (part of) solid angle.

The set consists of approximately equally spaced directions.
In this, the spacing in azimuth angle is zenith angle dependent.

Definition at line 66 of file JOmega3D.hh.

Member Typedef Documentation

typedef std::pair<double, double> JGEOMETRY3D::JOmega3D_t::range_type

inherited

Definition at line 32 of file JOmega3D.hh.

Constructor & Destructor Documentation

JGEOMETRY3D::JOmega3D::JOmega3D ( )

inline

Default constructor.

Definition at line 73 of file JOmega3D.hh.

               :
      JOmega3D_t()
    {}

JGEOMETRY3D::JOmega3D::JOmega3D ( const JAngle3D &  dir )

inline

Constructor.

Parameters

dir

direction

Definition at line 83 of file JOmega3D.hh.

                                  :
      JOmega3D_t()
    {
      push_back(dir);
    }

JGEOMETRY3D::JOmega3D::JOmega3D ( const double  grid )

inline

Constructor.

Parameters

grid	angular spacing [rad]

Definition at line 95 of file JOmega3D.hh.

                                :
      JOmega3D_t()
    {
      configure(JAngle3D(0.0,0.0), range_type(0.0,PI), grid);
    }

JGEOMETRY3D::JOmega3D::JOmega3D ( const JAngle3D &  dir, const range_type &  theta, const double  grid  )

inline

Constructor.

Parameters

dir	principal direction
theta	polar angle range [rad]
grid	angular spacing [rad]

Definition at line 110 of file JOmega3D.hh.

                                     :
      JOmega3D_t()
    {
      configure(dir, theta, grid);
    }

Member Function Documentation

void JGEOMETRY3D::JOmega3D::configure ( const JAngle3D &  dir, const range_type &  theta, const double  grid  )

inline

Configure direction set.

Parameters

dir	principal direction
theta	polar angle range [rad]
grid	angular spacing [rad]

Definition at line 126 of file JOmega3D.hh.

    {
      clear();

      // sanity checks

      double thetaMin = theta.first;
      double thetaMax = theta.second;

      if (thetaMin < 0.0) { thetaMin = 0.0; }
      if (thetaMin > PI)  { thetaMin = PI; }
      if (thetaMax < 0.0) { thetaMax = 0.0; }
      if (thetaMax > PI)  { thetaMax = PI; }

      if (thetaMax > thetaMin) {
        
        const JRotation3D R(dir);

        const double rad  = thetaMax - thetaMin;
        const double bin  = rad / floor(rad/(1.4*grid) + 0.5);        // polar angle step size
        const double unit = 2.0 * sqrt(1.0 - cos(grid));              // azimuth angle unit step size
        
        for (double theta = thetaMin; theta < thetaMax + 0.5*bin; theta += bin) {
          
          double step;
          
          if (theta > 0.5*bin && PI - theta > 0.5*bin)
            step = PI / floor(PI*sin(theta)/unit + 0.5);              // polar angle dependent step size
          else
            step = 2.0*PI;                                            // pole

          for (double phi = 0.0; phi < 2.0*PI - 0.5*step; phi += step) {
            push_back(JDirection3D(JAngle3D(theta,phi)).rotate_back(R));
          }
        }
      }
    }

int JGEOMETRY3D::JOmega3D_t::find ( const JAngle3D &  dir ) const

inlineinherited

Get index of direction closest to given direction.

Parameters

dir

direction

Returns

index (-1 if error)

Definition at line 41 of file JOmega3D.hh.

    {
      double dot   = -1.0;
      int    index = -1;

      for (const_iterator i = this->begin(); i != this->end(); ++i) {

        const double x = dir.getDot(*i);
        
        if (x > dot) {
          dot   = x;
          index = std::distance(this->begin(), i);
        }
      }
     
      return index;
    }

---

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

• JOmega3D.hh