Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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::JLocation 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...
 
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 and directions 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...
 
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 JLocation &location) const
 Less than method. More...
 
bool more (const int id) const
 More than method. More...
 
const JLocationgetLocation () const
 Get location. More...
 
JLocationgetLocation ()
 Get location. More...
 
void setLocation (const JLocation &location)
 Set location. More...
 
std::string toString () const
 Convert module location to string. More...
 
std::string toString (const std::string &fmt, const std::string target="%") const
 Convert module loation to string. More...
 
int getString () const
 Get string number. More...
 
int getFloor () const
 Get floor number. More...
 
const JPosition3DgetPosition () const
 Get position. More...
 
JPosition3DgetPosition ()
 Get position. More...
 
void setPosition (const JVector3D &pos)
 Set position. More...
 
 operator JAngle3D () const
 Type conversion operator. More...
 
 operator JVersor3D () const
 Type conversion operator. More...
 
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...
 
 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 JNullType &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

Data structure for a pair of addresses.

Definition at line 297 of file JFitK40.hh.

Constructor & Destructor Documentation

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  }
void configure(const int numberOfIndices)
Configure.
Double_t sigmaK40_ns
intrinsic K40 arrival time spread [ns]
Definition: JFitK40.hh:650
bool disable[NUMBER_OF_PMTS]
disable PMT from fit
Definition: JFitK40.hh:652
void sort(JComparator_t comparator)
Sort address pairs.
JModule()
Default constructor.
Definition: JModule.hh:60
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
Auxiliary class to sort pairs of PMT addresses within optical module.
int index_of_average_t0
index of t0 used for average time offset
Definition: JFitK40.hh:651
static Int_t getNumberOfModelParameters()
Get number of model parameters.
Definition: JFitK40.hh:99
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition: JDAQ.hh:26
JCALIBRATE::JFitK40::JFitK40 ( )
inlineprivate

Default constructor.

Definition at line 657 of file JFitK40.hh.

658  {}

