#include <SpectrumAnalysis.h>

Inheritance diagram for SpectrumAnalysis:

Collaboration diagram for SpectrumAnalysis:

Public Types
using	EffPeakSet = std::map< double, std::pair< double, double > >
	Container type for efficiency outputs. More...

Public Member Functions
	SpectrumAnalysis ()
	default constructor More...

void	loadHists (string a_fname, int a_modID, double a_countTime=0)

void	loadBackgroundHists (string a_fname, int a_modID, double a_countTime)
	loads GEANT4 output processed with GenAnalysis.generateSpectrumVeto() More...

void	addHists (string a_fname, int a_modID, double a_countTime=0)
	adds spectra to previously loaded spectra More...

void	addBackgroundHists (string a_fname, int a_modID, double a_countTime)

void	subtractBackgroundSpectra ()
	subtract background spectra from loaded spectra More...

void	setResolutionParameters (vector< double > a_params)

void	fixResolution (bool a_fixRes)

void	fixMean (bool a_fixMean)

void	setVerbosity (int a_verbosity)

void	rebinSpectrum (DetType a_detType, int a_detID, int a_elemID, int a_rebinFactor=2)
	rebins spectra More...

vector< PeakInfo >	FitPeaksMan (DetType a_detType, int a_detID, int a_elemID, double a_lower, double a_higher, vector< double > a_peakMeans, bool a_addBack=false, bool a_bgoVeto=false, bool a_checkPeakRes=true)

vector< PeakInfo >	FitPeaksMan (TH1F *a_hist, double a_lower, double a_higher, vector< double > a_peakMeans, int a_polyDegree=0)

vector< PeakInfo >	FitPeaksManPol1 (DetType a_detType, int a_detID, int a_elemID, double a_lower, double a_higher, vector< double > a_peakMeans, bool a_addBack=false, bool a_bgoVeto=false, bool a_checkPeakRes=true)
	same as FitPeaksMan but fit peaks with 1st degree polynomial background More...

void	printPeakFits ()
	prints information about fitted peaks (mean, area, resolution, chi2/ndf) More...

vector< PeakInfo >	FitPeaksErf (DetType a_detType, int a_detID, int a_elemID, double a_lower, double a_higher, double a_peakMean, bool a_addBack=false, bool a_bgoVeto=false)

vector< PeakInfo >	FitPeaksErf (TH1F *a_hist, double a_lower, double a_higher, double a_peakMean)

vector< PeakInfo >	FitPeaks6111Double (TH1F *a_hist)

EffPeakSet	Eu152Efficiency (DetType a_detType, int a_detID, int a_elemID, double a_activity, double a_countTime, double a_simBrTot=.2101, bool a_res=false, bool a_bgoVeto=false, bool a_addBack=false, bool a_draw=false)

EffPeakSet	Eu152Efficiency (TH1F *a_hist, double a_activity, double a_countTime, double a_simBrTot=.2101, bool a_res=false, bool a_draw=false)

EffPeakSet	Co56Efficiency (TH1F *a_hist, double a_activity, double a_countTime, double a_simBrTot=.9993999, bool a_res=false, bool a_draw=false)

EffPeakSet	Ga66Efficiency (TH1F *a_hist, double a_activity, double a_countTime, int multi_peaking=0, double a_simBrTot=.37, bool a_res=false, bool a_draw=false)

EffPeakSet	Cl35Efficiency (TH1F *a_hist, double a_activity, double a_countTime, bool double_peak=false, double a_simBrTot=1.0, bool a_res=false, bool a_draw=false)

EffPeakSet	createEffPeakSet (string infilename)

vector< double >	calcEffCurve (EffPeakSet const &a_ps, string a_fitType) const

void	calcAddBackFactor ()

map< int, TH1 * >	getHists ()
	returns m_hists, useful for drawing More...

map< int, vector< TH1 * > >	getCloverHists ()

