Jpp
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Protected Attributes | Private Member Functions | Static Private Member Functions | Private Attributes | List of all members
JCALIBRATE::JFitK40 Struct Reference

Parametrisation of two-fold coincidence rate due to K40 and other radioactive decays. More...

#include <JFitK40.hh>

Inheritance diagram for JCALIBRATE::JFitK40:
JCALIBRATE::JFitK40Parameters JDETECTOR::JModule JTOOLS::JCombinatorics JLANG::JObjectID JDETECTOR::JModuleLocation JGEOMETRY3D::JPosition3D std::vector< JPMT > JLANG::JComparable< JObjectID > JLANG::JComparable< JObjectID, int > JLANG::JComparable< JFirst_t, JSecond_t > JGEOMETRY3D::JVector3D JMATH::JMath< JVector3D >

Public Types

typedef JCombinatorics::pair_type pair_type
 Data structure for a pair of addresses. More...
 

Public Member Functions

 JFitK40 (const JModule &module, const Double_t xmin, const Double_t xmax, const Double_t ymin, const Double_t ymax, const bool option)
 Constructor. More...
 
Double_t getSigmaK40 () const
 Get intrinsic K40 arrival time spread. More...
 
void setSigmaK40 (const Double_t sigma)
 Set intrinsic K40 arrival time spread. More...
 
bool is_disabled (int pmt) const
 Test PMT status. More...
 
bool is_average_t0 (int pmt) const
 Test PMT status. More...
 
Double_t getT0 (const int pmt) const
 Get time offset of given PMT. More...
 
void enablePMT (const int pmt)
 Enable PMT. More...
 
void disablePMT (const int pmt)
 Disable PMT. More...
 
double getDot (const pair_type &pair) const
 Get cosine of space angle between PMT axes. More...
 
Double_t getSigma (const pair_type &pair) const
 Get time resolution of given PMT pair. More...
 
Double_t getT0 (const pair_type &pair) const
 Get time offset of given PMT pair. More...
 
Double_t getValue (const pair_type &pair) const
 Get K40 coincidence rate. More...
 
Double_t getValue (const pair_type &pair, const Double_t dt_ns) const
 Get K40 coincidence rate. More...
 
TFitResultPtr operator() (TH2 &h2, const std::string &option)
 Fit histogram. More...
 
Double_t getValue (const Double_t ct) const
 Get K40 coincidence rate as a function of cosine angle between PMT axes. More...
 
const Double_t * getModelParameters () const
 Get model parameters. More...
 
Double_t * getModelParameters ()
 Get model parameters. More...
 
void setModelParameters (const Double_t *data)
 Set model parameters. More...
 
Double_t getModelParameter (Int_t i) const
 Get model parameter. More...
 
JFitParameter_t getModelParameter (Double_t JFitK40Parameters::*p) const
 Get model parameter. More...
 
JFitParameter_t getModelParameter (Int_t pmt, Double_t JPMTParameters_t::*p) const
 Get model parameter. More...
 
Double_t getQE (const int pmt) const
 Get QE of given PMT. More...
 
void setQE (const int pmt, const Double_t QE)
 Set QE of given PMT. More...
 
Double_t getTTS (const int pmt) const
 Get time resolution of given PMT. More...
 
void setTTS (const int pmt, const Double_t TTS)
 Set time resolution of given PMT. More...
 
void setT0 (const int pmt, const Double_t t0)
 Set time offset of given PMT. More...
 
Double_t getValue (const Double_t ct) const
 Get K40 coincidence rate as a function of cosine angle between PMT axes. 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...
 
void compile ()
 Compile position of module from the positions of the PMTs. 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_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 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...
 
JModuleadd (const JVector3D &pos)
 Add position. More...
 
JModuleadd (const double t0)
 Add time offset. More...
 
JModulesub (const JVector3D &pos)
 Subtract position. More...
 
JModulesub (const double t0)
 Subtract time offset. More...
 
JModuleoperator+= (const JVector3D &pos)
 Add position. More...
 
JModuleoperator-= (const JVector3D &pos)
 Subtract position. More...
 
int getID () const
 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 JModuleLocation &location) const
 Less than method. More...
 
bool more (const int id) const
 More than method. More...
 
const JModuleLocationgetLocation () const
 Get location. More...
 
JModuleLocationgetLocation ()
 Get location. More...
 
