Vertex locator. More...

Inheritance diagram for JACOUSTICS::JVelo:

Public Member Functions
	JVelo (const JAbstractSoundVelocity &V, const JVector3D &center, const double RMax_m, const double Xv_m, const double factor=1.0)
	Constructor.

template<class T >
int	operator() (const JVertex3D &vx, T __begin, T __end) const
	Check vertex.

template<class T >
vertex_type	operator() (const JPosition3D &position, const hit_type &root, T __begin, T __end) const
	Update vertex at given position.

template<class T >
vertex_type	operator() (const hit_type &root, T __begin, T __end, const int numberOfHits=0) const
	Locate vertex.

Public Attributes
const JAbstractSoundVelocity &	V

const double	Tmax_s

Detailed Description

Vertex locator.

Definition at line 275 of file JAcousticsTriggerProcessor.cc.

Constructor & Destructor Documentation

◆ JVelo()

JACOUSTICS::JVelo::JVelo	(	const JAbstractSoundVelocity &	V,
		const JVector3D &	center,
		const double	RMax_m,
		const double	Xv_m,
		const double	factor = 1.0 )

inline

Constructor.

Parameters

V	sound velocity
center	center for generation of vertices
RMax_m	radius for generation of vertices [m]
Xv_m	step size for generation of vertices [m]
factor	multiplication factor for corresponding time window

Definition at line 287 of file JAcousticsTriggerProcessor.cc.

                                         :
      V(V),
      Tmax_s(factor * Xv_m / V())
    {
      this->push_back(center);
 
      if (RMax_m > 0.0 && Xv_m > 0.0) {
        for (double x = 0.5*Xv_m; x <= RMax_m; x += Xv_m) {
          for (double y = 0.5*Xv_m; y <= RMax_m; y += Xv_m) {
            if (x*x + y*y <= RMax_m*RMax_m) {
              this->push_back(JVector3D(center.getX() - x, center.getY() + y, center.getZ()));
              this->push_back(JVector3D(center.getX() + x, center.getY() + y, center.getZ()));
              this->push_back(JVector3D(center.getX() - x, center.getY() - y, center.getZ()));
              this->push_back(JVector3D(center.getX() + x, center.getY() - y, center.getZ()));
            }
          }
        }
      }
    }

Member Function Documentation

◆ operator()() [1/3]

template<class T >

int JACOUSTICS::JVelo::operator()	(	const JVertex3D &	vx,
		T	__begin,
		T	__end ) const

inline

Check vertex.

Parameters

vx	vertex
__begin	begin of data
__end	end of data

Returns: number of hits

Definition at line 321 of file JAcousticsTriggerProcessor.cc.

    {
      const JPosition3D& p0 = vx.getPosition();
      const double       t0 = vx.getT();
 
      int    n0 = 0;
      double q0 = 0.0;
 
      for (T p = __begin; p != __end; ++p) {
 
        const double t1 = p->getToA() - V.getTime(p0.getDistance(p->getPosition()), p0.getZ(), p->getZ());
 
        if (fabs(t1 - t0) <= Tmax_s) {
          n0 += 1;
          q0 += p->getQ();
        }
      }
 
      return n0;
    }

◆ operator()() [2/3]

template<class T >

vertex_type JACOUSTICS::JVelo::operator()	(	const JPosition3D &	position,
		const hit_type &	root,
		T	__begin,
		T	__end ) const

inline

Update vertex at given position.

Parameters

position	position
root	root hit
__begin	begin of data
__end	end of data

Returns: vertex

Definition at line 353 of file JAcousticsTriggerProcessor.cc.

    {
      const JPosition3D& p0 = position;
      const double       t0 = root.getToA() - V.getTime(p0.getDistance(root.getPosition()), p0.getZ(), root.getZ());
 
      int    n0 = 1;
      double q0 = root.getQ();
 
      for (T p = __begin; p != __end; ++p) {
 
        const double t1 = p->getToA() - V.getTime(p0.getDistance(p->getPosition()), p0.getZ(), p->getZ());
 
        if (fabs(t1 - t0) <= Tmax_s) {
          n0 += 1;
          q0 += p->getQ();
        }
      }
 
      return vertex_type(p0, t0, n0, q0);
    }

◆ operator()() [3/3]

template<class T >

vertex_type JACOUSTICS::JVelo::operator()	(	const hit_type &	root,
		T	__begin,
		T	__end,
		const int	numberOfHits = 0 ) const

inline

Locate vertex.

Parameters

root	root hit
__begin	begin of data
__end	end of data
numberOfHits	number of hits

Returns: vertex

Definition at line 385 of file JAcousticsTriggerProcessor.cc.

    {
      vertex_type vertex;
 
      for (const_iterator i = this->cbegin(); i != this->cend(); ++i) {
 
        const JPosition3D& p0 = *i;
        const double       t0 = root.getToA() - V.getTime(p0.getDistance(root.getPosition()), p0.getZ(), root.getZ());
 
        int    n0 = 1;
        double q0 = root.getQ();
 
        for (T p = __begin; p != __end && n0 + distance(p, __end) >= numberOfHits && n0 + distance(p, __end) >= vertex.N; ++p) {
 
          const double t1 = p->getToA() - V.getTime(p0.getDistance(p->getPosition()), p0.getZ(), p->getZ());
 
          if (fabs(t1 - t0) <= Tmax_s) {
            n0 += 1;
            q0 += p->getQ();
          }
        }
 
        if (q0 > vertex.Q) {
          vertex = vertex_type(p0, t0, n0, q0);
        }
      }
 
      return vertex;
    }