TH1 *	getHist (DetType a_detType, int a_detID, int a_elemID=0, bool a_addBack=false, bool a_bgoVeto=false)

map< int, map< double, PeakInfo > >	getPeakFits ()

vector< double >	ChannelResolution (DetType a_detType, int a_detID, int a_elemID=0)
	creates plot of resolution vs energy More...

int	autoEu152EnergyCalibration (int a_channel, int a_guess=5000)

int	autoEu152EnergyCalibration (TH1 *a_hist, int a_guess=5000)

int	autoEu152EnergyCalibration1408Guess (int a_channel, int a_guess)

int	autoEu152EnergyCalibration1408Guess (TH1 *a_hist, int a_guess)

void	removePeakFit (DetType a_detType, int a_detID, int a_elemID, double a_energy)
	removes last peak in m_peakFits More...

void	setUpMinimizer ()

double	fitPeakConstrained (const double *a_params)

std::vector< double >	fitSTOFBandProjection (TH1F *a_histogram, double a_minLight, double a_maxLight, vector< double > a_peakCentroid, int a_polOrder)

	SpectrumAnalysis ()
	default constructor More...

void	loadHists (string a_fname, int a_modID, double a_countTime=0)

void	loadSimHists (string a_fname)
	loads GEANT4 output processed with GenAnalysis.generateSpectrumVeto() More...

void	loadBackgroundHists (string a_fname, int a_modID, double a_countTime)
	loads a spectrum for background subtraction More...

void	addHists (string a_fname, int a_modID, double a_countTime=0)
	adds spectra to previously loaded spectra More...

void	addBackgroundHists (string a_fname, int a_modID, double a_countTime)

void	subtractBackgroundSpectra ()
	subtract background spectra from loaded spectra More...

void	rebinSpectrum (int a_channel, int a_rebinFactor=2, bool a_fullClover=false)
	rebins spectra More...

vector< Peak >	FitPeaksMan (int a_channel, double a_lower, double a_higher, vector< double > a_peakMeans, bool a_fullClover=false, bool a_addBack=false, bool a_checkPeakRes=true)

vector< Peak >	FitPeaksManPol1 (int a_channel, int a_lower, int a_higher, vector< double > a_peakMeans, bool a_checkPeakRes=true, bool a_fullClover=false, bool a_addBack=false)
	same as aboe but fit peaks with 1st degree polynomial background More...

void	printPeakFits ()
	prints information about fitted peaks (mean, area, resolution, chi2/ndf) More...

vector< Peak >	FitPeaksErf (int a_channel, double a_lower, double a_higher, int a_numPeaks=1, bool a_fullClover=false, bool a_addBack=false)

map< double, std::pair< double, double > >	Eu152Efficiency (int a_channel, double a_activity, double a_countTime, double a_simBrTot=.2101, bool a_res=false, bool a_fullClover=false, bool a_addBack=false, bool a_draw=false)

void	calcAddBackFactor ()

map< double, std::pair< double, double > >	U235Efficiency (int a_channel, double a_activity, double a_activityErr, double a_countTime)

map< int, TH1 * >	getHists ()
	returns m_hists, useful for drawing More...

map< int, vector< TH1 * > >	getCloverHists ()

map< int, vector< Peak > >	getPeakFits ()

void	ChannelResolution (int a_channel)
	creates plot of resolution vs energy More...

int	autoEu152EnergyCalibration (int a_channel)

void	removeLastPeak (int a_channel)
	removes last peak in m_peakFits More...

double	fitPeakConstrained (const double *a_params)

Static Public Member Functions
static Double_t	totalFit (Double_t x, Double_t par)
	https://doi.org/10.1016/S0168-9002(98)00778-5 More...

static Double_t	fullGeFit (Double_t x, Double_t par)

Public Attributes
map< DetType, map< int, map< int, TH1 * > > >	m_histsD

Member Typedef Documentation

◆ EffPeakSet

using SpectrumAnalysis::EffPeakSet = std::map<double, std::pair<double, double> >

Container type for efficiency outputs.

Constructor & Destructor Documentation

◆ SpectrumAnalysis() [1/2]

SpectrumAnalysis::SpectrumAnalysis ( )

default constructor

◆ SpectrumAnalysis() [2/2]

SpectrumAnalysis::SpectrumAnalysis ( )

default constructor

Member Function Documentation

◆ addBackgroundHists() [1/2]

void SpectrumAnalysis::addBackgroundHists	(	string	a_fname,
		int	a_modID,
		double	a_countTime
	)

Here is the caller graph for this function:

◆ addBackgroundHists() [2/2]

void SpectrumAnalysis::addBackgroundHists	(	string	a_fname,
		int	a_modID,
		double	a_countTime
	)

◆ addHists() [1/2]

void SpectrumAnalysis::addHists	(	string	a_fname,
		int	a_modID,
		double	a_countTime = `0`
	)

adds spectra to previously loaded spectra

Here is the caller graph for this function:

◆ addHists() [2/2]

void SpectrumAnalysis::addHists	(	string	a_fname,
		int	a_modID,
		double	a_countTime = `0`
	)

adds spectra to previously loaded spectra

◆ autoEu152EnergyCalibration() [1/3]

int SpectrumAnalysis::autoEu152EnergyCalibration ( int a_channel )

attempts to generate an energy calibration from Eu152 data when rootFileBasicHistogrammer file is loaded assumes biggest peak is either 121.8, 344.3, or 1408 keV very much a W.I.P. use with caution

Here is the call graph for this function:

◆ autoEu152EnergyCalibration() [2/3]

int SpectrumAnalysis::autoEu152EnergyCalibration	(	int	a_channel,
		int	a_guess = `5000`
	)

attempts to generate an energy calibration from Eu152 data when rootFileBasicHistogrammer file is loaded assumes biggest peak is either 121.8, 344.3, or 1408 keV very much a W.I.P. use with caution

◆ autoEu152EnergyCalibration() [3/3]

int SpectrumAnalysis::autoEu152EnergyCalibration	(	TH1 *	a_hist,
		int	a_guess = `5000`
	)

Here is the call graph for this function:

◆ autoEu152EnergyCalibration1408Guess() [1/2]

int SpectrumAnalysis::autoEu152EnergyCalibration1408Guess	(	int	a_channel,
		int	a_guess
	)

◆ autoEu152EnergyCalibration1408Guess() [2/2]

int SpectrumAnalysis::autoEu152EnergyCalibration1408Guess	(	TH1 *	a_hist,
		int	a_guess
	)

Here is the call graph for this function:

◆ calcAddBackFactor() [1/2]

void SpectrumAnalysis::calcAddBackFactor ( )

produces a plot of Add-back factor vs. gamma energy for Eu-152 spectra currently hard-coded to do Clover 0 and Clover 1

since the add back factor is a ratio, don't need accurate source activity, count time, etc. clover 0

Here is the call graph for this function:

◆ calcAddBackFactor() [2/2]

void SpectrumAnalysis::calcAddBackFactor ( )

produces a plot of Add-back factor vs. gamma energy for Eu-152 spectra currently hard-coded to do Clover 0 and Clover 1

◆ calcEffCurve()

vector< double > SpectrumAnalysis::calcEffCurve	(	EffPeakSet const &	a_ps,
		string	a_fitType
	)		const

Calculates efficiency curve on a given peak set from SpectrumAnalysis::Eu152Efficiency returns efficiency parameteres fitType: Which fit function to use, debertin or lin (see Kis paper cited in the totalFit and fullGeFit functions)

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ChannelResolution() [1/2]

vector< double > SpectrumAnalysis::ChannelResolution	(	DetType	a_detType,
		int	a_detID,
		int	a_elemID = `0`
	)

creates plot of resolution vs energy

◆ ChannelResolution() [2/2]

void SpectrumAnalysis::ChannelResolution ( int a_channel )

