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

#include <JFitK40.hh>

Inheritance diagram for JCALIBRATE::JFitK40:

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 Double_t ct) const
	Get K40 coincidence rate as a function of cosine angle between PMT axes. 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...

const JPMT &	getPMT (const int index) const
	Get PMT. More...

JPMT &	getPMT (const int index)
	Get PMT. More...

void	compile ()
	Compile position of module from the positions of the PMTs. More...

void	rotate (const JRotation3D &R)
	Rotate module. More...

JPosition3D &	rotate (const JRotation3X &R)
	Rotate around X-axis. More...

JPosition3D &	rotate (const JRotation3Y &R)
	Rotate around Y-axis. More...

JPosition3D &	rotate (const JRotation3Z &R)
	Rotate around Z-axis. More...

JPosition3D &	rotate (const JQuaternion3D &Q)
	Rotate. More...

void	rotate_back (const JRotation3D &R)
	Rotate back module. More...

JPosition3D &	rotate_back (const JRotation3X &R)
	Rotate back around X-axis. More...

JPosition3D &	rotate_back (const JRotation3Y &R)
	Rotate back around Y-axis. More...

JPosition3D &	rotate_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...

JVector3D &	transform (const JMatrix3D &T)
	Transform. More...

JModule &	set (const JVector3D &pos)
	Set position. More...

JModule &	set (const double t0)
	Set time offset. More...

JModule &	add (const JVector3D &pos)
	Add position. More...

JModule &	add (const double t0)
	Add time offset. More...

JModule &	sub (const JVector3D &pos)
	Subtract position. More...

JModule &	sub (const double t0)
	Subtract time offset. More...

JModule &	operator+= (const JVector3D &pos)
	Add position. More...

JModule &	operator-= (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 JModuleLocation &	getLocation () const
	Get location. More...

JModuleLocation &	getLocation ()
	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 JPosition3D &	getPosition () const
	Get position. More...

JPosition3D &	getPosition ()
	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...

JVector3D &	negate ()
	Negate vector. More...

JVector3D &	mul (const double factor)
	Scale vector. More...

JVector3D &	mul (const JNullType &object)
	Multiply with object. More...

JVector3D &	div (const double factor)
	Scale vector. More...

bool	equals (const JVector3D &vector) 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...

JVector3D &	cross (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_type &	getPair (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 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_type >	zbuf1D

std::vector< std::vector< int > >	zbuf2D

Private Member Functions
	JFitK40 ()
	Default constructor. More...

Static Private Member Functions
static JFitK40 &	getInstance ()
	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 273 of file JFitK40.hh.

Member Typedef Documentation

typedef JCombinatorics::pair_type JCALIBRATE::JFitK40::pair_type

Data structure for a pair of addresses.

Definition at line 282 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

module	detector module
xmin	minimal x
xmax	maximal x
ymin	minimal y
ymax	maximal y
option	true = fix average t0; false = free average t0

Definition at line 295 of file JFitK40.hh.

                                    :
 
       TF2("fitk40",
           JFitK40::getRate,
           xmin, xmax, ymin, ymax,
           JFitK40::getNumberOfModelParameters()),
 
       sigmaK40_ns(0.54)
     {
       static_cast<JModule&>(*this) = module;
     
       this->configure(module.size());
       
       this->sort(JPairwiseComparator(module));
 
       index_of_average_t0 = (option ? 0 : -1);
 
       for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
         disable[i] = false;
       }
     }

JCALIBRATE::JFitK40::JFitK40 ( )

inlineprivate

Default constructor.

Definition at line 654 of file JFitK40.hh.

655 {}

Member Function Documentation

Double_t JCALIBRATE::JFitK40::getSigmaK40 ( ) const

inline

Get intrinsic K40 arrival time spread.

Returns: sigma [ns]

Definition at line 328 of file JFitK40.hh.

     {
       return this->sigmaK40_ns;
     }

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

inline

Set intrinsic K40 arrival time spread.

Parameters

sigma sigma [ns]

Definition at line 339 of file JFitK40.hh.

     {
       this->sigmaK40_ns = sigma;
     }

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

inline

Test PMT status.

Parameters

pmt	pmt address

Returns: true if given pmt is disabled; else false

Definition at line 351 of file JFitK40.hh.

     {
       return disable[pmt];
     }

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

inline

Test PMT status.

Parameters

pmt	pmt address

Returns: true if given pmt is used for average t0; else false

Definition at line 363 of file JFitK40.hh.

     {
       return pmt == index_of_average_t0;
     }

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

pmt	pmt address

Returns: time offset [ns]

Definition at line 378 of file JFitK40.hh.

     {
       if (index_of_average_t0 == -1) {
 
         return JFitK40Parameters::getT0(pmt);
 
       } else {
 
         if (disable[pmt]) {
 
           return parameters[pmt].t0;
 
         } else {        
 
           Int_t    n0  =  0;       // number  of non-fixed t0's
           Double_t t0  =  0.0;     // average of non-fixed t0's
           
           for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
             if (!disable[i] && i != index_of_average_t0) {
               n0  +=  1;
               t0  +=  parameters[i].t0;
             }
           }
           
           t0  /=  (n0 + 1);
 
           if (pmt != index_of_average_t0)
             return parameters[pmt].t0 - t0;
           else
             return -t0;
         }
       }
     }

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

inline

Enable PMT.

Parameters

pmt	pmt address

Definition at line 418 of file JFitK40.hh.

     {
       disable[pmt] = false;
     }

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

pmt	pmt address

Definition at line 432 of file JFitK40.hh.

     {
       disable[pmt] = true;
 
       if (index_of_average_t0 == pmt) {
 
         for (index_of_average_t0 = 0; index_of_average_t0 != NUMBER_OF_PMTS; ++index_of_average_t0) {
           if (!disable[index_of_average_t0]) {
             break;
           }
         }
           
         if (index_of_average_t0 == NUMBER_OF_PMTS) {
           THROW(JIndexOutOfRange, "JFitK40::disable PMT: No free index for average t0.");
         }
       }
     }

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

inline

Get cosine of space angle between PMT axes.

Parameters

pair	pair of addresses

Returns: cosine

Definition at line 457 of file JFitK40.hh.

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

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

inline

Get time resolution of given PMT pair.

Parameters

pair	pair of addresses

Returns: sigma [ns]

Definition at line 470 of file JFitK40.hh.

     {
       return TMath::Sqrt(getTTS(pair.first)  * getTTS(pair.first)  +
                          getTTS(pair.second) * getTTS(pair.second) +
                          getSigmaK40()       * getSigmaK40());
     }

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

inline

Get time offset of given PMT pair.

Parameters

pair	pair of addresses

Returns: time offset [ns]

Definition at line 484 of file JFitK40.hh.

     {
       if (index_of_average_t0 == -1) {
 
         return     parameters[pair.first].t0 - parameters[pair.second].t0;
 
       } else {
 
         if        (pair.first  == index_of_average_t0) {
           
           return  -parameters[pair.second].t0;
           
         } else if (pair.second == index_of_average_t0) {
           
           return   parameters[pair.first].t0;
         
         } else {
           
           return   parameters[pair.first].t0 - parameters[pair.second].t0;
         }
       }
     }

Double_t JCALIBRATE::JFitK40::getValue ( const Double_t ct ) const

inline

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

Parameters

ct	cosine angle between PMT axes

Returns: rate [Hz]

Definition at line 514 of file JFitK40.hh.

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

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

inline

Get K40 coincidence rate.

Parameters

pair	pair of addresses

Returns: rate [Hz]

Definition at line 526 of file JFitK40.hh.

     {
       const Double_t ct = getDot(pair);
 
       return (getValue(ct)         *
               getQE(pair.first)    *
               getQE(pair.second));
     }

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

inline

Get K40 coincidence rate.

Parameters

pair	pair of addresses
dt_ns	time difference [ns]

Returns: rate [Hz/ns]

Definition at line 543 of file JFitK40.hh.

     {
       const Double_t t0    = getT0   (pair);
       const Double_t sigma = getSigma(pair);
 
       return getValue(pair) * TMath::Gaus(dt_ns, t0, sigma, kTRUE);
     }

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

h2	ROOT 2D-histogram
option	fit option

Definition at line 561 of file JFitK40.hh.

     {
       for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {
 
         parameters[pmt] = JPMTParameters_t();
 
         if (disable[pmt]) {
 
           parameters[pmt].QE = 0.0;
 
           fixParameter(*this, this->getModelParameter(pmt, &JPMTParameters_t::QE));
           fixParameter(*this, this->getModelParameter(pmt, &JPMTParameters_t::TTS));
           fixParameter(*this, this->getModelParameter(pmt, &JPMTParameters_t::t0));
 
         } else {
 
           // Note that setting limits to a parameter will also release it
 
           if (!isParameterFixed(*this, this->getModelParameter(pmt, &JPMTParameters_t::QE)))  {
             setParLimits(*this, this->getModelParameter(pmt, &JPMTParameters_t::QE),  FITK40_QE_MIN,            FITK40_QE_MAX);
           }
           if (!isParameterFixed(*this, this->getModelParameter(pmt, &JPMTParameters_t::TTS))) {
             setParLimits(*this, this->getModelParameter(pmt, &JPMTParameters_t::TTS), FITK40_TTS_MIN_NS,        FITK40_TTS_MAX_NS);
           }
           if (!isParameterFixed(*this, this->getModelParameter(pmt, &JPMTParameters_t::t0)))  {
             setParLimits(*this, this->getModelParameter(pmt, &JPMTParameters_t::t0),  h2.GetYaxis()->GetXmin(), h2.GetYaxis()->GetXmax());
           }
         }
       }
 
       if (index_of_average_t0 != -1) {
 
         this->setT0(index_of_average_t0, 0.0);
 
         fixParameter(*this, this->getModelParameter(index_of_average_t0, &JPMTParameters_t::t0));
       }
 
       this->SetParameters(this->getModelParameters());
 
       getInstance() = *this;
 
       const TFitResultPtr result = h2.Fit(this, option.c_str());
 
       this->setModelParameters(this->GetParameters());
 
       return result;
     }

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

x	pointer to data
data	pointer to parameter values

Returns: rate [Hz/ns]

Definition at line 620 of file JFitK40.hh.

     {
       static int       id    =  -1;
       static Double_t  t0;
       static Double_t  sigma;
       static Double_t  rate;
       
       const int      ix    = (int) x[0];
       const Double_t dt_ns = x[1];
 
       if (ix != id) {
         
         getInstance().setModelParameters(data);
         
         const pair_type& pair = getInstance().getPair(ix);
 
         t0     =  getInstance().getT0   (pair);
         sigma  =  getInstance().getSigma(pair);
         rate   =  getInstance().getValue(pair);
         id     =  ix;
       }
 
       return getInstance().bg  +  rate * (getInstance().cc + TMath::Gaus(dt_ns, t0, sigma, kTRUE));
     }

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

inlinestaticprivate

Get unique instance of fit object.

Returns: reference to fit object.

Definition at line 663 of file JFitK40.hh.

     {
       static JFitK40 object;
 
       return object;
     }

static Int_t JCALIBRATE::JFitK40Parameters::getNumberOfModelParameters ( )

inlinestaticinherited

Get number of model parameters.

Returns: number of parameters

Definition at line 98 of file JFitK40.hh.

     {
       return sizeof(JFitK40Parameters) / sizeof(Double_t);
     }

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

inlineinherited

Get model parameters.

Returns: pointer to parameters

Definition at line 109 of file JFitK40.hh.

     {
       return &Rate_Hz;
     }

Double_t* JCALIBRATE::JFitK40Parameters::getModelParameters ( )

inlineinherited

Get model parameters.

Returns: pointer to parameters

Definition at line 120 of file JFitK40.hh.

     {
       return &Rate_Hz;
     }

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

inlineinherited

Set model parameters.

Parameters

data	pointer to parameters

Definition at line 131 of file JFitK40.hh.

     {
       for (Int_t i = 0; i != getNumberOfModelParameters(); ++i) {
         getModelParameters()[i] = data[i];
       }
     }

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

inlineinherited

Get model parameter.

Parameters

i	parameter index

Returns: parameter value

Definition at line 145 of file JFitK40.hh.

     {
       return getModelParameters()[i];
     }

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

inlineinherited

Get model parameter.

Parameters

p	pointer to data member

Returns: parameter index and value

Definition at line 157 of file JFitK40.hh.

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

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

inlineinherited

Get model parameter.

Parameters

pmt	PMT number
p	pointer to data member of PMT parameters

Returns: parameter index and value

Definition at line 172 of file JFitK40.hh.

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

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

inlineinherited

Get QE of given PMT.

Parameters

pmt	pmt address

Returns: QE

Definition at line 186 of file JFitK40.hh.

     {
       return parameters[pmt].QE;
     }

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

inlineinherited

Set QE of given PMT.

Parameters

pmt	pmt address
QE	QE

Definition at line 198 of file JFitK40.hh.

     {
       parameters[pmt].QE = QE;
     }

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

inlineinherited

Get time resolution of given PMT.

Parameters

pmt	pmt address

Returns: TTS [ns]

Definition at line 210 of file JFitK40.hh.

     {
       return parameters[pmt].TTS;
     }

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

inlineinherited

Set time resolution of given PMT.

Parameters

pmt	pmt address
TTS	TTS [ns]

Definition at line 222 of file JFitK40.hh.

     {
       parameters[pmt].TTS = TTS;
     }

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

inlineinherited

Set time offset of given PMT.

Parameters

pmt	pmt address
t0	time offset [ns]

Definition at line 248 of file JFitK40.hh.

     {
       parameters[pmt].t0 = t0;
     }

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

inlineinherited

Get PMT.

Parameters

index readout channel (TDC)

Returns: PMT at given index

Definition at line 141 of file JModule.hh.

     {
       return at(index);
     }

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

inlineinherited

Get PMT.

Parameters

index readout channel (TDC)

Returns: PMT at given index

Definition at line 153 of file JModule.hh.

     {
       return at(index);
     }

void JDETECTOR::JModule::compile ( )

inlineinherited

Compile position of module from the positions of the PMTs.

Definition at line 162 of file JModule.hh.

     {
       if (!empty()) {
 
         JPosition3D pos;
 
         for (iterator i = begin(); i != end(); ++i) {
           pos.add(i->getPosition());
         }
 
         static_cast<JPosition3D&>(*this) = pos.div(size());
       }
     }

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

inlineinherited

Rotate module.

Parameters

R	rotation matrix

Definition at line 182 of file JModule.hh.

     {
       static_cast<JPosition3D&>(*this).rotate(R);
 
       for (iterator i = this->begin(); i != this->end(); ++i) {
         i->rotate(R);
       }
     }

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

inlineinherited

Rotate around X-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 213 of file JPosition3D.hh.

     {
       R.rotate(__y, __z);
 
       return *this;
     }

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

inlineinherited

Rotate around Y-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 241 of file JPosition3D.hh.

     {
       R.rotate(__x, __z);
 
       return *this;
     }

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

inlineinherited

Rotate around Z-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 269 of file JPosition3D.hh.

     {
       R.rotate(__x, __y);
 
       return *this;
     }

JPosition3D& JGEOMETRY3D::JPosition3D::rotate ( const JQuaternion3D & Q )

inlineinherited

Rotate.

Parameters

Q	quaternion

Returns: this position

Definition at line 297 of file JPosition3D.hh.

     {
       JQuaternion3D b(Q);
       JQuaternion3D a;
 
       a.mul(getPosition(), b.conjugate());
 
       setPosition(b.mul(Q,a));
       
       return *this;
     }

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

inlineinherited

Rotate back module.

Parameters

R	rotation matrix

Definition at line 197 of file JModule.hh.

     {
       static_cast<JPosition3D&>(*this).rotate_back(R);
 
       for (iterator i = this->begin(); i != this->end(); ++i) {
         i->rotate_back(R);
       }
     }

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

inlineinherited

Rotate back around X-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 227 of file JPosition3D.hh.

     {
       R.rotate_back(__y, __z);
 
       return *this;
     }

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

inlineinherited

Rotate back around Y-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 255 of file JPosition3D.hh.

     {
       R.rotate_back(__x, __z);
 
       return *this;
     }

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

inlineinherited

Rotate back around Z-axis.

Parameters

R	rotation matrix

Returns: this position

Definition at line 283 of file JPosition3D.hh.

     {
       R.rotate_back(__x, __y);
 
       return *this;
     }

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

inlineinherited

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

Parameters

R	rotation matrix
pos	position of origin (after rotation)

Definition at line 213 of file JModule.hh.

     {
       static_cast<JPosition3D&>(*this).transform(R, pos);
 
       for (iterator i = this->begin(); i != this->end(); ++i) {
         i->transform(R, pos);
       }
     }

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

inlineinherited

Transformation of geometry.

Parameters

T	transformation

Definition at line 229 of file JModule.hh.

     {
       static_cast<JPosition3D&>(*this).transform(T.getRotation(), T.getPosition());
 
       for (iterator i = this->begin(); i != this->end(); ++i) {
         i->transform(T);
       }
     }

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

inlineinherited

Transform.

Parameters

T matrix

Returns: this vector

Definition at line 204 of file JVector3D.hh.

     {
       T.transform(__x, __y, __z);
 
       return *this;
     }

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

inlineinherited

Set position.

Parameters

pos position

Returns: this module

Definition at line 245 of file JModule.hh.

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

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

inlineinherited

Set time offset.

Parameters

t0	time offset [ns]

Returns: this module

Definition at line 293 of file JModule.hh.

     {
       for (iterator i = begin(); i != end(); ++i) {
         i->setT0(t0);
       }
 
       return *this;
     }

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

inlineinherited

Add position.

Parameters

pos position

Returns: this module

Definition at line 257 of file JModule.hh.

     {
       for (iterator i = begin(); i != end(); ++i) {
         i->add(pos);
       }
 
       JVector3D::add(pos);
 
       return *this;
     }

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

inlineinherited

Add time offset.

Parameters

t0	time offset [ns]

Returns: this module

Definition at line 309 of file JModule.hh.

     {
       for (iterator i = begin(); i != end(); ++i) {
         i->addT0(t0);
       }
 
       return *this;
     }

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

inlineinherited

Subtract position.

Parameters

pos position

Returns: this module

Definition at line 275 of file JModule.hh.

     {
       for (iterator i = begin(); i != end(); ++i) {
         i->sub(pos);
       }
 
       JVector3D::sub(pos);
 
       return *this;
     }

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

inlineinherited

Subtract time offset.

Parameters

t0	time offset [ns]

Returns: this module

Definition at line 325 of file JModule.hh.

     {
       for (iterator i = begin(); i != end(); ++i) {
         i->subT0(t0);
       }
 
       return *this;
     }

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

inlineinherited

Add position.

Parameters

pos position

Returns: this module

Definition at line 341 of file JModule.hh.

     {
       return this->add(pos);
     }

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

inlineinherited

Subtract position.

Parameters

pos position

Returns: this module

Definition at line 353 of file JModule.hh.

     {
       return this->sub(pos);
     }

int JLANG::JObjectID::getID ( ) const

inlineinherited

Get identifier.

Returns: identifier

Definition at line 54 of file JObjectID.hh.

     { 
       return __id;
     }

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

inlineinherited

Set identifier.

Parameters

id	identifier

Definition at line 65 of file JObjectID.hh.

     { 
       this->__id = id;
     }

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

inlineinherited

Less than method.

Parameters

object object identifier

Returns: true if this identifier less than given identifier; else false

Definition at line 77 of file JObjectID.hh.

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

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

inlineinherited

Less than method.

Parameters

id	identifier

Returns: true if this identifier less than given identifier; else false

Definition at line 89 of file JObjectID.hh.

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

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

inlineinherited

Less than method.

Parameters

location module location

Returns: true if first location before second location; else false

Definition at line 156 of file JModuleLocation.hh.

     {
       if (this->getString() == location.getString())
         return this->getFloor()  < location.getFloor();
       else
         return this->getString() < location.getString();
     }

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

inlineinherited

More than method.

Parameters

id	identifier

Returns: true if this identifier greater than given identifier; else false

Definition at line 101 of file JObjectID.hh.

     { 
       return this->getID() > id;
     }

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

inlineinherited

Get location.

Returns: location

Definition at line 68 of file JModuleLocation.hh.

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

JModuleLocation& JDETECTOR::JModuleLocation::getLocation ( )

inlineinherited

Get location.

Returns: location

Definition at line 79 of file JModuleLocation.hh.

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

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

inlineinherited

Set location.

Parameters

location location

Definition at line 90 of file JModuleLocation.hh.

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

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

inlineinherited

Convert module location to string.

Returns: string

Definition at line 101 of file JModuleLocation.hh.

     {
       return toString("% %");
     }

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

fmt	format
target	target

Returns: string

Definition at line 117 of file JModuleLocation.hh.

     {
       JLANG::JString buffer(fmt);
 
       buffer.replace(target, string, 1);
       buffer.replace(target, floor,  1);
 
       return buffer;
     }

int JDETECTOR::JModuleLocation::getString ( ) const

inlineinherited

Get string number.

Returns: string number

Definition at line 133 of file JModuleLocation.hh.

     {
       return string;
     }

int JDETECTOR::JModuleLocation::getFloor ( ) const

inlineinherited

Get floor number.

Returns: floor number

Definition at line 144 of file JModuleLocation.hh.

     {
       return floor;
     }

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

inlineinherited

Get position.

Returns: position

Definition at line 129 of file JPosition3D.hh.

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

JPosition3D& JGEOMETRY3D::JPosition3D::getPosition ( )

inlineinherited

Get position.

Returns: position

Definition at line 140 of file JPosition3D.hh.

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

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

inlineinherited

Set position.

Parameters

pos position

Definition at line 151 of file JPosition3D.hh.

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

JGEOMETRY3D::JPosition3D::operator JAngle3D ( ) const

inlineinherited

Type conversion operator.

Returns: angle

Definition at line 162 of file JPosition3D.hh.

     {
       return JAngle3D(getX(), getY(), getZ());
     }

JGEOMETRY3D::JPosition3D::operator JVersor3D ( ) const

inlineinherited

Type conversion operator.

Returns: direction

Definition at line 173 of file JPosition3D.hh.

     {
       return JVersor3D(getX(), getY(), getZ());
     }

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

inlineinherited

Transform back position.

The final position is obtained as follows:

offset position with position pos;
rotation of postion according matrix R;

Parameters

R	rotation matrix
pos	position of origin (before rotation)

Definition at line 349 of file JPosition3D.hh.

     {
       // offset with respect to origin
 
       add(pos);
 
       // rotate back geometry to system with particle direction along z-axis
 
       rotate_back(R);
     }

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

inlineinherited

Get dot product.

Parameters

angle angle

Returns: dot product

Definition at line 368 of file JPosition3D.hh.

     {
       return
         getX() * angle.getDX() +
         getY() * angle.getDY() +
         getZ() * angle.getDZ();
     }

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

inlineinherited

Get dot product.

Parameters

dir direction

Returns: dot product

Definition at line 383 of file JPosition3D.hh.

     {
       return
         getX() * dir.getDX() +
         getY() * dir.getDY() +
         getZ() * dir.getDZ();
     }

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

inlineinherited

Get dot product.

Parameters

dir direction

Returns: dot product

Definition at line 398 of file JPosition3D.hh.

     {
       return
         getX() * dir.getDX() +
         getY() * dir.getDY() +
         getZ() * dir.getDZ();
     }

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

inlineinherited

Get dot product.

Parameters

vector vector

Returns: dot product

Definition at line 278 of file JVector3D.hh.

     {
       return
         getX() * vector.getX() +
         getY() * vector.getY() +
         getZ() * vector.getZ();
     }

JGEOMETRY3D::JVector3D::operator JVector2D ( ) const

inlineinherited

Type conversion operator.

Returns: JVector2D

Definition at line 81 of file JVector3D.hh.

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

double JGEOMETRY3D::JVector3D::getX ( ) const

inlineinherited

Get x position.

Returns: x position

Definition at line 92 of file JVector3D.hh.

     { 
       return __x;
     }

double JGEOMETRY3D::JVector3D::getY ( ) const

inlineinherited

Get y position.

Returns: y position

Definition at line 102 of file JVector3D.hh.

     { 
       return __y;
     }

double JGEOMETRY3D::JVector3D::getZ ( ) const

inlineinherited

Get z position.

Returns: z position

Definition at line 113 of file JVector3D.hh.

     { 
       return __z;
     }

JVector3D& JGEOMETRY3D::JVector3D::negate ( )

inlineinherited

Negate vector.

Returns: this vector

Definition at line 124 of file JVector3D.hh.

     {
       __x = -__x;
       __y = -__y;
       __z = -__z;
 
       return *this;
     }

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

inlineinherited

Scale vector.

Parameters

factor multiplication factor

Returns: this vector

Definition at line 172 of file JVector3D.hh.

     {
       __x *= factor;
       __y *= factor;
       __z *= factor;
 
       return *this;
     }

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

inlineinherited

Multiply with object.

Parameters

object object

Returns: result object

Definition at line 273 of file JMath.hh.

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

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

inlineinherited

Scale vector.

Parameters

factor division factor

Returns: this vector

Definition at line 188 of file JVector3D.hh.

     {
       __x /= factor;
       __y /= factor;
       __z /= factor;
 
       return *this;
     }

bool JGEOMETRY3D::JVector3D::equals ( const JVector3D & vector ) const

inlineinherited

Check equality.

Parameters

vector vector

Returns: true if vectors are equal; else false

Definition at line 218 of file JVector3D.hh.

     {
       return (getX() == vector.getX() && 
               getY() == vector.getY() &&
               getZ() == vector.getZ());
     }

double JGEOMETRY3D::JVector3D::getLengthSquared ( ) const

inlineinherited

Get length squared.

Returns: square of length

Definition at line 231 of file JVector3D.hh.

     { 
       return getX()*getX() + getY()*getY() + getZ()*getZ();
     }

double JGEOMETRY3D::JVector3D::getLength ( ) const

inlineinherited

Get length.

Returns: length

Definition at line 242 of file JVector3D.hh.

     { 
       return sqrt(getLengthSquared());
     }

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

inlineinherited

Get squared of distance to point.

Parameters

pos position

Returns: square of distance

Definition at line 254 of file JVector3D.hh.

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

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

inlineinherited

Get distance to point.

Parameters

pos position

Returns: distance

Definition at line 266 of file JVector3D.hh.

     {
       return sqrt(getDistanceSquared(pos));
     }

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

first	first vector
second	second vector

Returns: this vector

Definition at line 295 of file JVector3D.hh.

     {
       __x = first .getY() * second.getZ()  -  second.getY() * first .getZ();
       __y = second.getX() * first .getZ()  -  first .getX() * second.getZ();
       __z = first .getX() * second.getY()  -  second.getX() * first .getY();
 
       return *this;
     }

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

inlineinherited

Configure.

Parameters

numberOfIndices number of indices

Definition at line 54 of file JCombinatorics.hh.

     {
       zbuf1D.clear();
 
       zbuf2D.resize(numberOfIndices);
 
       for (int i = 0; i != numberOfIndices; ++i) {
         zbuf2D[i].resize(numberOfIndices);
       }
 
       for (int i = 0; i != numberOfIndices; ++i) {
 
         zbuf2D[i][i] = -1;
 
         for (int j = i; ++j != numberOfIndices; ) {
 
           zbuf2D[i][j] = zbuf1D.size();
           zbuf2D[j][i] = zbuf1D.size();
 
           zbuf1D.push_back(pair_type(i,j));
         }
       }
     }

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

inlineinherited

Get number of indices.

Returns: number of indices

Definition at line 84 of file JCombinatorics.hh.

     {
       return zbuf2D.size();
     }

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

inlineinherited

Get number of pairs.

Returns: number of pairs

Definition at line 95 of file JCombinatorics.hh.

     {
       return zbuf1D.size();
     }

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

inlineinherited

Get index of pair of indices.

Parameters

first	first address
second	second address

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

Definition at line 108 of file JCombinatorics.hh.

     {
       return zbuf2D[first][second];
     }

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

inlineinherited

Get pair of indices for given index.

Parameters

index index

Returns: pair of indices

Definition at line 120 of file JCombinatorics.hh.

     {
       return zbuf1D[index];
     }

template<class JComparator_t >

void JTOOLS::JCombinatorics::sort ( JComparator_t comparator )

inlineinherited

Sort address pairs.

Parameters

comparator comparator for pairs

Definition at line 132 of file JCombinatorics.hh.

     {
       std::stable_sort(zbuf1D.begin(), zbuf1D.end(), comparator);
 
       for (int i = 0; i != (int) zbuf1D.size(); ++i) {
         
         const pair_type pair = zbuf1D[i];
 
         zbuf2D[pair.first][pair.second] = i;
         zbuf2D[pair.second][pair.first] = i;
       }
 
       for (int i = 0; i != (int) zbuf2D.size(); ++i) {
         zbuf2D[i][i] = -1;
       }
     }

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

inlinestaticinherited

Sign of pair of indices.

Parameters

first	first address
second	second address

Returns: +1 if second >= first; else -1

Definition at line 157 of file JCombinatorics.hh.

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

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

inlinestaticinherited

Sign of pair of indices.

Parameters

pair	pair of indices

Returns: +1 if second >= first; else -1

Definition at line 169 of file JCombinatorics.hh.

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

Member Data Documentation

Double_t JCALIBRATE::JFitK40::sigmaK40_ns

private

intrinsic K40 arrival time spread [ns]

Definition at line 647 of file JFitK40.hh.

int JCALIBRATE::JFitK40::index_of_average_t0

private

index of t0 used for average time offset

Definition at line 648 of file JFitK40.hh.

bool JCALIBRATE::JFitK40::disable[NUMBER_OF_PMTS]

private

disable PMT from fit

Definition at line 649 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::Rate_Hz

inherited

maximal coincidence rate [Hz]

Definition at line 256 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::p1

inherited

angle dependence coincidence rate

Definition at line 257 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::p2

inherited

angle dependence coincidence rate

Definition at line 258 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::p3

inherited

angle dependence coincidence rate

Definition at line 259 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::p4

inherited

angle dependence coincidence rate

Definition at line 260 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::bg

inherited

remaining constant background

Definition at line 261 of file JFitK40.hh.

Double_t JCALIBRATE::JFitK40Parameters::cc

inherited

fraction of signal correlated background

Definition at line 262 of file JFitK40.hh.

JPMTParameters_t JCALIBRATE::JFitK40Parameters::parameters[NUMBER_OF_PMTS]

inherited

Definition at line 264 of file JFitK40.hh.

int JLANG::JObjectID::__id

protectedinherited

Definition at line 159 of file JObjectID.hh.

int JDETECTOR::JModuleLocation::string

protectedinherited

Definition at line 233 of file JModuleLocation.hh.

int JDETECTOR::JModuleLocation::floor

protectedinherited

Definition at line 234 of file JModuleLocation.hh.

double JGEOMETRY3D::JVector3D::__x

protectedinherited

Definition at line 306 of file JVector3D.hh.

double JGEOMETRY3D::JVector3D::__y

protectedinherited

Definition at line 307 of file JVector3D.hh.

double JGEOMETRY3D::JVector3D::__z

protectedinherited

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

JFitK40.hh

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Protected Attributes

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation