Parameter error calculation tool. More...

#include <cterror.hpp>

Inheritance diagram for cterror:

Public Member Functions
	cterror (void)
	Void constructor.

	cterror (const GObservations &obs)
	Observations constructor.

	cterror (int argc, char *argv[])
	Command line constructor.

	cterror (const cterror &app)
	Copy constructor.

virtual	~cterror (void)
	Destructor.

cterror &	operator= (const cterror &app)
	Assignment operator.

void	clear (void)
	Clear cterror tool.

void	process (void)
	Compute parameter errors using a likelihood profile method.

void	save (void)
	Save model.

Public Member Functions inherited from ctlikelihood
	ctlikelihood (const std::string &name, const std::string &version)
	Name constructor.

	ctlikelihood (const std::string &name, const std::string &version, const GApplicationPars &pars)
	Application parameters constructor.

	ctlikelihood (const std::string &name, const std::string &version, int argc, char *argv[])
	Command line constructor.

	ctlikelihood (const std::string &name, const std::string &version, const GObservations &obs)
	Observations constructor.

	ctlikelihood (const ctlikelihood &app)
	Copy constructor.

virtual	~ctlikelihood (void)
	Destructor.

ctlikelihood &	operator= (const ctlikelihood &app)
	Assignment operator.

void	opt (const GOptimizer *opt)
	Set optimizer.

const GOptimizer *	opt (void) const
	Return optimizer.

Public Member Functions inherited from ctobservation
	ctobservation (const std::string &name, const std::string &version)
	Name constructor.

	ctobservation (const std::string &name, const std::string &version, const GApplicationPars &pars)
	Application parameters constructor.

	ctobservation (const std::string &name, const std::string &version, int argc, char *argv[])
	Command line constructor.

	ctobservation (const std::string &name, const std::string &version, const GObservations &obs)
	Observations constructor.

	ctobservation (const ctobservation &app)
	Copy constructor.

virtual	~ctobservation (void)
	Destructor.

ctobservation &	operator= (const ctobservation &app)
	Assignment operator.

void	obs (const GObservations &obs)
	Set observation container.

const GObservations &	obs (void) const
	Return observation container.

Public Member Functions inherited from ctool
	ctool (const std::string &name, const std::string &version)
	Name constructor.

	ctool (const std::string &name, const std::string &version, const GApplicationPars &pars)
	Application parameter constructor.

	ctool (const std::string &name, const std::string &version, int argc, char *argv[])
	Command line constructor.

	ctool (const ctool &app)
	Copy constructor.

virtual	~ctool (void)
	Destructor.

ctool &	operator= (const ctool &app)
	Assignment operator.

virtual void	run (void)
	Run ctool.

virtual void	execute (void)
	Execute ctool.

Protected Member Functions
void	init_members (void)
	Initialise class members.

void	copy_members (const cterror &app)
	Copy class members.

void	free_members (void)
	Delete class members.

void	get_parameters (void)
	Get application parameters.

double	error_bisection (const double &min, const double &max)
	Calculate error using a bisection method.

void	update_intensities (void)
	Update intensities.

void	create_fits (void)
	Set result FITS file.

Protected Member Functions inherited from ctlikelihood
double	evaluate (GModelPar &par, const double &value)
	Evaluates the log-likelihood function.

void	init_members (void)
	Initialise class members.

void	copy_members (const ctlikelihood &app)
	Copy class members.

void	free_members (void)
	Delete class members.

Protected Member Functions inherited from ctobservation
GCTAObservation *	first_unbinned_observation (void)
	Return first unbinned CTA observation.

GCTAObservation *	next_unbinned_observation (void)
	Return next unbinned CTA observation.

const GCTAObservation *	first_unbinned_observation (void) const
	Return first unbinned CTA observation (const version)

const GCTAObservation *	next_unbinned_observation (void) const
	Return next unbinned CTA observation (const version)

