Jpp  master_rocky-37-gf0c5bc59d
the software that should make you happy
Classes | Public Member Functions | Static Public Member Functions | Public Attributes | Protected Attributes | Private Attributes | Friends | List of all members
JCALIBRATE::JModel Struct Reference

Fit model. More...

#include <JFitK40.hh>

Inheritance diagram for JCALIBRATE::JModel:
JCALIBRATE::JK40Parameters JDETECTOR::JModule JCALIBRATE::JCombinatorics_t JTOOLS::JCombinatorics std::vector< JPMT > JEEP::JStatus JDETECTOR::JCalibration JGEOMETRY3D::JQuaternion3D JGEOMETRY3D::JPosition3D JDETECTOR::JLocation JLANG::JObjectID JCALIBRATE::JK40Parameters_t

Classes

struct  real_type
 Auxiliary data structure for derived quantities of a given PMT pair. More...
 

Public Member Functions

 JModel ()
 Default constructor. More...
 
 JModel (const JModule &module, const JK40Parameters &parameters, const JTDC_t::range_type &TDC, const int option)
 Constructor. More...
 
 JModel (const JModule &module, const JK40Parameters &parameters)
 Constructor. More...
 
JOption_t getOption () const
 Get fit option. More...
 
void setOption (const int option)
 Set fit option. More...
 
bool hasFixedTimeOffset () const
 Check if time offset is fixed. More...
 
double getFixedTimeOffset () const
 Get time offset. More...
 
int getIndex () const
 Get index of PMT used for fixed time offset. More...
 
void setIndex ()
 Set index of PMT used for fixed time offset. More...
 
size_t getN () const
 Get number of fit parameters. More...
 
int getIndex (int pmt) const
 Get index of parameter. More...
 
int getIndex (int pmt, JParameter_t JPMTParameters_t::*p) const
 Get index of parameter. More...
 
double getSigmaK40 () const
 Get intrinsic K40 arrival time spread. More...
 
void setSigmaK40 (const double sigma)
 Set intrinsic K40 arrival time spread. More...
 
const real_typegetReal (const pair_type &pair) const
 Get derived parameters. More...
 
double getValue (const pair_type &pair, const double dt_ns) const
 Get K40 coincidence rate. More...
 
double getValue (const pair_type &pair) const
 Get K40 coincidence rate. More...
 
int getIndex (JParameter_t JK40Parameters::*p) const
 Get index of parameter. More...
 
double getValue (const double ct) const
 Get K40 coincidence rate as a function of cosine angle between PMT axes. More...
 
const JK40ParametersgetK40Parameters () const
 Get K40 parameters. More...
 
void setK40Parameters (const JK40Parameters_t &parameters)
 Set K40 parameters. More...
 
int getIndex (JParameter_t JK40Parameters::*p) const
 Get index of parameter. More...
 
double getValue (const double ct) const
 Get K40 coincidence rate as a function of cosine angle between PMT axes. More...
 
const JK40Parameters_tgetGradient (const double ct) const
 Get gradient. More...
 
void reset ()
 Reset. More...
 
const JPMTgetPMT (const int index) const
 Get PMT. More...
 
JPMTgetPMT (const int index)
 Get PMT. More...
 
void setPMT (const int index, const JPMT &pmt)
 Set PMT. More...
 
JVector3D getCenter () const
 Get center of module based on crossing point of PMT axes. More...
 
void compile ()
 Compile module data. More...
 
void rotate (const JRotation3D &R)
 Rotate module. More...
 
void rotate (const JQuaternion3D &Q)
 Rotate module. More...
 
JPosition3Drotate (const JRotation3X &R)
 Rotate around X-axis. More...
 
JPosition3Drotate (const JRotation3Y &R)
 Rotate around Y-axis. More...
 
JPosition3Drotate (const JRotation3Z &R)
 Rotate around Z-axis. More...
 
void rotate (double &__x, double &__y, double &__z) const
 Rotate. More...
 
void rotate_back (const JRotation3D &R)
 Rotate back module. More...
 
void rotate_back (const JQuaternion3D &Q)
 Rotate back module. More...
 
JPosition3Drotate_back (const JRotation3X &R)
 Rotate back around X-axis. More...
 
JPosition3Drotate_back (const JRotation3Y &R)
 Rotate back around Y-axis. More...
 
JPosition3Drotate_back (const JRotation3Z &R)
 Rotate back around Z-axis. More...
 
void rotate_back (double &__x, double &__y, double &__z) const
 Rotate back. More...
 
void transform (const JRotation3D &R, const JVector3D &pos)
 Transformation of geometry (see method JGEOMETRY3D::JPosition3D::transform(const JRotation3D&, const JVector3D&)). More...
 
void transform (const JTransformation3D &T)
 Transformation of geometry. More...
 
JVector3Dtransform (const JMatrix3D &T)
 Transform. More...
 
JVector3Dtransform (const JMatrix3D &T)
 Transform. More...
 
JModuleset (const JVector3D &pos)
 Set position. More...
 
JModuleset (const double t0)
 Set time offset. More...
 
void set (const int bit)
 Set PMT status. More...
 
void set (const int bit)
 Set PMT status. More...
 
JModuleadd (const JVector3D &pos)
 Add position. More...
 
JModuleadd (const double t0)
 Add time offset. More...
 
JQuaternion3Dadd (const JQuaternion3D &quaternion)
 Add quaternion. More...
 
JModulesub (const JVector3D &pos)
 Subtract position. More...
 
JModulesub (const double t0)
 Subtract time offset. More...
 
JQuaternion3Dsub (const JQuaternion3D &quaternion)
 Subtract quaternion. More...
 
JModuleoperator+= (const JVector3D &pos)
 Add position. More...
 
JModuleoperator-= (const JVector3D &pos)
 Subtract position. More...
 
bool has (const int bit) const
 Test PMT status. More...
 
bool has (const int bit) const
 Test PMT status. More...
 
void reset (const int bit)
 Reset PMT status. More...
 
void reset (const int bit)
 Reset PMT status. More...
 
int getID () const
 Get identifier. More...
 
int & getID ()
 Get identifier. More...
 
void setID (const int id)
 Set identifier. More...
 
bool less (const JObjectID &object) const
 Less than method. More...
 
bool less (const int id) const
 Less than method. More...
 
bool less (const JLocation &location) const
 Less than method. More...
 
bool more (const int id) const
 More than method. More...
 
const JLocationgetLocation () const
 Get location. More...
 
JLocationgetLocation ()
 Get location. More...
 
void setLocation (const JLocation &location)
 Set location. More...
 
std::string toString () const
 Convert module location to string. More...
 
std::string toString (const std::string &fmt, const std::string target="%") const
 Convert module loation to string. More...
 
int getString () const
 Get string number. More...
 
int getFloor () const
 Get floor number. More...
 
const JPosition3DgetPosition () const
 Get position. More...
 
JPosition3DgetPosition ()
 Get position. More...
 
void setPosition (const JVector3D &pos)
 Set position. More...
 
 operator JAngle3D () const
 Type conversion operator. More...
 
 operator JVersor3D () const
 Type conversion operator. More...
 
 operator JVersor3D () const
 Type conversion operator. More...
 
void transform_back (const JRotation3D &R, const JVector3D &pos)
 Transform back position. More...
 
double getDot (const JAngle3D &angle) const
 Get dot product. More...
 
double getDot (const JVersor3D &dir) const
 Get dot product. More...
 
double getDot (const JVersor3Z &dir) const
 Get dot product. More...
 
double getDot (const JVector3D &vector) const
 Get dot product. More...
 
double getDot (const JVector3D &vector) const
 Get dot product. More...
 
double getDot (const JQuaternion3D &quaternion) const
 Get dot product. More...
 
 operator JVector2D () const
 Type conversion operator. More...
 
double getX () const
 Get x position. More...
 
double getY () const
 Get y position. More...
 
double getZ () const
 Get z position. More...
 
JVector3Dnegate ()
 Negate vector. More...
 
JQuaternion3Dnegate ()
 Negate quaternion. More...
 
JVector3Dmul (const double factor)
 Scale vector. More...
 
JFirst_t & mul (const JSecond_t &object)
 Multiply with object. More...
 
JQuaternion3Dmul (const double factor)
 Scale quaternion. More...
 
JQuaternion3Dmul (const JQuaternion3D &first, const JQuaternion3X &second)
 Quaternion multiplication. More...
 
JQuaternion3Dmul (const JQuaternion3D &first, const JQuaternion3Y &second)
 Quaternion multiplication. More...
 
JQuaternion3Dmul (const JQuaternion3D &first, const JQuaternion3Z &second)
 Quaternion multiplication. More...
 
JQuaternion3Dmul (const JQuaternion3D &first, const JQuaternion3D &second)
 Quaternion multiplication. More...
 
JQuaternion3Dmul (const JQuaternion3X &qx, const JQuaternion3Y &qy, const JQuaternion3Z &qz)
 Quaternion multiplication. More...
 
JFirst_t & mul (const JSecond_t &object)
 Multiply with object. More...
 
JVector3Ddiv (const double factor)
 Scale vector. More...
 
JQuaternion3Ddiv (const double factor)
 Scale quaternion. More...
 
bool equals (const JVector3D &vector, const double precision=std::numeric_limits< double >::min()) const
 Check equality. More...
 
bool equals (const JQuaternion3D &quaternion, const double precision=std::numeric_limits< double >::min()) const
 Check equality. More...
 
bool equals (const JStatus &status) const
 Equal method. More...
 
double getLengthSquared () const
 Get length squared. More...
 
double getLengthSquared () const
 Get length squared. More...
 
double getLength () const
 Get length. More...
 
double getLength () const
 Get length. More...
 
double getDistanceSquared (const JVector3D &pos) const
 Get squared of distance to point. More...
 
double getDistanceSquared (const JQuaternion3D &quaternion) const
 Get squared of distance to quaternion. More...
 
double getDistance (const JVector3D &pos) const
 Get distance to point. More...
 
double getDistance (const JQuaternion3D &quaternion) const
 Get distance to quaternion. More...
 
JVector3DgetCross (const JVector3D &first, const JVector3D &second)
 Get cross product. More...
 
const JQuaternion3DgetQuaternion () const
 Get quaternion. More...
 
JQuaternion3DgetQuaternion ()
 Get quaternion. More...
 
void setQuaternion (const JQuaternion3D &quaternion)
 Set quaternion. More...
 
 operator JVector3D () const
 Type conversion operator. More...
 
double getAngle () const
 Get rotation angle. More...
 
double getA () const
 Get a value. More...
 
double getB () const
 Get b value. More...
 
double getC () const
 Get c value. More...
 
double getD () const
 Get d value. More...
 
JQuaternion3DsetIdentity ()
 Set to identity quaternion. More...
 
JQuaternion3Dconjugate ()
 Conjugate quaternion. More...
 
bool isIdentity (const double precision=std::numeric_limits< double >::min()) const
 Test identity. More...
 
JQuaternion3D getConjugate () const
 Get conjugate of this quaternion. More...
 
JQuaternion3Dnormalise ()
 Normalise quaternion. More...
 
JQuaternion3Dpow (const double y)
 Raise quaternion to given power. More...
 
JQuaternion3Dinterpolate (const JQuaternion3D &object, const double alpha)
 Interpolation between quaternions. More...
 
const JCalibrationgetCalibration () const
 Get calibration. More...
 
JCalibrationgetCalibration ()
 Get calibration. More...
 
void setCalibration (const JCalibration &cal)
 Set calibration. More...
 
double getT0 () const
 Get time offset. More...
 
void setT0 (const double t0)
 Set time offset. More...
 
void addT0 (const double t0)
 Add time offset. More...
 
void subT0 (const double t0)
 Subtract time offset. More...
 
int getStatus () const
 Get status. More...
 
JStatusgetStatus ()
 Get status. More...
 
int getStatus (const int mask) const
 Get status. More...
 
void setStatus (const JStatus &status)
 Set status. More...
 
