Jpp test-rotations-new
the software that should make you happy
Loading...
Searching...
No Matches
JCALIBRATE::JModel Struct Reference

Fit model. More...

#include <JFitK40.hh>

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

Classes

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

Public Member Functions

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

Static Public Member Functions

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

Public Attributes

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

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
 
double sigmaK40_ns = 0.54
 intrinsic K40 arrival time spread [ns]
 
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.
 

Detailed Description

Fit model.

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

Definition at line 832 of file JFitK40.hh.

Constructor & Destructor Documentation

◆ JModel() [1/3]

JCALIBRATE::JModel::JModel ( )
inline

Default constructor.

Definition at line 858 of file JFitK40.hh.

859 {}

◆ 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 870 of file JFitK40.hh.

873 :
874 JModule (module),
875 JCombinatorics_t(module)
876 {
878
879 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
881 }
882
883 for (JTDC_t::const_iterator i = TDC.first; i != TDC.second; ++i) {
884
885 if (i->second != JTDC_t::WILDCARD) {
886
887 this->parameters[i->second].t0.fix();
888
889 } else {
890
891 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
892 this->parameters[i].t0.fix();
893 }
894 }
895 }
896
897 this->index = (TDC.first == TDC.second ? 0 : INVALID_INDEX);
898
900 }
void fix()
Fix current value.
Definition JFitK40.hh:280
JModule()
Default constructor.
Definition JModule.hh:85
static const int INVALID_INDEX
invalid index
Definition JFitK40.hh:60
JCombinatorics_t()
Default constructor.
void setK40Parameters(const JK40Parameters_t &parameters)
Set K40 parameters.
Definition JFitK40.hh:644
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition JFitK40.hh:797
void setOption(const int option)
Set fit option.
Definition JFitK40.hh:938
int index
index of PMT used for fixed time offset
Definition JFitK40.hh:1241
JParameter_t t0
time offset [ns]
Definition JFitK40.hh:610
static const JPMTParameters_t & getInstance()
Get default values.
Definition JFitK40.hh:477

◆ JModel() [3/3]

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

Constructor.

Parameters
moduledetector module
parametersK40 parameters

Definition at line 909 of file JFitK40.hh.

910 :
911 JModule (module),
912 JCombinatorics_t(module)
913 {
915
916 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
918 }
919 }

Member Function Documentation

◆ getOption()

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

Get fit option.

Returns
option

Definition at line 927 of file JFitK40.hh.

928 {
929 return option;
930 }

◆ setOption()

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

Set fit option.

Parameters
optionoption

Definition at line 938 of file JFitK40.hh.

939 {
940 switch (option) {
941
942 case FIT_PMTS_t:
943
944 R .fix();
945 p1.fix();
946 p2.fix();
947 p3.fix();
948 p4.fix();
949
950 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
951 parameters[i].bg.fix();
952 }
953
954 break;
955
957
958 R .fix();
959
960 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
961 parameters[i].bg.fix();
962 }
963
964 break;
965
967
968 R .fix();
969 p1.fix();
970 p2.fix();
971 p3.fix();
972 p4.fix();
973
974 break;
975
977
978 R .fix();
979 p1.fix();
980 p2.fix();
981 p3.fix();
982 p4.fix();
983
984 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
985 parameters[i].QE.fix();
986 parameters[i].bg.fix();
987 }
988
989 break;
990
991 case FIT_MODEL_t:
992
993 cc.fix();
994
995 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
996 parameters[i].QE.fix();
997 parameters[i].t0.fix();
998 parameters[i].bg.fix();
999 }
1000
1001 break;
1002
1003 default:
1004
1005 THROW(JValueOutOfRange, "Invalid option " << option);
1006 }
1007
1008 this->option = static_cast<JOption_t>(option);
1009 }
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
@ 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:664
JParameter_t p1
1st order angle dependence coincidence rate
Definition JFitK40.hh:665
JParameter_t p2
2nd order angle dependence coincidence rate
Definition JFitK40.hh:666
JParameter_t p3
3rd order angle dependence coincidence rate
Definition JFitK40.hh:667
JParameter_t p4
4th order angle dependence coincidence rate
Definition JFitK40.hh:668
JParameter_t cc
fraction of signal correlated background
Definition JFitK40.hh:669
JParameter_t bg
background [Hz/ns]
Definition JFitK40.hh:611
JParameter_t QE
relative quantum efficiency [unit]
Definition JFitK40.hh:608

◆ hasFixedTimeOffset()

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

Check if time offset is fixed.

Returns
true if time offset fixed; else false

Definition at line 1017 of file JFitK40.hh.

1018 {
1019 return index != INVALID_INDEX;
1020 }

◆ getFixedTimeOffset()

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

Get time offset.

Returns
time offset

Definition at line 1028 of file JFitK40.hh.

1029 {
1030 double t0 = 0.0;
1031 size_t N = 0;
1032
1033 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
1034 if (parameters[i].t0.isFree()) {
1035 t0 += parameters[i].t0;
1036 N += 1;
1037 }
1038 }
1039
1040 return t0 /= N;
1041 }

◆ getIndex() [1/6]

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

Get index of PMT used for fixed time offset.

Returns
index

Definition at line 1049 of file JFitK40.hh.

1050 {
1051 return index;
1052 }

◆ setIndex()

void JCALIBRATE::JModel::setIndex ( )
inline

Set index of PMT used for fixed time offset.

Definition at line 1058 of file JFitK40.hh.

1059 {
1060 if (index != INVALID_INDEX) {
1061
1062 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
1063
1064 if (parameters[i].status) {
1065
1066 index = i;
1067
1069
1070 return;
1071 }
1072 }
1073
1074 THROW(JValueOutOfRange, "No valid index.");
1075 }
1076 }

◆ getN()

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

Get number of fit parameters.

Returns
number of parameters

Definition at line 1084 of file JFitK40.hh.

1085 {
1086 size_t N = JK40Parameters::getN();
1087
1088 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
1089 N += parameters[i].getN();
1090 }
1091
1092 return N;
1093 }
size_t getN() const
Get number of fit parameters.
Definition JFitK40.hh:715
size_t getN() const
Get number of fit parameters.
Definition JFitK40.hh:527

◆ getIndex() [2/6]

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

Get index of parameter.

Parameters
pmtPMT
Returns
index

Definition at line 1102 of file JFitK40.hh.

1103 {
1104 int N = JK40Parameters::getN();
1105
1106 for (int i = 0; i != pmt; ++i) {
1107 N += parameters[i].getN();
1108 }
1109
1110 return N;
1111 }

◆ getIndex() [3/6]

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 1121 of file JFitK40.hh.

1122 {
1123 return getIndex(pmt) + parameters[pmt].getIndex(p);
1124 }
int getIndex() const
Get index of PMT used for fixed time offset.
Definition JFitK40.hh:1049
int getIndex(JParameter_t JPMTParameters_t::*p) const
Get index of parameter.
Definition JFitK40.hh:542

◆ getSigmaK40()

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

Get intrinsic K40 arrival time spread.

Returns
sigma [ns]

Definition at line 1132 of file JFitK40.hh.

1133 {
1134 return this->sigmaK40_ns;
1135 }
double sigmaK40_ns
intrinsic K40 arrival time spread [ns]
Definition JFitK40.hh:1242

◆ setSigmaK40()

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

Set intrinsic K40 arrival time spread.

Parameters
sigmasigma [ns]

Definition at line 1143 of file JFitK40.hh.

1144 {
1145 this->sigmaK40_ns = sigma;
1146 }
const double sigma[]

◆ getReal()

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

Get derived parameters.

Parameters
pairPMT pair
Returns
parameters

Definition at line 1155 of file JFitK40.hh.

1156 {
1157 real.ct = JPP::getDot((*this)[pair.first].getDirection(), (*this)[pair.second].getDirection());
1158
1159 real.t0 = (pair.first == this->index ? -this->parameters[pair.second].t0() :
1160 pair.second == this->index ? +this->parameters[pair.first ].t0() :
1161 this->parameters[pair.first].t0() - this->parameters[pair.second].t0());
1162
1163 real.sigma = sqrt(this->parameters[pair.first ].TTS() * this->parameters[pair.first ].TTS() +
1164 this->parameters[pair.second].TTS() * this->parameters[pair.second].TTS() +
1165 this->getSigmaK40() * this->getSigmaK40());
1166
1167 real.signal = this->parameters[pair.first].QE() * this->parameters[pair.second].QE();
1168
1169 real.background = this->parameters[pair.first].bg() + this->parameters[pair.second].bg();
1170
1171 real.cc = real.signal * this->cc();
1172
1173 return real;
1174 }
double getDot(const JNeutrinoDirection &first, const JNeutrinoDirection &second)
Dot product.
double signal
combined signal
Definition JFitK40.hh:849
double sigma
total width [ns]
Definition JFitK40.hh:848
double cc
correlated background
Definition JFitK40.hh:851
double background
combined background
Definition JFitK40.hh:850
double t0
time offset [ns]
Definition JFitK40.hh:847
double ct
cosine angle between PMT axes
Definition JFitK40.hh:846
JParameter_t TTS
transition-time spread [ns]
Definition JFitK40.hh:609

◆ getValue() [1/3]

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 1184 of file JFitK40.hh.

1185 {
1186 using namespace std;
1187 using namespace JPP;
1188
1189 const real_type& real = getReal(pair);
1190
1191 const JGauss gauss(real.t0, real.sigma, real.signal);
1192
1193 const double R1 = this->getValue(real.ct);
1194 const double R2 = gauss.getValue(dt_ns);
1195
1196 return real.background + R1 * (real.cc + R2);
1197 }
#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).
const real_type & getReal(const pair_type &pair) const
Get derived parameters.
Definition JFitK40.hh:1155
double getValue(const pair_type &pair, const double dt_ns) const
Get K40 coincidence rate.
Definition JFitK40.hh:1184

◆ getValue() [2/3]

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

Get K40 coincidence rate.

Parameters
pairPMT pair
Returns
rate [Hz]

Definition at line 1206 of file JFitK40.hh.

1207 {
1208 using namespace std;
1209 using namespace JPP;
1210
1211 const real_type& real = getReal(pair);
1212
1213 const double R1 = this->getValue(real.ct);
1214 const double R2 = real.signal;
1215
1216 return real.background + R1 * (cc() + R2);
1217 }

◆ getIndex() [4/6]

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

Get index of parameter.

Parameters
ppointer to data member
Returns
index

Definition at line 732 of file JFitK40.hh.

733 {
734 if (!(this->*p).isFree()) {
735 return INVALID_INDEX;
736 }
737
738 int N = 0;
739
740 if (p == &JK40Parameters::R) { return N; } if (R .isFree()) { ++N; }
741 if (p == &JK40Parameters::p1) { return N; } if (p1.isFree()) { ++N; }
742 if (p == &JK40Parameters::p2) { return N; } if (p2.isFree()) { ++N; }
743 if (p == &JK40Parameters::p3) { return N; } if (p3.isFree()) { ++N; }
744 if (p == &JK40Parameters::p4) { return N; } if (p4.isFree()) { ++N; }
745 if (p == &JK40Parameters::cc) { return N; } if (cc.isFree()) { ++N; }
746
747 return INVALID_INDEX;
748 }
bool isFree() const
Check if parameter is free.
Definition JFitK40.hh:240

◆ getValue() [3/3]

double JCALIBRATE::JK40Parameters::getValue ( const double ct) const
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 757 of file JFitK40.hh.

758 {
759 return R * exp(ct*(p1+ct*(p2+ct*(p3+ct*p4))) - (p1+p2+p3+p4));
760 }

◆ 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 695 of file JFitK40.hh.

696 {
697 static JK40Parameters parameters;
698
699 parameters.R .set(18.460546);
700 parameters.p1.set( 3.0767);
701 parameters.p2.set(-1.2078);
702 parameters.p3.set( 0.9905);
703 parameters.p4.set( 0.9379);
704 parameters.cc.set( 0.0);
705
706 return parameters;
707 }
JK40Parameters()
Default constructor.
Definition JFitK40.hh:682

◆ 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 769 of file JFitK40.hh.

770 {
771 gradient.reset();
772
773 const double rate = getValue(ct);
774 const double ct2 = ct * ct;
775
776 if (R .isFree()) { gradient.R = rate / R; }
777 if (p1.isFree()) { gradient.p1 = rate * ct - rate; }
778 if (p2.isFree()) { gradient.p2 = rate * ct2 - rate; }
779 if (p3.isFree()) { gradient.p3 = rate * ct2 * ct - rate; }
780 if (p4.isFree()) { gradient.p4 = rate * ct2 * ct2 - rate; }
781 if (cc.isFree()) { gradient.cc = rate; }
782
783 return gradient;
784 }
JK40Parameters_t gradient
Definition JFitK40.hh:787
double getValue(const double ct) const
Get K40 coincidence rate as a function of cosine angle between PMT axes.
Definition JFitK40.hh:757

◆ getK40Parameters()

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

Get K40 parameters.

Returns
K40 parameters

Definition at line 633 of file JFitK40.hh.

634 {
635 return static_cast<const JK40Parameters_t&>(*this);
636 }
JK40Parameters_t()
Default constructor.
Definition JFitK40.hh:622

◆ setK40Parameters()

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

Set K40 parameters.

Parameters
parametersK40 parameters

Definition at line 644 of file JFitK40.hh.

645 {
646 static_cast<JK40Parameters_t&>(*this) = parameters;
647 }

◆ reset() [1/2]

void JCALIBRATE::JK40Parameters_t::reset ( )
inlineinherited

Reset.

Definition at line 653 of file JFitK40.hh.

654 {
655 R .set(0.0);
656 p1.set(0.0);
657 p2.set(0.0);
658 p3.set(0.0);
659 p4.set(0.0);
660 cc.set(0.0);
661 }
void set()
Set current value.
Definition JFitK40.hh:271

◆ 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 JAngle3D &angle) const
Get dot product.
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
Exception for accessing a value in a collection that is outside of its range.
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.
JPosition3D()
Default constructor.
const JPosition3D & getPosition() const
Get position.
JVector3D & add(const JVector3D &vector)
Add vector.
Definition JVector3D.hh:142
JVector3D & div(const double factor)
Scale vector.
Definition JVector3D.hh:190
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.

◆ 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 }

◆ 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.

◆ 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 }

◆ 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/3]

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 {
349
350 for (iterator i = this->begin(); i != this->end(); ++i) {
351 i->transform(R, pos);
352 }
353 }
void transform(const JRotation3D &R, const JVector3D &pos)
Transform position.

◆ transform() [2/3]

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

Transformation of geometry.

Parameters
Ttransformation

Definition at line 361 of file JModule.hh.

362 {
363 JPosition3D::transform(T.getRotation(), T.getPosition());
364
365 for (iterator i = this->begin(); i != this->end(); ++i) {
366 i->transform(T);
367 }
368 }

◆ transform() [3/3]

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/3]

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/3]

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/3]

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 }

◆ 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()

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/2]

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 }
JLocation()
Default constructor.
Definition JLocation.hh:45

◆ 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

◆ 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 }

◆ 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 }
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 }

◆ 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 }
double getB() const
Get b value.
double getD() const
Get d value.
double getC() const
Get c value.

◆ 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/5]

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 }

◆ getDot() [2/5]

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 }

◆ getDot() [3/5]

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 }

◆ getDot() [4/5]

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() [5/5]

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 }
double getA() const
Get a value.

◆ 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 }

◆ 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 }

◆ 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 }

◆ 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 }

◆ 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.

◆ 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 }

◆ 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:801

◆ 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() [5/6]

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() [6/6]

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 Symbol 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 1227 of file JFitK40.hh.

1228 {
1229 using namespace std;
1230
1231 out << "Module " << setw(10) << object.getID() << endl;
1232 out << "option " << object.option << endl;
1233 out << "index " << object.index << endl;
1234
1235 out << static_cast<const JModel_t&>(object);
1236
1237 return out;
1238 }

Member Data Documentation

◆ index

int JCALIBRATE::JModel::index
private

index of PMT used for fixed time offset

Definition at line 1241 of file JFitK40.hh.

◆ sigmaK40_ns

double JCALIBRATE::JModel::sigmaK40_ns = 0.54
private

intrinsic K40 arrival time spread [ns]

Definition at line 1242 of file JFitK40.hh.

◆ option

JOption_t JCALIBRATE::JModel::option
private

Definition at line 1243 of file JFitK40.hh.

◆ real

real_type JCALIBRATE::JModel::real
mutableprivate

Definition at line 1244 of file JFitK40.hh.

◆ parameters

JPMTParameters_t JCALIBRATE::JModel_t::parameters[NUMBER_OF_PMTS]
inherited

Definition at line 797 of file JFitK40.hh.

◆ gradient

JK40Parameters_t JCALIBRATE::JK40Parameters::gradient
mutableprivateinherited

Definition at line 787 of file JFitK40.hh.

◆ R

JParameter_t JCALIBRATE::JK40Parameters_t::R
inherited

maximal coincidence rate [Hz]

Definition at line 664 of file JFitK40.hh.

◆ p1

JParameter_t JCALIBRATE::JK40Parameters_t::p1
inherited

1st order angle dependence coincidence rate

Definition at line 665 of file JFitK40.hh.

◆ p2

JParameter_t JCALIBRATE::JK40Parameters_t::p2
inherited

2nd order angle dependence coincidence rate

Definition at line 666 of file JFitK40.hh.

◆ p3

JParameter_t JCALIBRATE::JK40Parameters_t::p3
inherited

3rd order angle dependence coincidence rate

Definition at line 667 of file JFitK40.hh.

◆ p4

JParameter_t JCALIBRATE::JK40Parameters_t::p4
inherited

4th order angle dependence coincidence rate

Definition at line 668 of file JFitK40.hh.

◆ cc

JParameter_t JCALIBRATE::JK40Parameters_t::cc
inherited

fraction of signal correlated background

Definition at line 669 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: