Auxiliary class for map of PMT parameters. More...

#include <JPMTParametersMap.hh>

Inheritance diagram for JDETECTOR::JPMTParametersMap:

Classes
class	JPMTParametersMapHelper
	Auxiliary class for I/O of PMT parameters map. More...

Public Member Functions
	JPMTParametersMap (const JPMTParameters &parameters=JPMTParameters())
	Constructor. More...

	JPMTParametersMap (const char *const file_name)
	Constructor. More...

const JPMTParameters &	getPMTParameters (const JPMTIdentifier &id) const
	Get PMT parameters. More...

double	getQE (const JPMTIdentifier &id) const
	Get QE of given PMT. More...

double	getHitProbability (const JPMTIdentifier &id) const
	Get ratio of hit probabilities of given PMT. More...

const JPMTParameters &	getDefaultPMTParameters () const
	Get default PMT parameters. More...

void	setDefaultPMTParameters (const JPMTParameters &parameters)
	Set default PMT parameters. More...

double	getQE () const
	Get global QE. More...

void	setQE (const double QE)
	Set global QE. More...

double	getMu () const
	Get expectation value for number of photo-electrons given two-fold (or higher) coincidence rate. More...

bool	is_valid () const
	Check validity of PMT parameters. More...

void	convertHitProbabilityToQE (const double mu)
	Convert the hit probabilities to QEs for given expectation value. More...

void	load (const char *file_name)
	Load from input file. More...

void	store (const char *file_name) const
	Store to output file. More...

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

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

static void	Throw (const bool option)
	Enable/disable throw option. More...

static int	Throw (const JException &error, const int value=-1)
	Throw exception or return error. More...

Public Attributes
JComment	comment

Protected Attributes
JPMTParameters	defaultPMTParameters
	Default PMT parameters. More...

double	QE
	Global QE. More...

double	mu
	Expectation value for number of photo-electrons given two-fold (or higher) coincidence rate. More...

Static Protected Attributes
static bool	do_throw = true
	throw option More...

Friends
std::istream &	operator>> (std::istream &in, JPMTParametersMap &object)
	Stream input. More...

std::ostream &	operator<< (std::ostream &out, const JPMTParametersMap &object)
	Stream output. More...

Detailed Description

Auxiliary class for map of PMT parameters.

Input syntax as follows:

           QE=<value>
           mu=<value>

           %.<key>=<value>

           PMT=<module> <pmt> <PMT parameters>

           pmt=<module> <pmt> <key>=<value>[, <key>=<value>]

           file=<file name>

where

QE refers to the global quantum efficiency;
mu refers to the expectation value for number of photo-electrons given two-fold (or higher) coincidence rate;
'%' to the default values;
<module> to the module identifier;
<pmt> to the readout channel;
<key> to one of the data members of the JDETECTOR::JPMTParameters data structure; and
<value> to the corresponding value.

The global QE is applied as a multiplicative factor to the default and PMT specific QEs.

The default values are provided by the method getDefaultPMTParameters and can also be modified using method setDefaultPMTParameters. If no optional values are given, the default values are maintained.

Finally, <PMT parameters> corresponds to the list of PMT parameter values that should be compatible with the I/O methods of the JDETECTOR::JPMTParameters data structure.

Note that multiple input sequences should be seperated by a semicolumn ';' or be terminated by a newline. When multiple input sequences appear for the same target, the last sequence prevails.

Definition at line 88 of file JPMTParametersMap.hh.

Constructor & Destructor Documentation

◆ JPMTParametersMap() [1/2]

JDETECTOR::JPMTParametersMap::JPMTParametersMap ( const JPMTParameters & parameters = JPMTParameters() )

inline

Constructor.

Parameters

parameters PMT parameters

Definition at line 99 of file JPMTParametersMap.hh.

                                                                            :
       JPMTParametersMap_t(),
       QE(1.0),
       mu(0.0)
     {
       setDefaultPMTParameters(parameters);
     }