const JCombinatoricsgetCombinatorics () const
 Get combinatorics. More...
 
void configure (const int numberOfIndices)
 Configure. More...
 
size_t getNumberOfIndices () const
 Get number of indices. More...
 
size_t getNumberOfPairs () const
 Get number of pairs. More...
 
int getIndex (const int first, const int second) const
 Get index of pair of indices. More...
 
int getIndex (const pair_type &pair) const
 Get index of pair. More...
 
const pair_typegetPair (const int index) const
 Get pair of indices for given index. More...
 
template<class JComparator_t >
void sort (JComparator_t comparator)
 Sort address pairs. More...
 

Static Public Member Functions

static const JK40ParametersgetInstance ()
 Get default values. More...
 
static JDetectorVersiongetVersion ()
 Get detector version. More...
 
static void setVersion (const JVersion &version)
 Set detector version. More...
 
static bool compare (const JModule &first, const JModule &second, const double precision=1.0e-3)
 Compare modules. More...
 
static const JQuaternion3DgetIdentity ()
 Get identity quaternion. More...
 
static int getSign (const int first, const int second)
 Sign of pair of indices. More...
 
static int getSign (const pair_type &pair)
 Sign of pair of indices. More...
 

Public Attributes

JPMTParameters_t parameters [NUMBER_OF_PMTS]
 
JParameter_t R
 maximal coincidence rate [Hz] More...
 
JParameter_t p1
 1st order angle dependence coincidence rate More...
 
JParameter_t p2
 2nd order angle dependence coincidence rate More...
 
JParameter_t p3
 3rd order angle dependence coincidence rate More...
 
JParameter_t p4
 4th order angle dependence coincidence rate More...
 
JParameter_t cc
 fraction of signal correlated background More...
 

Protected Attributes

int __id
 
int string
 
int floor
 
double __x
 
double __y
 
double __z
 
double __a
 
double __b
 
double __c
 
double __d
 
double t0
 
int status
 
std::vector< pair_typezbuf1D
 
std::vector< std::vector< int > > zbuf2D
 

Private Attributes

int index
 index of PMT used for fixed time offset More...
 
double sigmaK40_ns = 0.54
 intrinsic K40 arrival time spread [ns] More...
 
JOption_t option
 
real_type real
 
JK40Parameters_t gradient
 

Friends

std::ostream & operator<< (std::ostream &out, const JModel &object)
 Write model parameters to output stream. More...
 

Detailed Description

Fit model.

In the absence of TDC constraints, the average time offset is fixed to zero.

Definition at line 750 of file JFitK40.hh.

Constructor & Destructor Documentation

◆ JModel() [1/3]

JCALIBRATE::JModel::JModel ( )
inline

Default constructor.

Definition at line 775 of file JFitK40.hh.

776  {}

◆ JModel() [2/3]

JCALIBRATE::JModel::JModel ( const JModule module,
const JK40Parameters parameters,
const JTDC_t::range_type TDC,
const int  option 
)
inline

Constructor.

Parameters
moduledetector module
parametersK40 parameters
TDCTDC constraints
optionoption

Definition at line 787 of file JFitK40.hh.

790  :
791  JModule (module),
792  JCombinatorics_t(module)
793  {
795 
796  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
797  this->parameters[i] = JPMTParameters_t::getInstance();
798  }
799 
800  for (JTDC_t::const_iterator i = TDC.first; i != TDC.second; ++i) {
801 
802  if (i->second != JTDC_t::WILDCARD) {
803 
804  this->parameters[i->second].t0.fix();
805 
806  } else {
807 
808  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
809  this->parameters[i].t0.fix();
810  }
811  }
812  }
813 
814  this->index = (TDC.first == TDC.second ? 0 : INVALID_INDEX);
815 
816  setOption(option);
817  }
JModule()
Default constructor.
Definition: JModule.hh:85
static const int INVALID_INDEX
invalid index
Definition: JFitK40.hh:60
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition: JDAQ.hh:26
JCombinatorics_t()
Default constructor.
void setK40Parameters(const JK40Parameters_t &parameters)
Set K40 parameters.
Definition: JFitK40.hh:658
void setOption(const int option)
Set fit option.
Definition: JFitK40.hh:855
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:1163
JOption_t option
Definition: JFitK40.hh:1168
int index
index of PMT used for fixed time offset
Definition: JFitK40.hh:1166
static const JPMTParameters_t & getInstance()
Get default values.
Definition: JFitK40.hh:459

◆ JModel() [3/3]

JCALIBRATE::JModel::JModel ( const JModule module,
const JK40Parameters parameters 
)
inline

Constructor.

Parameters
moduledetector module
parametersK40 parameters

Definition at line 826 of file JFitK40.hh.

827  :
828  JModule (module),
829  JCombinatorics_t(module)
830  {
832 
833  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
834  this->parameters[i] = JPMTParameters_t::getInstance();
835  }
836  }

Member Function Documentation

◆ getOption()

JOption_t JCALIBRATE::JModel::getOption ( ) const
inline

Get fit option.

Returns
option

Definition at line 844 of file JFitK40.hh.

845  {
846  return option;
847  }

◆ setOption()

void JCALIBRATE::JModel::setOption ( const int  option)
inline

Set fit option.

Parameters
optionoption

Definition at line 855 of file JFitK40.hh.

856  {
857  switch (option) {
858 
859  case FIT_PMTS_t:
860 
861  R .fix();
862  p1.fix();
863  p2.fix();
864  p3.fix();
865  p4.fix();
866 
867  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
868  parameters[i].bg.fix();
869  }
870 
871  break;
872 
874 
875  R .fix();
876 
877  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
878  parameters[i].bg.fix();
879  }
880 
881  break;
882 
884 
885  R .fix();
886  p1.fix();
887  p2.fix();
888  p3.fix();
889  p4.fix();
890 
891  break;
892 
893  case FIT_PMTS_QE_FIXED_t:
894 
895  R .fix();
896  p1.fix();
897  p2.fix();
898  p3.fix();
899  p4.fix();
900 
901  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
902  parameters[i].QE.fix();
903  parameters[i].bg.fix();
904  }
905 
906  break;
907 
908  case FIT_MODEL_t:
909 
910  cc.fix();
911 
912  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
913  parameters[i].QE.fix();
914  parameters[i].t0.fix();
915  parameters[i].bg.fix();
916  }
917 
918  break;
919 
920  default:
921 
922  THROW(JValueOutOfRange, "Invalid option " << option);
923  }
924 
925  this->option = static_cast<JOption_t>(option);
926  }
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
Definition: JException.hh:712
void fix()
Fix current value.
Definition: JFitK40.hh:280
Exception for accessing a value in a collection that is outside of its range.
Definition: JException.hh:180
JOption_t
Fit options.
Definition: JFitK40.hh:52
@ FIT_PMTS_QE_FIXED_t
fit parameters of PMTs with QE fixed
Definition: JFitK40.hh:56
@ FIT_PMTS_AND_ANGULAR_DEPENDENCE_t
fit parameters of PMTs and angular dependence of K40 rate
Definition: JFitK40.hh:54
@ FIT_MODEL_t
fit parameters of K40 rate and TTSs of PMTs
Definition: JFitK40.hh:57
@ FIT_PMTS_AND_BACKGROUND_t
fit parameters of PMTs and background
Definition: JFitK40.hh:55
@ FIT_PMTS_t
fit parameters of PMTs
Definition: JFitK40.hh:53
JParameter_t R
maximal coincidence rate [Hz]
Definition: JFitK40.hh:596
JParameter_t p1
1st order angle dependence coincidence rate
Definition: JFitK40.hh:597
JParameter_t p2
2nd order angle dependence coincidence rate
Definition: JFitK40.hh:598
JParameter_t p3
3rd order angle dependence coincidence rate
Definition: JFitK40.hh:599
JParameter_t p4
4th order angle dependence coincidence rate
Definition: JFitK40.hh:600
JParameter_t cc
fraction of signal correlated background
Definition: JFitK40.hh:601
JParameter_t t0
time offset [ns]
Definition: JFitK40.hh:564
JParameter_t bg
background [Hz/ns]
Definition: JFitK40.hh:565
JParameter_t QE
relative quantum efficiency [unit]
Definition: JFitK40.hh:562

◆ hasFixedTimeOffset()

bool JCALIBRATE::JModel::hasFixedTimeOffset ( ) const
inline

Check if time offset is fixed.

Returns
true if time offset fixed; else false

Definition at line 934 of file JFitK40.hh.

935  {
936  return index != INVALID_INDEX;
937  }

◆ getFixedTimeOffset()

double JCALIBRATE::JModel::getFixedTimeOffset ( ) const
inline

Get time offset.

Returns
time offset

Definition at line 945 of file JFitK40.hh.

946  {
947  double t0 = 0.0;
948  size_t N = 0;
949 
950  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
951  if (parameters[i].t0.isFree()) {
952  t0 += parameters[i].t0;
953  N += 1;
954  }
955  }
956 
957  return t0 /= N;
958  }

◆ getIndex() [1/7]

int JCALIBRATE::JModel::getIndex ( ) const
inline

Get index of PMT used for fixed time offset.

Returns
index

Definition at line 966 of file JFitK40.hh.

967  {
968  return index;
969  }

◆ setIndex()

void JCALIBRATE::JModel::setIndex ( )
inline

Set index of PMT used for fixed time offset.

Definition at line 975 of file JFitK40.hh.

976  {
977  if (index != INVALID_INDEX) {
978 
979  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
980 
981  if (parameters[i].status) {
982 
983  index = i;
984 
985  parameters[index].t0.fix();
986 
987  return;
988  }
989  }
990 
991  THROW(JValueOutOfRange, "No valid index.");
992  }
993  }

◆ getN()

size_t JCALIBRATE::JModel::getN ( ) const
inline

Get number of fit parameters.

Returns
number of parameters

Definition at line 1001 of file JFitK40.hh.

1002  {
1003  size_t N = JK40Parameters::getN();
1004 
1005  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
1006  N += parameters[i].getN();
1007  }
1008 
1009  return N;
1010  }
size_t getN() const
Get number of fit parameters.
Definition: JFitK40.hh:669
size_t getN() const
Get number of fit parameters.
Definition: JFitK40.hh:481

◆ getIndex() [2/7]

int JCALIBRATE::JModel::getIndex ( int  pmt) const
inline

Get index of parameter.

Parameters
pmtPMT
Returns
index

Definition at line 1019 of file JFitK40.hh.

1020  {
1021  int N = JK40Parameters::getN();
1022 
1023  for (int i = 0; i != pmt; ++i) {
1024  N += parameters[i].getN();
1025  }
1026 
1027  return N;
1028  }

◆ getIndex() [3/7]

int JCALIBRATE::JModel::getIndex ( int  pmt,
JParameter_t JPMTParameters_t::*  p 
) const
inline

Get index of parameter.

Parameters
pmtPMT
ppointer to data member
Returns
index

Definition at line 1038 of file JFitK40.hh.

1039  {
1040  return getIndex(pmt) + parameters[pmt].getIndex(p);
1041  }
int getIndex() const
Get index of PMT used for fixed time offset.
Definition: JFitK40.hh:966
int getIndex(JParameter_t JPMTParameters_t::*p) const
Get index of parameter.
Definition: JFitK40.hh:496

◆ getSigmaK40()

double JCALIBRATE::JModel::getSigmaK40 ( ) const
inline

Get intrinsic K40 arrival time spread.

Returns
sigma [ns]

Definition at line 1049 of file JFitK40.hh.

1050  {
1051  return this->sigmaK40_ns;
1052  }
double sigmaK40_ns
intrinsic K40 arrival time spread [ns]
Definition: JFitK40.hh:1167

◆ setSigmaK40()

void JCALIBRATE::JModel::setSigmaK40 ( const double  sigma)
inline

Set intrinsic K40 arrival time spread.

Parameters
sigmasigma [ns]

Definition at line 1060 of file JFitK40.hh.

1061  {
1062  this->sigmaK40_ns = sigma;
1063  }
const double sigma[]
Definition: JQuadrature.cc:74

◆ getReal()

const real_type& JCALIBRATE::JModel::getReal ( const pair_type pair) const
inline

Get derived parameters.

Parameters
pairPMT pair
Returns
parameters

Definition at line 1072 of file JFitK40.hh.

1073  {
1074  real.ct = JPP::getDot((*this)[pair.first].getDirection(), (*this)[pair.second].getDirection());
1075 
1076  real.t0 = (pair.first == this->index ? -this->parameters[pair.second].t0() :
1077  pair.second == this->index ? +this->parameters[pair.first ].t0() :
1078  this->parameters[pair.first].t0() - this->parameters[pair.second].t0());
1079 
1080  real.sigma = sqrt(this->parameters[pair.first ].TTS() * this->parameters[pair.first ].TTS() +
1081  this->parameters[pair.second].TTS() * this->parameters[pair.second].TTS() +
1082  this->getSigmaK40() * this->getSigmaK40());
1083 
1084  real.signal = this->parameters[pair.first].QE() * this->parameters[pair.second].QE();
1085 
1086  real.background = this->parameters[pair.first].bg() + this->parameters[pair.second].bg();
1087 
1088  return real;
1089  }
double getDot(const JNeutrinoDirection &first, const JNeutrinoDirection &second)
Dot product.
Definition: JAstronomy.hh:676
double signal
combined signal
Definition: JFitK40.hh:767
double sigma
total width [ns]
Definition: JFitK40.hh:766
double background
combined background
Definition: JFitK40.hh:768
double t0
time offset [ns]
Definition: JFitK40.hh:765
double ct
cosine angle between PMT axes
Definition: JFitK40.hh:764
real_type real
Definition: JFitK40.hh:1169
JParameter_t TTS
transition-time spread [ns]
Definition: JFitK40.hh:563

◆ getValue() [1/4]

double JCALIBRATE::JModel::getValue ( const pair_type pair,
const double  dt_ns 
) const
inline

Get K40 coincidence rate.

Parameters
pairPMT pair
dt_nstime difference [ns]
Returns
rate [Hz/ns]

Definition at line 1099 of file JFitK40.hh.

1100  {
1101  using namespace std;
1102  using namespace JPP;
1103 
1104  const real_type& real = getReal(pair);
1105 
1106  const JGauss gauss(real.t0, real.sigma, real.signal);
1107 
1108  const double R1 = this->getValue(real.ct);
1109  const double R2 = gauss.getValue(dt_ns);
1110 
1111  return real.background + R1 * (cc() + R2);
1112  }
#define R1(x)
double gauss(const double x, const double sigma)
Gauss function (normalised to 1 at x = 0).
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition: JSTDTypes.hh:14
const real_type & getReal(const pair_type &pair) const
Get derived parameters.
Definition: JFitK40.hh:1072
double getValue(const double ct) const
Get K40 coincidence rate as a function of cosine angle between PMT axes.
Definition: JFitK40.hh:711
Gauss function object.
Definition: JGauss.hh:175

◆ getValue() [2/4]

double JCALIBRATE::JModel::getValue ( const pair_type pair) const
inline

Get K40 coincidence rate.

Parameters
pairPMT pair
Returns
rate [Hz]

Definition at line 1121 of file JFitK40.hh.

1122  {
1123  using namespace std;
1124  using namespace JPP;
1125 
1126  const real_type& real = getReal(pair);
1127 
1128  const double R1 = this->getValue(real.ct);
1129  const double R2 = real.signal;
1130 
1131  return real.background + R1 * (cc() + R2);
1132  }

◆ getIndex() [4/7]

int JCALIBRATE::JK40Parameters::getIndex
inline

Get index of parameter.

Parameters
ppointer to data member
Returns
index

Definition at line 686 of file JFitK40.hh.

687  {
688  if (!(this->*p).isFree()) {
689  return INVALID_INDEX;
690  }
691 
692  int N = 0;
693 
694  if (p == &JK40Parameters::R) { return N; } if (R .isFree()) { ++N; }
695  if (p == &JK40Parameters::p1) { return N; } if (p1.isFree()) { ++N; }
696  if (p == &JK40Parameters::p2) { return N; } if (p2.isFree()) { ++N; }
697  if (p == &JK40Parameters::p3) { return N; } if (p3.isFree()) { ++N; }
698  if (p == &JK40Parameters::p4) { return N; } if (p4.isFree()) { ++N; }
699  if (p == &JK40Parameters::cc) { return N; } if (cc.isFree()) { ++N; }
700 
701  return INVALID_INDEX;
702  }
bool isFree() const
Check if parameter is free.
Definition: JFitK40.hh:240

◆ getValue() [3/4]

double JCALIBRATE::JK40Parameters::getValue
inline

Get K40 coincidence rate as a function of cosine angle between PMT axes.

Parameters
ctcosine angle between PMT axes
Returns
rate [Hz]

Definition at line 711 of file JFitK40.hh.

712  {
713  return R * exp(ct*(p1+ct*(p2+ct*(p3+ct*p4))) - (p1+p2+p3+p4));
714  }

◆ getInstance()

static const JK40Parameters& JCALIBRATE::JK40Parameters::getInstance ( )
inlinestaticinherited

Get default values.

The default parameter values are set to those obtained from a designated simulation of K40 decays (see http://wiki.km3net.de/index.php/OMGsim_simulations_for_K40_fit).
If you change these values, you may also want to change the corresponding values in JK40DefaultSimulator.hh.

Returns
parameters

Definition at line 627 of file JFitK40.hh.

628  {
629  static JK40Parameters parameters;
630 
631  parameters.R .set(18.460546);
632  parameters.p1.set( 3.0767);
633  parameters.p2.set(-1.2078);
634  parameters.p3.set( 0.9905);
635  parameters.p4.set( 0.9379);
636  parameters.cc.set( 0.0);
637 
638  return parameters;
639  }
void set()
Set current value.
Definition: JFitK40.hh:271
Fit parameters for two-fold coincidence rate due to K40.
Definition: JFitK40.hh:610

◆ getK40Parameters()

const JK40Parameters& JCALIBRATE::JK40Parameters::getK40Parameters ( ) const
inlineinherited

Get K40 parameters.

Returns
K40 parameters

Definition at line 647 of file JFitK40.hh.

648  {
649  return static_cast<const JK40Parameters&>(*this);
650  }

◆ setK40Parameters()

void JCALIBRATE::JK40Parameters::setK40Parameters ( const JK40Parameters_t parameters)
inlineinherited

Set K40 parameters.

Parameters
parametersK40 parameters

Definition at line 658 of file JFitK40.hh.

659  {
660  static_cast<JK40Parameters_t&>(*this) = parameters;
661  }
Fit parameters for two-fold coincidence rate due to K40.
Definition: JFitK40.hh:572

◆ getIndex() [5/7]

int JCALIBRATE::JK40Parameters::getIndex ( JParameter_t JK40Parameters::*  p) const
inlineinherited

Get index of parameter.

Parameters
ppointer to data member
Returns
index

Definition at line 686 of file JFitK40.hh.

687  {
688  if (!(this->*p).isFree()) {
689  return INVALID_INDEX;
690  }
691 
692  int N = 0;
693 
694  if (p == &JK40Parameters::R) { return N; } if (R .isFree()) { ++N; }
695  if (p == &JK40Parameters::p1) { return N; } if (p1.isFree()) { ++N; }
696  if (p == &JK40Parameters::p2) { return N; } if (p2.isFree()) { ++N; }
697  if (p == &JK40Parameters::p3) { return N; } if (p3.isFree()) { ++N; }
698  if (p == &JK40Parameters::p4) { return N; } if (p4.isFree()) { ++N; }
699  if (p == &JK40Parameters::cc) { return N; } if (cc.isFree()) { ++N; }
700 
701  return INVALID_INDEX;
702  }

◆ getValue() [4/4]

double JCALIBRATE::JK40Parameters::getValue ( const double  ct) const
inlineinherited

Get K40 coincidence rate as a function of cosine angle between PMT axes.

Parameters
ctcosine angle between PMT axes
Returns
rate [Hz]

Definition at line 711 of file JFitK40.hh.

712  {
713  return R * exp(ct*(p1+ct*(p2+ct*(p3+ct*p4))) - (p1+p2+p3+p4));
714  }

◆ getGradient()

const JK40Parameters_t& JCALIBRATE::JK40Parameters::getGradient ( const double  ct) const
inlineinherited

Get gradient.

Parameters
ctcosine angle between PMT axes
Returns
gradient

Definition at line 723 of file JFitK40.hh.

724  {
725  gradient.reset();
726 
727  const double rate = getValue(ct);
728  const double ct2 = ct * ct;
729 
730  if (R .isFree()) { gradient.R = rate / R; }
731  if (p1.isFree()) { gradient.p1 = rate * ct - rate; }
732  if (p2.isFree()) { gradient.p2 = rate * ct2 - rate; }
733  if (p3.isFree()) { gradient.p3 = rate * ct2 * ct - rate; }
734  if (p4.isFree()) { gradient.p4 = rate * ct2 * ct2 - rate; }
735  if (cc.isFree()) { gradient.cc = rate; }
736 
737  return gradient;
738  }
JK40Parameters_t gradient
Definition: JFitK40.hh:741
double getValue(const double ct) const
Get K40 coincidence rate as a function of cosine angle between PMT axes.
Definition: JFitK40.hh:711

◆ reset() [1/3]

void JCALIBRATE::JK40Parameters_t::reset ( )
inlineinherited

Reset.

Definition at line 585 of file JFitK40.hh.

586  {
587  R .set(0.0);
588  p1.set(0.0);
589  p2.set(0.0);
590  p3.set(0.0);
591  p4.set(0.0);
592  cc.set(0.0);
593  }

◆ getVersion()

static JDetectorVersion& JDETECTOR::JModule::getVersion ( )
inlinestaticinherited

Get detector version.

Definition at line 117 of file JModule.hh.

118  {
119  static JDetectorVersion version;
120 
121  return version;
122  }

◆ setVersion()

static void JDETECTOR::JModule::setVersion ( const JVersion version)
inlinestaticinherited

Set detector version.

Parameters
versionversion

Definition at line 130 of file JModule.hh.

131  {
132  getVersion() = JDetectorVersion(version);
133  }
static JDetectorVersion & getVersion()
Get detector version.
Definition: JModule.hh:117

◆ compare()

static bool JDETECTOR::JModule::compare ( const JModule first,
const JModule second,
const double  precision = 1.0e-3 
)
inlinestaticinherited

Compare modules.

The comparison only covers the orientations of the modules.

Parameters
firstfirst module
secondsecond module
precisionprecision
Returns
true if two modules are equal; else false

Definition at line 146 of file JModule.hh.

149  {
150  if (first.size() == second.size()) {
151 
152  for (size_t i = 0; i != first.size(); ++i) {
153  if (first[i].getDirection().getDot(second[i].getDirection()) < 1.0 - precision) {
154  return false;
155  }
156  }
157 
158  return true;
159  }
160 
161  return false;
162  }
double getDot(const JVector3D &vector) const
Get dot product.
Definition: JVector3D.hh:282
JDirection3D getDirection(const Vec &dir)
Get direction.

◆ getPMT() [1/2]

const JPMT& JDETECTOR::JModule::getPMT ( const int  index) const
inlineinherited

Get PMT.

Parameters
indexreadout channel (TDC)
Returns
PMT at given index

Definition at line 172 of file JModule.hh.

173  {
174  return at(index);
175  }

◆ getPMT() [2/2]

JPMT& JDETECTOR::JModule::getPMT ( const int  index)
inlineinherited

Get PMT.

Parameters
indexreadout channel (TDC)
Returns
PMT at given index

Definition at line 184 of file JModule.hh.

185  {
186  return at(index);
187  }

◆ setPMT()

void JDETECTOR::JModule::setPMT ( const int  index,
const JPMT pmt 
)
inlineinherited

Set PMT.

Parameters
indexreadout channel (TDC)
pmtPMT

Definition at line 196 of file JModule.hh.

197  {
198  if (index >= (int) size()) {
199  resize(index + 1);
200  }
201 
202  (*this)[index] = pmt;
203  }

◆ getCenter()

JVector3D JDETECTOR::JModule::getCenter ( ) const
inlineinherited

Get center of module based on crossing point of PMT axes.

This method perform a fit of the crossing point of the PMT axes.
A general purpose implementation of such a fit is available in JFIT::JEstimator<JPoint3D>.

Returns
center

Definition at line 214 of file JModule.hh.

215  {
216  using namespace JPP;
217 
218  if (this->size() > 1u) {
219 
220  double x = 0;
221  double y = 0;
222  double z = 0;
223 
224  JMatrix3S V;
225 
226  for (const_iterator pmt = this->begin(); pmt != this->end(); ++pmt) {
227 
228  const double xx = 1.0 - pmt->getDX() * pmt->getDX();
229  const double yy = 1.0 - pmt->getDY() * pmt->getDY();
230  const double zz = 1.0 - pmt->getDZ() * pmt->getDZ();
231 
232  const double xy = -pmt->getDX() * pmt->getDY();
233  const double xz = -pmt->getDX() * pmt->getDZ();
234  const double yz = -pmt->getDY() * pmt->getDZ();
235 
236  V.a00 += xx;
237  V.a01 += xy;
238  V.a02 += xz;
239 
240  V.a11 += yy;
241  V.a12 += yz;
242 
243  V.a22 += zz;
244 
245  x += xx * pmt->getX() + xy * pmt->getY() + xz * pmt->getZ();
246  y += xy * pmt->getX() + yy * pmt->getY() + yz * pmt->getZ();
247  z += xz * pmt->getX() + yz * pmt->getY() + zz * pmt->getZ();
248  }
249 
250  V.a10 = V.a01;
251  V.a20 = V.a02;
252  V.a21 = V.a12;
253 
254  V.invert();
255 
256  return JVector3D(V.a00 * x + V.a01 * y + V.a02 * z,
257  V.a10 * x + V.a11 * y + V.a12 * z,
258  V.a20 * x + V.a21 * y + V.a22 * z);
259 
260  } else {
261  throw JValueOutOfRange("JModule::getCenter(): Not enough PMTs.");
262  }
263  }
JVector3D()
Default constructor.
Definition: JVector3D.hh:41
3 x 3 symmetric matrix
Definition: JMatrix3S.hh:31
void invert()
Invert matrix.
Definition: JMatrix3S.hh:80

◆ compile()

void JDETECTOR::JModule::compile ( )
inlineinherited

Compile module data.

For detector versions before JDetectorVersion::V4, the position, quaternion data and time offset of the module should be set.
The position is set to the intersection point of the PMT axes (or their average position if this is not possible).
The quaternion data are maintained.
For an optical module (i.e. floor > 0), the time offset is set to the average time offset of the PMTs and for a base module (i.e. floor = 0), the time offset is set to zero.
These time offsets should be corrected for delay in the piezo sensor and hydrophone, respectively.

Definition at line 281 of file JModule.hh.

282  {
283  using namespace std;
284  using namespace JPP;
285 
286  if (!this->empty()) {
287 
288  JPosition3D& pos = this->getPosition();
289 
290  try {
291  pos = this->getCenter();
292  }
293  catch(const exception&) {
294 
295  pos = JPosition3D(0.0, 0.0, 0.0);
296 
297  for (iterator i = this->begin(); i != this->end(); ++i) {
298  pos.add(i->getPosition());
299  }
300 
301  pos.div(size());
302  }
303  }
304 
305  this->setCalibration(getAverage(make_array(this->begin(), this->end(), &JModule::getT0), 0.0));
306  }
double getT0() const
Get time offset.
void setCalibration(const JCalibration &cal)
Set calibration.
JVector3D getCenter() const
Get center of module based on crossing point of PMT axes.
Definition: JModule.hh:214
Data structure for position in three dimensions.
Definition: JPosition3D.hh:38
const JPosition3D & getPosition() const
Get position.
Definition: JPosition3D.hh:130
JPosition3D()
Default constructor.
Definition: JPosition3D.hh:48
JVector3D & add(const JVector3D &vector)
Add vector.
Definition: JVector3D.hh:142
JVector3D & div(const double factor)
Scale vector.
Definition: JVector3D.hh:190
const array_type< JValue_t > & make_array(const JValue_t(&array)[N])
Method to create array of values.
Definition: JVectorize.hh:54
std::iterator_traits< T >::value_type getAverage(T __begin, T __end)
Get average.
Definition: JMath.hh:494

◆ rotate() [1/6]

void JDETECTOR::JModule::rotate ( const JRotation3D R)
inlineinherited

Rotate module.

Parameters
Rrotation matrix

Definition at line 314 of file JModule.hh.

315  {
317 
318  for (iterator i = this->begin(); i != this->end(); ++i) {
319  i->rotate(R);
320  }
321  }
JPosition3D & rotate(const JRotation3D &R)
Rotate.
Definition: JPosition3D.hh:186

◆ rotate() [2/6]

void JDETECTOR::JModule::rotate ( const JQuaternion3D Q)
inlineinherited

Rotate module.

Parameters
Qquaternion

Definition at line 376 of file JModule.hh.

377  {
379 
380  for (iterator i = this->begin(); i != this->end(); ++i) {
381  i->rotate(Q);
382  }
383  }

◆ rotate() [3/6]

JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3X R)
inlineinherited

Rotate around X-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 214 of file JPosition3D.hh.

215  {
216  R.rotate(__y, __z);
217 
218  return *this;
219  }
void rotate(double &__x, double &__y) const
Rotate.
Definition: JRotation2D.hh:96

◆ rotate() [4/6]

JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3Y R)
inlineinherited

Rotate around Y-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 242 of file JPosition3D.hh.

243  {
244  R.rotate(__x, __z);
245 
246  return *this;
247  }

◆ rotate() [5/6]

JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JRotation3Z R)
inlineinherited

Rotate around Z-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 270 of file JPosition3D.hh.

271  {
272  R.rotate(__x, __y);
273 
274  return *this;
275  }

◆ rotate() [6/6]

void JGEOMETRY3D::JQuaternion3D::rotate ( double &  __x,
double &  __y,
double &  __z 
) const
inlineinherited

Rotate.

Parameters
__xx value
__yy value
__zz value

Definition at line 805 of file JQuaternion3D.hh.

806  {
807  const double qx = 2.0 * (__c*__z - __d*__y);
808  const double qy = 2.0 * (__d*__x - __b*__z);
809  const double qz = 2.0 * (__b*__y - __c*__x);
810 
811  __x = __x + __c*qz - __d*qy + __a*qx;
812  __y = __y - __b*qz + __a*qy + __d*qx;
813  __z = __z + __a*qz + __b*qy - __c*qx;
814  }

◆ rotate_back() [1/6]

void JDETECTOR::JModule::rotate_back ( const JRotation3D R)
inlineinherited

Rotate back module.

Parameters
Rrotation matrix

Definition at line 329 of file JModule.hh.

330  {
332 
333  for (iterator i = this->begin(); i != this->end(); ++i) {
334  i->rotate_back(R);
335  }
336  }
JPosition3D & rotate_back(const JRotation3D &R)
Rotate back.
Definition: JPosition3D.hh:200

◆ rotate_back() [2/6]

void JDETECTOR::JModule::rotate_back ( const JQuaternion3D Q)
inlineinherited

Rotate back module.

Parameters
Qquaternion

Definition at line 391 of file JModule.hh.

392  {
394 
395  for (iterator i = this->begin(); i != this->end(); ++i) {
396  i->rotate_back(Q);
397  }
398  }

◆ rotate_back() [3/6]

JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3X R)
inlineinherited

Rotate back around X-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 228 of file JPosition3D.hh.

229  {
230  R.rotate_back(__y, __z);
231 
232  return *this;
233  }
void rotate_back(double &__x, double &__y) const
Rotate back.
Definition: JRotation2D.hh:112

◆ rotate_back() [4/6]

JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3Y R)
inlineinherited

Rotate back around Y-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 256 of file JPosition3D.hh.

257  {
258  R.rotate_back(__x, __z);
259 
260  return *this;
261  }

◆ rotate_back() [5/6]

JPosition3D& JGEOMETRY3D::JPosition3D::rotate_back ( const JRotation3Z R)
inlineinherited

Rotate back around Z-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 284 of file JPosition3D.hh.

285  {
286  R.rotate_back(__x, __y);
287 
288  return *this;
289  }

◆ rotate_back() [6/6]

void JGEOMETRY3D::JQuaternion3D::rotate_back ( double &  __x,
double &  __y,
double &  __z 
) const
inlineinherited

Rotate back.

Parameters
__xx value
__yy value
__zz value

Definition at line 824 of file JQuaternion3D.hh.

825  {
826  const double qx = 2.0 * (__d*__y - __c*__z);
827  const double qy = 2.0 * (__b*__z - __d*__x);
828  const double qz = 2.0 * (__c*__x - __b*__y);
829 
830  __x = __x - __c*qz + __d*qy + __a*qx;
831  __y = __y + __b*qz + __a*qy - __d*qx;
832  __z = __z + __a*qz - __b*qy + __c*qx;
833  }

◆ transform() [1/4]

void JDETECTOR::JModule::transform ( const JRotation3D R,
const JVector3D pos 
)
inlineinherited

Transformation of geometry (see method JGEOMETRY3D::JPosition3D::transform(const JRotation3D&, const JVector3D&)).

Parameters
Rrotation matrix
posposition of origin (after rotation)

Definition at line 345 of file JModule.hh.

347  {
348  JPosition3D::transform(R, pos);
349 
350  for (iterator i = this->begin(); i != this->end(); ++i) {
351  i->transform(R, pos);
352  }
353  }
JVector3D & transform(const JMatrix3D &T)
Transform.
Definition: JVector3D.hh:206

◆ transform() [2/4]

void JDETECTOR::JModule::transform ( const JTransformation3D T)
inlineinherited

Transformation of geometry.

Parameters
Ttransformation

Definition at line 361 of file JModule.hh.

362  {
364 
365  for (iterator i = this->begin(); i != this->end(); ++i) {
366  i->transform(T);
367  }
368  }
const JRotation3D & getRotation() const
Get rotation.
Definition: JRotation3D.hh:272

◆ transform() [3/4]

JVector3D& JGEOMETRY3D::JVector3D::transform
inlineinherited

Transform.

Parameters
Tmatrix
Returns
this vector

Definition at line 206 of file JVector3D.hh.

207  {
208  T.transform(__x, __y, __z);
209 
210  return *this;
211  }

◆ transform() [4/4]

JVector3D& JGEOMETRY3D::JVector3D::transform ( const JMatrix3D T)
inlineinherited

Transform.

Parameters
Tmatrix
Returns
this vector

Definition at line 206 of file JVector3D.hh.

207  {
208  T.transform(__x, __y, __z);
209 
210  return *this;
211  }

◆ set() [1/4]

JModule& JDETECTOR::JModule::set ( const JVector3D pos)
inlineinherited

Set position.

Parameters
posposition
Returns
this module

Definition at line 407 of file JModule.hh.

408  {
409  return add(pos - this->getPosition());
410  }
JModule & add(const JVector3D &pos)
Add position.
Definition: JModule.hh:419

◆ set() [2/4]

JModule& JDETECTOR::JModule::set ( const double  t0)
inlineinherited

Set time offset.

Parameters
t0time offset [ns]
Returns
this module

Definition at line 455 of file JModule.hh.

456  {
457  for (iterator i = begin(); i != end(); ++i) {
458  i->setT0(t0);
459  }
460 
461  return *this;
462  }

◆ set() [3/4]

void JEEP::JStatus::set
inlineinherited

Set PMT status.

Parameters
bitbit

Definition at line 131 of file JStatus.hh.

132  {
133  this->status |= (1<<bit);
134  }

◆ set() [4/4]

void JEEP::JStatus::set ( const int  bit)
inlineinherited

Set PMT status.

Parameters
bitbit

Definition at line 131 of file JStatus.hh.

132  {
133  this->status |= (1<<bit);
134  }

◆ add() [1/3]

JModule& JDETECTOR::JModule::add ( const JVector3D pos)
inlineinherited

Add position.

Parameters
posposition
Returns
this module

Definition at line 419 of file JModule.hh.

420  {
421  for (iterator i = begin(); i != end(); ++i) {
422  i->add(pos);
423  }
424 
425  JPosition3D::add(pos);
426 
427  return *this;
428  }

◆ add() [2/3]

JModule& JDETECTOR::JModule::add ( const double  t0)
inlineinherited

Add time offset.

Parameters
t0time offset [ns]
Returns
this module

Definition at line 471 of file JModule.hh.

472  {
473  for (iterator i = begin(); i != end(); ++i) {
474  i->addT0(t0);
475  }
476 
477  return *this;
478  }

◆ add() [3/3]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::add ( const JQuaternion3D quaternion)
inlineinherited

Add quaternion.

Parameters
quaternionquaternion
Returns
this quaternion

Definition at line 637 of file JQuaternion3D.hh.

638  {
639  __a += quaternion.getA();
640  __b += quaternion.getB();
641  __c += quaternion.getC();
642  __d += quaternion.getD();
643 
644  return *this;
645  }
double getB() const
Get b value.
double getD() const
Get d value.
double getC() const
Get c value.
double getA() const
Get a value.

◆ sub() [1/3]

JModule& JDETECTOR::JModule::sub ( const JVector3D pos)
inlineinherited

Subtract position.

Parameters
posposition
Returns
this module

Definition at line 437 of file JModule.hh.

438  {
439  for (iterator i = begin(); i != end(); ++i) {
440  i->sub(pos);
441  }
442 
443  JPosition3D::sub(pos);
444 
445  return *this;
446  }
JVector3D & sub(const JVector3D &vector)
Subtract vector.
Definition: JVector3D.hh:158

◆ sub() [2/3]

JModule& JDETECTOR::JModule::sub ( const double  t0)
inlineinherited

Subtract time offset.

Parameters
t0time offset [ns]
Returns
this module

Definition at line 487 of file JModule.hh.

488  {
489  for (iterator i = begin(); i != end(); ++i) {
490  i->subT0(t0);
491  }
492 
493  return *this;
494  }

◆ sub() [3/3]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::sub ( const JQuaternion3D quaternion)
inlineinherited

Subtract quaternion.

Parameters
quaternionquaternion
Returns
this quaternion

Definition at line 654 of file JQuaternion3D.hh.

655  {
656  __a -= quaternion.getA();
657  __b -= quaternion.getB();
658  __c -= quaternion.getC();
659  __d -= quaternion.getD();
660 
661  return *this;
662  }

◆ operator+=()

JModule& JDETECTOR::JModule::operator+= ( const JVector3D pos)
inlineinherited

Add position.

Parameters
posposition
Returns
this module

Definition at line 503 of file JModule.hh.

504  {
505  return this->add(pos);
506  }

◆ operator-=()

JModule& JDETECTOR::JModule::operator-= ( const JVector3D pos)
inlineinherited

Subtract position.

Parameters
posposition
Returns
this module

Definition at line 515 of file JModule.hh.

516  {
517  return this->sub(pos);
518  }
JModule & sub(const JVector3D &pos)
Subtract position.
Definition: JModule.hh:437

◆ has() [1/2]

bool JEEP::JStatus::has
inlineinherited

Test PMT status.

Parameters
bitbit

Definition at line 120 of file JStatus.hh.

121  {
122  return (this->status & (1<<bit)) != 0;
123  }

◆ has() [2/2]

bool JEEP::JStatus::has ( const int  bit) const
inlineinherited

Test PMT status.

Parameters
bitbit

Definition at line 120 of file JStatus.hh.

121  {
122  return (this->status & (1<<bit)) != 0;
123  }

◆ reset() [2/3]

void JEEP::JStatus::reset
inlineinherited

Reset PMT status.

Parameters
bitbit

Definition at line 142 of file JStatus.hh.

143  {
144  this->status &= ~(1<<bit);
145  }

◆ reset() [3/3]

void JEEP::JStatus::reset ( const int  bit)
inlineinherited

Reset PMT status.

Parameters
bitbit

Definition at line 142 of file JStatus.hh.

143  {
144  this->status &= ~(1<<bit);
145  }

◆ getID() [1/2]

int JLANG::JObjectID::getID ( ) const
inlineinherited

Get identifier.

Returns
identifier

Definition at line 50 of file JObjectID.hh.

51  {
52  return __id;
53  }

◆ getID() [2/2]

int& JLANG::JObjectID::getID ( )
inlineinherited

Get identifier.

Returns
identifier

Definition at line 61 of file JObjectID.hh.

62  {
63  return __id;
64  }

◆ setID()

void JLANG::JObjectID::setID ( const int  id)
inlineinherited

Set identifier.

Parameters
ididentifier

Definition at line 72 of file JObjectID.hh.

73  {
74  this->__id = id;
75  }

◆ less() [1/3]

bool JLANG::JObjectID::less ( const JObjectID object) const
inlineinherited

Less than method.

Parameters
objectobject identifier
Returns
true if this identifier less than given identifier; else false

Definition at line 84 of file JObjectID.hh.

85  {
86  return this->getID() < object.getID();
87  }
int getID() const
Get identifier.
Definition: JObjectID.hh:50

◆ less() [2/3]

bool JLANG::JObjectID::less ( const int  id) const
inlineinherited

Less than method.

Parameters
ididentifier
Returns
true if this identifier less than given identifier; else false

Definition at line 96 of file JObjectID.hh.

97  {
98  return this->getID() < id;
99  }

◆ less() [3/3]

bool JDETECTOR::JLocation::less ( const JLocation location) const
inlineinherited

Less than method.

Parameters
locationmodule location
Returns
true if first location before second location; else false

Definition at line 158 of file JLocation.hh.

159  {
160  if (this->getString() == location.getString())
161  return this->getFloor() < location.getFloor();
162  else
163  return this->getString() < location.getString();
164  }
int getFloor() const
Get floor number.
Definition: JLocation.hh:146
int getString() const
Get string number.
Definition: JLocation.hh:135

◆ more()

bool JLANG::JObjectID::more ( const int  id) const
inlineinherited

More than method.

Parameters
ididentifier
Returns
true if this identifier greater than given identifier; else false

Definition at line 108 of file JObjectID.hh.

109  {
110  return this->getID() > id;
111  }

◆ getLocation() [1/2]

const JLocation& JDETECTOR::JLocation::getLocation ( ) const
inlineinherited

Get location.

Returns
location

Definition at line 70 of file JLocation.hh.

71  {
72  return static_cast<const JLocation&>(*this);
73  }
Logical location of module.
Definition: JLocation.hh:40

◆ getLocation() [2/2]

JLocation& JDETECTOR::JLocation::getLocation ( )
inlineinherited

Get location.

Returns
location

Definition at line 81 of file JLocation.hh.

82  {
83  return static_cast<JLocation&>(*this);
84  }

◆ setLocation()

void JDETECTOR::JLocation::setLocation ( const JLocation location)
inlineinherited

Set location.

Parameters
locationlocation

Definition at line 92 of file JLocation.hh.

93  {
94  static_cast<JLocation&>(*this) = location;
95  }

◆ toString() [1/2]

std::string JDETECTOR::JLocation::toString ( ) const
inlineinherited

Convert module location to string.

Returns
string

Definition at line 103 of file JLocation.hh.

104  {
105  return toString("% %");
106  }
std::string toString() const
Convert module location to string.
Definition: JLocation.hh:103

◆ toString() [2/2]

std::string JDETECTOR::JLocation::toString ( const std::string &  fmt,
const std::string  target = "%" 
) const
inlineinherited

Convert module loation to string.

The targets target in the format string fmt are consecutively replaced by floor and string.

Parameters
fmtformat
targettarget
Returns
string

Definition at line 119 of file JLocation.hh.

120  {
121  JLANG::JString buffer(fmt);
122 
123  buffer.replace(target, string, 1);
124  buffer.replace(target, floor, 1);
125 
126  return buffer;
127  }
Wrapper class around STL string class.
Definition: JString.hh:29
Target.
Definition: JHead.hh:300

◆ getString()

int JDETECTOR::JLocation::getString ( ) const
inlineinherited

Get string number.

Returns
string number

Definition at line 135 of file JLocation.hh.

136  {
137  return string;
138  }

◆ getFloor()

int JDETECTOR::JLocation::getFloor ( ) const
inlineinherited

Get floor number.

Returns
floor number

Definition at line 146 of file JLocation.hh.

147  {
148  return floor;
149  }

◆ getPosition() [1/2]

const JPosition3D& JGEOMETRY3D::JPosition3D::getPosition ( ) const
inlineinherited

Get position.

Returns
position

Definition at line 130 of file JPosition3D.hh.

131  {
132  return static_cast<const JPosition3D&>(*this);
133  }

◆ getPosition() [2/2]

JPosition3D& JGEOMETRY3D::JPosition3D::getPosition ( )
inlineinherited

Get position.

Returns
position

Definition at line 141 of file JPosition3D.hh.

142  {
143  return static_cast<JPosition3D&>(*this);
144  }

◆ setPosition()

void JGEOMETRY3D::JPosition3D::setPosition ( const JVector3D pos)
inlineinherited

Set position.

Parameters
posposition

Definition at line 152 of file JPosition3D.hh.

153  {
154  static_cast<JVector3D&>(*this) = pos;
155  }
Data structure for vector in three dimensions.
Definition: JVector3D.hh:36

◆ operator JAngle3D()

JGEOMETRY3D::JPosition3D::operator JAngle3D ( ) const
inlineinherited

Type conversion operator.

Returns
angle

Definition at line 163 of file JPosition3D.hh.

164  {
165  return JAngle3D(getX(), getY(), getZ());
166  }
Data structure for angles in three dimensions.
Definition: JAngle3D.hh:35
double getY() const
Get y position.
Definition: JVector3D.hh:104
double getZ() const
Get z position.
Definition: JVector3D.hh:115
double getX() const
Get x position.
Definition: JVector3D.hh:94

◆ operator JVersor3D() [1/2]

JGEOMETRY3D::JPosition3D::operator JVersor3D ( ) const
inlineinherited

Type conversion operator.

Returns
direction

Definition at line 174 of file JPosition3D.hh.

175  {
176  return JVersor3D(getX(), getY(), getZ());
177  }
Data structure for normalised vector in three dimensions.
Definition: JVersor3D.hh:28

◆ operator JVersor3D() [2/2]

JGEOMETRY3D::JQuaternion3D::operator JVersor3D ( ) const
inlineinherited

Type conversion operator.

Returns
direction

Definition at line 520 of file JQuaternion3D.hh.

521  {
522  return JVersor3D(this->getB(), this->getC(), this->getD());
523  }

◆ transform_back()

void JGEOMETRY3D::JPosition3D::transform_back ( const JRotation3D R,
const JVector3D pos 
)
inlineinherited

Transform back position.

The final position is obtained as follows:

  1. offset position with position pos;
  2. rotation of postion according matrix R;
Parameters
Rrotation matrix
posposition of origin (before rotation)

Definition at line 359 of file JPosition3D.hh.

361  {
362  // offset with respect to origin
363 
364  add(pos);
365 
366  // rotate back geometry to system with particle direction along z-axis
367 
368  rotate_back(R);
369  }

◆ getDot() [1/6]

double JGEOMETRY3D::JPosition3D::getDot ( const JAngle3D angle) const
inlineinherited

Get dot product.

Parameters
angleangle
Returns
dot product

Definition at line 378 of file JPosition3D.hh.

379  {
380  return
381  getX() * angle.getDX() +
382  getY() * angle.getDY() +
383  getZ() * angle.getDZ();
384  }
double getDY() const
Get y direction.
Definition: JAngle3D.hh:119
double getDZ() const
Get z direction.
Definition: JAngle3D.hh:130
double getDX() const
Get x direction.
Definition: JAngle3D.hh:108

◆ getDot() [2/6]

double JGEOMETRY3D::JPosition3D::getDot ( const JVersor3D dir) const
inlineinherited

Get dot product.

Parameters
dirdirection
Returns
dot product

Definition at line 393 of file JPosition3D.hh.

394  {
395  return
396  getX() * dir.getDX() +
397  getY() * dir.getDY() +
398  getZ() * dir.getDZ();
399  }
double getDY() const
Get y direction.
Definition: JVersor3D.hh:106
double getDX() const
Get x direction.
Definition: JVersor3D.hh:95
double getDZ() const
Get z direction.
Definition: JVersor3D.hh:117

◆ getDot() [3/6]

double JGEOMETRY3D::JPosition3D::getDot ( const JVersor3Z dir) const
inlineinherited

Get dot product.

Parameters
dirdirection
Returns
dot product

Definition at line 408 of file JPosition3D.hh.

409  {
410  return
411  getX() * dir.getDX() +
412  getY() * dir.getDY() +
413  getZ() * dir.getDZ();
414  }
double getDZ() const
Get z direction.
Definition: JVersor3Z.hh:169
double getDY() const
Get y direction.
Definition: JVersor3Z.hh:158
double getDX() const
Get x direction.
Definition: JVersor3Z.hh:147

◆ getDot() [4/6]

double JGEOMETRY3D::JVector3D::getDot
inlineinherited

Get dot product.

Parameters
vectorvector
Returns
dot product

Definition at line 282 of file JVector3D.hh.

283  {
284  return
285  getX() * vector.getX() +
286  getY() * vector.getY() +
287  getZ() * vector.getZ();
288  }

◆ getDot() [5/6]

double JGEOMETRY3D::JVector3D::getDot ( const JVector3D vector) const
inlineinherited

Get dot product.

Parameters
vectorvector
Returns
dot product

Definition at line 282 of file JVector3D.hh.

283  {
284  return
285  getX() * vector.getX() +
286  getY() * vector.getY() +
287  getZ() * vector.getZ();
288  }

◆ getDot() [6/6]

double JGEOMETRY3D::JQuaternion3D::getDot ( const JQuaternion3D quaternion) const
inlineinherited

Get dot product.

Parameters
quaternionquaternion
Returns
dot product

Definition at line 933 of file JQuaternion3D.hh.

934  {
935  return
936  getA() * quaternion.getA() -
937  getB() * quaternion.getB() -
938  getC() * quaternion.getC() -
939  getD() * quaternion.getD();
940  }

◆ operator JVector2D()

JGEOMETRY3D::JVector3D::operator JVector2D ( ) const
inlineinherited

Type conversion operator.

Returns
JVector2D

Definition at line 83 of file JVector3D.hh.

84  {
85  return JVector2D(this->getX(), this->getY());
86  }

◆ getX()

double JGEOMETRY3D::JVector3D::getX ( ) const
inlineinherited

Get x position.

Returns
x position

Definition at line 94 of file JVector3D.hh.

95  {
96  return __x;
97  }

◆ getY()

double JGEOMETRY3D::JVector3D::getY ( ) const
inlineinherited

Get y position.

Returns
y position

Definition at line 104 of file JVector3D.hh.

105  {
106  return __y;
107  }

◆ getZ()

double JGEOMETRY3D::JVector3D::getZ ( ) const
inlineinherited

Get z position.

Returns
z position

Definition at line 115 of file JVector3D.hh.

116  {
117  return __z;
118  }

◆ negate() [1/2]

JVector3D& JGEOMETRY3D::JVector3D::negate ( )
inlineinherited

Negate vector.

Returns
this vector

Definition at line 126 of file JVector3D.hh.

127  {
128  __x = -__x;
129  __y = -__y;
130  __z = -__z;
131 
132  return *this;
133  }

◆ negate() [2/2]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::negate ( )
inlineinherited

Negate quaternion.

Returns
this quaternion

Definition at line 620 of file JQuaternion3D.hh.

621  {
622  __a = -__a;
623  __b = -__b;
624  __c = -__c;
625  __d = -__d;
626 
627  return *this;
628  }

◆ mul() [1/9]

JVector3D& JGEOMETRY3D::JVector3D::mul ( const double  factor)
inlineinherited

Scale vector.

Parameters
factormultiplication factor
Returns
this vector

Definition at line 174 of file JVector3D.hh.

175  {
176  __x *= factor;
177  __y *= factor;
178  __z *= factor;
179 
180  return *this;
181  }

◆ mul() [2/9]

template<class JFirst_t , class JSecond_t >
JFirst_t& JMATH::JMath< JFirst_t, JSecond_t >::mul ( const JSecond_t &  object)
inlineinherited

Multiply with object.

Parameters
objectobject
Returns
result object

Definition at line 354 of file JMath.hh.

355  {
356  return static_cast<JFirst_t&>(*this) = JFirst_t().mul(static_cast<const JFirst_t&>(*this), object);
357  }

◆ mul() [3/9]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::mul ( const double  factor)
inlineinherited

Scale quaternion.

Parameters
factormultiplication factor
Returns
this quaternion

Definition at line 671 of file JQuaternion3D.hh.

672  {
673  __a *= factor;
674  __b *= factor;
675  __c *= factor;
676  __d *= factor;
677 
678  return *this;
679  }

◆ mul() [4/9]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::mul ( const JQuaternion3D first,
const JQuaternion3X second 
)
inlineinherited

Quaternion multiplication.

This method evaluates the Hamilton product (also called cross product).

Parameters
firstfirst quaternion
secondsecond quaternion
Returns
this quaternion

Definition at line 708 of file JQuaternion3D.hh.

709  {
710  __a = first.getA() * second.getA() - first.getB() * second.getB();
711  __b = first.getA() * second.getB() + first.getB() * second.getA();
712  __c = first.getC() * second.getA() + first.getD() * second.getB();
713  __d = -first.getC() * second.getB() + first.getD() * second.getA();
714 
715  return *this;
716  }
double getA() const
Get a value.
double getB() const
Get b value.

◆ mul() [5/9]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::mul ( const JQuaternion3D first,
const JQuaternion3Y second 
)
inlineinherited

Quaternion multiplication.

This method evaluates the Hamilton product (or cross product).

Parameters
firstfirst quaternion
secondsecond quaternion
Returns
this quaternion

Definition at line 728 of file JQuaternion3D.hh.

730  {
731  __a = first.getA() * second.getA() - first.getC() * second.getC();
732  __b = first.getB() * second.getA() - first.getD() * second.getC();
733  __c = first.getA() * second.getC() + first.getC() * second.getA();
734  __d = first.getB() * second.getC() + first.getD() * second.getA();
735 
736  return *this;
737  }
double getC() const
Get c value.

◆ mul() [6/9]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::mul ( const JQuaternion3D first,
const JQuaternion3Z second 
)
inlineinherited

Quaternion multiplication.

This method evaluates the Hamilton product (or cross product).

Parameters
firstfirst quaternion
secondsecond quaternion
Returns
this quaternion

Definition at line 749 of file JQuaternion3D.hh.

751  {
752  __a = first.getA() * second.getA() - first.getD() * second.getD();
753  __b = first.getB() * second.getA() + first.getC() * second.getD();
754  __c = -first.getB() * second.getD() + first.getC() * second.getA();
755  __d = first.getA() * second.getD() + first.getD() * second.getA();
756 
757  return *this;
758  }
double getD() const
Get d value.

◆ mul() [7/9]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::mul ( const JQuaternion3D first,
const JQuaternion3D second 
)
inlineinherited

Quaternion multiplication.

This method evaluates the Hamilton product (or cross product).

Parameters
firstfirst quaternion
secondsecond quaternion
Returns
this quaternion

Definition at line 770 of file JQuaternion3D.hh.

772  {
773  __a = first.getA() * second.getA() - first.getB() * second.getB() - first.getC() * second.getC() - first.getD() * second.getD();
774  __b = first.getA() * second.getB() + first.getB() * second.getA() + first.getC() * second.getD() - first.getD() * second.getC();
775  __c = first.getA() * second.getC() - first.getB() * second.getD() + first.getC() * second.getA() + first.getD() * second.getB();
776  __d = first.getA() * second.getD() + first.getB() * second.getC() - first.getC() * second.getB() + first.getD() * second.getA();
777 
778  return *this;
779  }

◆ mul() [8/9]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::mul ( const JQuaternion3X qx,
const JQuaternion3Y qy,
const JQuaternion3Z qz 
)
inlineinherited

Quaternion multiplication.

Parameters
qxrotation around x-axis
qyrotation around y-axis
qzrotation around x-axis
Returns
this quaternion

Definition at line 790 of file JQuaternion3D.hh.

793  {
794  return *this = JQuaternion3D().setIdentity().mul(qz).mul(qy).mul(qx);
795  }
JQuaternion3D()
Default constructor.

◆ mul() [9/9]

template<class JFirst_t , class JSecond_t >
JFirst_t& JMATH::JMath< JFirst_t, JSecond_t >::mul ( const JSecond_t &  object)
inlineinherited

Multiply with object.

Parameters
objectobject
Returns
result object

Definition at line 354 of file JMath.hh.

355  {
356  return static_cast<JFirst_t&>(*this) = JFirst_t().mul(static_cast<const JFirst_t&>(*this), object);
357  }

◆ div() [1/2]

JVector3D& JGEOMETRY3D::JVector3D::div ( const double  factor)
inlineinherited

Scale vector.

Parameters
factordivision factor
Returns
this vector

Definition at line 190 of file JVector3D.hh.

191  {
192  __x /= factor;
193  __y /= factor;
194  __z /= factor;
195 
196  return *this;
197  }

◆ div() [2/2]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::div ( const double  factor)
inlineinherited

Scale quaternion.

Parameters
factordivision factor
Returns
this quaternion

Definition at line 688 of file JQuaternion3D.hh.

689  {
690  __a /= factor;
691  __b /= factor;
692  __c /= factor;
693  __d /= factor;
694 
695  return *this;
696  }

◆ equals() [1/3]

bool JGEOMETRY3D::JVector3D::equals ( const JVector3D vector,
const double  precision = std::numeric_limits<double>::min() 
) const
inlineinherited

Check equality.

Parameters
vectorvector
precisionprecision
Returns
true if vectors are equal; else false

Definition at line 221 of file JVector3D.hh.

223  {
224  return (fabs(getX() - vector.getX()) <= precision &&
225  fabs(getY() - vector.getY()) <= precision &&
226  fabs(getZ() - vector.getZ()) <= precision);
227  }

◆ equals() [2/3]

bool JGEOMETRY3D::JQuaternion3D::equals ( const JQuaternion3D quaternion,
const double  precision = std::numeric_limits<double>::min() 
) const
inlineinherited

Check equality.

Parameters
quaternionquaternion
precisionnumerical precision
Returns
true if quaternions are equal; else false

Definition at line 843 of file JQuaternion3D.hh.

845  {
846  return (fabs(getA() - quaternion.getA()) <= precision &&
847  fabs(getB() - quaternion.getB()) <= precision &&
848  fabs(getC() - quaternion.getC()) <= precision &&
849  fabs(getD() - quaternion.getD()) <= precision);
850  }

◆ equals() [3/3]

bool JEEP::JStatus::equals ( const JStatus status) const
inlineinherited

Equal method.

Parameters
statusstatus
Returns
true if this status equal to given status; else false

Definition at line 109 of file JStatus.hh.

110  {
111  return this->status == status.status;
112  }

◆ getLengthSquared() [1/2]

double JGEOMETRY3D::JVector3D::getLengthSquared ( ) const
inlineinherited

Get length squared.

Returns
square of length

Definition at line 235 of file JVector3D.hh.

236  {
237  return getX()*getX() + getY()*getY() + getZ()*getZ();
238  }

◆ getLengthSquared() [2/2]

double JGEOMETRY3D::JQuaternion3D::getLengthSquared ( ) const
inlineinherited

Get length squared.

Returns
square of length

Definition at line 886 of file JQuaternion3D.hh.

887  {
888  return getA()*getA() + getB()*getB() + getC()*getC() + getD()*getD();
889  }

◆ getLength() [1/2]

double JGEOMETRY3D::JVector3D::getLength ( ) const
inlineinherited

Get length.

Returns
length

Definition at line 246 of file JVector3D.hh.

247  {
248  return sqrt(getLengthSquared());
249  }
double getLengthSquared() const
Get length squared.
Definition: JVector3D.hh:235

◆ getLength() [2/2]

double JGEOMETRY3D::JQuaternion3D::getLength ( ) const
inlineinherited

Get length.

Returns
length

Definition at line 897 of file JQuaternion3D.hh.

898  {
899  return sqrt(getLengthSquared());
900  }
double getLengthSquared() const
Get length squared.

◆ getDistanceSquared() [1/2]

double JGEOMETRY3D::JVector3D::getDistanceSquared ( const JVector3D pos) const
inlineinherited

Get squared of distance to point.

Parameters
posposition
Returns
square of distance

Definition at line 258 of file JVector3D.hh.

259  {
260  return JVector3D(pos).sub(*this).getLengthSquared();
261  }

◆ getDistanceSquared() [2/2]

double JGEOMETRY3D::JQuaternion3D::getDistanceSquared ( const JQuaternion3D quaternion) const
inlineinherited

Get squared of distance to quaternion.

Parameters
quaternionquaternion
Returns
square of distance

Definition at line 909 of file JQuaternion3D.hh.

910  {
911  return JQuaternion3D(quaternion).sub(*this).getLengthSquared();
912  }

◆ getDistance() [1/2]

double JGEOMETRY3D::JVector3D::getDistance ( const JVector3D pos) const
inlineinherited

Get distance to point.

Parameters
posposition
Returns
distance

Definition at line 270 of file JVector3D.hh.

271  {
272  return sqrt(getDistanceSquared(pos));
273  }
double getDistanceSquared(const JVector3D &pos) const
Get squared of distance to point.
Definition: JVector3D.hh:258

◆ getDistance() [2/2]

double JGEOMETRY3D::JQuaternion3D::getDistance ( const JQuaternion3D quaternion) const
inlineinherited

Get distance to quaternion.

Parameters
quaternionquaternion
Returns
distance

Definition at line 921 of file JQuaternion3D.hh.

922  {
923  return sqrt(getDistanceSquared(quaternion));
924  }
double getDistanceSquared(const JQuaternion3D &quaternion) const
Get squared of distance to quaternion.

◆ getCross()

JVector3D& JGEOMETRY3D::JVector3D::getCross ( const JVector3D first,
const JVector3D second 
)
inlineinherited

Get cross product.

Note that this vector should not overlap with the first or second vector,

Parameters
firstfirst vector
secondsecond vector
Returns
this vector

Definition at line 299 of file JVector3D.hh.

301  {
302  __x = first .getY() * second.getZ() - second.getY() * first .getZ();
303  __y = second.getX() * first .getZ() - first .getX() * second.getZ();
304  __z = first .getX() * second.getY() - second.getX() * first .getY();
305 
306  return *this;
307  }

◆ getIdentity()

static const JQuaternion3D& JGEOMETRY3D::JQuaternion3D::getIdentity ( )
inlinestaticinherited

Get identity quaternion.

Returns
identity quaternion

Definition at line 463 of file JQuaternion3D.hh.

464  {
465  static const JQuaternion3D Q = JQuaternion3D().setIdentity();
466 
467  return Q;
468  }
Data structure for unit quaternion in three dimensions.

◆ getQuaternion() [1/2]

const JQuaternion3D& JGEOMETRY3D::JQuaternion3D::getQuaternion ( ) const
inlineinherited

Get quaternion.

Returns
quaternion

Definition at line 476 of file JQuaternion3D.hh.

477  {
478  return static_cast<const JQuaternion3D&>(*this);
479  }

◆ getQuaternion() [2/2]

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::getQuaternion ( )
inlineinherited

Get quaternion.

Returns
quaternion

Definition at line 487 of file JQuaternion3D.hh.

488  {
489  return static_cast<JQuaternion3D&>(*this);
490  }

◆ setQuaternion()

void JGEOMETRY3D::JQuaternion3D::setQuaternion ( const JQuaternion3D quaternion)
inlineinherited

Set quaternion.

Parameters
quaternionquaternion

Definition at line 498 of file JQuaternion3D.hh.

499  {
500  static_cast<JQuaternion3D&>(*this) = quaternion;
501  }

◆ operator JVector3D()

JGEOMETRY3D::JQuaternion3D::operator JVector3D ( ) const
inlineinherited

Type conversion operator.

Returns
position

Definition at line 509 of file JQuaternion3D.hh.

510  {
511  return JVector3D(this->getB(), this->getC(), this->getD());
512  }

◆ getAngle()

double JGEOMETRY3D::JQuaternion3D::getAngle ( ) const
inlineinherited

Get rotation angle.

Returns
angle [rad]

Definition at line 531 of file JQuaternion3D.hh.

532  {
533  return atan2(sqrt(this->getB()*this->getB() +
534  this->getC()*this->getC() +
535  this->getD()*this->getD()),
536  this->getA()) * 2.0;
537  }

◆ getA()

double JGEOMETRY3D::JQuaternion3D::getA ( ) const
inlineinherited

Get a value.

Returns
a value

Definition at line 545 of file JQuaternion3D.hh.

546  {
547  return __a;
548  }

◆ getB()

double JGEOMETRY3D::JQuaternion3D::getB ( ) const
inlineinherited

Get b value.

Returns
b value

Definition at line 556 of file JQuaternion3D.hh.

557  {
558  return __b;
559  }

◆ getC()

double JGEOMETRY3D::JQuaternion3D::getC ( ) const
inlineinherited

Get c value.

Returns
c value

Definition at line 567 of file JQuaternion3D.hh.

568  {
569  return __c;
570  }

◆ getD()

double JGEOMETRY3D::JQuaternion3D::getD ( ) const
inlineinherited

Get d value.

Returns
d value

Definition at line 578 of file JQuaternion3D.hh.

579  {
580  return __d;
581  }

◆ setIdentity()

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::setIdentity ( )
inlineinherited

Set to identity quaternion.

Returns
this quaternion

Definition at line 589 of file JQuaternion3D.hh.

590  {
591  __a = 1.0;
592  __b = 0.0;
593  __c = 0.0;
594  __d = 0.0;
595 
596  return *this;
597  }

◆ conjugate()

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::conjugate ( )
inlineinherited

Conjugate quaternion.

Returns
this quaternion

Definition at line 605 of file JQuaternion3D.hh.

606  {
607  __b = -__b;
608  __c = -__c;
609  __d = -__d;
610 
611  return *this;
612  }

◆ isIdentity()

bool JGEOMETRY3D::JQuaternion3D::isIdentity ( const double  precision = std::numeric_limits<double>::min()) const
inlineinherited

Test identity.

Parameters
precisionprecision
Returns
true if identity quaternion; else false

Definition at line 859 of file JQuaternion3D.hh.

860  {
861  if (fabs(getA()) <= precision) {
862 
863  if (fabs(getB()) <= precision)
864  return ((fabs(getC()) <= precision && fabs(fabs(getD()) - 1.0) <= precision) ||
865  (fabs(fabs(getC()) - 1.0) <= precision && fabs(getD()) <= precision));
866  else
867  return (fabs(fabs(getB()) - 1.0) <= precision &&
868  fabs(getC()) <= precision &&
869  fabs(getD()) <= precision);
870 
871  } else {
872 
873  return (fabs(fabs(getA()) - 1.0) <= precision &&
874  fabs(getB()) <= precision &&
875  fabs(getC()) <= precision &&
876  fabs(getD()) <= precision);
877  }
878  }

◆ getConjugate()

JQuaternion3D JGEOMETRY3D::JQuaternion3D::getConjugate ( ) const
inlineinherited

Get conjugate of this quaternion.

Returns
conjugate quaternion

Definition at line 948 of file JQuaternion3D.hh.

949  {
950  return JQuaternion3D(*this).conjugate();
951  }

◆ normalise()

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::normalise ( )
inlineinherited

Normalise quaternion.

Returns
this quaternion

Definition at line 959 of file JQuaternion3D.hh.

960  {
961  const double v = getLength();
962 
963  if (v != 0.0) {
964  mul(1.0 / v);
965  }
966 
967  return *this;
968  }
double getLength() const
Get length.
JQuaternion3D & mul(const double factor)
Scale quaternion.
data_type v[N+1][M+1]
Definition: JPolint.hh:866

◆ pow()

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::pow ( const double  y)
inlineinherited

Raise quaternion to given power.

Parameters
ypower
Returns
this object

Definition at line 977 of file JQuaternion3D.hh.

978  {
979  const double v = sqrt(getB() * getB() +
980  getC() * getC() +
981  getD() * getD());
982 
983  if (v != 0.0) {
984 
985  const JVersor3D u(getB(), getC(), getD());
986 
987  const double theta = atan2(v, getA());
988 
989  *this = JQuaternion3D(2.0 * theta * y, u);
990  }
991 
992  return *this;
993  }
double u[N+1]
Definition: JPolint.hh:865

◆ interpolate()

JQuaternion3D& JGEOMETRY3D::JQuaternion3D::interpolate ( const JQuaternion3D object,
const double  alpha 
)
inlineinherited

Interpolation between quaternions.

The result is equal to *this = (1 - alpha) * (*this) + (alpha) * (object).

Parameters
objectobject
alphainterpolation factor [0, 1]
Returns
this object

Definition at line 1004 of file JQuaternion3D.hh.

1006  {
1007  const double MAXIMAL_DOT_PRODUCT = 0.9995;
1008 
1009  JQuaternion3D v0(*this);
1010  JQuaternion3D v1(object);
1011 
1012  v0.normalise();
1013  v1.normalise();
1014 
1015  double dot = v0.getDot(v1);
1016 
1017  if (dot < 0.0) {
1018  v1 = -v1;
1019  dot = -dot;
1020  }
1021 
1022  double s1 = alpha;
1023  double s0 = 1.0 - alpha;
1024 
1025  if (dot <= MAXIMAL_DOT_PRODUCT) {
1026 
1027  const double theta_0 = acos(dot);
1028  const double theta_1 = theta_0 * alpha;
1029 
1030  s1 = sin(theta_1) / sin(theta_0);
1031  s0 = cos(theta_1) - dot * s1;
1032  }
1033 
1034  *this = (s0 * v0) + (s1 * v1);
1035 
1036  return normalise();
1037  }
JQuaternion3D & normalise()
Normalise quaternion.

◆ getCalibration() [1/2]

const JCalibration& JDETECTOR::JCalibration::getCalibration ( ) const
inlineinherited

Get calibration.

Returns
calibration

Definition at line 67 of file JDetector/JCalibration.hh.

68  {
69  return *this;
70  }

◆ getCalibration() [2/2]

JCalibration& JDETECTOR::JCalibration::getCalibration ( )
inlineinherited

Get calibration.

Returns
calibration

Definition at line 78 of file JDetector/JCalibration.hh.

79  {
80  return *this;
81  }

◆ setCalibration()

void JDETECTOR::JCalibration::setCalibration ( const JCalibration cal)
inlineinherited

Set calibration.

Parameters
calcalibration

Definition at line 89 of file JDetector/JCalibration.hh.

90  {
91  *this = cal;
92  }

◆ getT0()

double JDETECTOR::JCalibration::getT0 ( ) const
inlineinherited

Get time offset.

Returns
time offset [ns]

Definition at line 100 of file JDetector/JCalibration.hh.

101  {
102  return t0;
103  }

◆ setT0()

void JDETECTOR::JCalibration::setT0 ( const double  t0)
inlineinherited

Set time offset.

Parameters
t0time offset [ns]

Definition at line 111 of file JDetector/JCalibration.hh.

112  {
113  this->t0 = t0;
114  }

◆ addT0()

void JDETECTOR::JCalibration::addT0 ( const double  t0)
inlineinherited

Add time offset.

Parameters
t0time offset [ns]

Definition at line 122 of file JDetector/JCalibration.hh.

123  {
124  this->t0 += t0;
125  }

◆ subT0()

void JDETECTOR::JCalibration::subT0 ( const double  t0)
inlineinherited

Subtract time offset.

Parameters
t0time offset [ns]

Definition at line 133 of file JDetector/JCalibration.hh.

134  {
135  this->t0 -= t0;
136  }

◆ getStatus() [1/3]

int JEEP::JStatus::getStatus ( ) const
inlineinherited

Get status.

Returns
status

Definition at line 63 of file JStatus.hh.

64  {
65  return this->status;
66  }

◆ getStatus() [2/3]

JStatus& JEEP::JStatus::getStatus ( )
inlineinherited

Get status.

Returns
status

Definition at line 74 of file JStatus.hh.

75  {
76  return *this;
77  }

◆ getStatus() [3/3]

int JEEP::JStatus::getStatus ( const int  mask) const
inlineinherited

Get status.

Parameters
maskmask
Returns
status

Definition at line 86 of file JStatus.hh.

87  {
88  return (this->status & mask);
89  }

◆ setStatus()

void JEEP::JStatus::setStatus ( const JStatus status)
inlineinherited

Set status.

Parameters
statusstatus

Definition at line 97 of file JStatus.hh.

98  {
99  this->status = status.status;
100  }

◆ getCombinatorics()

const JCombinatorics& JTOOLS::JCombinatorics::getCombinatorics ( ) const
inlineinherited

Get combinatorics.

Returns
combinatorics

Definition at line 89 of file JCombinatorics.hh.

90  {
91  return *this;
92  }

◆ configure()

void JTOOLS::JCombinatorics::configure ( const int  numberOfIndices)
inlineinherited

Configure.

Parameters
numberOfIndicesnumber of indices

Definition at line 100 of file JCombinatorics.hh.

101  {
102  zbuf1D.clear();
103 
104  zbuf2D.resize(numberOfIndices);
105 
106  for (int i = 0; i != numberOfIndices; ++i) {
107  zbuf2D[i].resize(numberOfIndices);
108  }
109 
110  for (int i = 0; i != numberOfIndices; ++i) {
111 
112  zbuf2D[i][i] = -1;
113 
114  for (int j = i; ++j != numberOfIndices; ) {
115 
116  zbuf2D[i][j] = zbuf1D.size();
117  zbuf2D[j][i] = zbuf1D.size();
118 
119  zbuf1D.push_back(pair_type(i,j));
120  }
121  }
122  }
std::vector< pair_type > zbuf1D
std::vector< std::vector< int > > zbuf2D
JCombinatorics::pair_type pair_type
int j
Definition: JPolint.hh:792

◆ getNumberOfIndices()

size_t JTOOLS::JCombinatorics::getNumberOfIndices ( ) const
inlineinherited

Get number of indices.

Returns
number of indices

Definition at line 130 of file JCombinatorics.hh.

131  {
132  return zbuf2D.size();
133  }

◆ getNumberOfPairs()

size_t JTOOLS::JCombinatorics::getNumberOfPairs ( ) const
inlineinherited

Get number of pairs.

Returns
number of pairs

Definition at line 141 of file JCombinatorics.hh.

142  {
143  return zbuf1D.size();
144  }

◆ getIndex() [6/7]

int JTOOLS::JCombinatorics::getIndex ( const int  first,
const int  second 
) const
inlineinherited

Get index of pair of indices.

Parameters
firstfirst address
secondsecond address
Returns
index (-1 if first and second address are equal)

Definition at line 154 of file JCombinatorics.hh.

155  {
156  return zbuf2D[first][second];
157  }

◆ getIndex() [7/7]

int JTOOLS::JCombinatorics::getIndex ( const pair_type pair) const
inlineinherited

Get index of pair.

Parameters
pairpair
Returns
index (-1 if first and second address are equal)

Definition at line 166 of file JCombinatorics.hh.

167  {
168  return getIndex(pair.first, pair.second);
169  }
int getIndex(const int first, const int second) const
Get index of pair of indices.

◆ getPair()

const pair_type& JTOOLS::JCombinatorics::getPair ( const int  index) const
inlineinherited

Get pair of indices for given index.

Parameters
indexindex
Returns
pair of indices

Definition at line 178 of file JCombinatorics.hh.

179  {
180  return zbuf1D[index];
181  }

◆ sort()

template<class JComparator_t >
void JTOOLS::JCombinatorics::sort ( JComparator_t  comparator)
inlineinherited

Sort address pairs.

Parameters
comparatorcomparator for pairs

Definition at line 190 of file JCombinatorics.hh.

191  {
192  std::stable_sort(zbuf1D.begin(), zbuf1D.end(), comparator);
193 
194  for (int i = 0; i != (int) zbuf1D.size(); ++i) {
195 
196  const pair_type pair = zbuf1D[i];
197 
198  zbuf2D[pair.first][pair.second] = i;
199  zbuf2D[pair.second][pair.first] = i;
200  }
201 
202  for (int i = 0; i != (int) zbuf2D.size(); ++i) {
203  zbuf2D[i][i] = -1;
204  }
205  }

◆ getSign() [1/2]

static int JTOOLS::JCombinatorics::getSign ( const int  first,
const int  second 
)
inlinestaticinherited

Sign of pair of indices.

Parameters
firstfirst address
secondsecond address
Returns
+1 if second >= first; else -1

Definition at line 215 of file JCombinatorics.hh.

216  {
217  return (second >= first ? +1 : -1);
218  }

◆ getSign() [2/2]

static int JTOOLS::JCombinatorics::getSign ( const pair_type pair)
inlinestaticinherited

Sign of pair of indices.

Parameters
pairpair of indices
Returns
+1 if second >= first; else -1

Definition at line 227 of file JCombinatorics.hh.

228  {
229  return getSign(pair.first, pair.second);
230  }
static int getSign(const int first, const int second)
Sign of pair of indices.

Friends And Related Function Documentation

◆ operator<<

std::ostream& operator<< ( std::ostream &  out,
const JModel object 
)
friend

Write model parameters to output stream.

Parameters
outoutput stream
objectmodel parameters
Returns
output stream

Definition at line 1142 of file JFitK40.hh.

1143  {
1144  using namespace std;
1145 
1146  out << "Module " << setw(10) << object.getID() << endl;
1147  out << "option " << object.option << endl;
1148  out << "index " << object.index << endl;
1149  out << "Rate [Hz] " << FIXED(12,6) << object.R << endl;
1150  out << "p1 " << FIXED(12,6) << object.p1 << endl;
1151  out << "p2 " << FIXED(12,6) << object.p2 << endl;
1152  out << "p3 " << FIXED(12,6) << object.p3 << endl;
1153  out << "p4 " << FIXED(12,6) << object.p4 << endl;
1154  out << "cc " << FIXED(12,6) << object.cc << endl;
1155 
1156  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
1157  out << "PMT[" << FILL(2,'0') << i << FILL() << "]." << object.parameters[i].status << endl << object.parameters[i];
1158  }
1159 
1160  return out;
1161  }
Auxiliary data structure for sequence of same character.
Definition: JManip.hh:330
Auxiliary data structure for floating point format specification.
Definition: JManip.hh:448

Member Data Documentation

◆ parameters

JPMTParameters_t JCALIBRATE::JModel::parameters[NUMBER_OF_PMTS]

Definition at line 1163 of file JFitK40.hh.

◆ index

int JCALIBRATE::JModel::index
private

index of PMT used for fixed time offset

Definition at line 1166 of file JFitK40.hh.

◆ sigmaK40_ns

double JCALIBRATE::JModel::sigmaK40_ns = 0.54
private

intrinsic K40 arrival time spread [ns]

Definition at line 1167 of file JFitK40.hh.

◆ option

JOption_t JCALIBRATE::JModel::option
private

Definition at line 1168 of file JFitK40.hh.

◆ real

real_type JCALIBRATE::JModel::real
mutableprivate

Definition at line 1169 of file JFitK40.hh.

◆ gradient

JK40Parameters_t JCALIBRATE::JK40Parameters::gradient
mutableprivateinherited

Definition at line 741 of file JFitK40.hh.

◆ R

JParameter_t JCALIBRATE::JK40Parameters_t::R
inherited

maximal coincidence rate [Hz]

Definition at line 596 of file JFitK40.hh.

◆ p1

JParameter_t JCALIBRATE::JK40Parameters_t::p1
inherited

1st order angle dependence coincidence rate

Definition at line 597 of file JFitK40.hh.

◆ p2

JParameter_t JCALIBRATE::JK40Parameters_t::p2
inherited

2nd order angle dependence coincidence rate

Definition at line 598 of file JFitK40.hh.

◆ p3

JParameter_t JCALIBRATE::JK40Parameters_t::p3
inherited

3rd order angle dependence coincidence rate

Definition at line 599 of file JFitK40.hh.

◆ p4

JParameter_t JCALIBRATE::JK40Parameters_t::p4
inherited

4th order angle dependence coincidence rate

Definition at line 600 of file JFitK40.hh.

◆ cc

JParameter_t JCALIBRATE::JK40Parameters_t::cc
inherited

fraction of signal correlated background

Definition at line 601 of file JFitK40.hh.

◆ __id

int JLANG::JObjectID::__id
protectedinherited

Definition at line 140 of file JObjectID.hh.

◆ string

int JDETECTOR::JLocation::string
protectedinherited

Definition at line 235 of file JLocation.hh.

◆ floor

int JDETECTOR::JLocation::floor
protectedinherited

Definition at line 236 of file JLocation.hh.

◆ __x

double JGEOMETRY3D::JVector3D::__x
protectedinherited

Definition at line 310 of file JVector3D.hh.

◆ __y

double JGEOMETRY3D::JVector3D::__y
protectedinherited

Definition at line 311 of file JVector3D.hh.

◆ __z

double JGEOMETRY3D::JVector3D::__z
protectedinherited

Definition at line 312 of file JVector3D.hh.

◆ __a

double JGEOMETRY3D::JQuaternion3D::__a
protectedinherited

Definition at line 1117 of file JQuaternion3D.hh.

◆ __b

double JGEOMETRY3D::JQuaternion3D::__b
protectedinherited

Definition at line 1118 of file JQuaternion3D.hh.

◆ __c

double JGEOMETRY3D::JQuaternion3D::__c
protectedinherited

Definition at line 1119 of file JQuaternion3D.hh.

◆ __d

double JGEOMETRY3D::JQuaternion3D::__d
protectedinherited

Definition at line 1120 of file JQuaternion3D.hh.

◆ t0

double JDETECTOR::JCalibration::t0
protectedinherited

Definition at line 202 of file JDetector/JCalibration.hh.

◆ status

int JEEP::JStatus::status
protectedinherited

Definition at line 200 of file JStatus.hh.

◆ zbuf1D

std::vector<pair_type> JTOOLS::JCombinatorics::zbuf1D
protectedinherited

Definition at line 234 of file JCombinatorics.hh.

◆ zbuf2D

std::vector<std::vector<int> > JTOOLS::JCombinatorics::zbuf2D
protectedinherited

Definition at line 235 of file JCombinatorics.hh.


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