Auxiliary classes and methods for PMT calibration. More...

Classes
struct	data_type
	Data structure for measured coincidence rates of all pairs of PMTs in optical module. More...

struct	JCombinatorics_t
	PMT combinatorics for optical module. More...

class	JFit
	Fit. More...

struct	JFitToT
	Parametrisation of time-over-threshold distribution. More...

struct	JFitToTParameters
	Fit parameters for two-fold coincidence rate due to K40. More...

struct	JHVInterpolator
	Auxiliary data structure to store high-voltage versus gain data and interpolate the nominal high-voltage. More...

struct	JK40Parameters
	Fit parameters for two-fold coincidence rate due to K40. More...

struct	JK40Parameters_t
	Fit parameters for two-fold coincidence rate due to K40. More...

struct	JModel
	Fit model. More...

struct	JModel_t
	Fit model. More...

struct	JModuleRouter_t
	Module router. More...

struct	JPairwiseComparator
	Auxiliary class to sort pairs of PMT addresses within optical module. More...

class	JParameter_t
	Auxiliary class for fit parameter with optional limits. More...

struct	JParameters
	Auxiliary class for PMT parameters including threshold. More...

struct	JPMTParameters_t
	Fit parameters for single PMT. More...

struct	JRouter_t
	Interface for routing module identifier to index and vice versa. More...

struct	JStringRouter_t
	String router. More...

struct	JTDC_t
	Auxiliary class for TDC constraints. More...

struct	rate_type
	Data structure for measured coincidence rate of pair of PMTs. More...

Typedefs
typedef JCombinatorics::pair_type	pair_type

Enumerations
enum	JOption_t { FIT_PMTS_t = 1 , FIT_PMTS_AND_ANGULAR_DEPENDENCE_t = 2 , FIT_PMTS_AND_BACKGROUND_t = 3 , FIT_PMTS_QE_FIXED_t = 4 , FIT_MODEL_t = 5 }
	Fit options. More...

Functions
double	getP (const double E1, const double E2, const double ED, const int M_min, const int M_max)
	Get coincidence probability of two PMTs within one module due to random background.

Variables
static const char *const	weights_hist_t = "weights_hist"
	Histogram naming.

static const char *const	W1_overall_t = "W1_overall"
	Named bin for duration of the run.

static const char *const	WS_t = "WS"
	Named bin for time residual bin width.

static const char *const	WB_t = "WB"
	Named bin for value of TMax_ns in JCalibrateK40.

static const char *const	_2S = ".2S"
	Name extension for 2D counts.

static const char *const	_1B = ".1B"
	Name extension for 1D background.

static const char *const	_1L = ".1L"
	Name extension for 1D live time.

static const char *const	_2R = ".2R"
	Name extension for 2D rate measured.

static const char *const	_2F = ".2F"
	Name extension for 2F rate fitted.

static const char *const	_1D = ".1D"
	Name extension for 1D time offset.

static const char *const	_1F = ".1F"
	Name extension for 1D time offset fit.

const char *const	module_t = "module"
	routing by module

const char *const	string_t = "string"
	routing by string

static const char *const	_2SToT = ".2SToT"
	Histogram naming.

static const int	INVALID_INDEX = -1
	invalid index

static double	TEROSTAT_DZ = 0.40
	maximal PMT inclination

static double	TEROSTAT_R1 = 1.00
	scaling factor

static double	BELL_SHAPE = 1.55
	Bell shape.

static const std::string	FITTOT_SUFFIX = ".1ToT"

static const std::string	FITTOT_FNAME = "fittot"

static const char *const	FITTOT_GAIN_PARNAME = "gain"

static const char *const	FITTOT_GAINSPREAD_PARNAME = "gainSpread"

static const char *const	FITTOT_NORMALIZATION_PARNAME = "normalization"

static const double	FITTOT_GAIN_MIN = 0.25
	Default minimal gain.

static const double	FITTOT_GAIN_MAX = 3.00
	Default maximal gain.

static const double	FITTOT_GAINSPREAD_MIN = 0.05
	Default minimal gain spread.

static const double	FITTOT_GAINSPREAD_MAX = 2.00
	Default maximal gain spread.

static const std::string	HVINTERPOLATOR_SUFFIX = ".HVxG"

static const char *const	HVINTERPOLATOR_DATA = "HVINTERPOLATOR_DATA"

static const char *const	HVINTERPOLATOR_RESULT = "HVINTERPOLATOR_RESULT"

Detailed Description

Auxiliary classes and methods for PMT calibration.

Author: mdejong; mkarel, bjung; bjung; acreusot

Typedef Documentation

◆ pair_type

typedef JCombinatorics::pair_type JCALIBRATE::pair_type