creates plot of resolution vs energy

◆ Cl35Efficiency()

SpectrumAnalysis::EffPeakSet SpectrumAnalysis::Cl35Efficiency	(	TH1F *	a_hist,
		double	a_activity,
		double	a_countTime,
		bool	double_peak = `false`,
		double	a_simBrTot = `1.0`,
		bool	a_res = `false`,
		bool	a_draw = `false`
	)

calculates efficiency from an Cl-35(n,gamma) calibrated spectrum produces a graph of efficiency activity in mciroCurie, countTime in seconds Intensities (in the form of cross section, aka absolute intensity per 100 neutron captures) from A=36 Nuclear Data Sheets returns map of energy: area, area uncertainty If using this for relative efficiency to scale to an absolute source, just use an arbitrary value for the activity

gamma energies, branching ratios, BR uncertainties

Here is the call graph for this function:

◆ Co56Efficiency()

SpectrumAnalysis::EffPeakSet SpectrumAnalysis::Co56Efficiency	(	TH1F *	a_hist,
		double	a_activity,
		double	a_countTime,
		double	a_simBrTot = `.9993999`,
		bool	a_res = `false`,
		bool	a_draw = `false`
	)

calculates efficiency from an Co-56 calibrated spectrum produces a graph of efficiency activity in mciroCurie, countTime in seconds branching ratios (intensities) from A=56 Nuclear Data Sheets returns map of energy: area, area uncertainty doesn't fit all Co-56 peaks (but could) If using this for relative efficiency to scale to an absolute source, just use an arbitrary value for the activity

gamma energies, branching ratios, BR uncertainties

Here is the call graph for this function:

◆ createEffPeakSet()

SpectrumAnalysis::EffPeakSet SpectrumAnalysis::createEffPeakSet ( string infilename )

Here is the call graph for this function:

◆ Eu152Efficiency() [1/3]

SpectrumAnalysis::EffPeakSet SpectrumAnalysis::Eu152Efficiency	(	DetType	a_detType,
		int	a_detID,
		int	a_elemID,
		double	a_activity,
		double	a_countTime,
		double	a_simBrTot = `.2101`,
		bool	a_res = `false`,
		bool	a_bgoVeto = `false`,
		bool	a_addBack = `false`,
		bool	a_draw = `false`
	)

calculates efficiency from an Eu-152 calibrated spectrum produces a graph of efficiency activity in mciroCurie, countTime in seconds branching ratios from Knoll 3rd edition page 447 returns map of energy: area, area uncertainty doesn't fit all Eu-152 peaks (but could)

Here is the caller graph for this function:

◆ Eu152Efficiency() [2/3]

map< double, std::pair< double, double > > SpectrumAnalysis::Eu152Efficiency	(	int	a_channel,
		double	a_activity,
		double	a_countTime,
		double	a_simBrTot = `.2101`,
		bool	a_res = `false`,
		bool	a_fullClover = `false`,
		bool	a_addBack = `false`,
		bool	a_draw = `false`
	)

calculates efficiency from an Eu-152 calibrated spectrum produces a graph of efficiency activity in mciroCurie, countTime in seconds branching ratios from Knoll 3rd edition page 447 returns map of energy: area, area uncertainty doesn't fit all Eu-152 peaks (but could)

gamma energies, branching ratios, BR uncertainties

Here is the call graph for this function:

◆ Eu152Efficiency() [3/3]

SpectrumAnalysis::EffPeakSet SpectrumAnalysis::Eu152Efficiency	(	TH1F *	a_hist,
		double	a_activity,
		double	a_countTime,
		double	a_simBrTot = `.2101`,
		bool	a_res = `false`,
		bool	a_draw = `false`
	)

gamma energies, branching ratios, BR uncertainties

Here is the call graph for this function:

◆ fitPeakConstrained() [1/2]

double SpectrumAnalysis::fitPeakConstrained ( const double * a_params )