◆ JPMTParametersMap() [2/2]

JDETECTOR::JPMTParametersMap::JPMTParametersMap ( const char *const file_name )

inline

Constructor.

Parameters

file_name file name

Definition at line 113 of file JPMTParametersMap.hh.

                                                    : 
       JPMTParametersMap_t(),
       QE(1.0),
       mu(0.0)
     {
       using namespace std;
  
       ifstream in(file_name);
  
       in >> *this;
  
       in.close();
     }

Member Function Documentation

◆ getPMTParameters()

const JPMTParameters& JDETECTOR::JPMTParametersMap::getPMTParameters ( const JPMTIdentifier & id ) const

inline

Get PMT parameters.

This method returns the default PMT parameters if the parameters corresponding to the given PMT identifier have not been defined.

Note that the value of QE as part of the return value has been scaled with the global QE of the PMT parameters map.

Parameters

id	PMT identifier

Returns: PMT parameters

Definition at line 139 of file JPMTParametersMap.hh.

     {
       static JPMTParameters parameters;
       
       JPMTParametersMap_t::const_iterator i = find(id);
  
       if (i != end())
         parameters = i->second;
       else
         parameters = getDefaultPMTParameters();
  
       parameters.QE *= this->QE;
       
       return parameters;
     }

◆ getQE() [1/2]

double JDETECTOR::JPMTParametersMap::getQE ( const JPMTIdentifier & id ) const

inline

Get QE of given PMT.

Parameters

id	PMT identifier

Returns: QE

Definition at line 162 of file JPMTParametersMap.hh.

     {
       return getPMTParameters(id).QE;
     }

◆ getHitProbability()

double JDETECTOR::JPMTParametersMap::getHitProbability ( const JPMTIdentifier & id ) const

inline

Get ratio of hit probabilities of given PMT.

Parameters

id	PMT identifier

Returns: ratio

Definition at line 174 of file JPMTParametersMap.hh.

     {
       using JDETECTOR::getHitProbability;
  
       return getHitProbability(getQE(id), this->mu);
     }

◆ getDefaultPMTParameters()

const JPMTParameters& JDETECTOR::JPMTParametersMap::getDefaultPMTParameters ( ) const

inline

Get default PMT parameters.

Returns: PMT parameters

Definition at line 187 of file JPMTParametersMap.hh.

     {
       return defaultPMTParameters;
     }

◆ setDefaultPMTParameters()

void JDETECTOR::JPMTParametersMap::setDefaultPMTParameters ( const JPMTParameters & parameters )

inline

Set default PMT parameters.

Parameters

parameters PMT parameters

Definition at line 198 of file JPMTParametersMap.hh.

     {
       defaultPMTParameters = parameters;
     }

◆ getQE() [2/2]

double JDETECTOR::JPMTParametersMap::getQE ( ) const

inline

Get global QE.

Returns: QE

Definition at line 209 of file JPMTParametersMap.hh.

     {
       return this->QE;
     }

◆ setQE()

void JDETECTOR::JPMTParametersMap::setQE ( const double QE )

inline

Set global QE.

Parameters

QE QE

Definition at line 220 of file JPMTParametersMap.hh.

     {
       this->QE = QE;
     }

◆ getMu()

double JDETECTOR::JPMTParametersMap::getMu ( ) const

inline

Get expectation value for number of photo-electrons given two-fold (or higher) coincidence rate.

Returns: expectation value

Definition at line 231 of file JPMTParametersMap.hh.

     {
       return this->mu;
     }

◆ is_valid()

bool JDETECTOR::JPMTParametersMap::is_valid ( ) const

inline

Check validity of PMT parameters.

Returns: true if valid; else false

Definition at line 242 of file JPMTParametersMap.hh.

     {
       try {
  
         if (this->QE < 0) {
           THROW(JException, "Invalid global QE " << QE);
         }
  
         if (this->mu < 0) {
           THROW(JException, "Invalid expectation value number of photo-electrons " << mu);
         }
         
         if (!defaultPMTParameters.is_valid()) {
           THROW(JException, "Invalid default PMT parameters");
         }
         
         for (const_iterator i = this->begin(); i != this->end(); ++i) {
           
           if (!i->second.is_valid()) {
             THROW(JException, "Invalid PMT parameters at " << i->first);
           }
  
           if (i->second.QE * this->QE > 1.0 ||
               i->second.QE * this->QE < 0.0) {
             THROW(JException, "Invalid PMT QE at " << i->first << ' ' << i->second.QE * this->QE);
           }
         }
  
         return true;
       }
       catch(const JException& error) {
  
         Throw(error);
  
         return false;
       }
     }

◆ convertHitProbabilityToQE()

void JDETECTOR::JPMTParametersMap::convertHitProbabilityToQE ( const double mu )

inline

Convert the hit probabilities to QEs for given expectation value.

The expectation value corresponds to the number of photo-electrons given the multiplicity range of coincidences as specified at JCalibrateK40.cc.
The survival probability is taken into account in the conversion (see method JDETECTOR::getSurvivalProbability).

Parameters

mu	expectation value

Definition at line 291 of file JPMTParametersMap.hh.

     {
       using namespace std;
       using JDETECTOR::getHitProbability;
       using JDETECTOR::getQE;
  
       this->mu = mu;
  
       const double Pmax = getMaximalHitProbability(this->mu);
  
       for (JPMTParametersMap::const_iterator i = this->begin(); i != this->end(); ++i) {
         if (i->second.QE > Pmax) {
           THROW(JValueOutOfRange, "Hit probability PMT " << i->first << ' ' << i->second.QE << " > maximum probability given expectation value " << this->mu << endl);
         }
       }
       
       for (JPMTParametersMap::iterator i = this->begin(); i != this->end(); ++i) {
         
         const double P = getSurvivalProbability(i->second);
         
         if (i->second.QE > 0.0 && P > 0.0)
           i->second.QE = getQE(i->second.QE, this->mu) / P; 
         else
           i->second.QE = 0.0;
       }
     }

◆ getEquationParameters()

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

inlinestatic

Get equation parameters.

Returns: equation parameters

Definition at line 387 of file JPMTParametersMap.hh.

     {
       static JEquationParameters parameters;
  
       return parameters;
     }

◆ setEquationParameters()

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

inlinestatic

Set equation parameters.

Parameters

equation equation parameters

Definition at line 400 of file JPMTParametersMap.hh.

     {
       getEquationParameters() = equation;
     }

◆ load()

template<class T>

void JLANG::JObjectStreamIO< T >::load ( const char * file_name )

inlineinherited

Load from input file.

Parameters

file_name file name

Definition at line 30 of file JObjectStreamIO.hh.

     {
       JLANG::load<std::ifstream>(file_name, static_cast<T&>(*this));
     }

◆ store()

template<class T>

void JLANG::JObjectStreamIO< T >::store ( const char * file_name ) const

inlineinherited

Store to output file.

Parameters

file_name file name

Definition at line 41 of file JObjectStreamIO.hh.

     {
       JLANG::store<std::ofstream>(file_name, static_cast<const T&>(*this));
     }

◆ Throw() [1/2]

template<class T>

static void JLANG::JThrow< T >::Throw ( const bool option )

inlinestaticinherited

Enable/disable throw option.

Parameters

option true enable; false disable

Definition at line 37 of file JThrow.hh.

     {
       do_throw = option;
     }

◆ Throw() [2/2]

template<class T>

static int JLANG::JThrow< T >::Throw	(	const JException &	error,
		const int	value = `-1`
	)

inlinestaticinherited

Throw exception or return error.

Parameters

error	exception
value	return code

Returns: return code

Definition at line 50 of file JThrow.hh.

     {
       using namespace std;
       
       if (do_throw) {
         throw error;
       }
  
       cerr << error.what() << endl;
       
       return value;
     }

Friends And Related Function Documentation

◆ operator>>

std::istream& operator>>	(	std::istream &	in,
		JPMTParametersMap &	object
	)

friend

Stream input.

Parameters

in	input stream
object	PMT parameters map

Returns: input stream

Definition at line 326 of file JPMTParametersMap.hh.

     {
       using namespace std;
       using namespace JPP;
       
       JStringStream is(in);
  
       if (getFileStatus(is.str().c_str())) {
         is.load();
       }
  
       JProperties properties(getEquationParameters(), 1);
         
       JPMTParametersMapHelper  helper(object, getEquationParameters());
       JProperties              demo = object.defaultPMTParameters.getProperties();
       
       properties["PMT"] = static_cast<JPMTParametersMap_t&>(object);
       properties["pmt"] = helper;
       properties["%"]   = demo;
       properties["QE"]  = object.QE;
       properties["mu"]  = object.mu;
  
       is >> object.comment;
       is >> properties;
       
       return in;
     }

◆ operator<<

std::ostream& operator<<	(	std::ostream &	out,
		const JPMTParametersMap &	object
	)

friend

Stream output.

Parameters

out	output stream
object	PMT parameters map

Returns: output stream

Definition at line 362 of file JPMTParametersMap.hh.

     {
       using namespace JLANG;
       
       JProperties properties(getEquationParameters());
  
       JProperties demo = object.defaultPMTParameters.getProperties();
  
       properties["PMT"] = static_cast<const JPMTParametersMap_t&>(object);
       properties["%"]   = demo;
  
       out << object.comment;
       out << setequation(getEquationParameters()) << JEquation::make_equation("QE", object.QE);
       out << setequation(getEquationParameters()) << JEquation::make_equation("mu", object.mu);
       out << properties;
  
       return out;
     }

Member Data Documentation

◆ comment

JComment JDETECTOR::JPMTParametersMap::comment

Definition at line 406 of file JPMTParametersMap.hh.

◆ defaultPMTParameters

JPMTParameters JDETECTOR::JPMTParametersMap::defaultPMTParameters

protected

Default PMT parameters.

Definition at line 412 of file JPMTParametersMap.hh.

◆ QE

double JDETECTOR::JPMTParametersMap::QE

protected

Global QE.

Definition at line 417 of file JPMTParametersMap.hh.

◆ mu

double JDETECTOR::JPMTParametersMap::mu

protected

Expectation value for number of photo-electrons given two-fold (or higher) coincidence rate.

Definition at line 423 of file JPMTParametersMap.hh.

◆ do_throw

template<class T>

bool JLANG::JThrow< T >::do_throw = true

staticprotectedinherited

throw option

Set default throw option to true.

Definition at line 28 of file JThrow.hh.

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

JPMTParametersMap.hh

Classes

Public Member Functions

Static Public Member Functions

Public Attributes

Protected Attributes

Static Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ JPMTParametersMap() [1/2]

◆ JPMTParametersMap() [2/2]

Member Function Documentation

◆ getPMTParameters()

◆ getQE() [1/2]

◆ getHitProbability()

◆ getDefaultPMTParameters()

◆ setDefaultPMTParameters()

◆ getQE() [2/2]

◆ setQE()

◆ getMu()

◆ is_valid()

◆ convertHitProbabilityToQE()

◆ getEquationParameters()

◆ setEquationParameters()

◆ load()

◆ store()

◆ Throw() [1/2]

◆ Throw() [2/2]

Friends And Related Function Documentation

◆ operator>>

◆ operator<<

Member Data Documentation

◆ comment

◆ defaultPMTParameters

◆ QE

◆ mu

◆ do_throw