void	read_ogip_keywords (GFitsHDU *hdu) const
	Read OGIP keywords from FITS HDU.

void	write_ogip_keywords (GFitsHDU *hdu) const
	Write OGIP keywords in FITS HDU.

void	set_obs_statistic (const std::string &statistic)
	Set fit statistic for CTA observations.

void	set_obs_bounds ()
	Set observation boundaries for CTA observations.

void	save_events_fits (void)
	Save event list in FITS format.

void	save_events_xml (void)
	Save event list(s) in XML format.

void	init_members (void)
	Initialise class members.

void	copy_members (const ctobservation &app)
	Copy class members.

void	free_members (void)
	Delete class members.

Protected Member Functions inherited from ctool
void	init_members (void)
	Initialise class members.

void	copy_members (const ctool &app)
	Copy class members.

void	free_members (void)
	Delete class members.

void	sync_pfiles (void)
	Synchronise parameter files.

const bool &	read_ahead (void) const
	Signal whether parameters should be read ahead.

void	setup_observations (GObservations &obs, const bool &response=true, const bool &list=true, const bool &cube=true)
	Setup observation container.

void	setup_models (GObservations &obs, const std::string &name="")
	Setup model container.

GEbounds	create_ebounds (void)
	Create energy boundaries from user parameters.

GEnergies	create_energies (void)
	Create energies from user parameters.

GSkyMap	create_map (const GObservations &obs)
	Create a skymap from user parameters.

GCTAEventCube	create_cube (const GObservations &obs)
	Create a CTA event cube from user parameters.

GCTAObservation	create_cta_obs (void)
	Create a CTA observation from User parameters.

void	require_inobs (const std::string &method)
	Throws exception if inobs parameter is not valid.

void	require_inobs_nolist (const std::string &method)

void	require_inobs_nocube (const std::string &method)
	Throws exception if inobs parameter is a counts cube.

GCTARoi	get_roi (const GCTAPointing &pnt=GCTAPointing())
	Return RoI from User parameters.

GEbounds	get_ebounds (void)
	Return energy boundaries from User parameters.

GGti	get_gti (const GTimeReference &ref)
	Return Good Time Intervals from User parameter.

GCTAPointing	get_pointing (void)
	Return CTA pointing from User parameters.

GSkyDir	get_skydir (void)
	Return sky direction from User parameters.

std::string	set_outfile_name (const std::string &filename)
	Set output file name.

bool	is_stacked (void)
	Query user parameters for stacked analysis.

bool	is_onoff (void)
	Query user parameters for On/Off analysis.

void	log_parameters (const GChatter &chatter)
	Log application parameters.

void	log_observations (const GChatter &chatter, const GObservations &obs, const std::string &what="Observation")
	Log observation container.

void	log_models (const GChatter &chatter, const GModels &models, const std::string &what="Model")
	Log model container.

void	set_response (GObservations &obs)
	Set response for all CTA observations in container.

std::vector< bool >	set_edisp (GObservations &obs, const bool &edisp) const
	Set energy dispersion to CTA observations.

void	restore_edisp (GObservations &obs, const std::vector< bool > &edisp) const
	Restore energy dispersion flags of CTA observations.

void	set_obs_response (GCTAObservation *obs)
	Set response for CTA observation.

GObservations	get_observations (const bool &get_response=true)
	Get observation container.

GSkyDir	get_mean_pointing (const GObservations &obs)
	Derives mean pointing from CTA observations.

size_t	get_current_rss (void)
	Get current resident set size (physical memory use) in Bytes.

std::string	get_obs_header (const GObservation *obs) const
	Return observation header string.

GEnergies	insert_energy_boundaries (const GEnergies &energies, const GCTAObservation &obs)
	Insert observation energy boundaries into list of energies.

std::vector< bool >	cube_layer_usage (const GEbounds &cube_ebounds, const GEbounds &list_ebounds) const
	Determine the counts cube layer usage.

std::string	get_gtiname (const std::string &filename, const std::string &evtname) const
	Get Good Time Intervals extension name.

void	save_event_list (const GCTAObservation *obs, const std::string &infile, const std::string &evtname, const std::string &gtiname, const std::string &outfile) const
	Save event list into FITS file.

std::string	warn_too_few_energies (const GEnergies &energies) const
	Set warning string if there are too few energies.

std::string	warn_xml_suffix (const GFilename &filename) const
	Set warning string if file has no .xml suffix.

void	provide_help (void) const
	Dump help text in the console.

Protected Attributes
std::string	m_srcname
	Name of source.

double	m_confidence
	Confidence level.

double	m_tol
	Tolerance for limit determination.

int	m_max_iter
	Maximum number of iterations.

bool	m_apply_edisp
	Apply energy dispersion?

GChatter	m_chatter
	Chattiness.

double	m_value
	Parameter value.

double	m_dlogL
	Likelihood difference for upper limit computation.

GModelPar *	m_model_par
	Pointer to model parameter.

double	m_best_logL
	Best fit log likelihood of given model.

GEnergies	m_energies
	Energy values.

std::vector< double >	m_intensities
	Model intensity.

std::vector< double >	m_min_intensities
	Minimum intensities.

std::vector< double >	m_max_intensities
	Maximum intensities.

GFits	m_fits
	FITS file holding butterfly.

Protected Attributes inherited from ctlikelihood
GOptimizerLM	m_opt
	Optimizer.

Protected Attributes inherited from ctobservation
GObservations	m_obs
	Observation container.

Protected Attributes inherited from ctool
bool	m_read_ahead
	Read ahead output parameters.

bool	m_use_xml
	Use XML file instead of FITS file for observations.

Detailed Description

Parameter error calculation tool.

Definition at line 42 of file cterror.hpp.

Constructor & Destructor Documentation

◆ cterror() [1/4]

cterror::cterror ( void )

Void constructor.

Constructs an empty cterror tool.

Definition at line 60 of file cterror.cpp.

References init_members().

◆ cterror() [2/4]

cterror::cterror ( const GObservations & obs )

explicit

Observations constructor.

param[in] obs Observation container.

Constructs cterror tool from an observations container.

Definition at line 77 of file cterror.cpp.

References init_members().

◆ cterror() [3/4]

cterror::cterror	(	int	argc,
		char *	argv[] )

Command line constructor.

Parameters

[in]	argc	Number of arguments in command line.
[in]	argv	Array of command line arguments.

Constructs an instance of the cterror tool that will parse user parameters that are provided as command line arguments.

Definition at line 97 of file cterror.cpp.

References init_members().

◆ cterror() [4/4]

cterror::cterror ( const cterror & app )

Copy constructor.

Parameters

[in] app Application.

Constructs an instance of the cterror tool by copying information from another ctulimit tool.

Definition at line 116 of file cterror.cpp.

References copy_members(), and init_members().

◆ ~cterror()

cterror::~cterror ( void )

virtual

Destructor.

Destructs the cterror tool.

Definition at line 134 of file cterror.cpp.

References free_members().

Member Function Documentation

◆ clear()

void cterror::clear ( void )

virtual

Clear cterror tool.

Clears cterror tool.

Implements ctlikelihood.

Definition at line 193 of file cterror.cpp.

References free_members(), ctlikelihood::free_members(), ctobservation::free_members(), ctool::free_members(), init_members(), ctlikelihood::init_members(), ctobservation::init_members(), and ctool::init_members().

◆ copy_members()

void cterror::copy_members ( const cterror & app )

protected

Copy class members.

Parameters

[in] app Application.

Definition at line 471 of file cterror.cpp.

References m_apply_edisp, m_best_logL, m_chatter, m_confidence, m_dlogL, m_energies, m_fits, m_intensities, m_max_intensities, m_max_iter, m_min_intensities, m_model_par, m_srcname, m_tol, and m_value.