Definition at line 41 of file JCalibrateK40.hh.

Enumeration Type Documentation

◆ JOption_t

enum JCALIBRATE::JOption_t

Fit options.

Enumerator
FIT_PMTS_t	fit parameters of PMTs
FIT_PMTS_AND_ANGULAR_DEPENDENCE_t	fit parameters of PMTs and angular dependence of K40 rate
FIT_PMTS_AND_BACKGROUND_t	fit parameters of PMTs and background
FIT_PMTS_QE_FIXED_t	fit parameters of PMTs with QE fixed
FIT_MODEL_t	fit parameters of K40 rate and TTSs of PMTs

Definition at line 52 of file JFitK40.hh.

                 {
    FIT_PMTS_t                         =   1,     //!< fit parameters of PMTs
    FIT_PMTS_AND_ANGULAR_DEPENDENCE_t  =   2,     //!< fit parameters of PMTs and angular dependence of K40 rate
    FIT_PMTS_AND_BACKGROUND_t          =   3,     //!< fit parameters of PMTs and background
    FIT_PMTS_QE_FIXED_t                =   4,     //!< fit parameters of PMTs with QE fixed
    FIT_MODEL_t                        =   5      //!< fit parameters of K40 rate and TTSs of PMTs
  };

Function Documentation

◆ getP()

double JCALIBRATE::getP	(	const double	E1,
		const double	E2,
		const double	ED,
		const int	M_min,
		const int	M_max )

inline

Get coincidence probability of two PMTs within one module due to random background.

Parameters

E1	expected counts of 1st PMT
E2	expected counts of 2nd PMT
ED	expected counts of module, excluding PMT 1 and PMT 2
M_min	multiplicity lower limit
M_max	multiplicity upper limit

Returns: probability

Definition at line 136 of file JCalibrateK40.hh.

  {
    // 1st, always calculate the probability for twofold multiplicity
 
    double P = E1 * E2 * exp(-(E1 + E2 + ED));
      
    // 2nd, calculate the probability of multiplicity at lower limit
 
    int MD = 1;
 
    for ( ; MD + 2 <= M_min; ++MD) {
      P *= ED/MD;
    }
  
    // 3rd, add the probabilities of multiplicities up to upper limit
 
    double P_sum = P;
 
    for ( ; MD + 2 <= M_max; ++MD) {
      P     *= ED/MD;
      P_sum += P;
    }
  
    return P_sum;
  } 

Variable Documentation

◆ weights_hist_t

const char* const JCALIBRATE::weights_hist_t = "weights_hist"

static

Histogram naming.

Name of histogram with normalisation constants.

Definition at line 27 of file JCalibrateK40.hh.

◆ W1_overall_t

const char* const JCALIBRATE::W1_overall_t = "W1_overall"

static

Named bin for duration of the run.

Definition at line 29 of file JCalibrateK40.hh.

◆ WS_t

const char* const JCALIBRATE::WS_t = "WS"

static

Named bin for time residual bin width.

Definition at line 30 of file JCalibrateK40.hh.

◆ WB_t

const char* const JCALIBRATE::WB_t = "WB"

static

Named bin for value of TMax_ns in JCalibrateK40.

Definition at line 31 of file JCalibrateK40.hh.

◆ _2S

const char* const JCALIBRATE::_2S = ".2S"

static

Name extension for 2D counts.

Definition at line 33 of file JCalibrateK40.hh.

◆ _1B

const char* const JCALIBRATE::_1B = ".1B"

static

Name extension for 1D background.

Definition at line 34 of file JCalibrateK40.hh.

◆ _1L

const char* const JCALIBRATE::_1L = ".1L"

static

Name extension for 1D live time.

Definition at line 35 of file JCalibrateK40.hh.

◆ _2R

const char* const JCALIBRATE::_2R = ".2R"

static

Name extension for 2D rate measured.

Definition at line 36 of file JCalibrateK40.hh.

◆ _2F

const char* const JCALIBRATE::_2F = ".2F"

static

Name extension for 2F rate fitted.

Definition at line 37 of file JCalibrateK40.hh.

◆ _1D

const char* const JCALIBRATE::_1D = ".1D"

static

Name extension for 1D time offset.

Definition at line 38 of file JCalibrateK40.hh.

◆ _1F

const char* const JCALIBRATE::_1F = ".1F"

static

Name extension for 1D time offset fit.

Definition at line 39 of file JCalibrateK40.hh.

◆ module_t

const char* const JCALIBRATE::module_t = "module"

routing by module

Definition at line 172 of file JCalibrateMuon.hh.

◆ string_t

const char* const JCALIBRATE::string_t = "string"

routing by string

Definition at line 173 of file JCalibrateMuon.hh.

◆ _2SToT

const char* const JCALIBRATE::_2SToT = ".2SToT"

static

Histogram naming.

Name extension for 2D counts

Definition at line 16 of file JCalibrateToT.hh.

◆ INVALID_INDEX

const int JCALIBRATE::INVALID_INDEX = -1

static

invalid index

Definition at line 60 of file JFitK40.hh.

◆ TEROSTAT_DZ

double JCALIBRATE::TEROSTAT_DZ = 0.40

static

maximal PMT inclination

Definition at line 876 of file JFitK40.hh.

◆ TEROSTAT_R1

double JCALIBRATE::TEROSTAT_R1 = 1.00

static

scaling factor

Definition at line 877 of file JFitK40.hh.

◆ BELL_SHAPE

double JCALIBRATE::BELL_SHAPE = 1.55

static

Bell shape.

Definition at line 878 of file JFitK40.hh.

◆ FITTOT_SUFFIX

const std::string JCALIBRATE::FITTOT_SUFFIX = ".1ToT"

static

Definition at line 37 of file JFitToT.hh.

◆ FITTOT_FNAME

const std::string JCALIBRATE::FITTOT_FNAME = "fittot"

static

Definition at line 38 of file JFitToT.hh.

◆ FITTOT_GAIN_PARNAME

const char* const JCALIBRATE::FITTOT_GAIN_PARNAME = "gain"

static

Definition at line 40 of file JFitToT.hh.

◆ FITTOT_GAINSPREAD_PARNAME

const char* const JCALIBRATE::FITTOT_GAINSPREAD_PARNAME = "gainSpread"

static

Definition at line 41 of file JFitToT.hh.

◆ FITTOT_NORMALIZATION_PARNAME

const char* const JCALIBRATE::FITTOT_NORMALIZATION_PARNAME = "normalization"

static

Definition at line 42 of file JFitToT.hh.

◆ FITTOT_GAIN_MIN

const double JCALIBRATE::FITTOT_GAIN_MIN = 0.25

static

Default minimal gain.

Definition at line 44 of file JFitToT.hh.

◆ FITTOT_GAIN_MAX

const double JCALIBRATE::FITTOT_GAIN_MAX = 3.00

static

Default maximal gain.

Definition at line 45 of file JFitToT.hh.

◆ FITTOT_GAINSPREAD_MIN

const double JCALIBRATE::FITTOT_GAINSPREAD_MIN = 0.05

static

Default minimal gain spread.

Definition at line 47 of file JFitToT.hh.

◆ FITTOT_GAINSPREAD_MAX

const double JCALIBRATE::FITTOT_GAINSPREAD_MAX = 2.00

static

Default maximal gain spread.

Definition at line 48 of file JFitToT.hh.

◆ HVINTERPOLATOR_SUFFIX

const std::string JCALIBRATE::HVINTERPOLATOR_SUFFIX = ".HVxG"

static

Definition at line 29 of file JHVInterpolator.hh.

◆ HVINTERPOLATOR_DATA

const char* const JCALIBRATE::HVINTERPOLATOR_DATA = "HVINTERPOLATOR_DATA"

static

Definition at line 31 of file JHVInterpolator.hh.

◆ HVINTERPOLATOR_RESULT

const char* const JCALIBRATE::HVINTERPOLATOR_RESULT = "HVINTERPOLATOR_RESULT"

static

Definition at line 32 of file JHVInterpolator.hh.

Classes

Typedefs

Enumerations

Functions

Variables

Detailed Description

Typedef Documentation

◆ pair_type

Enumeration Type Documentation

◆ JOption_t

Function Documentation

◆ getP()

Variable Documentation

◆ weights_hist_t

◆ W1_overall_t

◆ WS_t

◆ WB_t

◆ _2S

◆ _1B

◆ _1L

◆ _2R

◆ _2F

◆ _1D

◆ _1F

◆ module_t

◆ string_t

◆ _2SToT

◆ INVALID_INDEX

◆ TEROSTAT_DZ

◆ TEROSTAT_R1

◆ BELL_SHAPE

◆ FITTOT_SUFFIX

◆ FITTOT_FNAME

◆ FITTOT_GAIN_PARNAME

◆ FITTOT_GAINSPREAD_PARNAME

◆ FITTOT_NORMALIZATION_PARNAME

◆ FITTOT_GAIN_MIN

◆ FITTOT_GAIN_MAX

◆ FITTOT_GAINSPREAD_MIN

◆ FITTOT_GAINSPREAD_MAX

◆ HVINTERPOLATOR_SUFFIX

◆ HVINTERPOLATOR_DATA

◆ HVINTERPOLATOR_RESULT