void setLocation (const JModuleLocation &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...
 
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...
 
 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...
 
JVector3Dmul (const double factor)
 Scale vector. More...
 
JVector3Dmul (const JSecond_t &object)
 Multiply with object. More...
 
JVector3Ddiv (const double factor)
 Scale vector. More...
 
bool equals (const JVector3D &vector, const double precision=std::numeric_limits< double >::min()) const
 Check equality. More...
 
double getLengthSquared () const
 Get length squared. More...
 
double getLength () const
 Get length. More...
 
double getDistanceSquared (const JVector3D &pos) const
 Get squared of distance to point. More...
 
double getDistance (const JVector3D &pos) const
 Get distance to point. More...
 
JVector3Dcross (const JVector3D &first, const JVector3D &second)
 Get cross product. More...
 
void configure (const int numberOfIndices)
 Configure. More...
 
unsigned int getNumberOfIndices () const
 Get number of indices. More...
 
unsigned int getNumberOfPairs () const
 Get number of pairs. More...
 
int getIndex (const int first, const int second) const
 Get index of pair of indices. 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 Double_t getRate (const Double_t *x, const Double_t *data)
 Get K40 coincidence rate as a function of the fit parameters. More...
 
static Int_t getNumberOfModelParameters ()
 Get number of model parameters. More...
 
static bool compare (const JModule &first, const JModule &second, const double precision=1.0e-3)
 Compare modules. 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

Double_t Rate_Hz
 maximal coincidence rate [Hz] More...
 
Double_t p1
 angle dependence coincidence rate More...
 
Double_t p2
 angle dependence coincidence rate More...
 
Double_t p3
 angle dependence coincidence rate More...
 
Double_t p4
 angle dependence coincidence rate More...
 
Double_t bg
 remaining constant background More...
 
Double_t cc
 fraction of signal correlated background More...
 
JPMTParameters_t parameters [NUMBER_OF_PMTS]
 

Protected Attributes

int __id
 
int string
 
int floor
 
double __x
 
double __y
 
double __z
 
std::vector< pair_typezbuf1D
 
std::vector< std::vector< int > > zbuf2D
 

Private Member Functions

 JFitK40 ()
 Default constructor. More...
 

Static Private Member Functions

static JFitK40getInstance ()
 Get unique instance of fit object. More...
 

Private Attributes

Double_t sigmaK40_ns
 intrinsic K40 arrival time spread [ns] More...
 
int index_of_average_t0
 index of t0 used for average time offset More...
 
bool disable [NUMBER_OF_PMTS]
 disable PMT from fit More...
 

Detailed Description

Parametrisation of two-fold coincidence rate due to K40 and other radioactive decays.

Note that for use in ROOT fit operations, the member method JFitK40::getRate is static.

Definition at line 286 of file JFitK40.hh.

Member Typedef Documentation

◆ pair_type

Data structure for a pair of addresses.

Definition at line 297 of file JFitK40.hh.

Constructor & Destructor Documentation

◆ JFitK40() [1/2]

JCALIBRATE::JFitK40::JFitK40 ( const JModule module,
const Double_t  xmin,
const Double_t  xmax,
const Double_t  ymin,
const Double_t  ymax,
const bool  option 
)
inline

Constructor.

Parameters
moduledetector module
xminminimal x
xmaxmaximal x
yminminimal y
ymaxmaximal y
optiontrue = fix average t0; false = free average t0

Definition at line 310 of file JFitK40.hh.

315  :
316 
317  TF2("fitk40",
319  xmin, xmax, ymin, ymax,
321 
322  sigmaK40_ns(0.54)
323  {
324  static_cast<JModule&>(*this) = module;
325 
326  this->configure(module.size());
327 
328  this->sort(JPairwiseComparator(module));
329 
330  index_of_average_t0 = (option ? 0 : -1);
331 
332  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
333  disable[i] = false;
334  }
335  }

◆ JFitK40() [2/2]

JCALIBRATE::JFitK40::JFitK40 ( )
inlineprivate

Default constructor.

Definition at line 657 of file JFitK40.hh.

658  {}

Member Function Documentation

◆ getSigmaK40()

Double_t JCALIBRATE::JFitK40::getSigmaK40 ( ) const
inline

Get intrinsic K40 arrival time spread.

Returns
sigma [ns]

Definition at line 343 of file JFitK40.hh.

344  {
345  return this->sigmaK40_ns;
346  }

◆ setSigmaK40()

void JCALIBRATE::JFitK40::setSigmaK40 ( const Double_t  sigma)
inline

Set intrinsic K40 arrival time spread.

Parameters
sigmasigma [ns]

Definition at line 354 of file JFitK40.hh.

355  {
356  this->sigmaK40_ns = sigma;
357  }

◆ is_disabled()

bool JCALIBRATE::JFitK40::is_disabled ( int  pmt) const
inline

Test PMT status.

Parameters
pmtpmt address
Returns
true if given pmt is disabled; else false

Definition at line 366 of file JFitK40.hh.

367  {
368  return disable[pmt];
369  }

◆ is_average_t0()

bool JCALIBRATE::JFitK40::is_average_t0 ( int  pmt) const
inline

Test PMT status.

Parameters
pmtpmt address
Returns
true if given pmt is used for average t0; else false

Definition at line 378 of file JFitK40.hh.

379  {
380  return pmt == index_of_average_t0;
381  }

◆ getT0() [1/2]

Double_t JCALIBRATE::JFitK40::getT0 ( const int  pmt) const
inline

Get time offset of given PMT.

Note that if the average time offset is constraint, the time offset of each PMT is corrected by the average time offset of all PMTs.

Parameters
pmtpmt address
Returns
time offset [ns]

Definition at line 393 of file JFitK40.hh.

394  {
395  if (index_of_average_t0 == -1) {
396 
397  return JFitK40Parameters::getT0(pmt);
398 
399  } else {
400 
401  if (disable[pmt]) {
402 
403  return parameters[pmt].t0;
404 
405  } else {
406 
407  Int_t n0 = 0; // number of non-fixed t0's
408  Double_t t0 = 0.0; // average of non-fixed t0's
409 
410  for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
411  if (!disable[i] && i != index_of_average_t0) {
412  n0 += 1;
413  t0 += parameters[i].t0;
414  }
415  }
416 
417  t0 /= (n0 + 1);
418 
419  if (pmt != index_of_average_t0)
420  return parameters[pmt].t0 - t0;
421  else
422  return -t0;
423  }
424  }
425  }

◆ enablePMT()

void JCALIBRATE::JFitK40::enablePMT ( const int  pmt)
inline

Enable PMT.

Parameters
pmtpmt address

Definition at line 433 of file JFitK40.hh.

434  {
435  disable[pmt] = false;
436  }

◆ disablePMT()

void JCALIBRATE::JFitK40::disablePMT ( const int  pmt)
inline

Disable PMT.

Note that if the average time offset is constraint to zero, the corresponding index may change.

Parameters
pmtpmt address

Definition at line 447 of file JFitK40.hh.

448  {
449  disable[pmt] = true;
450 
451  if (index_of_average_t0 == pmt) {
452 
455  break;
456  }
457  }
458 
460  THROW(JIndexOutOfRange, "JFitK40::disable PMT: No free index for average t0.");
461  }
462  }
463  }

◆ getDot() [1/6]

double JCALIBRATE::JFitK40::getDot ( const pair_type pair) const
inline

Get cosine of space angle between PMT axes.

Parameters
pairpair of addresses
Returns
cosine

Definition at line 472 of file JFitK40.hh.

473  {
474  return JMATH::getDot(this->getPMT(pair.first),
475  this->getPMT(pair.second));
476  }

◆ getSigma()

Double_t JCALIBRATE::JFitK40::getSigma ( const pair_type pair) const
inline

Get time resolution of given PMT pair.

Parameters
pairpair of addresses
Returns
sigma [ns]

Definition at line 485 of file JFitK40.hh.

486  {
487  return TMath::Sqrt(getTTS(pair.first) * getTTS(pair.first) +
488  getTTS(pair.second) * getTTS(pair.second) +
489  getSigmaK40() * getSigmaK40());
490  }

◆ getT0() [2/2]

Double_t JCALIBRATE::JFitK40::getT0 ( const pair_type pair) const
inline

Get time offset of given PMT pair.

Parameters
pairpair of addresses
Returns
time offset [ns]

Definition at line 499 of file JFitK40.hh.

500  {
501  if (index_of_average_t0 == -1) {
502 
503  return parameters[pair.first].t0 - parameters[pair.second].t0;
504 
505  } else {
506 
507  if (pair.first == index_of_average_t0) {
508 
509  return -parameters[pair.second].t0;
510 
511  } else if (pair.second == index_of_average_t0) {
512 
513  return parameters[pair.first].t0;
514 
515  } else {
516 
517  return parameters[pair.first].t0 - parameters[pair.second].t0;
518  }
519  }
520  }

◆ getValue() [1/4]

Double_t JCALIBRATE::JFitK40::getValue ( const pair_type pair) const
inline

Get K40 coincidence rate.

Parameters
pairpair of addresses
Returns
rate [Hz]