set-up the fit function

calculate chi2 between fit and histogram

penalize for distance to expected width(s)

Here is the caller graph for this function:

◆ fitPeakConstrained() [2/2]

double SpectrumAnalysis::fitPeakConstrained ( const double * a_params )

◆ FitPeaks6111Double()

vector< PeakInfo > SpectrumAnalysis::FitPeaks6111Double ( TH1F * a_hist )

estimate amplitude and background

containers for fit results

set up the minimizer with initial conditions

do the minimization

get the results – need to check if minimizer actually succeeds?

save fit results

Here is the call graph for this function:

Here is the caller graph for this function:

◆ FitPeaksErf() [1/3]

vector< PeakInfo > SpectrumAnalysis::FitPeaksErf	(	DetType	a_detType,
		int	a_detID,
		int	a_elemID,
		double	a_lower,
		double	a_higher,
		double	a_peakMean,
		bool	a_addBack = `false`,
		bool	a_bgoVeto = `false`
	)

fits peaks using the full RadWare fit function currently only works for single peaks (a_numPeaks is meaningless)

Here is the caller graph for this function:

◆ FitPeaksErf() [2/3]

vector< Peak > SpectrumAnalysis::FitPeaksErf	(	int	a_channel,
		double	a_lower,
		double	a_higher,
		int	a_numPeaks = `1`,
		bool	a_fullClover = `false`,
		bool	a_addBack = `false`
	)

fits peaks using the full RadWare fit function currently only works for single peaks (a_numPeaks is meaningless)

Here is the call graph for this function:

◆ FitPeaksErf() [3/3]

vector< PeakInfo > SpectrumAnalysis::FitPeaksErf	(	TH1F *	a_hist,
		double	a_lower,
		double	a_higher,
		double	a_peakMean
	)

Here is the call graph for this function:

◆ FitPeaksMan() [1/3]

vector< PeakInfo > SpectrumAnalysis::FitPeaksMan	(	DetType	a_detType,
		int	a_detID,
		int	a_elemID,
		double	a_lower,
		double	a_higher,
		vector< double >	a_peakMeans,
		bool	a_addBack = `false`,
		bool	a_bgoVeto = `false`,
		bool	a_checkPeakRes = `true`
	)

takes in channel, upper and lower bounds over which to fit a peak fits with linear (0 degree polynomial) background stores fit in m_peakFits a_peakMeans -> used to seed peak means for fitting doublets put larger peak first in list a_fullClover -> will fit clover spectra (a_channel then is Clover ID) a_addBack -> will fit either added-back or non clover spectra a_checkPeakRes -> only stores fit is resolution is reasonable passes relevant hists to FitPeaksMan and saves fits to m_peakFits

Here is the caller graph for this function:

◆ FitPeaksMan() [2/3]

vector< Peak > SpectrumAnalysis::FitPeaksMan	(	int	a_channel,
		double	a_lower,
		double	a_higher,
		vector< double >	a_peakMeans,
		bool	a_fullClover = `false`,
		bool	a_addBack = `false`,
		bool	a_checkPeakRes = `true`
	)

takes in channel, upper and lower bounds over which to fit a peak fits with linear (0 degree polynomial) background stores fit in m_peakFits a_peakMeans -> used to seed peak means for fitting doublets a_fullClover -> will fit clover spectra (a_channel then is Clover ID) a_addBack -> will fit either added-back or non clover spectra a_checkPeakRes -> only stores fit is resolution is reasonable

get parameters for seeding the fit

seed fit

save the fitted peak info

r1->Ndf();

Here is the call graph for this function:

◆ FitPeaksMan() [3/3]

vector< PeakInfo > SpectrumAnalysis::FitPeaksMan	(	TH1F *	a_hist,
		double	a_lower,
		double	a_higher,
		vector< double >	a_peakMeans,
		int	a_polyDegree = `0`
	)

performs fit on a_hist a_polyDegree is degree of polynomial for background (=0, 1) can fit singlets, doublets, triplets

estimate amplitude and background

containers for fit results

expected width

gauss + constant background

seed fit with educated guesses

number of parameters to fit

location of sigma in parameter list (for penalty)

set up the minimizer with initial conditions

do the minimization

get the results – need to check if minimizer actually succeeds?

save fit results

Here is the call graph for this function:

◆ FitPeaksManPol1() [1/2]

vector< PeakInfo > SpectrumAnalysis::FitPeaksManPol1	(	DetType	a_detType,
		int	a_detID,
		int	a_elemID,
		double	a_lower,
		double	a_higher,
		vector< double >	a_peakMeans,
		bool	a_addBack = `false`,
		bool	a_bgoVeto = `false`,
		bool	a_checkPeakRes = `true`
	)

same as FitPeaksMan but fit peaks with 1st degree polynomial background

Here is the call graph for this function:

Here is the caller graph for this function:

◆ FitPeaksManPol1() [2/2]

vector< Peak > SpectrumAnalysis::FitPeaksManPol1	(	int	a_channel,
		int	a_lower,
		int	a_higher,
		vector< double >	a_peakMeans,
		bool	a_checkPeakRes = `true`,
		bool	a_fullClover = `false`,
		bool	a_addBack = `false`
	)

same as aboe but fit peaks with 1st degree polynomial background

Here is the call graph for this function:

◆ fitSTOFBandProjection()

std::vector< double > SpectrumAnalysis::fitSTOFBandProjection	(	TH1F *	a_histogram,
		double	a_minLight,
		double	a_maxLight,
		vector< double >	a_peakCentroid,
		int	a_polOrder
	)

Here is the call graph for this function:

Here is the caller graph for this function:

◆ fixMean()

void SpectrumAnalysis::fixMean ( bool a_fixMean )

Here is the caller graph for this function:

◆ fixResolution()

void SpectrumAnalysis::fixResolution ( bool a_fixRes )

fixes the resolution based on resolution params for peak fits sets m_fixRes, default is false

◆ fullGeFit()

Double_t SpectrumAnalysis::fullGeFit	(	Double_t *	x,
		Double_t *	par
	)

static

Fitting function for efficiency curve Lin fit function [see z. kis NIMA (1998)] Good for 300 - 9700 keV https://doi.org/10.1016/S0168-9002(98)00778-5

Here is the caller graph for this function:

◆ Ga66Efficiency()

SpectrumAnalysis::EffPeakSet SpectrumAnalysis::Ga66Efficiency	(	TH1F *	a_hist,
		double	a_activity,
		double	a_countTime,
		int	multi_peaking = `0`,
		double	a_simBrTot = `.37`,
		bool	a_res = `false`,
		bool	a_draw = `false`
	)

calculates efficiency from an Ga-66 (from 66Zn(d,2n)) calibrated spectrum produces a graph of efficiency activity in mciroCurie, countTime in seconds branching ratios (intensities) from A=66 Nuclear Data Sheets returns map of energy: area, area uncertainty If using this for relative efficiency to scale to an absolute source, just use an arbitrary value for the activity

gamma energies, branching ratios, BR uncertainties

Here is the call graph for this function:

◆ getCloverHists() [1/2]

map< int, vector< TH1 * > > SpectrumAnalysis::getCloverHists ( )

◆ getCloverHists() [2/2]

map<int,vector<TH1*> > SpectrumAnalysis::getCloverHists ( )

◆ getHist()

TH1 * SpectrumAnalysis::getHist	(	DetType	a_detType,
		int	a_detID,
		int	a_elemID = `0`,
		bool	a_addBack = `false`,
		bool	a_bgoVeto = `false`
	)

Here is the caller graph for this function:

◆ getHists() [1/2]

map< int, TH1 * > SpectrumAnalysis::getHists ( )

returns m_hists, useful for drawing

low energy efficiency from U235 decays activity in Bq

◆ getHists() [2/2]

map<int, TH1*> SpectrumAnalysis::getHists ( )

returns m_hists, useful for drawing

◆ getPeakFits() [1/2]

map< int, vector< Peak > > SpectrumAnalysis::getPeakFits ( )

◆ getPeakFits() [2/2]

map<int, vector<Peak> > SpectrumAnalysis::getPeakFits ( )

◆ loadBackgroundHists() [1/2]

void SpectrumAnalysis::loadBackgroundHists	(	string	a_fname,
		int	a_modID,
		double	a_countTime
	)

loads GEANT4 output processed with GenAnalysis.generateSpectrumVeto()

loads a spectrum for background subtraction

Here is the caller graph for this function:

◆ loadBackgroundHists() [2/2]

void SpectrumAnalysis::loadBackgroundHists	(	string	a_fname,
		int	a_modID,
		double	a_countTime
	)

loads a spectrum for background subtraction

◆ loadHists() [1/2]

void SpectrumAnalysis::loadHists	(	string	a_fname,
		int	a_modID,
		double	a_countTime = `0`
	)

loads histograms (with >10 counts) from file containing histograms currently only loads histograms generated from rootFileBasicHistogrammer or "calibrated hists" (GenesisAnalysis.buildGammaSpectra()) stores hists in m_hists --> this map is cleared at the top of the function m_hists is indexed by channel - hist name and channel are printed if the file if CalibratedSpectra with full clover / add-back/BGO etc. those hists are stored in m_cloverHists

Here is the caller graph for this function:

◆ loadHists() [2/2]

void SpectrumAnalysis::loadHists	(	string	a_fname,
		int	a_modID,
		double	a_countTime = `0`
	)

loads histograms (with >10 counts) from file containing histograms currently only loads histograms generated from "basic histogrammer" or "calibrated hists" (GenCoinAnalysis.buildGammaSpectra())

◆ loadSimHists()

void SpectrumAnalysis::loadSimHists ( string a_fname )

loads GEANT4 output processed with GenAnalysis.generateSpectrumVeto()

◆ printPeakFits() [1/2]

void SpectrumAnalysis::printPeakFits ( )

prints information about fitted peaks (mean, area, resolution, chi2/ndf)

◆ printPeakFits() [2/2]

void SpectrumAnalysis::printPeakFits ( )

prints information about fitted peaks (mean, area, resolution, chi2/ndf)

◆ rebinSpectrum() [1/2]

void SpectrumAnalysis::rebinSpectrum	(	DetType	a_detType,
		int	a_detID,
		int	a_elemID,
		int	a_rebinFactor = `2`
	)

rebins spectra

◆ rebinSpectrum() [2/2]

void SpectrumAnalysis::rebinSpectrum	(	int	a_channel,
		int	a_rebinFactor = `2`,
		bool	a_fullClover = `false`
	)

rebins spectra

◆ removeLastPeak()

void SpectrumAnalysis::removeLastPeak ( int a_channel )

removes last peak in m_peakFits

◆ removePeakFit()

void SpectrumAnalysis::removePeakFit	(	DetType	a_detType,
		int	a_detID,
		int	a_elemID,
		double	a_energy
	)

removes last peak in m_peakFits

◆ setResolutionParameters()

void SpectrumAnalysis::setResolutionParameters ( vector< double > a_params )

set the expected peak resolution for constraining peak fits parameterizing resolution as FWHM = sqrt([0]*E_gamma + [1])

Here is the caller graph for this function:

◆ setUpMinimizer()

void SpectrumAnalysis::setUpMinimizer ( )

need to add check for Minuit2

Here is the call graph for this function:

Here is the caller graph for this function:

◆ setVerbosity()

void SpectrumAnalysis::setVerbosity ( int a_verbosity )

Here is the caller graph for this function:

◆ subtractBackgroundSpectra() [1/2]

void SpectrumAnalysis::subtractBackgroundSpectra ( )

subtract background spectra from loaded spectra

Here is the call graph for this function:

Here is the caller graph for this function:

◆ subtractBackgroundSpectra() [2/2]

void SpectrumAnalysis::subtractBackgroundSpectra ( )

subtract background spectra from loaded spectra

◆ totalFit()

Double_t SpectrumAnalysis::totalFit	(	Double_t *	x,
		Double_t *	par
	)

static

https://doi.org/10.1016/S0168-9002(98)00778-5

Fitting function for efficiency curve debertin fit function [see z. kis NIMA (1998)] Good for 120-1400 keV

Here is the caller graph for this function:

◆ U235Efficiency()

map< double, std::pair< double, double > > SpectrumAnalysis::U235Efficiency	(	int	a_channel,
		double	a_activity,
		double	a_activityErr,
		double	a_countTime
	)

low energy efficiency from U235 decays activity in Bq

Here is the call graph for this function:

Member Data Documentation

◆ m_histsD

map<DetType, map<int, map<int,TH1*> > > SpectrumAnalysis::m_histsD

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

analysis/include/SpectrumAnalysis.h
analysis/src/SpectrumAnalysis.cpp

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Member Typedef Documentation

◆ EffPeakSet

Constructor & Destructor Documentation

◆ SpectrumAnalysis() [1/2]

◆ SpectrumAnalysis() [2/2]

Member Function Documentation

◆ addBackgroundHists() [1/2]

◆ addBackgroundHists() [2/2]

◆ addHists() [1/2]

◆ addHists() [2/2]

◆ autoEu152EnergyCalibration() [1/3]

◆ autoEu152EnergyCalibration() [2/3]

◆ autoEu152EnergyCalibration() [3/3]

◆ autoEu152EnergyCalibration1408Guess() [1/2]

◆ autoEu152EnergyCalibration1408Guess() [2/2]

◆ calcAddBackFactor() [1/2]

◆ calcAddBackFactor() [2/2]

◆ calcEffCurve()

◆ ChannelResolution() [1/2]

◆ ChannelResolution() [2/2]

◆ Cl35Efficiency()

◆ Co56Efficiency()

◆ createEffPeakSet()

◆ Eu152Efficiency() [1/3]

◆ Eu152Efficiency() [2/3]

◆ Eu152Efficiency() [3/3]

◆ fitPeakConstrained() [1/2]

◆ fitPeakConstrained() [2/2]

◆ FitPeaks6111Double()

◆ FitPeaksErf() [1/3]

◆ FitPeaksErf() [2/3]

◆ FitPeaksErf() [3/3]

◆ FitPeaksMan() [1/3]

◆ FitPeaksMan() [2/3]

◆ FitPeaksMan() [3/3]

◆ FitPeaksManPol1() [1/2]

◆ FitPeaksManPol1() [2/2]

◆ fitSTOFBandProjection()

◆ fixMean()

◆ fixResolution()

◆ fullGeFit()

◆ Ga66Efficiency()

◆ getCloverHists() [1/2]

◆ getCloverHists() [2/2]

◆ getHist()

◆ getHists() [1/2]

◆ getHists() [2/2]

◆ getPeakFits() [1/2]

◆ getPeakFits() [2/2]

◆ loadBackgroundHists() [1/2]

◆ loadBackgroundHists() [2/2]

◆ loadHists() [1/2]

◆ loadHists() [2/2]

◆ loadSimHists()

◆ printPeakFits() [1/2]

◆ printPeakFits() [2/2]

◆ rebinSpectrum() [1/2]

◆ rebinSpectrum() [2/2]

◆ removeLastPeak()

◆ removePeakFit()

◆ setResolutionParameters()

◆ setUpMinimizer()

◆ setVerbosity()

◆ subtractBackgroundSpectra() [1/2]

◆ subtractBackgroundSpectra() [2/2]

◆ totalFit()

◆ U235Efficiency()

Member Data Documentation

◆ m_histsD