Member Function Documentation

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  }
Double_t sigmaK40_ns
intrinsic K40 arrival time spread [ns]
Definition: JFitK40.hh:650
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  }
Double_t sigmaK40_ns
intrinsic K40 arrival time spread [ns]
Definition: JFitK40.hh:650
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  }
bool disable[NUMBER_OF_PMTS]
disable PMT from fit
Definition: JFitK40.hh:652
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
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  }
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
int index_of_average_t0
index of t0 used for average time offset
Definition: JFitK40.hh:651
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 
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  }
Double_t t0
time offset [ns]
Definition: JFitK40.hh:55
bool disable[NUMBER_OF_PMTS]
disable PMT from fit
Definition: JFitK40.hh:652
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
Double_t getT0(const int pmt) const
Get time offset of given PMT.
Definition: JFitK40.hh:237
int index_of_average_t0
index of t0 used for average time offset
Definition: JFitK40.hh:651
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition: JDAQ.hh:26
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  }
bool disable[NUMBER_OF_PMTS]
disable PMT from fit
Definition: JFitK40.hh:652
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
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  }
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
Definition: JException.hh:670
bool disable[NUMBER_OF_PMTS]
disable PMT from fit
Definition: JFitK40.hh:652
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
int index_of_average_t0
index of t0 used for average time offset
Definition: JFitK40.hh:651
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition: JDAQ.hh:26
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  }
double getDot(const JFirst_t &first, const JSecond_t &second)
Get dot product of objects.
const JPMT & getPMT(const int index) const
Get PMT.
Definition: JModule.hh:174
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  }
Double_t getSigmaK40() const
Get intrinsic K40 arrival time spread.
Definition: JFitK40.hh:343
Double_t getTTS(const int pmt) const
Get time resolution of given PMT.
Definition: JFitK40.hh:211
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  }
Double_t t0
time offset [ns]
Definition: JFitK40.hh:55
int index_of_average_t0
index of t0 used for average time offset
Definition: JFitK40.hh:651
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
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  }
Double_t getQE(const int pmt) const
Get QE of given PMT.
Definition: JFitK40.hh:187
Double_t getValue(const pair_type &pair) const
Get K40 coincidence rate.
Definition: JFitK40.hh:529
double getDot(const pair_type &pair) const
Get cosine of space angle between PMT axes.
Definition: JFitK40.hh:472
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  }
Double_t getT0(const int pmt) const
Get time offset of given PMT.
Definition: JFitK40.hh:393
Double_t getValue(const pair_type &pair) const
Get K40 coincidence rate.
Definition: JFitK40.hh:529
Double_t getSigma(const pair_type &pair) const
Get time resolution of given PMT pair.
Definition: JFitK40.hh:485
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 
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  }
Double_t TTS
transition-time spread [ns]
Definition: JFitK40.hh:54
void setT0(const int pmt, const Double_t t0)
Set time offset of given PMT.
Definition: JFitK40.hh:249
Double_t t0
time offset [ns]
Definition: JFitK40.hh:55
Double_t QE
quantum efficiency [unit]
Definition: JFitK40.hh:53
bool disable[NUMBER_OF_PMTS]
disable PMT from fit
Definition: JFitK40.hh:652
static const double FITK40_QE_MAX
Maximal quantum efficiency [unit].
Definition: JFitK40.hh:34
bool isParameterFixed(TF1 &f1, const Int_t index)
Check if fit parameter is fixed.
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
const Double_t * getModelParameters() const
Get model parameters.
Definition: JFitK40.hh:110
void setModelParameters(const Double_t *data)
Set model parameters.
Definition: JFitK40.hh:132
return result
Definition: JPolint.hh:695
Double_t getModelParameter(Int_t i) const
Get model parameter.
Definition: JFitK40.hh:146
int index_of_average_t0
index of t0 used for average time offset
Definition: JFitK40.hh:651
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
static JFitK40 & getInstance()
Get unique instance of fit object.
Definition: JFitK40.hh:666
bool setParLimits(TF1 &f1, const Int_t index, Double_t xmin, Double_t xmax)
Set fit parameter limits.
Fit parameters for single PMT.
Definition: JFitK40.hh:42
static const double FITK40_QE_MIN
Minimal quantum efficiency [unit].
Definition: JFitK40.hh:33
static const double FITK40_TTS_MIN_NS
Minimal transition-time spread [ns].
Definition: JFitK40.hh:35
bool fixParameter(TF1 &f1, const JFitParameter_t &parameter)
Fix fit parameter.
static const double FITK40_TTS_MAX_NS
Maximal transition-time spread [ns].
Definition: JFitK40.hh:36
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition: JDAQ.hh:26
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  }
Double_t bg
remaining constant background
Definition: JFitK40.hh:274
const pair_type & getPair(const int index) const
Get pair of indices for given index.
Double_t getT0(const int pmt) const
Get time offset of given PMT.
Definition: JFitK40.hh:393
then set_variable singlesRate set_variable doublesRate set_variable numberOfSlices echo Generating random background echo Singles rate
void setModelParameters(const Double_t *data)
Set model parameters.
Definition: JFitK40.hh:132
static JFitK40 & getInstance()
Get unique instance of fit object.
Definition: JFitK40.hh:666
Double_t getValue(const pair_type &pair) const
Get K40 coincidence rate.
Definition: JFitK40.hh:529
Double_t getSigma(const pair_type &pair) const
Get time resolution of given PMT pair.
Definition: JFitK40.hh:485
Double_t cc
fraction of signal correlated background
Definition: JFitK40.hh:275
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  }
Parametrisation of two-fold coincidence rate due to K40 and other radioactive decays.
Definition: JFitK40.hh:286
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  }
JFitK40Parameters()
Default constructor.
Definition: JFitK40.hh:70
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  }
Double_t Rate_Hz
maximal coincidence rate [Hz]
Definition: JFitK40.hh:269
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  }
Double_t Rate_Hz
maximal coincidence rate [Hz]
Definition: JFitK40.hh:269
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  }
const Double_t * getModelParameters() const
Get model parameters.
Definition: JFitK40.hh:110
static Int_t getNumberOfModelParameters()
Get number of model parameters.
Definition: JFitK40.hh:99
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  }
const Double_t * getModelParameters() const
Get model parameters.
Definition: JFitK40.hh:110
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  }
const Double_t * getModelParameters() const
Get model parameters.
Definition: JFitK40.hh:110
Double_t getModelParameter(Int_t i) const
Get model parameter.
Definition: JFitK40.hh:146
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  }
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
const Double_t * getModelParameters() const
Get model parameters.
Definition: JFitK40.hh:110
Double_t getModelParameter(Int_t i) const
Get model parameter.
Definition: JFitK40.hh:146
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
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  }
Double_t QE
quantum efficiency [unit]
Definition: JFitK40.hh:53
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
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  }
Double_t QE
quantum efficiency [unit]
Definition: JFitK40.hh:53
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
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  }
Double_t TTS
transition-time spread [ns]
Definition: JFitK40.hh:54
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
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  }
Double_t TTS
transition-time spread [ns]
Definition: JFitK40.hh:54
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
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  }
Double_t t0
time offset [ns]
Definition: JFitK40.hh:55
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
JPMTParameters_t parameters[NUMBER_OF_PMTS]
Definition: JFitK40.hh:277
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  }
Double_t p3
angle dependence coincidence rate
Definition: JFitK40.hh:272
Double_t p4
angle dependence coincidence rate
Definition: JFitK40.hh:273
Double_t Rate_Hz
maximal coincidence rate [Hz]
Definition: JFitK40.hh:269
Double_t p1
angle dependence coincidence rate
Definition: JFitK40.hh:270
Double_t p2
angle dependence coincidence rate
Definition: JFitK40.hh:271
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 148 of file JModule.hh.