Referenced by cterror(), and operator=().

◆ create_fits()

void cterror::create_fits ( void )

protected

Set result FITS file.

Set the FITS file with the butterfly results.

Definition at line 751 of file cterror.cpp.

References m_energies, m_fits, m_intensities, m_max_intensities, and m_min_intensities.

Referenced by process().

◆ error_bisection()

double cterror::error_bisection	(	const double &	min,
		const double &	max )

protected

Calculate error using a bisection method.

Parameters

[in]	min	Minimum parameter value
[in]	max	Maximum parameter value

Calculates the error using a bisection method.

Definition at line 578 of file cterror.cpp.

References ctlikelihood::evaluate(), G_ERR_BISECTION, m_best_logL, m_dlogL, m_max_iter, m_model_par, m_tol, and m_value.

Referenced by process().

◆ free_members()

void cterror::free_members ( void )

protected

Delete class members.

Definition at line 502 of file cterror.cpp.

Referenced by clear(), operator=(), and ~cterror().

◆ get_parameters()

void cterror::get_parameters ( void )

protected

Get application parameters.

Get all task parameters from parameter file.

Definition at line 514 of file cterror.cpp.

References ctool::create_energies(), ctool::log_parameters(), m_apply_edisp, m_chatter, m_confidence, m_dlogL, m_energies, m_max_iter, ctobservation::m_obs, ctlikelihood::m_opt, m_srcname, m_tol, ctool::read_ahead(), ctobservation::set_obs_statistic(), ctool::setup_models(), and ctool::setup_observations().

Referenced by process().

◆ init_members()

void cterror::init_members ( void )

protected

Initialise class members.

Definition at line 438 of file cterror.cpp.

References m_apply_edisp, m_best_logL, m_chatter, m_confidence, m_dlogL, m_energies, m_fits, m_intensities, m_max_intensities, m_max_iter, m_min_intensities, m_model_par, m_srcname, m_tol, and m_value.

Referenced by clear(), cterror(), cterror(), cterror(), cterror(), and operator=().

◆ operator=()

cterror & cterror::operator= ( const cterror & app )

Assignment operator.

Parameters

[in] app Application.

Returns: Application.

Assigns a cterror tool.

Definition at line 158 of file cterror.cpp.

References copy_members(), free_members(), init_members(), and ctlikelihood::operator=().

◆ process()

void cterror::process ( void )

virtual

Compute parameter errors using a likelihood profile method.

Computes the parameter errors using a likelihood profile method.

Implements ctlikelihood.

Definition at line 223 of file cterror.cpp.

References create_fits(), error_bisection(), get_parameters(), ctool::log_observations(), m_apply_edisp, m_best_logL, m_chatter, m_confidence, m_dlogL, m_energies, m_intensities, m_max_intensities, m_min_intensities, m_model_par, ctobservation::m_obs, ctlikelihood::m_opt, m_srcname, m_value, ctool::restore_edisp(), ctool::set_edisp(), and update_intensities().

◆ save()

void cterror::save ( void )

virtual

Save model.

Saves the model into an XML file.

Implements ctlikelihood.

Definition at line 381 of file cterror.cpp.

References m_fits, and ctobservation::m_obs.

◆ update_intensities()

void cterror::update_intensities ( void )

protected

Update intensities.

Update the intensities using the current model in the observation container. If the method is called with empty vectors, intensity vectors are allocated.

Definition at line 687 of file cterror.cpp.

References m_energies, m_intensities, m_max_intensities, m_min_intensities, ctobservation::m_obs, and m_srcname.

Referenced by process().

Member Data Documentation

◆ m_apply_edisp

bool cterror::m_apply_edisp

protected

Apply energy dispersion?

Definition at line 75 of file cterror.hpp.

Referenced by copy_members(), get_parameters(), init_members(), and process().

◆ m_best_logL

double cterror::m_best_logL

protected

Best fit log likelihood of given model.

Definition at line 82 of file cterror.hpp.

Referenced by copy_members(), error_bisection(), init_members(), and process().

◆ m_chatter

GChatter cterror::m_chatter

protected

Chattiness.

Definition at line 76 of file cterror.hpp.

Referenced by copy_members(), get_parameters(), init_members(), and process().

◆ m_confidence

double cterror::m_confidence

protected

Confidence level.

Definition at line 72 of file cterror.hpp.

Referenced by copy_members(), get_parameters(), init_members(), and process().

◆ m_dlogL

double cterror::m_dlogL

protected

Likelihood difference for upper limit computation.

Definition at line 80 of file cterror.hpp.

Referenced by copy_members(), error_bisection(), get_parameters(), init_members(), and process().

◆ m_energies

GEnergies cterror::m_energies

protected

Energy values.

Definition at line 85 of file cterror.hpp.

Referenced by copy_members(), create_fits(), get_parameters(), init_members(), process(), and update_intensities().

◆ m_fits

GFits cterror::m_fits

protected

FITS file holding butterfly.

Definition at line 89 of file cterror.hpp.

Referenced by copy_members(), create_fits(), init_members(), and save().

◆ m_intensities

std::vector<double> cterror::m_intensities

protected

Model intensity.

Definition at line 86 of file cterror.hpp.

Referenced by copy_members(), create_fits(), init_members(), process(), and update_intensities().

◆ m_max_intensities

std::vector<double> cterror::m_max_intensities

protected

Maximum intensities.

Definition at line 88 of file cterror.hpp.

Referenced by copy_members(), create_fits(), init_members(), process(), and update_intensities().

◆ m_max_iter

int cterror::m_max_iter

protected

Maximum number of iterations.

Definition at line 74 of file cterror.hpp.

Referenced by copy_members(), error_bisection(), get_parameters(), and init_members().

◆ m_min_intensities

std::vector<double> cterror::m_min_intensities

protected

Minimum intensities.

Definition at line 87 of file cterror.hpp.

Referenced by copy_members(), create_fits(), init_members(), process(), and update_intensities().

◆ m_model_par

GModelPar* cterror::m_model_par

protected

Pointer to model parameter.

Definition at line 81 of file cterror.hpp.

Referenced by copy_members(), error_bisection(), init_members(), and process().

◆ m_srcname

std::string cterror::m_srcname

protected

Name of source.

Definition at line 71 of file cterror.hpp.

Referenced by copy_members(), get_parameters(), init_members(), process(), and update_intensities().

◆ m_tol

double cterror::m_tol

protected

Tolerance for limit determination.

Definition at line 73 of file cterror.hpp.

Referenced by copy_members(), error_bisection(), get_parameters(), and init_members().

◆ m_value

double cterror::m_value

protected

Parameter value.

Definition at line 79 of file cterror.hpp.

Referenced by copy_members(), error_bisection(), init_members(), and process().

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

src/cterror/cterror.hpp
src/cterror/cterror.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ cterror() [1/4]

◆ cterror() [2/4]

◆ cterror() [3/4]

◆ cterror() [4/4]

◆ ~cterror()

Member Function Documentation

◆ clear()

◆ copy_members()

◆ create_fits()

◆ error_bisection()

◆ free_members()

◆ get_parameters()

◆ init_members()

◆ operator=()

◆ process()

◆ save()

◆ update_intensities()

Member Data Documentation

◆ m_apply_edisp

◆ m_best_logL

◆ m_chatter

◆ m_confidence

◆ m_dlogL

◆ m_energies

◆ m_fits

◆ m_intensities

◆ m_max_intensities

◆ m_max_iter

◆ m_min_intensities

◆ m_model_par

◆ m_srcname

◆ m_tol

◆ m_value