Definition at line 529 of file JFitK40.hh.

530  {
531  const Double_t ct = getDot(pair);
532 
533  return (getValue(ct) *
534  getQE(pair.first) *
535  getQE(pair.second));
536  }

◆ getValue() [2/4]

Double_t JCALIBRATE::JFitK40::getValue ( const pair_type pair,
const Double_t  dt_ns 
) const
inline

Get K40 coincidence rate.

Parameters
pairpair of addresses
dt_nstime difference [ns]
Returns
rate [Hz/ns]

Definition at line 546 of file JFitK40.hh.

547  {
548  const Double_t t0 = getT0 (pair);
549  const Double_t sigma = getSigma(pair);
550 
551  return getValue(pair) * TMath::Gaus(dt_ns, t0, sigma, kTRUE);
552  }

◆ operator()()

TFitResultPtr JCALIBRATE::JFitK40::operator() ( TH2 &  h2,
const std::string &  option 
)
inline

Fit histogram.

Note that the PMT parameters which are part of the model are reset before the fit according the status of each PMT and the obtained fit parameters are copied back to the model parameters after the fit.

Parameters
h2ROOT 2D-histogram
optionfit option

Definition at line 564 of file JFitK40.hh.

565  {
566  for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
567 
568  parameters[pmt] = JPMTParameters_t();
569 
570  if (disable[pmt]) {
571 
572  parameters[pmt].QE = 0.0;
573 
577 
578  } else {
579 
580  // Note that setting limits to a parameter will also release it
581 
582  if (!isParameterFixed(*this, this->getModelParameter(pmt, &JPMTParameters_t::QE))) {
584  }
585  if (!isParameterFixed(*this, this->getModelParameter(pmt, &JPMTParameters_t::TTS))) {
587  }
588  if (!isParameterFixed(*this, this->getModelParameter(pmt, &JPMTParameters_t::t0))) {
589  setParLimits(*this, this->getModelParameter(pmt, &JPMTParameters_t::t0), h2.GetYaxis()->GetXmin(), h2.GetYaxis()->GetXmax());
590  }
591  }
592  }
593 
594  if (index_of_average_t0 != -1) {
595 
596  this->setT0(index_of_average_t0, 0.0);
597 
599  }
600 
601  this->SetParameters(this->getModelParameters());
602 
603  getInstance() = *this;
604 
605  const TFitResultPtr result = h2.Fit(this, option.c_str());
606 
607  this->setModelParameters(this->GetParameters());
608 
609  return result;
610  }

◆ getRate()

static Double_t JCALIBRATE::JFitK40::getRate ( const Double_t *  x,
const Double_t *  data 
)
inlinestatic

Get K40 coincidence rate as a function of the fit parameters.

To speed up the calculation, it is assumed that the parameter values do not change unless also the x[0] value changes (i.e. the index of the PMT pair).

Parameters
xpointer to data
datapointer to parameter values
Returns
rate [Hz/ns]

Definition at line 623 of file JFitK40.hh.

624  {
625  static int id = -1;
626  static Double_t t0;
627  static Double_t sigma;
628  static Double_t rate;
629 
630  const int ix = (int) x[0];
631  const Double_t dt_ns = x[1];
632 
633  if (ix != id) {
634 
636 
637  const pair_type& pair = getInstance().getPair(ix);
638 
639  t0 = getInstance().getT0 (pair);
640  sigma = getInstance().getSigma(pair);
641  rate = getInstance().getValue(pair);
642  id = ix;
643  }
644 
645  return getInstance().bg + rate * (getInstance().cc + TMath::Gaus(dt_ns, t0, sigma, kTRUE));
646  }

◆ getInstance()

static JFitK40& JCALIBRATE::JFitK40::getInstance ( )
inlinestaticprivate

Get unique instance of fit object.

Returns
reference to fit object.

Definition at line 666 of file JFitK40.hh.

667  {
668  static JFitK40 object;
669 
670  return object;
671  }

◆ getValue() [3/4]

Double_t JCALIBRATE::JFitK40Parameters::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 261 of file JFitK40.hh.

262  {
263  return Rate_Hz * TMath::Exp(-(p1+p2+p3+p4)) * TMath::Exp(ct*(p1+ct*(p2+ct*(p3+ct*p4))));
264  }

◆ getNumberOfModelParameters()

static Int_t JCALIBRATE::JFitK40Parameters::getNumberOfModelParameters ( )
inlinestaticinherited

Get number of model parameters.

Returns
number of parameters

Definition at line 99 of file JFitK40.hh.

100  {
101  return sizeof(JFitK40Parameters) / sizeof(Double_t);
102  }

◆ getModelParameters() [1/2]

const Double_t* JCALIBRATE::JFitK40Parameters::getModelParameters ( ) const
inlineinherited

Get model parameters.

Returns
pointer to parameters

Definition at line 110 of file JFitK40.hh.

111  {
112  return &Rate_Hz;
113  }

◆ getModelParameters() [2/2]

Double_t* JCALIBRATE::JFitK40Parameters::getModelParameters ( )
inlineinherited

Get model parameters.

Returns
pointer to parameters

Definition at line 121 of file JFitK40.hh.

122  {
123  return &Rate_Hz;
124  }

◆ setModelParameters()

void JCALIBRATE::JFitK40Parameters::setModelParameters ( const Double_t *  data)
inlineinherited

Set model parameters.

Parameters
datapointer to parameters

Definition at line 132 of file JFitK40.hh.

133  {
134  for (Int_t i = 0; i != getNumberOfModelParameters(); ++i) {
135  getModelParameters()[i] = data[i];
136  }
137  }

◆ getModelParameter() [1/3]

Double_t JCALIBRATE::JFitK40Parameters::getModelParameter ( Int_t  i) const
inlineinherited

Get model parameter.

Parameters
iparameter index
Returns
parameter value

Definition at line 146 of file JFitK40.hh.

147  {
148  return getModelParameters()[i];
149  }

◆ getModelParameter() [2/3]

JFitParameter_t JCALIBRATE::JFitK40Parameters::getModelParameter ( Double_t JFitK40Parameters::*  p) const
inlineinherited

Get model parameter.

Parameters
ppointer to data member
Returns
parameter index and value

Definition at line 158 of file JFitK40.hh.

159  {
160  const Int_t i = &(this->*p) - getModelParameters();
161 
162  return JFitParameter_t(i, getModelParameter(i));
163  }

◆ getModelParameter() [3/3]

JFitParameter_t JCALIBRATE::JFitK40Parameters::getModelParameter ( Int_t  pmt,
Double_t JPMTParameters_t::*  p 
) const
inlineinherited

Get model parameter.

Parameters
pmtPMT number
ppointer to data member of PMT parameters
Returns
parameter index and value

Definition at line 173 of file JFitK40.hh.

174  {
175  const Int_t i = &(parameters[pmt].*p) - getModelParameters();
176 
177  return JFitParameter_t(i, getModelParameter(i));
178  }

◆ getQE()

Double_t JCALIBRATE::JFitK40Parameters::getQE ( const int  pmt) const
inlineinherited

Get QE of given PMT.

Parameters
pmtpmt address
Returns
QE

Definition at line 187 of file JFitK40.hh.

188  {
189  return parameters[pmt].QE;
190  }

◆ setQE()

void JCALIBRATE::JFitK40Parameters::setQE ( const int  pmt,
const Double_t  QE 
)
inlineinherited

Set QE of given PMT.

Parameters
pmtpmt address
QEQE

Definition at line 199 of file JFitK40.hh.

200  {
201  parameters[pmt].QE = QE;
202  }

◆ getTTS()

Double_t JCALIBRATE::JFitK40Parameters::getTTS ( const int  pmt) const
inlineinherited

Get time resolution of given PMT.

Parameters
pmtpmt address
Returns
TTS [ns]

Definition at line 211 of file JFitK40.hh.

212  {
213  return parameters[pmt].TTS;
214  }

◆ setTTS()

void JCALIBRATE::JFitK40Parameters::setTTS ( const int  pmt,
const Double_t  TTS 
)
inlineinherited

Set time resolution of given PMT.

Parameters
pmtpmt address
TTSTTS [ns]

Definition at line 223 of file JFitK40.hh.

224  {
225  parameters[pmt].TTS = TTS;
226  }

◆ setT0()

void JCALIBRATE::JFitK40Parameters::setT0 ( const int  pmt,
const Double_t  t0 
)
inlineinherited

Set time offset of given PMT.

Parameters
pmtpmt address
t0time offset [ns]

Definition at line 249 of file JFitK40.hh.

250  {
251  parameters[pmt].t0 = t0;
252  }

◆ getValue() [4/4]

Double_t JCALIBRATE::JFitK40Parameters::getValue ( const Double_t  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 261 of file JFitK40.hh.

262  {
263  return Rate_Hz * TMath::Exp(-(p1+p2+p3+p4)) * TMath::Exp(ct*(p1+ct*(p2+ct*(p3+ct*p4))));
264  }

◆ 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 147 of file JModule.hh.

150  {
151  if (first.size() == second.size()) {
152 
153  for (size_t i = 0; i != first.size(); ++i) {
154  if (first[i].getDirection().getDot(second[i].getDirection()) < 1.0 - precision) {
155  return false;
156  }
157  }
158 
159  return true;
160  }
161 
162  return false;
163  }

◆ 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 173 of file JModule.hh.

174  {
175  return at(index);
176  }

◆ 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 185 of file JModule.hh.

186  {
187  return at(index);
188  }

◆ setPMT()

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

Set PMT.

Parameters
indexreadout channel (TDC)
pmtPMT

Definition at line 197 of file JModule.hh.

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

◆ compile()

void JDETECTOR::JModule::compile ( )
inlineinherited

Compile position of module from the positions of the PMTs.

Definition at line 210 of file JModule.hh.

211  {
212  if (!empty()) {
213 
214  JPosition3D pos;
215 
216  for (iterator i = begin(); i != end(); ++i) {
217  pos.add(i->getPosition());
218  }
219 
220  static_cast<JPosition3D&>(*this) = pos.div(size());
221  }
222  }

◆ rotate() [1/5]

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

Rotate module.

Parameters
Rrotation matrix

Definition at line 230 of file JModule.hh.

231  {
232  static_cast<JPosition3D&>(*this).rotate(R);
233 
234  for (iterator i = this->begin(); i != this->end(); ++i) {
235  i->rotate(R);
236  }
237  }

◆ rotate() [2/5]

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

Rotate module.

Parameters
Qquaternion

Definition at line 292 of file JModule.hh.

293  {
294  static_cast<JPosition3D&>(*this).rotate(Q);
295 
296  for (iterator i = this->begin(); i != this->end(); ++i) {
297  i->rotate(Q);
298  }
299  }

◆ rotate() [3/5]

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

Rotate around X-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 213 of file JPosition3D.hh.

214  {
215  R.rotate(__y, __z);
216 
217  return *this;
218  }

◆ rotate() [4/5]

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

Rotate around Y-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 241 of file JPosition3D.hh.

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

◆ rotate() [5/5]

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

Rotate around Z-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 269 of file JPosition3D.hh.

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

◆ rotate_back() [1/5]

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

Rotate back module.

Parameters
Rrotation matrix

Definition at line 245 of file JModule.hh.

246  {
247  static_cast<JPosition3D&>(*this).rotate_back(R);
248 
249  for (iterator i = this->begin(); i != this->end(); ++i) {
250  i->rotate_back(R);
251  }
252  }

◆ rotate_back() [2/5]

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

Rotate back module.

Parameters
Qquaternion

Definition at line 307 of file JModule.hh.

308  {
309  static_cast<JPosition3D&>(*this).rotate_back(Q);
310 
311  for (iterator i = this->begin(); i != this->end(); ++i) {
312  i->rotate_back(Q);
313  }
314  }

◆ rotate_back() [3/5]

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

Rotate back around X-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 227 of file JPosition3D.hh.

228  {
229  R.rotate_back(__y, __z);
230 
231  return *this;
232  }

◆ rotate_back() [4/5]

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

Rotate back around Y-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 255 of file JPosition3D.hh.

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

◆ rotate_back() [5/5]

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

Rotate back around Z-axis.

Parameters
Rrotation matrix
Returns
this position

Definition at line 283 of file JPosition3D.hh.

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

◆ 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 261 of file JModule.hh.

263  {
264  static_cast<JPosition3D&>(*this).transform(R, pos);
265 
266  for (iterator i = this->begin(); i != this->end(); ++i) {
267  i->transform(R, pos);
268  }
269  }

◆ transform() [2/4]

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

Transformation of geometry.

Parameters
Ttransformation

Definition at line 277 of file JModule.hh.

278  {
279  static_cast<JPosition3D&>(*this).transform(T.getRotation(), T.getPosition());
280 
281  for (iterator i = this->begin(); i != this->end(); ++i) {
282  i->transform(T);
283  }
284  }

◆ transform() [3/4]

JVector3D& JGEOMETRY3D::JVector3D::transform
inlineinherited

Transform.

Parameters
Tmatrix
Returns
this vector

Definition at line 205 of file JVector3D.hh.

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

◆ transform() [4/4]

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

Transform.

Parameters
Tmatrix
Returns
this vector

Definition at line 205 of file JVector3D.hh.

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

◆ set() [1/2]

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

Set position.

Parameters
posposition
Returns
this module

Definition at line 323 of file JModule.hh.

324  {
325  return add(pos - static_cast<JPosition3D&>(*this));
326  }

◆ set() [2/2]

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

Set time offset.

Parameters
t0time offset [ns]
Returns
this module

Definition at line 371 of file JModule.hh.

372  {
373  for (iterator i = begin(); i != end(); ++i) {
374  i->setT0(t0);
375  }
376 
377  return *this;
378  }

◆ add() [1/2]

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

Add position.

Parameters
posposition
Returns
this module

Definition at line 335 of file JModule.hh.

336  {
337  for (iterator i = begin(); i != end(); ++i) {
338  i->add(pos);
339  }
340 
341  JVector3D::add(pos);
342 
343  return *this;
344  }

◆ add() [2/2]

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

Add time offset.

Parameters
t0time offset [ns]
Returns
this module

Definition at line 387 of file JModule.hh.

388  {
389  for (iterator i = begin(); i != end(); ++i) {
390  i->addT0(t0);
391  }
392 
393  return *this;
394  }

◆ sub() [1/2]

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

Subtract position.

Parameters
posposition
Returns
this module

Definition at line 353 of file JModule.hh.

354  {
355  for (iterator i = begin(); i != end(); ++i) {
356  i->sub(pos);
357  }
358 
359  JVector3D::sub(pos);
360 
361  return *this;
362  }

◆ sub() [2/2]

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

Subtract time offset.

Parameters
t0time offset [ns]
Returns
this module

Definition at line 403 of file JModule.hh.

404  {
405  for (iterator i = begin(); i != end(); ++i) {
406  i->subT0(t0);
407  }
408 
409  return *this;
410  }

◆ operator+=()

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

Add position.

Parameters
posposition
Returns
this module

Definition at line 419 of file JModule.hh.

420  {
421  return this->add(pos);
422  }

◆ operator-=()

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

Subtract position.

Parameters
posposition
Returns
this module

Definition at line 431 of file JModule.hh.

432  {
433  return this->sub(pos);
434  }

◆ getID()

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

Get identifier.

Returns
identifier

Definition at line 55 of file JObjectID.hh.

56  {
57  return __id;
58  }

◆ setID()

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

Set identifier.

Parameters
ididentifier

Definition at line 66 of file JObjectID.hh.

67  {
68  this->__id = id;
69  }

◆ 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 78 of file JObjectID.hh.

79  {
80  return this->getID() < object.getID();
81  }

◆ 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 90 of file JObjectID.hh.

91  {
92  return this->getID() < id;
93  }

◆ less() [3/3]

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

Less than method.

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

Definition at line 156 of file JModuleLocation.hh.

157  {
158  if (this->getString() == location.getString())
159  return this->getFloor() < location.getFloor();
160  else
161  return this->getString() < location.getString();
162  }

◆ 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 102 of file JObjectID.hh.

103  {
104  return this->getID() > id;
105  }

◆ getLocation() [1/2]

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

Get location.

Returns
location

Definition at line 68 of file JModuleLocation.hh.

69  {
70  return static_cast<const JModuleLocation&>(*this);
71  }

◆ getLocation() [2/2]

JModuleLocation& JDETECTOR::JModuleLocation::getLocation ( )
inlineinherited

Get location.

Returns
location

Definition at line 79 of file JModuleLocation.hh.

80  {
81  return static_cast<JModuleLocation&>(*this);
82  }

◆ setLocation()

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

Set location.

Parameters
locationlocation

Definition at line 90 of file JModuleLocation.hh.

91  {
92  static_cast<JModuleLocation&>(*this) = location;
93  }

◆ toString() [1/2]

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

Convert module location to string.

Returns
string

Definition at line 101 of file JModuleLocation.hh.

102  {
103  return toString("% %");
104  }

◆ toString() [2/2]

std::string JDETECTOR::JModuleLocation::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 117 of file JModuleLocation.hh.

118  {
119  JLANG::JString buffer(fmt);
120 
121  buffer.replace(target, string, 1);
122  buffer.replace(target, floor, 1);
123 
124  return buffer;
125  }

◆ getString()

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

Get string number.

Returns
string number

Definition at line 133 of file JModuleLocation.hh.

134  {
135  return string;
136  }

◆ getFloor()

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

Get floor number.

Returns
floor number

Definition at line 144 of file JModuleLocation.hh.

145  {
146  return floor;
147  }

◆ getPosition() [1/2]

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

Get position.

Returns
position

Definition at line 129 of file JPosition3D.hh.

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

◆ getPosition() [2/2]

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

Get position.

Returns
position

Definition at line 140 of file JPosition3D.hh.

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

◆ setPosition()

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

Set position.

Parameters
posposition

Definition at line 151 of file JPosition3D.hh.

152  {
153  static_cast<JVector3D&>(*this) = pos;
154  }

◆ operator JAngle3D()

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

Type conversion operator.

Returns
angle

Definition at line 162 of file JPosition3D.hh.

163  {
164  return JAngle3D(getX(), getY(), getZ());
165  }

◆ operator JVersor3D()

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

Type conversion operator.

Returns
direction

Definition at line 173 of file JPosition3D.hh.

174  {
175  return JVersor3D(getX(), getY(), getZ());
176  }

◆ 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 358 of file JPosition3D.hh.

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

◆ getDot() [2/6]

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

Get dot product.

Parameters
angleangle
Returns
dot product

Definition at line 377 of file JPosition3D.hh.

378  {
379  return
380  getX() * angle.getDX() +
381  getY() * angle.getDY() +
382  getZ() * angle.getDZ();
383  }

◆ getDot() [3/6]

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

Get dot product.

Parameters
dirdirection
Returns
dot product

Definition at line 392 of file JPosition3D.hh.

393  {
394  return
395  getX() * dir.getDX() +
396  getY() * dir.getDY() +
397  getZ() * dir.getDZ();
398  }

◆ getDot() [4/6]

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

Get dot product.

Parameters
dirdirection
Returns
dot product

Definition at line 407 of file JPosition3D.hh.

408  {
409  return
410  getX() * dir.getDX() +
411  getY() * dir.getDY() +
412  getZ() * dir.getDZ();
413  }

◆ getDot() [5/6]

double JGEOMETRY3D::JVector3D::getDot
inlineinherited

Get dot product.

Parameters
vectorvector
Returns
dot product

Definition at line 281 of file JVector3D.hh.

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

◆ getDot() [6/6]

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

Get dot product.

Parameters
vectorvector
Returns
dot product

Definition at line 281 of file JVector3D.hh.

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

◆ operator JVector2D()

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

Type conversion operator.

Returns
JVector2D

Definition at line 82 of file JVector3D.hh.

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

◆ getX()

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

Get x position.

Returns
x position

Definition at line 93 of file JVector3D.hh.

94  {
95  return __x;
96  }

◆ getY()

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

Get y position.

Returns
y position

Definition at line 103 of file JVector3D.hh.

104  {
105  return __y;
106  }

◆ getZ()

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

Get z position.

Returns
z position

Definition at line 114 of file JVector3D.hh.

115  {
116  return __z;
117  }

◆ negate()

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

Negate vector.

Returns
this vector

Definition at line 125 of file JVector3D.hh.

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

◆ mul() [1/2]

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

Scale vector.

Parameters
factormultiplication factor
Returns
this vector

Definition at line 173 of file JVector3D.hh.

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

◆ mul() [2/2]

JVector3D & JMATH::JMath< JVector3D , JSecond_t >::mul ( const JSecond_t &  object)
inlineinherited

Multiply with object.

Parameters
objectobject
Returns
result object

Definition at line 273 of file JMath.hh.

274  {
275  return static_cast<JFirst_t&>(*this) = JCalculator<JFirst_t>::calculator.mul(static_cast<const JFirst_t&>(*this), object);
276  }

◆ div()

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

Scale vector.

Parameters
factordivision factor
Returns
this vector

Definition at line 189 of file JVector3D.hh.

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

◆ equals()

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 220 of file JVector3D.hh.

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

◆ getLengthSquared()

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

Get length squared.

Returns
square of length

Definition at line 234 of file JVector3D.hh.

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

◆ getLength()

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

Get length.

Returns
length

Definition at line 245 of file JVector3D.hh.

246  {
247  return sqrt(getLengthSquared());
248  }

◆ getDistanceSquared()

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

Get squared of distance to point.

Parameters
posposition
Returns
square of distance

Definition at line 257 of file JVector3D.hh.

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

◆ getDistance()

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

Get distance to point.

Parameters
posposition
Returns
distance

Definition at line 269 of file JVector3D.hh.

270  {
271  return sqrt(getDistanceSquared(pos));
272  }

◆ cross()

JVector3D& JGEOMETRY3D::JVector3D::cross ( 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 298 of file JVector3D.hh.

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

◆ configure()

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

Configure.

Parameters
numberOfIndicesnumber of indices

Definition at line 54 of file JCombinatorics.hh.

55  {
56  zbuf1D.clear();
57 
58  zbuf2D.resize(numberOfIndices);
59 
60  for (int i = 0; i != numberOfIndices; ++i) {
61  zbuf2D[i].resize(numberOfIndices);
62  }
63 
64  for (int i = 0; i != numberOfIndices; ++i) {
65 
66  zbuf2D[i][i] = -1;
67 
68  for (int j = i; ++j != numberOfIndices; ) {
69 
70  zbuf2D[i][j] = zbuf1D.size();
71  zbuf2D[j][i] = zbuf1D.size();
72 
73  zbuf1D.push_back(pair_type(i,j));
74  }
75  }
76  }

◆ getNumberOfIndices()

unsigned int JTOOLS::JCombinatorics::getNumberOfIndices ( ) const
inlineinherited

Get number of indices.

Returns
number of indices

Definition at line 84 of file JCombinatorics.hh.

85  {
86  return zbuf2D.size();
87  }

◆ getNumberOfPairs()

unsigned int JTOOLS::JCombinatorics::getNumberOfPairs ( ) const
inlineinherited

Get number of pairs.

Returns
number of pairs

Definition at line 95 of file JCombinatorics.hh.

96  {
97  return zbuf1D.size();
98  }

◆ getIndex()

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 108 of file JCombinatorics.hh.

109  {
110  return zbuf2D[first][second];
111  }

◆ 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 120 of file JCombinatorics.hh.

121  {
122  return zbuf1D[index];
123  }

◆ sort()

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

Sort address pairs.

Parameters
comparatorcomparator for pairs

Definition at line 132 of file JCombinatorics.hh.

133  {
134  std::stable_sort(zbuf1D.begin(), zbuf1D.end(), comparator);
135 
136  for (int i = 0; i != (int) zbuf1D.size(); ++i) {
137 
138  const pair_type pair = zbuf1D[i];
139 
140  zbuf2D[pair.first][pair.second] = i;
141  zbuf2D[pair.second][pair.first] = i;
142  }
143 
144  for (int i = 0; i != (int) zbuf2D.size(); ++i) {
145  zbuf2D[i][i] = -1;
146  }
147  }

◆ 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 157 of file JCombinatorics.hh.

158  {
159  return (second >= first ? +1 : -1);
160  }

◆ 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 169 of file JCombinatorics.hh.

170  {
171  return getSign(pair.first, pair.second);
172  }

Member Data Documentation

◆ sigmaK40_ns

Double_t JCALIBRATE::JFitK40::sigmaK40_ns
private

intrinsic K40 arrival time spread [ns]

Definition at line 650 of file JFitK40.hh.

◆ index_of_average_t0

int JCALIBRATE::JFitK40::index_of_average_t0
private

index of t0 used for average time offset

Definition at line 651 of file JFitK40.hh.

◆ disable

bool JCALIBRATE::JFitK40::disable[NUMBER_OF_PMTS]
private

disable PMT from fit

Definition at line 652 of file JFitK40.hh.

◆ Rate_Hz

Double_t JCALIBRATE::JFitK40Parameters::Rate_Hz
inherited

maximal coincidence rate [Hz]

Definition at line 269 of file JFitK40.hh.

◆ p1

Double_t JCALIBRATE::JFitK40Parameters::p1
inherited

angle dependence coincidence rate

Definition at line 270 of file JFitK40.hh.

◆ p2

Double_t JCALIBRATE::JFitK40Parameters::p2
inherited

angle dependence coincidence rate

Definition at line 271 of file JFitK40.hh.

◆ p3

Double_t JCALIBRATE::JFitK40Parameters::p3
inherited

angle dependence coincidence rate

Definition at line 272 of file JFitK40.hh.

◆ p4

Double_t JCALIBRATE::JFitK40Parameters::p4
inherited

angle dependence coincidence rate

Definition at line 273 of file JFitK40.hh.

◆ bg

Double_t JCALIBRATE::JFitK40Parameters::bg
inherited

remaining constant background

Definition at line 274 of file JFitK40.hh.

◆ cc

Double_t JCALIBRATE::JFitK40Parameters::cc
inherited

fraction of signal correlated background

Definition at line 275 of file JFitK40.hh.

◆ parameters

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

Definition at line 277 of file JFitK40.hh.

◆ __id

int JLANG::JObjectID::__id
protectedinherited

Definition at line 160 of file JObjectID.hh.

◆ string

int JDETECTOR::JModuleLocation::string
protectedinherited

Definition at line 233 of file JModuleLocation.hh.

◆ floor

int JDETECTOR::JModuleLocation::floor
protectedinherited

Definition at line 234 of file JModuleLocation.hh.

◆ __x

double JGEOMETRY3D::JVector3D::__x
protectedinherited

Definition at line 309 of file JVector3D.hh.

◆ __y

double JGEOMETRY3D::JVector3D::__y
protectedinherited

Definition at line 310 of file JVector3D.hh.

◆ __z

double JGEOMETRY3D::JVector3D::__z
protectedinherited

Definition at line 311 of file JVector3D.hh.

◆ zbuf1D

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

Definition at line 176 of file JCombinatorics.hh.

◆ zbuf2D

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

Definition at line 177 of file JCombinatorics.hh.


The documentation for this struct was generated from the following file:
JGEOMETRY3D::JVersor3D::getDZ
double getDZ() const
Get z direction.
Definition: JVersor3D.hh:114
JCALIBRATE::JPMTParameters_t::QE
Double_t QE
quantum efficiency [unit]
Definition: JFitK40.hh:53
JCALIBRATE::JFitK40::getDot
double getDot(const pair_type &pair) const
Get cosine of space angle between PMT axes.
Definition: JFitK40.hh:472
JGEOMETRY2D::JRotation2D::rotate_back
void rotate_back(double &__x, double &__y) const
Rotate back.
Definition: JRotation2D.hh:112
std::iterator
Definition: JSTDTypes.hh:18
JCALIBRATE::JFitK40Parameters::bg
Double_t bg
remaining constant background
Definition: JFitK40.hh:274
JCALIBRATE::JFitK40::getRate
static Double_t getRate(const Double_t *x, const Double_t *data)
Get K40 coincidence rate as a function of the fit parameters.
Definition: JFitK40.hh:623
JDETECTOR::JModule::sub
JModule & sub(const JVector3D &pos)
Subtract position.
Definition: JModule.hh:353
JROOT::isParameterFixed
bool isParameterFixed(TF1 &f1, const Int_t index)
Check if fit parameter is fixed.
Definition: JRootToolkit.hh:319
JGEOMETRY3D::JVector3D::JVector3D
JVector3D()
Default constructor.
Definition: JVector3D.hh:40
JTOOLS::JCombinatorics::getSign
static int getSign(const int first, const int second)
Sign of pair of indices.
Definition: JCombinatorics.hh:157
JCALIBRATE::JFitK40Parameters::p2
Double_t p2
angle dependence coincidence rate
Definition: JFitK40.hh:271
JCALIBRATE::JFitK40::disable
bool disable[NUMBER_OF_PMTS]
disable PMT from fit
Definition: JFitK40.hh:652
JGEOMETRY3D::JVector3D::__y
double __y
Definition: JVector3D.hh:310
JCALIBRATE::JFitK40Parameters::setT0
void setT0(const int pmt, const Double_t t0)
Set time offset of given PMT.
Definition: JFitK40.hh:249
JGEOMETRY3D::JVector3D::getLengthSquared
double getLengthSquared() const
Get length squared.
Definition: JVector3D.hh:234
JDETECTOR::JModuleLocation::floor
int floor
Definition: JModuleLocation.hh:234
JTOOLS::JCombinatorics::configure
void configure(const int numberOfIndices)
Configure.
Definition: JCombinatorics.hh:54
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition: JVector3D.hh:114
JCALIBRATE::JFitK40::getSigma
Double_t getSigma(const pair_type &pair) const
Get time resolution of given PMT pair.
Definition: JFitK40.hh:485
JCALIBRATE::FITK40_QE_MAX
static const double FITK40_QE_MAX
Maximal quantum efficiency [unit].
Definition: JFitK40.hh:34
JCALIBRATE::JFitK40Parameters::setModelParameters
void setModelParameters(const Double_t *data)
Set model parameters.
Definition: JFitK40.hh:132
JCALIBRATE::JFitK40Parameters::getQE
Double_t getQE(const int pmt) const
Get QE of given PMT.
Definition: JFitK40.hh:187
JGEOMETRY3D::JDirection3D::getDot
double getDot(const JAngle3D &angle) const
Get dot product.
Definition: JDirection3D.hh:333
JGEOMETRY3D::JPosition3D::JPosition3D
JPosition3D()
Default constructor.
Definition: JPosition3D.hh:47
JTOOLS::j
int j
Definition: JPolint.hh:634
JGEOMETRY3D::JVersor3D
Data structure for normalised vector in three dimensions.
Definition: JVersor3D.hh:23
KM3NETDAQ::NUMBER_OF_PMTS
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition: JDAQ.hh:26
JCALIBRATE::JPMTParameters_t::t0
Double_t t0
time offset [ns]
Definition: JFitK40.hh:55
JCALIBRATE::JFitK40::getValue
Double_t getValue(const pair_type &pair) const
Get K40 coincidence rate.
Definition: JFitK40.hh:529
JCALIBRATE::JFitK40Parameters::getModelParameter
Double_t getModelParameter(Int_t i) const
Get model parameter.
Definition: JFitK40.hh:146
JCALIBRATE::JFitK40Parameters::parameters
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
JCALIBRATE::JFitK40Parameters::p3
Double_t p3
angle dependence coincidence rate
Definition: JFitK40.hh:272
JCALIBRATE::FITK40_TTS_MAX_NS
static const double FITK40_TTS_MAX_NS
Maximal transition-time spread [ns].
Definition: JFitK40.hh:36
JTOOLS::JCombinatorics::getPair
const pair_type & getPair(const int index) const
Get pair of indices for given index.
Definition: JCombinatorics.hh:120
JTOOLS::JCombinatorics::zbuf2D
std::vector< std::vector< int > > zbuf2D
Definition: JCombinatorics.hh:177
JGEOMETRY2D::JRotation2D::rotate
void rotate(double &__x, double &__y) const
Rotate.
Definition: JRotation2D.hh:96
JDETECTOR::JModule::getPMT
const JPMT & getPMT(const int index) const
Get PMT.
Definition: JModule.hh:173
JDETECTOR::JModuleLocation::string
int string
Definition: JModuleLocation.hh:233
JLANG::JString
Wrapper class around STL string class.
Definition: JString.hh:28
JCALIBRATE::JPMTParameters_t::TTS
Double_t TTS
transition-time spread [ns]
Definition: JFitK40.hh:54
JLANG::JObjectID::getID
int getID() const
Get identifier.
Definition: JObjectID.hh:55
JTOOLS::result
return result
Definition: JPolint.hh:695
JGEOMETRY3D::JPosition3D::rotate_back
JPosition3D & rotate_back(const JRotation3D &R)
Rotate back.
Definition: JPosition3D.hh:199
JCALIBRATE::JFitK40Parameters::Rate_Hz
Double_t Rate_Hz
maximal coincidence rate [Hz]
Definition: JFitK40.hh:269
JCALIBRATE::JFitK40::getSigmaK40
Double_t getSigmaK40() const
Get intrinsic K40 arrival time spread.
Definition: JFitK40.hh:343
THROW
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
Definition: JException.hh:670
JGEOMETRY3D::JAngle3D::getDX
double getDX() const
Get x direction.
Definition: JAngle3D.hh:106
JGEOMETRY3D::JVersor3D::getDX
double getDX() const
Get x direction.
Definition: JVersor3D.hh:92
JCALIBRATE::JFitK40Parameters::p1
Double_t p1
angle dependence coincidence rate
Definition: JFitK40.hh:270
JCALIBRATE::JFitK40Parameters::getTTS
Double_t getTTS(const int pmt) const
Get time resolution of given PMT.
Definition: JFitK40.hh:211
JCALIBRATE::JFitK40Parameters::getNumberOfModelParameters
static Int_t getNumberOfModelParameters()
Get number of model parameters.
Definition: JFitK40.hh:99
JGEOMETRY3D::JVersor3D::getDY
double getDY() const
Get y direction.
Definition: JVersor3D.hh:103
std::pair
Definition: JSTDTypes.hh:15
JCALIBRATE::JFitK40
Parametrisation of two-fold coincidence rate due to K40 and other radioactive decays.
Definition: JFitK40.hh:286
JCALIBRATE::JPMTParameters_t
Fit parameters for single PMT.
Definition: JFitK40.hh:42
JDETECTOR::JModuleLocation::toString
std::string toString() const
Convert module location to string.
Definition: JModuleLocation.hh:101
JROOT::setParLimits
bool setParLimits(TF1 &f1, const Int_t index, Double_t xmin, Double_t xmax)
Set fit parameter limits.
Definition: JRootToolkit.hh:293
JAANET::getDirection
JDirection3D getDirection(const Vec &v)
Get direction.
Definition: JAAnetToolkit.hh:221
JTOOLS::JCombinatorics::zbuf1D
std::vector< pair_type > zbuf1D
Definition: JCombinatorics.hh:176
JGEOMETRY3D::JAngle3D
Data structure for angles in three dimensions.
Definition: JAngle3D.hh:31
JGEOMETRY3D::JAngle3D::getDZ
double getDZ() const
Get z direction.
Definition: JAngle3D.hh:128
JCALIBRATE::JFitK40Parameters::p4
Double_t p4
angle dependence coincidence rate
Definition: JFitK40.hh:273
JLANG::JObjectID::__id
int __id
Definition: JObjectID.hh:160
JGEOMETRY3D::JVector3D::getDistanceSquared
double getDistanceSquared(const JVector3D &pos) const
Get squared of distance to point.
Definition: JVector3D.hh:257
JROOT::fixParameter
bool fixParameter(TF1 &f1, const JFitParameter_t &parameter)
Fix fit parameter.
Definition: JRootToolkit.hh:249
JCALIBRATE::JFitK40::getInstance
static JFitK40 & getInstance()
Get unique instance of fit object.
Definition: JFitK40.hh:666
JCALIBRATE::JFitK40Parameters::JFitK40Parameters
JFitK40Parameters()
Default constructor.
Definition: JFitK40.hh:70
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition: JVector3D.hh:103
JGEOMETRY3D::JAngle3D::getDY
double getDY() const
Get y direction.
Definition: JAngle3D.hh:117
JCALIBRATE::JFitK40::index_of_average_t0
int index_of_average_t0
index of t0 used for average time offset
Definition: JFitK40.hh:651
JDETECTOR::JModule::add
JModule & add(const JVector3D &pos)
Add position.
Definition: JModule.hh:335
JGEOMETRY3D::JVersor3Z::getDY
double getDY() const
Get y direction.
Definition: JVersor3Z.hh:156
JTOOLS::JCombinatorics::pair_type
std::pair< int, int > pair_type
Data structure for a pair of indices.
Definition: JCombinatorics.hh:28
JGEOMETRY3D::JVector3D::__z
double __z
Definition: JVector3D.hh:311
JMATH::JCalculator
Auxiliary class for arithmetic operations on objects.
Definition: JCalculator.hh:18
JCALIBRATE::JFitK40Parameters::getModelParameters
const Double_t * getModelParameters() const
Get model parameters.
Definition: JFitK40.hh:110
JGEOMETRY3D::JVector3D::__x
double __x
Definition: JVector3D.hh:309
JGEOMETRY3D::JVector3D::add
JVector3D & add(const JVector3D &vector)
Add vector.
Definition: JVector3D.hh:141
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition: JVector3D.hh:93
JCALIBRATE::JFitK40::sigmaK40_ns
Double_t sigmaK40_ns
intrinsic K40 arrival time spread [ns]
Definition: JFitK40.hh:650
JCALIBRATE::JPairwiseComparator
Auxiliary class to sort pairs of PMT addresses within optical module.
Definition: JCalibrateK40.hh:45
JDETECTOR::JModuleLocation::getFloor
int getFloor() const
Get floor number.
Definition: JModuleLocation.hh:144
JGEOMETRY3D::JVersor3Z::getDZ
double getDZ() const
Get z direction.
Definition: JVersor3Z.hh:167
JCALIBRATE::JFitK40::getT0
Double_t getT0(const int pmt) const
Get time offset of given PMT.
Definition: JFitK40.hh:393
JGEOMETRY3D::JVector3D::sub
JVector3D & sub(const JVector3D &vector)
Subtract vector.
Definition: JVector3D.hh:157
JCALIBRATE::FITK40_QE_MIN
static const double FITK40_QE_MIN
Minimal quantum efficiency [unit].
Definition: JFitK40.hh:33
JMATH::getDot
double getDot(const JFirst_t &first, const JSecond_t &second)
Get dot product of objects.
Definition: JMathToolkit.hh:131
JTOOLS::JCombinatorics::sort
void sort(JComparator_t comparator)
Sort address pairs.
Definition: JCombinatorics.hh:132
JDETECTOR::JModuleLocation::getString
int getString() const
Get string number.
Definition: JModuleLocation.hh:133
JCALIBRATE::FITK40_TTS_MIN_NS
static const double FITK40_TTS_MIN_NS
Minimal transition-time spread [ns].
Definition: JFitK40.hh:35
JGEOMETRY3D::JVersor3Z::getDX
double getDX() const
Get x direction.
Definition: JVersor3Z.hh:145
JCALIBRATE::JFitK40Parameters::cc
Double_t cc
fraction of signal correlated background
Definition: JFitK40.hh:275
JCALIBRATE::JFitK40Parameters::getT0
Double_t getT0(const int pmt) const
Get time offset of given PMT.
Definition: JFitK40.hh:237