151  {
152  if (first.size() == second.size()) {
153 
154  for (size_t i = 0; i != first.size(); ++i) {
155  if (first[i].getDirection().getDot(second[i].getDirection()) < 1.0 - precision) {
156  return false;
157  }
158  }
159 
160  return true;
161  }
162 
163  return false;
164  }
double getDot(const JAngle3D &angle) const
Get dot product.
JDirection3D getDirection(const Vec &v)
Get direction.
const JPMT& JDETECTOR::JModule::getPMT ( const int  index) const
inlineinherited

Get PMT.

Parameters
indexreadout channel (TDC)
Returns
PMT at given index

Definition at line 174 of file JModule.hh.

175  {
176  return at(index);
177  }
JPMT& JDETECTOR::JModule::getPMT ( const int  index)
inlineinherited

Get PMT.

Parameters
indexreadout channel (TDC)
Returns
PMT at given index

Definition at line 186 of file JModule.hh.

187  {
188  return at(index);
189  }
void JDETECTOR::JModule::setPMT ( const int  index,
const JPMT pmt 
)
inlineinherited

Set PMT.

Parameters
indexreadout channel (TDC)
pmtPMT

Definition at line 198 of file JModule.hh.

199  {
200  if (index >= (int) size()) {
201  resize(index + 1);
202  }
203 
204  (*this)[index] = pmt;
205  }
esac $JPP_DIR examples JDetector JTransitTime o $OUTPUT_FILE n N $NPE T $TTS_NS d $DEBUG for HISTOGRAM in tts tt2 pmt
Definition: JTransitTime.sh:36
void JDETECTOR::JModule::compile ( )
inlineinherited

Compile position of module from the positions and directions of the PMTs.

Definition at line 211 of file JModule.hh.

212  {
213  using namespace std;
214  using namespace JPP;
215 
216  if (!empty()) {
217 
218  JPosition3D& pos = getPosition();
219 
220  try {
221  pos = JEstimator<JPoint3D>(this->begin(), this->end());
222  }
223  catch(const exception&) {
224 
225  pos = JPosition3D(0.0, 0.0, 0.0);
226 
227  for (iterator i = begin(); i != end(); ++i) {
228  pos.add(i->getPosition());
229  }
230 
231  pos.div(size());
232  }
233  }
234  }
const JPosition3D & getPosition() const
Get position.
Definition: JPosition3D.hh:129
JPosition3D()
Default constructor.
Definition: JPosition3D.hh:47
void JDETECTOR::JModule::rotate ( const JRotation3D R)
inlineinherited

Rotate module.

Parameters
Rrotation matrix

Definition at line 242 of file JModule.hh.

243  {
245 
246  for (iterator i = this->begin(); i != this->end(); ++i) {
247  i->rotate(R);
248  }
249  }
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:40
JPosition3D & rotate(const JRotation3D &R)
Rotate.
Definition: JPosition3D.hh:185
void JDETECTOR::JModule::rotate ( const JQuaternion3D Q)
inlineinherited

Rotate module.

Parameters
Qquaternion

Definition at line 304 of file JModule.hh.

305  {
307 
308  for (iterator i = this->begin(); i != this->end(); ++i) {
309  i->rotate(Q);
310  }
311  }
JPosition3D & rotate(const JRotation3D &R)
Rotate.
Definition: JPosition3D.hh:185
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  }
void rotate(double &__x, double &__y) const
Rotate.
Definition: JRotation2D.hh:96
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  }
void rotate(double &__x, double &__y) const
Rotate.
Definition: JRotation2D.hh:96
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  }
void rotate(double &__x, double &__y) const
Rotate.
Definition: JRotation2D.hh:96
void JDETECTOR::JModule::rotate_back ( const JRotation3D R)
inlineinherited

Rotate back module.

Parameters
Rrotation matrix

Definition at line 257 of file JModule.hh.

258  {
260 
261  for (iterator i = this->begin(); i != this->end(); ++i) {
262  i->rotate_back(R);
263  }
264  }
JPosition3D & rotate_back(const JRotation3D &R)
Rotate back.
Definition: JPosition3D.hh:199
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:40
void JDETECTOR::JModule::rotate_back ( const JQuaternion3D Q)
inlineinherited

Rotate back module.

Parameters
Qquaternion

Definition at line 319 of file JModule.hh.

320  {
322 
323  for (iterator i = this->begin(); i != this->end(); ++i) {
324  i->rotate_back(Q);
325  }
326  }
JPosition3D & rotate_back(const JRotation3D &R)
Rotate back.
Definition: JPosition3D.hh:199
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  }
void rotate_back(double &__x, double &__y) const
Rotate back.
Definition: JRotation2D.hh:112
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  }
void rotate_back(double &__x, double &__y) const
Rotate back.
Definition: JRotation2D.hh:112
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  }
void rotate_back(double &__x, double &__y) const
Rotate back.
Definition: JRotation2D.hh:112
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 273 of file JModule.hh.

275  {
277 
278  for (iterator i = this->begin(); i != this->end(); ++i) {
279  i->transform(R, pos);
280  }
281  }
void transform(const JRotation3D &R, const JVector3D &pos)
Transform position.
Definition: JPosition3D.hh:330
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:40
void JDETECTOR::JModule::transform ( const JTransformation3D T)
inlineinherited

Transformation of geometry.

Parameters
Ttransformation

Definition at line 289 of file JModule.hh.

290  {
291  JPosition3D::transform(T.getRotation(), T.getPosition());
292 
293  for (iterator i = this->begin(); i != this->end(); ++i) {
294  i->transform(T);
295  }
296  }
void transform(const JRotation3D &R, const JVector3D &pos)
Transform position.
Definition: JPosition3D.hh:330
do set_variable OUTPUT_DIRECTORY $WORKDIR T
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  }
do set_variable OUTPUT_DIRECTORY $WORKDIR T
JModule& JDETECTOR::JModule::set ( const JVector3D pos)
inlineinherited

Set position.

Parameters
posposition
Returns
this module

Definition at line 335 of file JModule.hh.

336  {
337  return add(pos - static_cast<JPosition3D&>(*this));
338  }
JModule & add(const JVector3D &pos)
Add position.
Definition: JModule.hh:347
JModule& JDETECTOR::JModule::set ( const double  t0)
inlineinherited

Set time offset.

Parameters
t0time offset [ns]
Returns
this module

Definition at line 383 of file JModule.hh.

384  {
385  for (iterator i = begin(); i != end(); ++i) {
386  i->setT0(t0);
387  }
388 
389  return *this;
390  }
JModule& JDETECTOR::JModule::add ( const JVector3D pos)
inlineinherited

Add position.

Parameters
posposition
Returns
this module

Definition at line 347 of file JModule.hh.

348  {
349  for (iterator i = begin(); i != end(); ++i) {
350  i->add(pos);
351  }
352 
353  JPosition3D::add(pos);
354 
355  return *this;
356  }
JVector3D & add(const JVector3D &vector)
Add vector.
Definition: JVector3D.hh:141
JModule& JDETECTOR::JModule::add ( const double  t0)
inlineinherited

Add time offset.

Parameters
t0time offset [ns]
Returns
this module

Definition at line 399 of file JModule.hh.

400  {
401  for (iterator i = begin(); i != end(); ++i) {
402  i->addT0(t0);
403  }
404 
405  return *this;
406  }
JModule& JDETECTOR::JModule::sub ( const JVector3D pos)
inlineinherited

Subtract position.

Parameters
posposition
Returns
this module

Definition at line 365 of file JModule.hh.

366  {
367  for (iterator i = begin(); i != end(); ++i) {
368  i->sub(pos);
369  }
370 
371  JPosition3D::sub(pos);
372 
373  return *this;
374  }
JVector3D & sub(const JVector3D &vector)
Subtract vector.
Definition: JVector3D.hh:157
JModule& JDETECTOR::JModule::sub ( const double  t0)
inlineinherited

Subtract time offset.

Parameters
t0time offset [ns]
Returns
this module

Definition at line 415 of file JModule.hh.

416  {
417  for (iterator i = begin(); i != end(); ++i) {
418  i->subT0(t0);
419  }
420 
421  return *this;
422  }
JModule& JDETECTOR::JModule::operator+= ( const JVector3D pos)
inlineinherited

Add position.

Parameters
posposition
Returns
this module

Definition at line 431 of file JModule.hh.

432  {
433  return this->add(pos);
434  }
JModule & add(const JVector3D &pos)
Add position.
Definition: JModule.hh:347
JModule& JDETECTOR::JModule::operator-= ( const JVector3D pos)
inlineinherited

Subtract position.

Parameters
posposition
Returns
this module

Definition at line 443 of file JModule.hh.

444  {
445  return this->sub(pos);
446  }
JModule & sub(const JVector3D &pos)
Subtract position.
Definition: JModule.hh:365
int JLANG::JObjectID::getID ( ) const
inlineinherited

Get identifier.

Returns
identifier

Definition at line 55 of file JObjectID.hh.

56  {
57  return __id;
58  }
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  }
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  }
int getID() const
Get identifier.
Definition: JObjectID.hh:55
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  }
int getID() const
Get identifier.
Definition: JObjectID.hh:55
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 157 of file JLocation.hh.

158  {
159  if (this->getString() == location.getString())
160  return this->getFloor() < location.getFloor();
161  else
162  return this->getString() < location.getString();
163  }
int getFloor() const
Get floor number.
Definition: JLocation.hh:145
int getString() const
Get string number.
Definition: JLocation.hh:134
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  }
int getID() const
Get identifier.
Definition: JObjectID.hh:55
const JLocation& JDETECTOR::JLocation::getLocation ( ) const
inlineinherited

Get location.

Returns
location

Definition at line 69 of file JLocation.hh.

70  {
71  return static_cast<const JLocation&>(*this);
72  }
Logical location of module.
Definition: JLocation.hh:37
JLocation& JDETECTOR::JLocation::getLocation ( )
inlineinherited

Get location.

Returns
location

Definition at line 80 of file JLocation.hh.

81  {
82  return static_cast<JLocation&>(*this);
83  }
Logical location of module.
Definition: JLocation.hh:37
void JDETECTOR::JLocation::setLocation ( const JLocation location)
inlineinherited

Set location.

Parameters
locationlocation

Definition at line 91 of file JLocation.hh.

92  {
93  static_cast<JLocation&>(*this) = location;
94  }
Logical location of module.
Definition: JLocation.hh:37
std::string JDETECTOR::JLocation::toString ( ) const
inlineinherited

Convert module location to string.

Returns
string

Definition at line 102 of file JLocation.hh.

103  {
104  return toString("% %");
105  }
std::string toString() const
Convert module location to string.
Definition: JLocation.hh:102
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 118 of file JLocation.hh.

119  {
120  JLANG::JString buffer(fmt);
121 
122  buffer.replace(target, string, 1);
123  buffer.replace(target, floor, 1);
124 
125  return buffer;
126  }
Wrapper class around STL string class.
Definition: JString.hh:27
int JDETECTOR::JLocation::getString ( ) const
inlineinherited

Get string number.

Returns
string number

Definition at line 134 of file JLocation.hh.

135  {
136  return string;
137  }
int JDETECTOR::JLocation::getFloor ( ) const
inlineinherited

Get floor number.

Returns
floor number

Definition at line 145 of file JLocation.hh.

146  {
147  return floor;
148  }
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  }
Data structure for position in three dimensions.
Definition: JPosition3D.hh:35
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  }
Data structure for position in three dimensions.
Definition: JPosition3D.hh:35
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  }
Data structure for vector in three dimensions.
Definition: JVector3D.hh:33
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  }
Data structure for angles in three dimensions.
Definition: JAngle3D.hh:31
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getX() const
Get x position.
Definition: JVector3D.hh:93
double getZ() const
Get z position.
Definition: JVector3D.hh:114
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getX() const
Get x position.
Definition: JVector3D.hh:93
Data structure for normalised vector in three dimensions.
Definition: JVersor3D.hh:23
double getZ() const
Get z position.
Definition: JVector3D.hh:114
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  }
JPosition3D & rotate_back(const JRotation3D &R)
Rotate back.
Definition: JPosition3D.hh:199
JVector3D & add(const JVector3D &vector)
Add vector.
Definition: JVector3D.hh:141
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getX() const
Get x position.
Definition: JVector3D.hh:93
double getZ() const
Get z position.
Definition: JVector3D.hh:114
double getDX() const
Get x direction.
Definition: JAngle3D.hh:106
double getDZ() const
Get z direction.
Definition: JAngle3D.hh:128
double getDY() const
Get y direction.
Definition: JAngle3D.hh:117
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  }
double getDY() const
Get y direction.
Definition: JVersor3D.hh:103
double getDX() const
Get x direction.
Definition: JVersor3D.hh:92
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getX() const
Get x position.
Definition: JVector3D.hh:93
double getZ() const
Get z position.
Definition: JVector3D.hh:114
double getDZ() const
Get z direction.
Definition: JVersor3D.hh:114
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  }
double getDY() const
Get y direction.
Definition: JVersor3Z.hh:156
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getDX() const
Get x direction.
Definition: JVersor3Z.hh:145
double getX() const
Get x position.
Definition: JVector3D.hh:93
double getDZ() const
Get z direction.
Definition: JVersor3Z.hh:167
double getZ() const
Get z position.
Definition: JVector3D.hh:114
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getX() const
Get x position.
Definition: JVector3D.hh:93
double getZ() const
Get z position.
Definition: JVector3D.hh:114
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getX() const
Get x position.
Definition: JVector3D.hh:93
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  }
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  }
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  }
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  }
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  }
JVector3D & JMATH::JMath< JVector3D , JNullType >::mul ( const JNullType 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  }
Auxiliary class for arithmetic operations on objects.
Definition: JCalculator.hh:18
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  }
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getX() const
Get x position.
Definition: JVector3D.hh:93
double getZ() const
Get z position.
Definition: JVector3D.hh:114
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getX() const
Get x position.
Definition: JVector3D.hh:93
double getZ() const
Get z position.
Definition: JVector3D.hh:114
double JGEOMETRY3D::JVector3D::getLength ( ) const
inlineinherited

Get length.

Returns
length

Definition at line 245 of file JVector3D.hh.

246  {
247  return sqrt(getLengthSquared());
248  }
double getLengthSquared() const
Get length squared.
Definition: JVector3D.hh:234
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  }
JVector3D()
Default constructor.
Definition: JVector3D.hh:40
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  }
double getDistanceSquared(const JVector3D &pos) const
Get squared of distance to point.
Definition: JVector3D.hh:257
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  }
double getY() const
Get y position.
Definition: JVector3D.hh:103
double getX() const
Get x position.
Definition: JVector3D.hh:93
double getZ() const
Get z position.
Definition: JVector3D.hh:114
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  }
std::vector< std::vector< int > > zbuf2D
std::pair< int, int > pair_type
Data structure for a pair of indices.
std::vector< pair_type > zbuf1D
int j
Definition: JPolint.hh:634
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  }
std::vector< std::vector< int > > zbuf2D
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  }
std::vector< pair_type > zbuf1D
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  }
then echo The file $DIR KM3NeT_00000001_00000000 root already please rename or remove it first
std::vector< std::vector< int > > zbuf2D
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  }
std::vector< pair_type > zbuf1D
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  }
std::vector< std::vector< int > > zbuf2D
std::pair< int, int > pair_type
Data structure for a pair of indices.
std::vector< pair_type > zbuf1D
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  }
then echo The file $DIR KM3NeT_00000001_00000000 root already please rename or remove it first
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  }
static int getSign(const int first, const int second)
Sign of pair of indices.

Member Data Documentation

Double_t JCALIBRATE::JFitK40::sigmaK40_ns
private

intrinsic K40 arrival time spread [ns]

Definition at line 650 of file JFitK40.hh.

int JCALIBRATE::JFitK40::index_of_average_t0
private

index of t0 used for average time offset

Definition at line 651 of file JFitK40.hh.

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

disable PMT from fit

Definition at line 652 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::Rate_Hz
inherited

maximal coincidence rate [Hz]

Definition at line 269 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::p1
inherited

angle dependence coincidence rate

Definition at line 270 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::p2
inherited

angle dependence coincidence rate

Definition at line 271 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::p3
inherited

angle dependence coincidence rate

Definition at line 272 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::p4
inherited

angle dependence coincidence rate

Definition at line 273 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::bg
inherited

remaining constant background

Definition at line 274 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::cc
inherited

fraction of signal correlated background

Definition at line 275 of file JFitK40.hh.

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

Definition at line 277 of file JFitK40.hh.

int JLANG::JObjectID::__id
protectedinherited

Definition at line 160 of file JObjectID.hh.

int JDETECTOR::JLocation::string
protectedinherited

Definition at line 234 of file JLocation.hh.

int JDETECTOR::JLocation::floor
protectedinherited

Definition at line 235 of file JLocation.hh.

double JGEOMETRY3D::JVector3D::__x
protectedinherited

Definition at line 309 of file JVector3D.hh.

double JGEOMETRY3D::JVector3D::__y
protectedinherited

Definition at line 310 of file JVector3D.hh.

double JGEOMETRY3D::JVector3D::__z
protectedinherited

Definition at line 311 of file JVector3D.hh.

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

Definition at line 176 of file JCombinatorics.hh.

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: