Monte Carlo detector (i.e. More...

#include <JMonteCarloDetector.hh>

Inheritance diagram for JDETECTOR::JMonteCarloDetector:

Classes
class	LCM_logic
	Auxiliary class for LCM logic parameters. More...

class	LCM_reverse_logic
	Auxiliary class for LCM logic parameters. More...

class	OM
	Auxiliary class for OM. More...

class	OM_cluster
	Auxiliary class for OM cluster parameters. More...

class	OM_cluster_data
	Auxiliary class for OM cluster data. More...

class	OM_position
	Auxiliary class for OM position. More...

class	String
	Auxiliary class for string parameters. More...

Public Types
enum	JVersion_t { V1 = 1, V2 = 2, V3 = 3 }
	Enumeration of version types. More...

typedef JRange< T, JComparator_t >	range_type

typedef JLANG::JClass< T > ::argument_type	argument_type

typedef T	key_type

typedef T	mapped_type

Public Member Functions
	JMonteCarloDetector ()
	Default constructor. More...

	JMonteCarloDetector (const bool useLogic)
	Constructor. More...

void	setUseLogic (const bool useLogic)
	Set usage of logic. More...

void	setVersion (const JVersion &version)
	Set version. More...

void	setVersion (const JDetectorVersion::JVersion_t &version)
	Set version. More...

JDetector &	operator+= (const JVector3D &pos)
	Move detector elements. More...

JDetector &	operator-= (const JVector3D &pos)
	Move detector elements. More...

const JModule &	getModule (const JModuleAddress &address) const
	Get module parameters. More...

JModule &	getModule (const JModuleAddress &address)
	Get module parameters. More...

const JModule &	getModule (const JModuleLocation &location) const
	Get module parameters. More...

JModule &	getModule (const JModuleLocation &location)
	Get module parameters. More...

const JPMT &	getPMT (const JPMTAddress &address) const
	Get PMT parameters. More...

JPMT &	getPMT (const JPMTAddress &address)
	Get PMT parameters. More...

virtual JReader &	read (JReader &in)
	Read from input. More...

virtual JWriter &	write (JWriter &out) const
	Write to output. 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	more (const int id) const
	More than method. More...

const JVersion &	getVersion () const
	Get version. More...

JVersion &	getVersion ()
	Get version. More...

const std::string &	getVariant () const
	Get variant. More...

void	setVariant (const std::string &variant)
	Set variant. More...

bool	equals (const JVersion &version) const
	Equal method. More...

bool	equals (const range_type &range) const
	Equal method. More...

JProperties	getProperties (const JEquationParameters &equation=JDetectorHeader::getEquationParameters())
	Get properties of this class. More...

JProperties	getProperties (const JEquationParameters &equation=JDetectorHeader::getEquationParameters()) const
	Get properties of this class. More...

const range_type &	getRange () const
	Get range. More...

void	setRange (const range_type &range)
	Set range. More...

void	setRange (argument_type x, argument_type y)
	Set lower and uppper limit. More...

template<class iterator_type , class value_type , class result_type >
void	setRange (iterator_type __begin, iterator_type __end, result_type value_type::*member)
	Set lower and uppper limit. More...

template<class iterator_type , class value_type , class result_type >
void	setRange (iterator_type __begin, iterator_type __end, result_type(value_type::*function)() const)
	Set lower and uppper limit. More...

T	getLowerLimit () const
	Get lower limit. More...

T	getUpperLimit () const
	Get upper limit. More...

void	setLowerLimit (argument_type x)
	Set lower limit. More...

void	setUpperLimit (argument_type y)
	Set upper limit. More...

void	fixLowerLimit (argument_type x)
	Fix lower limit. More...

void	fixUpperLimit (argument_type y)
	Fix upper limit. More...

T	getLength () const
	Get length (difference between upper and lower limit). More...

bool	is_valid () const
	Check validity of range. More...

bool	operator() (argument_type x) const
	Test whether value is inside range. More...

T	constrain (argument_type x) const
	Constrain value to range. More...

bool	overlap (const range_type &range) const
	Test overlap with given range. More...

range_type	include (argument_type x)
	Include given value to range. More...

range_type &	join (const range_type &range)
	Join ranges. More...

range_type &	combine (const range_type &range)
	Combine ranges. More...

range_type &	add (argument_type x)
	Add offset. More...

range_type &	add (const range_type &range)
	Add offset. More...

range_type &	sub (argument_type x)
	Subtract offset. More...

range_type &	sub (const range_type &range)
	Subtract offset. More...

range_type &	mul (const double factor)
	Multiply range. More...

JFirst_t &	mul (const JSecond_t &object)
	Multiply with object. More...

range_type &	div (const double factor)
	Divide range. More...

T	getN (const double R) const
	Get expected number of occurances of given rate within this interval. More...

const key_type &	getKey () const

const std::string &	getKey () const
	Get key. More...

const mapped_type &	getValue () const

mapped_type &	getValue ()

const std::string &	getWGS () const
	Get WGS. More...

const std::string &	getUTMZone () const
	Get UTM zone. More...

std::string	toString () const
	Convert UTM grid. More...

JPosition3D	getPosition () const
	Get position. More...

double	getUTMEast () const
	Get UTM east. More...

double	getUTMNorth () const
	Get UTM north. More...

double	getUTMZ () const
	Get UTM Z. More...

Static Public Member Functions
static JEquationParameters &	getEquationParameters ()
	Get equation parameters. More...

static void	setEquationParameters (const JEquationParameters &equation)
	Set equation parameters. More...

static T	getMinimum ()
	Get minimum possible value. More...

static T	getMaximum ()
	Get maximum possible value. More...

static JUTMGrid	valueOf (const std::string &buffer)
	Extract UTM grid. More...

Public Attributes
JComment	comment

JComparator_t	compare
	Function object. More...

T	first

T	second

Static Public Attributes
static const JRange< T, JComparator_t >	DEFAULT_RANGE
	Default range. More...

Protected Member Functions
JDetector &	put (const int moduleID, const int pmtAddress, const JModuleLocation &location, const JPMT &pmt)
	Set PMT. More...

Static Protected Member Functions
template<class T >
static T	find (T __begin, T __end, const int id)
	Binary search method. More...

static bool	compare (const JModule &first, const JModule &second)
	Module comparator. More...

template<class JHead_t , class JTail_t >
static bool	eq (const JClass_t &first, const JClass_t &second, const JType< JTypeList< JHead_t, JTail_t > > &type)
	Equals method for composite data types. More...

template<class JHead_t >
static bool	eq (const JClass_t &first, const JClass_t &second, const JType< JTypeList< JHead_t, JNullType > > &type)
	Equals method for composite data types. More...

Protected Attributes
int	__id

std::string	__variant

std::string	key

std::string	wgs

std::string	zone

double	east

double	north

double	z

Private Attributes
bool	use_logic

Friends
std::istream &	operator>> (std::istream &in, JMonteCarloDetector &detector)
	Read detector from input. More...

Detailed Description

Monte Carlo detector (i.e.

so-called .det file).

Definition at line 43 of file JMonteCarloDetector.hh.

Member Typedef Documentation

template<class T, class JComparator_t = std::less<T>>

typedef JRange<T, JComparator_t> JTOOLS::JRange< T, JComparator_t >::range_type

inherited

Definition at line 40 of file JRange.hh.

template<class T, class JComparator_t = std::less<T>>

typedef JLANG::JClass<T>::argument_type JTOOLS::JRange< T, JComparator_t >::argument_type

inherited

Definition at line 41 of file JRange.hh.

typedef T JTOOLS::JPair< T , T >::key_type

inherited

Definition at line 32 of file JPair.hh.

typedef T JTOOLS::JPair< T , T >::mapped_type

inherited

Definition at line 33 of file JPair.hh.

Member Enumeration Documentation

enum JDETECTOR::JDetectorVersion::JVersion_t

inherited

Enumeration of version types.

Additional version types should be included in JGetDetectorVersion::JGetDetectorVersion.

Enumerator
V1	First version.
V2	Version with UTC time and UTM position data.
V3	Version with PMT status field and comments.

Definition at line 41 of file JDetectorVersion.hh.

                     {
       V1 = 1,        //!< First version
       V2 = 2,        //!< Version with UTC time and UTM position data
       V3 = 3         //!< Version with PMT status field and comments
     };

Constructor & Destructor Documentation

JDETECTOR::JMonteCarloDetector::JMonteCarloDetector ( )

inline

Default constructor.

Definition at line 256 of file JMonteCarloDetector.hh.

                           :
       JDetector(),
       use_logic(false)
     {}

JDETECTOR::JMonteCarloDetector::JMonteCarloDetector ( const bool useLogic )

inline

Constructor.

Parameters

useLogic apply logic data in file

Definition at line 267 of file JMonteCarloDetector.hh.

                                              :
       JDetector(),
       use_logic(useLogic)
     {}

Member Function Documentation

void JDETECTOR::JMonteCarloDetector::setUseLogic ( const bool useLogic )

inline

Set usage of logic.

Parameters

useLogic apply logic data in file

Definition at line 527 of file JMonteCarloDetector.hh.

     {
       use_logic = useLogic;
     }

JDetector& JDETECTOR::JMonteCarloDetector::put	(	const int	moduleID,
		const int	pmtAddress,
		const JModuleLocation &	location,
		const JPMT &	pmt
	)

inlineprotected

Set PMT.

Parameters

moduleID	module identifier
pmtAddress	PMT address
location	module location
pmt	JPMT

Returns: this JDetector

Definition at line 543 of file JMonteCarloDetector.hh.

     {
       JDetector::iterator p = std::lower_bound(this->begin(), this->end(), moduleID);
 
       if (p == this->end() || p->getID() != moduleID) {
         
         JModule module(moduleID, location);
         
         module.resize(pmtAddress + 1);
 
         module[pmtAddress] = pmt;
 
         module.setPosition(pmt.getPosition());
         
         this->insert(p, module);
         
       } else {
         
         // average position
         
         int n = 0;
 
         for (JModule::const_iterator i = p->begin(); i != p->end(); ++i) {
           if (i->getID() != getUndefinedObjectID().getID()) {
             ++n;
           }
         }
 
         if (pmtAddress + 1 > (int) p->size()) {
           p->resize(pmtAddress + 1);
         }
 
         JPosition3D& P = static_cast<JPosition3D&>(*p);
 
         P.mul(n);
 
         if (p->at(pmtAddress).getID() != getUndefinedObjectID().getID())
           P.sub(p->at(pmtAddress).getPosition());
         else
           ++n;
 
         P.add(pmt.getPosition());
         P.div(n);
 
         p->at(pmtAddress) = pmt;
       }
 
       return *this;
     }

template<class T >

static T JDETECTOR::JMonteCarloDetector::find	(	T	__begin,
		T	__end,
		const int	id
	)

inlinestaticprotected

Binary search method.

Parameters

__begin	begin of data
__end	end of data
id	target value

Returns: pointer to corresponding element

Definition at line 606 of file JMonteCarloDetector.hh.

     {
       T i = std::lower_bound(__begin, __end, id);
 
       if (i != __end && *i != id)
         return __end;
       else
         return i;
     }

static bool JDETECTOR::JMonteCarloDetector::compare	(	const JModule &	first,
		const JModule &	second
	)

inlinestaticprotected

Module comparator.

The comparison is based on:

string identifier;
z-position.

Parameters

first	first module
second	second module

Returns: true if first module before second; else false

Definition at line 627 of file JMonteCarloDetector.hh.

     {
       if (first.getString() == second.getString())
         return first.getZ()      < second.getZ();
       else
         return first.getString() < second.getString();
     }

void JDETECTOR::JDetector::setVersion ( const JVersion & version )

inlineinherited

Set version.

Parameters

version version

Definition at line 123 of file JDetector.hh.

     {
       static_cast<JDetectorVersion&>(*this).setVersion(version);
 
       JPMT::setVersion(getVersion());
     }

void JDETECTOR::JDetector::setVersion ( const JDetectorVersion::JVersion_t & version )

inlineinherited

Set version.

Parameters

version version

Definition at line 136 of file JDetector.hh.

     {
       setVersion(putDetectorVersion(version));
     }

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

inlineinherited

Move detector elements.

Parameters

pos	offset position

Returns: this JDetector

Definition at line 148 of file JDetector.hh.

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

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

inlineinherited

Move detector elements.

Parameters

pos	offset position

Returns: this JDetector

Definition at line 164 of file JDetector.hh.

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

const JModule& JDETECTOR::JDetector::getModule ( const JModuleAddress & address ) const

inlineinherited

Get module parameters.

Parameters

address module address

Returns: module parameters

Definition at line 180 of file JDetector.hh.

     {
       return at(address.first);
     }

JModule& JDETECTOR::JDetector::getModule ( const JModuleAddress & address )

inlineinherited

Get module parameters.

Parameters

address module address

Returns: module parameters

Definition at line 192 of file JDetector.hh.

     {
       return at(address.first);
     }

const JModule& JDETECTOR::JDetector::getModule ( const JModuleLocation & location ) const

inlineinherited

Get module parameters.

Parameters

location module location

Returns: module parameters

Definition at line 204 of file JDetector.hh.

     {
       for (JDetector::const_iterator module = this->begin(); module != this->end(); ++module) {
         if (*module == location) {
           return *module;
         }
       }
 
       THROW(JIndexOutOfRange, "JDetector::getModule()");
     }

JModule& JDETECTOR::JDetector::getModule ( const JModuleLocation & location )

inlineinherited

Get module parameters.

Parameters

location module location

Returns: module parameters

Definition at line 222 of file JDetector.hh.

     {
       for (JDetector::iterator module = this->begin(); module != this->end(); ++module) {
         if (*module == location) {
           return *module;
         }
       }
 
       THROW(JIndexOutOfRange, "JDetector::getModule()");
     }

const JPMT& JDETECTOR::JDetector::getPMT ( const JPMTAddress & address ) const

inlineinherited

Get PMT parameters.

Parameters

address JPMTAddress

Returns: JPMT

Definition at line 240 of file JDetector.hh.

     {
       return at(address.first).at(address.second);
     }

JPMT& JDETECTOR::JDetector::getPMT ( const JPMTAddress & address )

inlineinherited

Get PMT parameters.

Parameters

address JPMTAddress

Returns: JPMT

Definition at line 252 of file JDetector.hh.

     {
       return at(address.first).at(address.second);
     }

virtual JReader& JDETECTOR::JDetector::read ( JReader & in )

inlinevirtualinherited

Read from input.

Parameters

in reader

Returns: reader

Implements JIO::JSerialisable.

Definition at line 367 of file JDetector.hh.

     {
       clear();
 
       int n;
 
       in >> static_cast<JObjectID&>       (*this);
       in >> static_cast<JDetectorVersion&>(*this);
       in >> static_cast<JDetectorHeader&> (*this);
       in >> n;
 
       JPMT::setVersion(getVersion());
 
       for (JModule module; n != 0; --n) {
 
         in >> module;
 
         push_back(module);
       }
 
       return in;
     }

virtual JWriter& JDETECTOR::JDetector::write ( JWriter & out ) const

inlinevirtualinherited

Write to output.

Parameters

out writer

Returns: writer

Implements JIO::JSerialisable.

Definition at line 397 of file JDetector.hh.

     {
       int n = size();
 
       out << static_cast<const JObjectID&>       (*this);
       out << static_cast<const JDetectorVersion&>(*this);
       out << static_cast<const JDetectorHeader&> (*this);
       out << n;
 
       for (const_iterator i = begin(); i != end(); ++i) {
         out << *i;
       }
 
       return out;
     }

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 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 JVersion& JEEP::JVersion::getVersion ( ) const

inlineinherited

Get version.

Returns: version

Definition at line 57 of file JVersion.hh.

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

JVersion& JEEP::JVersion::getVersion ( )

inlineinherited

Get version.

Returns: version

Definition at line 68 of file JVersion.hh.

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

const std::string& JEEP::JVersion::getVariant ( ) const

inlineinherited

Get variant.

Returns: variant

Definition at line 90 of file JVersion.hh.

     { 
       return this->__variant;
     }

void JEEP::JVersion::setVariant ( const std::string & variant )

inlineinherited

Set variant.

Parameters

variant variant

Definition at line 101 of file JVersion.hh.

     { 
       this->__variant = variant;
     }

bool JEEP::JVersion::equals ( const JVersion & version ) const

inlineinherited

Equal method.

Parameters

version version

Returns: true if this version equal to given version; else false

Definition at line 113 of file JVersion.hh.

     {
       return this->getVariant() == version.getVariant();
     }

template<class T, class JComparator_t = std::less<T>>

bool JTOOLS::JRange< T, JComparator_t >::equals ( const range_type & range ) const

inlineinherited

Equal method.

Parameters

range range

Returns: true if this module range equal to given module range; else false

Definition at line 253 of file JRange.hh.

     {
       return (!this->compare(this->getLowerLimit(), range.getLowerLimit()) &&
               !this->compare(range.getLowerLimit(), this->getLowerLimit()) &&
               !this->compare(this->getUpperLimit(), range.getUpperLimit()) &&
               !this->compare(range.getUpperLimit(), this->getUpperLimit()));
     }

static JEquationParameters& JDETECTOR::JDetectorHeader::getEquationParameters ( )

inlinestaticinherited

Get equation parameters.

Returns: equation parameters

Definition at line 144 of file JDetectorHeader.hh.

     {
       static JEquationParameters equation;
 
       return equation;
     }

static void JDETECTOR::JDetectorHeader::setEquationParameters ( const JEquationParameters & equation )

inlinestaticinherited

Set equation parameters.

Parameters

equation equation parameters

Definition at line 158 of file JDetectorHeader.hh.

     {
       getEquationParameters() = equation;
     }

JProperties JDETECTOR::JDetectorHeader::getProperties ( const JEquationParameters & equation = JDetectorHeader::getEquationParameters() )

inlineinherited

Get properties of this class.

Parameters

equation equation parameters

Definition at line 169 of file JDetectorHeader.hh.

     {
       return JDetectorHeaderHelper(*this, equation);
     }

JProperties JDETECTOR::JDetectorHeader::getProperties ( const JEquationParameters & equation = JDetectorHeader::getEquationParameters() ) const

inlineinherited

Get properties of this class.

Parameters

equation equation parameters

Definition at line 180 of file JDetectorHeader.hh.

     {
       return JDetectorHeaderHelper(*this, equation);
     }

template<class T, class JComparator_t = std::less<T>>

const range_type& JTOOLS::JRange< T, JComparator_t >::getRange ( ) const

inlineinherited

Get range.

Returns: range

Definition at line 109 of file JRange.hh.

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

template<class T, class JComparator_t = std::less<T>>

void JTOOLS::JRange< T, JComparator_t >::setRange ( const range_type & range )

inlineinherited

Set range.

Parameters

range range

Definition at line 120 of file JRange.hh.

     {
       static_cast<range_type&>(*this) = range;
     }

template<class T, class JComparator_t = std::less<T>>

void JTOOLS::JRange< T, JComparator_t >::setRange	(	argument_type	x,
		argument_type	y
	)

inlineinherited

Set lower and uppper limit.

Parameters

x	lower limit
y	upper limit

Definition at line 132 of file JRange.hh.

     {
       this->first  = x;
       this->second = y;
     }

template<class T, class JComparator_t = std::less<T>>

template<class iterator_type , class value_type , class result_type >

void JTOOLS::JRange< T, JComparator_t >::setRange	(	iterator_type	__begin,
		iterator_type	__end,
		result_type value_type::*	member
	)

inlineinherited

Set lower and uppper limit.

Parameters

__begin	begin of data
__end	end of data
member	pointer to data member

Definition at line 147 of file JRange.hh.

     {
       setRange(getMaximum(), getMinimum());
 
       for (iterator_type i = __begin; i != __end; ++i) {
         include((*i).*member);
       }
     }

template<class T, class JComparator_t = std::less<T>>

template<class iterator_type , class value_type , class result_type >

void JTOOLS::JRange< T, JComparator_t >::setRange	(	iterator_type	__begin,
		iterator_type	__end,
		result_type(value_type::*)() const	function
	)

inlineinherited

Set lower and uppper limit.

Parameters

__begin	begin of data
__end	end of data
function	pointer to member method

Definition at line 165 of file JRange.hh.

     {
       setRange(getMaximum(), getMinimum());
 
       for (iterator_type i = __begin; i != __end; ++i) {
         include(((*i).*function)());
       }
     }

template<class T, class JComparator_t = std::less<T>>

T JTOOLS::JRange< T, JComparator_t >::getLowerLimit ( ) const

inlineinherited

Get lower limit.

Returns: lower limit

Definition at line 180 of file JRange.hh.

     {
       return this->first;
     }

template<class T, class JComparator_t = std::less<T>>

T JTOOLS::JRange< T, JComparator_t >::getUpperLimit ( ) const

inlineinherited

Get upper limit.

Returns: upper limit

Definition at line 191 of file JRange.hh.

     {
       return this->second;
     }

template<class T, class JComparator_t = std::less<T>>

void JTOOLS::JRange< T, JComparator_t >::setLowerLimit ( argument_type x )

inlineinherited

Set lower limit.

Parameters

x	lower limit

Definition at line 202 of file JRange.hh.

     {
       this->first  = x;
     }

template<class T, class JComparator_t = std::less<T>>

void JTOOLS::JRange< T, JComparator_t >::setUpperLimit ( argument_type y )

inlineinherited

Set upper limit.

Parameters

y	upper limit

Definition at line 213 of file JRange.hh.

     {
       this->second = y;
     }

template<class T, class JComparator_t = std::less<T>>

void JTOOLS::JRange< T, JComparator_t >::fixLowerLimit ( argument_type x )

inlineinherited

Fix lower limit.

The range is shifted to the given lower limit.

Parameters

x	lower limit

Definition at line 226 of file JRange.hh.

     {
       this->second += x - this->first;
       this->first   = x;
     }

template<class T, class JComparator_t = std::less<T>>

void JTOOLS::JRange< T, JComparator_t >::fixUpperLimit ( argument_type y )

inlineinherited

Fix upper limit.

The range is shifted to the given upper limit.

Parameters

y	upper limit

Definition at line 240 of file JRange.hh.

     {
       this->first  += y - this->second;
       this->second  = y;
     }

template<class T, class JComparator_t = std::less<T>>

T JTOOLS::JRange< T, JComparator_t >::getLength ( ) const

inlineinherited

Get length (difference between upper and lower limit).

Returns: length

Definition at line 267 of file JRange.hh.

     {
       return getUpperLimit() - getLowerLimit();
     }

template<class T, class JComparator_t = std::less<T>>

bool JTOOLS::JRange< T, JComparator_t >::is_valid ( ) const

inlineinherited

Check validity of range.

Returns: true if lower limit less than or equal to upper limit; else false

Definition at line 278 of file JRange.hh.

     {
       return !compare(getUpperLimit(), getLowerLimit());
     }

template<class T, class JComparator_t = std::less<T>>

bool JTOOLS::JRange< T, JComparator_t >::operator() ( argument_type x ) const

inlineinherited

Test whether value is inside range.

Parameters

x value

Returns: true if lower limit <= value <= upper limit; else false

Definition at line 290 of file JRange.hh.

     {
       return (!compare(x, getLowerLimit()) && 
               !compare(getUpperLimit(), x));
     }

template<class T, class JComparator_t = std::less<T>>

T JTOOLS::JRange< T, JComparator_t >::constrain ( argument_type x ) const

inlineinherited

Constrain value to range.

This method returns the original value if it is in this range, else lower limit if value < lower limit or upper limit if value > upper limit.

Parameters

x value

Returns: lower limit <= x <= upper limit

Definition at line 305 of file JRange.hh.

     {
       if (compare(x, getLowerLimit())) return getLowerLimit();
       if (compare(getUpperLimit(), x)) return getUpperLimit();
 
       return x;
     }

template<class T, class JComparator_t = std::less<T>>

bool JTOOLS::JRange< T, JComparator_t >::overlap ( const range_type & range ) const

inlineinherited

Test overlap with given range.

The result is equivalent to join(range).is_valid().

Parameters

range range

Returns: true if there is a non-zero overlap; else false

Definition at line 321 of file JRange.hh.

     {  
       return (compare(getLowerLimit(), range.getUpperLimit()) &&
               compare(range.getLowerLimit(), getUpperLimit()));
     }

template<class T, class JComparator_t = std::less<T>>

range_type JTOOLS::JRange< T, JComparator_t >::include ( argument_type x )

inlineinherited

Include given value to range.

The new lower limit is the minimim of the original lower limit and given value and the new upper limit is the maximum of the original upper limit and given value;

Parameters

x value

Returns: range

Definition at line 336 of file JRange.hh.

     {
       if (compare(x, getLowerLimit())) setLowerLimit(x);
       if (compare(getUpperLimit(), x)) setUpperLimit(x);
 
       return *this;
     }

template<class T, class JComparator_t = std::less<T>>

range_type& JTOOLS::JRange< T, JComparator_t >::join ( const range_type & range )

inlineinherited

Join ranges.

The new lower limit is the maximim of the two lower limits and the new upper limit is the minimum of the two upper limits. This operation results in an equal or smaller range and may result in an unphysical range (i.e. lower limit > upper limit).

Parameters

range range

Definition at line 354 of file JRange.hh.

     {  
       if (compare(getLowerLimit(), range.getLowerLimit())) setLowerLimit(range.getLowerLimit());
       if (compare(range.getUpperLimit(), getUpperLimit())) setUpperLimit(range.getUpperLimit());
 
       return *this;      
     }

template<class T, class JComparator_t = std::less<T>>

range_type& JTOOLS::JRange< T, JComparator_t >::combine ( const range_type & range )

inlineinherited

Combine ranges.

The new lower limit is the minimim of the two lower limits and the new upper limit is the maximum of the two upper limits. This operation results in an equal or larger range.

Parameters

range range

Definition at line 371 of file JRange.hh.

     {  
       if (compare(range.getLowerLimit(), getLowerLimit())) setLowerLimit(range.getLowerLimit());
       if (compare(getUpperLimit(), range.getUpperLimit())) setUpperLimit(range.getUpperLimit());
 
       return *this;
     }

template<class T, class JComparator_t = std::less<T>>

range_type& JTOOLS::JRange< T, JComparator_t >::add ( argument_type x )

inlineinherited

Add offset.

Parameters

x offset

Definition at line 385 of file JRange.hh.

     { 
       this->first  += x;
       this->second += x;
 
       return *this;
     }

template<class T, class JComparator_t = std::less<T>>

range_type& JTOOLS::JRange< T, JComparator_t >::add ( const range_type & range )

inlineinherited

Add offset.

The new lower limit is the sum of the two lower limits and the new upper limit is the sum of the two upper limits.

Parameters

range offset

Definition at line 415 of file JRange.hh.

     { 
       this->first  += range.getLowerLimit();
       this->second += range.getUpperLimit();
 
       return *this;
     }

template<class T, class JComparator_t = std::less<T>>

range_type& JTOOLS::JRange< T, JComparator_t >::sub ( argument_type x )

inlineinherited

Subtract offset.

Parameters

x offset

Definition at line 399 of file JRange.hh.

     { 
       this->first  -= x;
       this->second -= x;
 
       return *this;
     }

template<class T, class JComparator_t = std::less<T>>

range_type& JTOOLS::JRange< T, JComparator_t >::sub ( const range_type & range )

inlineinherited

Subtract offset.

The new lower limit is the difference of the two lower limits and the new upper limit is the difference of the two upper limits.

Parameters

range offset

Definition at line 431 of file JRange.hh.

     { 
       this->first  -= range.getLowerLimit();
       this->second -= range.getUpperLimit();
 
       return *this;
     }

template<class T, class JComparator_t = std::less<T>>

range_type& JTOOLS::JRange< T, JComparator_t >::mul ( const double factor )

inlineinherited

Multiply range.

Parameters

factor factor

Definition at line 445 of file JRange.hh.

     { 
       this->first  *= factor;
       this->second *= factor;
 
       return *this;
     }

template<class JFirst_t, class JSecond_t = JNullType>

JFirst_t& JMATH::JMath< JFirst_t, JSecond_t >::mul ( const JSecond_t & 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);
     }

template<class T, class JComparator_t = std::less<T>>

range_type& JTOOLS::JRange< T, JComparator_t >::div ( const double factor )

inlineinherited

Divide range.

Parameters

factor factor

Definition at line 459 of file JRange.hh.

     { 
       this->first  /= factor;
       this->second /= factor;
 
       return *this;
     }

template<class T, class JComparator_t = std::less<T>>

T JTOOLS::JRange< T, JComparator_t >::getN ( const double R ) const

inlineinherited

Get expected number of occurances of given rate within this interval.

Parameters

R rate

Returns: expectation value

Definition at line 474 of file JRange.hh.

     {
       return R * (getUpperLimit() - getLowerLimit());
     }

template<class T, class JComparator_t = std::less<T>>

static T JTOOLS::JRange< T, JComparator_t >::getMinimum ( )

inlinestaticinherited

Get minimum possible value.

Returns: minimum possible value

Definition at line 485 of file JRange.hh.

     { 
       return JMATH::JLimits<T>::min();
     }

template<class T, class JComparator_t = std::less<T>>

static T JTOOLS::JRange< T, JComparator_t >::getMaximum ( )

inlinestaticinherited

Get maximum possible value.

Returns: maximum possible value

Definition at line 496 of file JRange.hh.

     { 
       return JMATH::JLimits<T>::max();
     }

const key_type& JTOOLS::JPair< T , T >::getKey ( ) const

inlineinherited

Definition at line 123 of file JPair.hh.

123 { return this->first; }

JTOOLS::JPair< T, T >::first

T first

Definition: JPair.hh:128

const std::string& JUTM::JUTMGrid::getKey ( ) const

inlineinherited

Get key.

Returns: key

Definition at line 68 of file JUTMGrid.hh.

     {
       return key;
     }

const mapped_type& JTOOLS::JPair< T , T >::getValue ( ) const

inlineinherited

Definition at line 124 of file JPair.hh.

124 { return this->second; }

JTOOLS::JPair< T, T >::second

T second

Definition: JPair.hh:129

mapped_type& JTOOLS::JPair< T , T >::getValue ( )

inlineinherited

Definition at line 126 of file JPair.hh.

126 { return this->second; }

JTOOLS::JPair< T, T >::second

T second

Definition: JPair.hh:129

const std::string& JUTM::JUTMGrid::getWGS ( ) const

inlineinherited

Get WGS.

Returns: WGS

Definition at line 79 of file JUTMGrid.hh.

     {
       return wgs;
     }

const std::string& JUTM::JUTMGrid::getUTMZone ( ) const

inlineinherited

Get UTM zone.

Returns: UTM zone

Definition at line 90 of file JUTMGrid.hh.

     {
       return zone;
     }

std::string JUTM::JUTMGrid::toString ( ) const

inlineinherited

Convert UTM grid.

Returns: UTM grid

Definition at line 101 of file JUTMGrid.hh.

     {
       return (key + " " + wgs + " " + zone);
     }

static JUTMGrid JUTM::JUTMGrid::valueOf ( const std::string & buffer )

inlinestaticinherited

Extract UTM grid.

Parameters

buffer WGS and UTM zone

Returns: UTM grid

Definition at line 113 of file JUTMGrid.hh.

     {
       JUTMGrid grid;
 
       std::istringstream is(buffer);
 
       if (is >> grid)
         return grid;
       else
         throw JParseError("JUTMGrid::valueOf()");
     }

JPosition3D JUTM::JUTMPosition::getPosition ( ) const

inlineinherited

Get position.

Returns: position

Definition at line 76 of file JUTMPosition.hh.

     {
       return JPosition3D(east, north, z);
     }

double JUTM::JUTMPosition::getUTMEast ( ) const

inlineinherited

Get UTM east.

Returns: UTM East

Definition at line 87 of file JUTMPosition.hh.

     {
       return this->east;
     }

double JUTM::JUTMPosition::getUTMNorth ( ) const

inlineinherited

Get UTM north.

Returns: UTM North

Definition at line 98 of file JUTMPosition.hh.

     {
       return this->north;
     }

double JUTM::JUTMPosition::getUTMZ ( ) const

inlineinherited

Get UTM Z.

Returns: UTM Z

Definition at line 109 of file JUTMPosition.hh.

     {
       return this->z;
     }

template<class JClass_t, class JTypelist_t>

template<class JHead_t , class JTail_t >

static bool JLANG::JMultiEquals< JClass_t, JTypelist_t >::eq	(	const JClass_t &	first,
		const JClass_t &	second,
		const JType< JTypeList< JHead_t, JTail_t > > &	type
	)

inlinestaticprotectedinherited

Equals method for composite data types.

Parameters

first	first object
second	second object
type	type

Returns: true if two objects are equal; else false

Definition at line 42 of file JMultiEquals.hh.

     {
       return (static_cast<const JHead_t&>(first)  ==
               static_cast<const JHead_t&>(second) &&
               eq(first, second, JType<JTail_t>()));
     }

template<class JClass_t, class JTypelist_t>

template<class JHead_t >

static bool JLANG::JMultiEquals< JClass_t, JTypelist_t >::eq	(	const JClass_t &	first,
		const JClass_t &	second,
		const JType< JTypeList< JHead_t, JNullType > > &	type
	)

inlinestaticprotectedinherited

Equals method for composite data types.

Parameters

first	first object
second	second object
type	type

Returns: true if two objects are equal; else false

Definition at line 61 of file JMultiEquals.hh.

     {
       return (static_cast<const JHead_t&>(first) ==
               static_cast<const JHead_t&>(second));
     }

Friends And Related Function Documentation

std::istream& operator>>	(	std::istream &	in,
		JMonteCarloDetector &	detector
	)

friend

Read detector from input.

Parameters

in	input stream
detector	detector

Returns: input stream

Definition at line 281 of file JMonteCarloDetector.hh.

     {
       using namespace std;
 
       detector.clear();
 
 
       vector<OM>              OMaddress;
       vector<OM_cluster>      OMcluster;
       vector<OM_position>     OMposition;
       vector<String>          DetectorString;
       vector<OM_cluster_data> OMclusterdata;
       vector<LCM_logic>       LCMlogic;
 
 
       for (string key, buffer; in >> key; ) {
 
         if      (key == "OM:")
 
           OMaddress.push_back(OM(in));
 
         else if (key == "OM_cluster:")
 
           OMcluster.push_back(OM_cluster(in));
 
         else if (key == "OM_position:")
 
           OMposition.push_back(OM_position(in));
 
         else if (key == "string:")
 
           DetectorString.push_back(String(in));
 
         else if (key == "OM_cluster_data:")
 
           OMclusterdata.push_back(OM_cluster_data(in));
 
         else if (key == "LCM_logic:")
 
           LCMlogic.push_back(LCM_logic(in));
 
         else
 
           getline(in, buffer);
       }
 
 
       // invert LCM_Logic data for fast lookup
 
       vector<LCM_reverse_logic> LCMreverselogic;
 
       for (vector<LCM_logic>::const_iterator i  = LCMlogic.begin(); i != LCMlogic.end(); ++i) {
         
         const int lcm_id = i->getID();
 
         for (LCM_logic::const_iterator j = i->begin(); j != i->end(); ++j) {
 
           const int id     = *j;
           const int pmt_id = distance(i->begin(),j) + 1;
 
           LCMreverselogic.push_back(LCM_reverse_logic(id, lcm_id, pmt_id));
         }
       }
 
 
       // sort data according to identifier
 
       sort(OMaddress      .begin(), OMaddress      .end());
       sort(OMcluster      .begin(), OMcluster      .end());
       sort(OMposition     .begin(), OMposition     .end());
       sort(DetectorString .begin(), DetectorString .end());
       sort(LCMreverselogic.begin(), LCMreverselogic.end());
 
 
       vector<OM_cluster>     ::const_iterator   OMclusterIterator;
       vector<String>         ::const_iterator   StringIterator;
       vector<int>            ::const_iterator   addressIterator;
       vector<OM>             ::const_iterator   OMIterator;
       vector<OM_position>    ::const_iterator   OMpositionIterator;
       vector<OM_cluster_data>::const_iterator   OMclusterdataIterator;
 
 
       for (OMclusterIterator  = OMcluster.begin(); 
            OMclusterIterator != OMcluster.end(); 
            OMclusterIterator++) {
         
         StringIterator = find(DetectorString.begin(),
                               DetectorString.end(),
                               OMclusterIterator->string_id);
         
         if (StringIterator != DetectorString.end()) {
           
           const JModuleLocation location(StringIterator->getID(), -1);
 
           for (addressIterator  = OMclusterIterator->begin(); 
                addressIterator != OMclusterIterator->end(); 
                addressIterator++) {
             
             OMIterator = find(OMaddress.begin(),
                               OMaddress.end(),
                               *addressIterator);
             
             if (OMIterator != OMaddress.end()) {
               
               OMpositionIterator = find(OMposition.begin(),
                                         OMposition.end(),
                                         OMIterator->address);
               
               if (OMpositionIterator != OMposition.end()) {
                 
                 // rotate string
                 
                 double ct = cos(StringIterator->twist * PI / 180.0);
                 double st = sin(StringIterator->twist * PI / 180.0);
                 
                 double x  =  ct * OMpositionIterator->x - st * OMpositionIterator->y;
                 double y  =  st * OMpositionIterator->x + ct * OMpositionIterator->y;
 
                 JPosition3D  pos(StringIterator->x + x,
                                  StringIterator->y + y,
                                  StringIterator->z + OMpositionIterator->z + OMclusterIterator->z);
                 
                 JDirection3D dir(JAngle3D(OMpositionIterator->theta, 
                                           StringIterator->twist * PI / 180.0 + OMpositionIterator->phi));
                 
                 int lcm_id  = OMclusterIterator->getID();
                 int pmt_id  = OMIterator->address;
                 
                 if (detector.use_logic) {
                   
                   vector<LCM_reverse_logic>::const_iterator p = find(LCMreverselogic.begin(), 
                                                                      LCMreverselogic.end(), 
                                                                      OMIterator->getID());
                   
                   if (p != LCMreverselogic.end()) {
                     lcm_id = p->lcm_id;
                     pmt_id = p->pmt_id;
                   }
                 }
                 
                 detector.put(lcm_id, pmt_id - 1, location, JPMT(OMIterator->getID(), JAxis3D(pos, dir)));
               }
             }
           }
         }
       }
 
 
       // 'event' data
 
       for (OMclusterdataIterator  = OMclusterdata.begin(); 
            OMclusterdataIterator != OMclusterdata.end(); 
            OMclusterdataIterator++) {
           
         const JPosition3D    P(OMclusterdataIterator->x,
                                OMclusterdataIterator->y,
                                OMclusterdataIterator->z);
         
         const JEulerMatrix3D R(JEulerAngle3D(OMclusterdataIterator->theta,
                                              OMclusterdataIterator->phi,
                                              OMclusterdataIterator->psi));
         
         OMclusterIterator = find(OMcluster.begin(),
                                  OMcluster.end(),
                                  OMclusterdataIterator->getID());
 
         if (OMclusterIterator != OMcluster.end()) {
 
           const JModuleLocation location(OMclusterIterator->string_id, -1);
 
           for (addressIterator  = OMclusterIterator->begin(); 
                addressIterator != OMclusterIterator->end(); 
                addressIterator++) {
 
             OMIterator = find(OMaddress.begin(),
                               OMaddress.end(),
                               *addressIterator);
               
             if (OMIterator != OMaddress.end()) {
                 
               OMpositionIterator = find(OMposition.begin(),
                                         OMposition.end(),
                                         OMIterator->address);
                 
               if (OMpositionIterator != OMposition.end()) {
                   
                 JPosition3D  pos(OMpositionIterator->x,
                                  OMpositionIterator->y,
                                  OMpositionIterator->z);
                   
                 JDirection3D dir(JAngle3D(OMpositionIterator->theta, 
                                           OMpositionIterator->phi));
                                   
                 pos.transform(R);
                 dir.transform(R);
                   
                 pos += P;
                   
                 int lcm_id  = OMclusterIterator->getID();
                 int pmt_id  = OMIterator->address;
 
                 if (detector.use_logic) {
 
                   vector<LCM_reverse_logic>::const_iterator p = find(LCMreverselogic.begin(), 
                                                                      LCMreverselogic.end(), 
                                                                      OMIterator->getID());
                   
                   if (p != LCMreverselogic.end()) {
                     lcm_id = p->lcm_id;
                     pmt_id = p->pmt_id;
                   }
                 }
 
                 detector.put(lcm_id, pmt_id - 1, location, JPMT(OMIterator->getID(), JAxis3D(pos, dir)));
               }
             }
           }
         }
       }
 
       // set floor identifier
 
       sort(detector.begin(), detector.end(), JMonteCarloDetector::compare);
 
       int string = -1;
       int floor  = -1;
       
       for (JDetector::iterator module = detector.begin(); module != detector.end(); ++module) {
         
         if (module->getString() != string) {
           string = module->getString();
           floor  = 0;
         }
 
         module->setLocation(JModuleLocation(string,++floor));
       }
       
       return in;
     }

Member Data Documentation

bool JDETECTOR::JMonteCarloDetector::use_logic

private

Definition at line 636 of file JMonteCarloDetector.hh.

JComment JDETECTOR::JDetector::comment

inherited

Definition at line 415 of file JDetector.hh.

int JLANG::JObjectID::__id

protectedinherited

Definition at line 159 of file JObjectID.hh.

std::string JEEP::JVersion::__variant

protectedinherited

Definition at line 171 of file JVersion.hh.

template<class T, class JComparator_t = std::less<T>>

const JRange< T, JComparator_t > JTOOLS::JRange< T, JComparator_t >::DEFAULT_RANGE

staticinherited

Default range.

This range corresponds to an unphysical range.

Definition at line 506 of file JRange.hh.

template<class T, class JComparator_t = std::less<T>>

JComparator_t JTOOLS::JRange< T, JComparator_t >::compare

inherited

Function object.

Parameters

first	first argument
second	second argument

Returns: true if first < second; else false

Definition at line 516 of file JRange.hh.

T JTOOLS::JPair< T , T >::first

inherited

Definition at line 128 of file JPair.hh.

T JTOOLS::JPair< T , T >::second

inherited

Definition at line 129 of file JPair.hh.

std::string JUTM::JUTMGrid::key

protectedinherited

Definition at line 183 of file JUTMGrid.hh.

std::string JUTM::JUTMGrid::wgs

protectedinherited

Definition at line 184 of file JUTMGrid.hh.

std::string JUTM::JUTMGrid::zone

protectedinherited

Definition at line 185 of file JUTMGrid.hh.

double JUTM::JUTMPosition::east

protectedinherited

Definition at line 169 of file JUTMPosition.hh.

double JUTM::JUTMPosition::north

protectedinherited

Definition at line 170 of file JUTMPosition.hh.

double JUTM::JUTMPosition::z

protectedinherited

Definition at line 171 of file JUTMPosition.hh.

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

JMonteCarloDetector.hh

Classes

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Private Attributes

Friends

Detailed Description

Member Typedef Documentation

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation