JACOUSTICS Namespace Reference

Auxiliary classes and methods for acoustic position calibration. More...

Namespaces
	JGEOMETRY
	Auxiliary namespace to encapsulate different geometries.
	JMODEL
	Auxiliary namespace to encapsulate different model parameters.

Classes
struct	JAbstractSoundVelocity
	Interface for depth dependend velocity of sound. More...
struct	JAttenuationLength
	Auxiliary data structure for calculation of attenuation length. More...
struct	JWaveform
	Utility class for emitter power and frequency. More...
struct	JCounter
	Acoustic counter. More...
struct	JEKey
	Emitter key. More...
struct	JEmitter
	Acoustic emitter. More...
struct	JEmitterID
	Auxiliary class for emitter identification. More...
struct	JEvent
	Acoustic event. More...
struct	JEventOverlap
	Match of two events considering overlap in time. More...
struct	JFit
	Acoustic single fit. More...
struct	JHead
	Acoustic event header. More...
struct	JEvt
	Acoustic event fit. More...
struct	getModel
	Auxiliary data structure to convert event to model. More...
struct	getEvt
	Auxiliary data structure to convert model to event. More...
struct	JFitParameters
	Fit parameters. More...
struct	JHit
	Acoustics hit. More...
struct	JKatoomba
	Template definition of fit function of acoustic model. More...
struct	JKatoomba< JType >
	Auxiliary base class for fit function of acoustic model. More...
struct	JKatoomba< JAbstractMinimiser >
	Template specialisation of fit function of acoustic model based on JAbstractMinimiser minimiser. More...
struct	TMatrixDS
	Auxiliary data structure for compatibility of symmetric matrix. More...
struct	JKatoomba< JEstimator >
	Template specialisation of fit function of acoustic model based on linear approximation. More...
struct	JKatoomba< JSimplex >
	Template specialisation of fit function of acoustic model based on JSimplex minimiser. More...
class	JGandalf
	Place holder for custom implementation. More...
struct	JKatoomba< JGandalf >
	Template specialisation of fit function of acoustic model based on JGandalf minimiser. More...
struct	JKatoomba_t
	Worker class for complete fit procedure of acoustic model. More...
struct	JMechanics
	Auxiliary data structure for parameters of mechanical model. More...
struct	JDetectorMechanics
	Auxiliary data structure for mechanical model parameters of strings in a given detector. More...
struct	JModel
	Model for fit to acoustics data. More...
struct	JNarrabri
	Dynamic position calibration. More...
struct	JReceiver
	Acoustic receiver. More...
struct	JSoundVelocity
	Implementation for depth dependend velocity of sound. More...
struct	JSuperEvt
	Acoustic super event fit. More...
struct	getSuperEvt
	Auxiliary data structure to convert model to super event. More...
struct	JTransceiver
	Acoustic transceiver. More...
struct	JTransmission
	Acoustic transmission. More...
struct	JTransmission_t
	Acoustic transmission identifier. More...
class	JTriggerOutput
	Set of acoustics events. More...
struct	JTriggerParameters
	Trigger parameters. More...
struct	hit_type
	Acoustic hit. More...
class	JMatch3D
	3D match criterion for acoustic signals. More...

Typedefs
typedef JGEOMETRY::JDetector	JGeometry
	Type definition of detector geometry. More...
typedef JTOOLS::JRange< double >	JTimeRange
	Type definition for time range (unit [s]). More...

Enumerations
enum	JMatrix_t { SVD_t = 1, LDU_t = 2 }
	Algorithm for matrix inversion. More...

Functions
double	getQ (const double D_m, const double f_kHz, const double d_m)
	Get relative quality for given frequency at given distance. More...
template<class T , class JTypename_t , class JComparator_t >
JVector3D	getPosition (T __begin, T __end, const JPredicate< JTypename_t, JComparator_t > &predicate)
	Get position from element in data which corresponds to given predicate. More...
template<class T , class JTypename_t , class JComparator_t >
JVector3D	getPosition (T __begin, T __end, const JPredicate< JTypename_t, JComparator_t > &predicate, const JVector3D &position)
	Get position from element in data which corresponds to given predicate. More...
double	getUNIXTime (const KM3NETDAQ::JDAQChronometer &chronometer)
	Get UNIX time of given DAQ object. More...
static bool	operator< (const JCounter &first, const JCounter &second)
	Less-than operator for two counters. More...
static bool	operator== (const JCounter &first, const JCounter &second)
	Equals operator for two counters. More...
bool	operator< (const JHead &first, const JHead &second)
	Less than operator for acoustics event headers. More...
double	getQuality (const double chi2, const int NDF)
	Get quality of fit. More...
template<class T >
size_t	getNumberOfEmitters (T __begin, T __end)
	Get number of emitters. More...
JMODEL::JString	getString (const JFit &fit)
	Get model parameters of string. More...
JFit	getFit (const int id, const JMODEL::JString &string)
	Get fit parameters of string. More...
template<class T >
double	getWeight (T __begin, T __end)
	Get total weight of data points. More...
JTreeParameters	getTreeParameters (const JType< JEvent > &type)
JTreeParameters	getTreeParameters (const JType< JEvt > &type)
JTreeParameters	getTreeParameters (const JType< JSuperEvt > &type)
static bool	operator< (const JTransmission &first, const JTransmission &second)
	Less-than operator for two transmissions. More...
static bool	operator== (const JTransmission &first, const JTransmission &second)
	Equals operator for two transmissions. More...
double	getVelocity (const double D, const double S, const double T)
	Get sound velocity according UNESCO equation. More...

Variables
static const double	TOAMIN_S = 0.0000000
	Minimal allowed time-of-arrival [s]. More...
static const double	TOAMAX_S = 0.6310912
	Maximal allowed time-of-arrival [s]. More...
static const JAttenuationLength	getAttenuationLength (13.2, 38.0, 8.0)
	Function object to calculate attenutation length. More...
static JEmitterID	getEmitterID
	Function object for emitter identification. More...
static JDetectorMechanics	getMechanics
	Function object to get string mechanics. More...
static const JSoundVelocity	getSoundVelocity (1541.0,-17.0e-3,-2000.0)
	Function object for velocity of sound. More...

Detailed Description

Auxiliary classes and methods for acoustic position calibration.

Author

mdejong

Typedef Documentation

typedef JGEOMETRY::JDetector JACOUSTICS::JGeometry

Type definition of detector geometry.

Definition at line 625 of file JGeometry.hh.

typedef JTOOLS::JRange<double> JACOUSTICS::JTimeRange

Type definition for time range (unit [s]).

Definition at line 19 of file JAcoustics/JTimeRange.hh.

Enumeration Type Documentation

enum JACOUSTICS::JMatrix_t

Algorithm for matrix inversion.

Enumerator
SVD_t	SVD.
LDU_t	LDU.

Definition at line 335 of file JKatoomba.hh.

                 {
    SVD_t   =   1,      //!< SVD
    LDU_t   =   2       //!< LDU
  };

Function Documentation

double JACOUSTICS::getQ ( const double  D_m, const double  f_kHz, const double  d_m  )

inline

Get relative quality for given frequency at given distance.

Parameters

D_m	depth [m]
f_kHz	frequency [kHz]
d_m	distance [m]

Returns

quality

Definition at line 92 of file JAcousticsSupportkit.hh.

  {
    return exp(-0.5 * d_m / getAttenuationLength(D_m, f_kHz)) / d_m;
  }

template<class T , class JTypename_t , class JComparator_t >

JVector3D JACOUSTICS::getPosition ( T  __begin, T  __end, const JPredicate< JTypename_t, JComparator_t > &  predicate  )

inline

Get position from element in data which corresponds to given predicate.

Parameters

__begin	begin of data
__end	end of data
predicate	predicate

Returns

position

Definition at line 37 of file JAcousticsToolkit.hh.

  {
    T p = std::find_if(__begin, __end, predicate);

    if (p != __end) 
      return p->getPosition();
    else
      THROW(JValueOutOfRange, "No element in container which corresponds to given predicate.");
  }

template<class T , class JTypename_t , class JComparator_t >

JVector3D JACOUSTICS::getPosition ( T  __begin, T  __end, const JPredicate< JTypename_t, JComparator_t > &  predicate, const JVector3D &  position  )

inline

Get position from element in data which corresponds to given predicate.

Parameters

__begin	begin of data
__end	end of data
predicate	predicate
position	default position

Returns

position

Definition at line 60 of file JAcousticsToolkit.hh.

  {
    try {
      return getPosition(__begin, __end, predicate);
    }
    catch(const std::exception&) {
      return position;
    }
  }

double JACOUSTICS::getUNIXTime ( const KM3NETDAQ::JDAQChronometer &  chronometer )

inline

Get UNIX time of given DAQ object.

Parameters

chronometer

chronometer

Returns

UNIX time [s]

Definition at line 80 of file JAcousticsToolkit.hh.

  {
    return chronometer.getTimesliceStart().getTimeNanoSecond() * 1.0e-9;
  }

static bool JACOUSTICS::operator< ( const JCounter &  first, const JCounter &  second  )

inlinestatic

Less-than operator for two counters.

Parameters

first	first counter
second	second counter

Returns

true if first counter less than second; else false

Definition at line 100 of file JAcoustics/JCounter.hh.

  {
    return first.getCounter() < second.getCounter();
  }

static bool JACOUSTICS::operator== ( const JCounter &  first, const JCounter &  second  )

inlinestatic

Equals operator for two counters.

Parameters

first	first counter
second	second counter

Returns

true if first counter equal to second; else false

Definition at line 113 of file JAcoustics/JCounter.hh.

  {
    return first.getCounter() == second.getCounter();
  }

bool JACOUSTICS::operator< ( const JHead &  first, const JHead &  second  )

inline

Less than operator for acoustics event headers.

The less than operator is applied to the object identifier, the UNIX start time and the UNIX stop time, respectively.

Parameters

first	first header
second	second header

Returns

true if first event header earliear than second; else false

Definition at line 276 of file JAcoustics/JEvt.hh.

  {
    if (first.oid == second.oid) {

      if (first.UNIXTimeStart == second.UNIXTimeStart) {

        return first.UNIXTimeStop  < second.UNIXTimeStop;

      } else {

        return first.UNIXTimeStart < second.UNIXTimeStart;
      }

    } else {

      return first.oid < second.oid;
    }
  }

double JACOUSTICS::getQuality ( const double  chi2, const int  NDF  )

inline

Get quality of fit.

The larger the quality, the better the fit.

Parameters

chi2	chi2
NDF	number of degrees of freedom

Returns

quality

Definition at line 29 of file JAcoustics/JEvtToolkit.hh.

  {
    return NDF  -  0.25 * chi2 / NDF;
  }

template<class T >

size_t JACOUSTICS::getNumberOfEmitters ( T  __begin, T  __end  )

inline

Get number of emitters.

Parameters

__begin	begin of events
__end	end of events

Returns

number of emitters

Definition at line 43 of file JAcoustics/JEvtToolkit.hh.

  {
    std::set<int> buffer;

    for (T i = __begin; i != __end; ++i) {
      if (!i->empty()) {
        buffer.insert(i->getID());
      }
    }

    return buffer.size();
  }

JMODEL::JString JACOUSTICS::getString ( const JFit &  fit )

inline

Get model parameters of string.

Parameters

fit	fit parameters

Returns

model parameters

Definition at line 63 of file JAcoustics/JEvtToolkit.hh.

  {
    return JMODEL::JString(fit.tx, fit.ty, fit.tx2, fit.ty2, fit.vs);
  }

JFit JACOUSTICS::getFit ( const int  id, const JMODEL::JString &  string  )

inline

Get fit parameters of string.

Parameters

id	identifier
string	model parameters

Returns

fit parameters

Definition at line 76 of file JAcoustics/JEvtToolkit.hh.

  {
    return JFit(id, string.tx, string.ty, string.tx2, string.ty2, string.vs);
  }

template<class T >

double JACOUSTICS::getWeight ( T  __begin, T  __end  )

inline

Get total weight of data points.

Parameters

__begin	begin of data
__end	end of data

Returns

total weight

Definition at line 60 of file JKatoomba.hh.

  {
    double W = 0;

    for (T i = __begin; i != __end; ++i) {
      W += i->getWeight();
    }

    return W;
  }

JTreeParameters JACOUSTICS::getTreeParameters ( const JType< JEvent > &  type )

inline

Parameters

type

data type

Returns

TTree parameters

Definition at line 35 of file JAcoustics/JSupport.hh.

  {
    return JTreeParameters("ACOUSTICS", "acoustics_event", "",  2, 65536, 2);
  }

JTreeParameters JACOUSTICS::getTreeParameters ( const JType< JEvt > &  type )

inline

Parameters

type

data type

Returns

TTree parameters

Definition at line 44 of file JAcoustics/JSupport.hh.

  {
    return JTreeParameters("ACOUSTICS_FIT", "acoustics_fit", "",  2, 65536, 2);
  }

JTreeParameters JACOUSTICS::getTreeParameters ( const JType< JSuperEvt > &  type )

inline

Parameters

type

data type

Returns

TTree parameters

Definition at line 53 of file JAcoustics/JSupport.hh.

  {
    return JTreeParameters("ACOUSTICS_SUPER_FIT", "acoustics_super_fit", "",  2, 65536, 2);
  }

static bool JACOUSTICS::operator< ( const JTransmission &  first, const JTransmission &  second  )

inlinestatic

Less-than operator for two transmissions.

The less-than operator is first applied to the time-of-emission and then to the identifier.

Parameters

first	first transmission
second	second transmission

Returns

true if first transmission earlier than second; else false

Definition at line 288 of file JTransmission.hh.

  {
    if (first.getToE() == second.getToE())
      return first.getID()  < second.getID();
    else
      return first.getToE() < second.getToE();
  }

static bool JACOUSTICS::operator== ( const JTransmission &  first, const JTransmission &  second  )

inlinestatic

Equals operator for two transmissions.

The equal operator is applied to the time-of-emission and the identifier.

Parameters

first	first transmission
second	second transmission

Returns

true if first transmission equal to second; else false

Definition at line 306 of file JTransmission.hh.

  {
    return (first.getID()  == second.getID()   &&
            first.getToE() == second.getToE());
  }

double JACOUSTICS::getVelocity ( const double  D, const double  S, const double  T  )

inline

Get sound velocity according UNESCO equation.

Parameters

D	Depth [m]
S	Salinity [ppk]
T	Temperature [C]

Returns

sound velocity [m/s]

Definition at line 23 of file JUNESCO.hh.

  {
    return (1448.96
            + 4.591 * T
            - 5.304e-2  * T*T
            + 2.374e-4  * T*T*T
            + 1.340 * (S - 35)
            + 1.630e-2  * D
            + 1.675e-7  * D*D
            - 1.025e-2  * T * (S - 35)
            - 7.139e-13 * T*D*D*D);
  }

Variable Documentation

const double JACOUSTICS::TOAMIN_S = 0.0000000

static

Minimal allowed time-of-arrival [s].

Definition at line 20 of file JAcousticsSupportkit.hh.

const double JACOUSTICS::TOAMAX_S = 0.6310912

static

Maximal allowed time-of-arrival [s].

Definition at line 21 of file JAcousticsSupportkit.hh.

const JAttenuationLength JACOUSTICS::getAttenuationLength(13.2,38.0,8.0)

static

Function object to calculate attenutation length.

JEmitterID JACOUSTICS::getEmitterID

static

Function object for emitter identification.

Definition at line 119 of file JEmitterID.hh.

JDetectorMechanics JACOUSTICS::getMechanics

static

Function object to get string mechanics.

Definition at line 243 of file JMechanics.hh.

const JSoundVelocity JACOUSTICS::getSoundVelocity(1541.0,-17.0e-3,-2000.0)

static

Function object for velocity of sound.