JDETECTOR::JPMTParametersMap Class Reference

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 which is reserved for the so-called light scaling;
• 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. It is set (read overwritten) in the JTriggerEfficiency[RunByRun].sh scripts.

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 89 of file JPMTParametersMap.hh.

Constructor & Destructor Documentation

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

inline

Constructor.

Parameters

parameters

PMT parameters

Definition at line 100 of file JPMTParametersMap.hh.

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

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

inline

Constructor.

Parameters

file_name

file name

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

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 140 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;
    }

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

inline

Get QE of given PMT.

Parameters

id	PMT identifier

Returns

QE

Definition at line 163 of file JPMTParametersMap.hh.

    {
      return getPMTParameters(id).QE;
    }

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 175 of file JPMTParametersMap.hh.

    {
      using JDETECTOR::getHitProbability;

      return getHitProbability(getQE(id), this->mu);
    }

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

inline

Get default PMT parameters.

Returns

PMT parameters

Definition at line 188 of file JPMTParametersMap.hh.

    {
      return defaultPMTParameters;
    }

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

inline

Set default PMT parameters.

Parameters

parameters

PMT parameters

Definition at line 199 of file JPMTParametersMap.hh.

    {
      defaultPMTParameters = parameters;
    }

double JDETECTOR::JPMTParametersMap::getQE ( ) const

inline

Get global QE.

Returns

QE

Definition at line 210 of file JPMTParametersMap.hh.

    {
      return this->QE;
    }

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

inline

Set global QE.

Parameters

QE

QE

Definition at line 221 of file JPMTParametersMap.hh.

    {
      this->QE = QE;
    }

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 232 of file JPMTParametersMap.hh.

    {
      return this->mu;
    }

bool JDETECTOR::JPMTParametersMap::is_valid ( ) const

inline

Check validity of PMT parameters.

Returns

true if valid; else false

Definition at line 243 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;
      }
    }

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 292 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;
      }
    }

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

inlinestatic

Get equation parameters.

Returns

equation parameters

Definition at line 388 of file JPMTParametersMap.hh.

    {
      static JEquationParameters parameters;

      return parameters;
    }

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

inlinestatic

Set equation parameters.

Parameters

equation

equation parameters

Definition at line 401 of file JPMTParametersMap.hh.

    {
      getEquationParameters() = equation;
    }

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));
    }

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));
    }

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;
    }

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

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 327 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;
    }

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 363 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

JComment JDETECTOR::JPMTParametersMap::comment

Definition at line 407 of file JPMTParametersMap.hh.

JPMTParameters JDETECTOR::JPMTParametersMap::defaultPMTParameters

protected

Default PMT parameters.

Definition at line 413 of file JPMTParametersMap.hh.

double JDETECTOR::JPMTParametersMap::QE

protected

Global QE.

Definition at line 418 of file JPMTParametersMap.hh.

double JDETECTOR::JPMTParametersMap::mu

protected

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

Definition at line 424 of file JPMTParametersMap.hh.

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