#include "GausFitCache.h"

Inheritance diagram for hit::details::CompiledGausFitCacheBaseStruct:

Classes
struct	InitializeFuncSumVector
	A helper class initializing the function vector. More...

struct	InitializeFuncSumVector< 0U, Func >

Public Types
template<unsigned int NGaus, unsigned int CutOff>
using	NGaussTruncClass = FuncSum< NGaus, gaus_trunc< CutOff >, 3U >
	Class around sum of NGaus Gaussian functions truncated at CutOff sigmas. More...

Public Member Functions
virtual TF1 *	GetClone (size_t nGaus)
	Throws an exception (ROOT does not support cloning compiled functions) More...

virtual unsigned int	MaxGaussians () const
	Returns the maximum number of Gaussians in a function that we support. More...

template<>
Double_t	ngaus (Double_t const x, Double_t const params)

std::string	GetName () const
	Return the name of this cache. More...

virtual TF1 *	Get (size_t nFunc)
	Returns a function sum of nFunc base functions. More...

virtual std::string	FunctionName (size_t nFunc) const
	Returns a name for the function with nFunc base functions. More...

Static Public Member Functions
static Double_t	gaus (Double_t const x, Double_t const params)
	Single Gaussian function. More...

template<unsigned int CutOff>
static Double_t	gaus_trunc (Double_t const x, Double_t const params)

template<unsigned int NGaus>
static Double_t	ngaus (Double_t const x, Double_t const params)

template<unsigned int NGaus, unsigned int CutOff>
static Double_t	ngaus_trunc (Double_t const x, Double_t const params)
	Sum of NGaus Gaussian functions truncated at CutOff sigmas. More...

Protected Member Functions
template<unsigned int NGaus>
void	InitializeCompiledGausFitVector ()
	Returns a vector initialized with multigaussians. More...

template<unsigned int NGaus>
void	AppendFunction ()
	Adds one function. More...

void	CannotCreateFunction (size_t nGaus) const
	Throws an error asserting compiled functions can't be cretead run-time. More...

template<>
void	InitializeCompiledGausFitVector ()

virtual TF1 *	CreateFunction (size_t nFunc) const
	Creates a new sum function. More...

Protected Attributes
std::string	name
	name of the cache More...

std::vector< TF1 * >	funcs

Detailed Description

Definition at line 152 of file GausFitCache.h.

Member Typedef Documentation

template<unsigned int NGaus, unsigned int CutOff>

using hit::details::CompiledGausFitCacheBaseStruct::NGaussTruncClass = FuncSum<NGaus, gaus_trunc<CutOff>, 3U>

Class around sum of NGaus Gaussian functions truncated at CutOff sigmas.

Definition at line 221 of file GausFitCache.h.

Member Function Documentation

template<unsigned int NGaus>

void hit::details::CompiledGausFitCacheBaseStruct::AppendFunction ( )

protected

Adds one function.

Definition at line 400 of file GausFitCache.h.

References hit::GausFitCache::funcs, and hit::GausFitCache::FunctionName().

                                                         {
       // create a function in the ficticious range [ 0, 1 ]:
       funcs.push_back
         (new TF1(FunctionName(NGaus).c_str(), &ngaus<NGaus>, 0., 1., 3*NGaus));
     } // CompiledGausFitCacheBaseStruct::AppendFunction()

void hit::details::CompiledGausFitCacheBaseStruct::CannotCreateFunction ( size_t nGaus ) const

protected

Throws an error asserting compiled functions can't be cretead run-time.

Definition at line 102 of file GausFitCache.cxx.

References art::errors::LogicError.

   {
     throw art::Exception(art::errors::LogicError)
       << name << " function cache can't create functions at run-time; " << nFunc
       << "-addend function was requested, but we have only functions with up to "
       << MaxGaussians() << " addends available\n";
   } // CompiledGausFitCache<>::CannotCreateFunction()

TF1 * hit::GausFitCache::CreateFunction ( size_t nFunc ) const

protectedvirtualinherited

Creates a new sum function.

Reimplemented in hit::CompiledTruncatedGausFitCache< MaxGaus, CutOff >, and hit::CompiledGausFitCache< MaxGaus >.

Definition at line 60 of file GausFitCache.cxx.

                                                       {
     
     std::string func_name = FunctionName(nFunc);
     
     // no gaussian, no nothing
     if (nFunc == 0) return new TF1(func_name.c_str(), "0");
     
     std::ostringstream sstr;
     sstr << "gaus(0)";
     for (size_t iGaus = 1; iGaus < nFunc; ++iGaus)
       sstr << " + gaus(" << (iGaus*3) << ")";
     
     // create and return the function
     return new TF1(func_name.c_str(), sstr.str().c_str());
   } // GausFitCache::CreateFunction()

std::string hit::GausFitCache::FunctionName ( size_t nFunc ) const

virtualinherited

Returns a name for the function with nFunc base functions.

Definition at line 78 of file GausFitCache.cxx.

Referenced by AppendFunction(), hit::details::CompiledGausFitCacheBaseStruct::InitializeFuncSumVector< NFunc, Func >::fill(), hit::details::CompiledGausFitCacheBaseStruct::InitializeFuncSumVector< 0U, Func >::fill(), and hit::GausFitCache::GetName().

                                                          {
     std::ostringstream sstr;
     sstr << name << "_" << nFunc;
     return sstr.str();
   } // GausFitCache::FunctionName()

Double_t hit::details::CompiledGausFitCacheBaseStruct::gaus	(	Double_t const *	x,
		Double_t const *	params
	)

static

Single Gaussian function.

Parameters

x	vector; [0] variable value
params	vector: [0] coefficient of the exponent [1] mean [2] sigma

Returns: the Gaussian evaluated at *x

This function is equivalent to ROOT's "gaus".

Definition at line 89 of file GausFitCache.cxx.

References geo::sqr().

   {
     return params[0] * std::exp(-0.5*sqr((x[0] - params[1])/params[2]));
   } // details::CompiledGausFitCacheBaseStruct::gaus()

template<unsigned int CutOff>

Double_t hit::details::CompiledGausFitCacheBaseStruct::gaus_trunc	(	Double_t const *	x,
		Double_t const *	params
	)

static

Definition at line 409 of file GausFitCache.h.

References hit::details::sqr(), and z.

     {
       const Double_t z = (x[0] - params[1])/params[2];
       return ((z > -((Double_t) CutOff)) && (z < (Double_t) CutOff))?
         params[0] * std::exp(-0.5*sqr(z)): 0.;
     } // CompiledGausFitCacheBaseStruct::gaus_trunc()

TF1 * hit::GausFitCache::Get ( size_t nFunc )

virtualinherited

Returns a function sum of nFunc base functions.

Parameters

nGaus the number of base functions in the function

Returns: a pointer to a TF1 with the required characteristics

The returned function must not be deleted at any time!

This implementation returns a function sum of nFunc Gaussians. The parameters are sorted by Gaussian: first normalization (not the area, but the coefficient), first mean, first sigma, second normalization etc.

Definition at line 41 of file GausFitCache.cxx.

Referenced by hit::GausFitCache::GetName().

                                      {
     
     // expand the list if needed
     if (nFunc >= funcs.size()) funcs.resize(nFunc + 1, nullptr);
     
     // get the pointer we need...
     TF1*& pFunc = funcs[nFunc];
     
     if (!pFunc) pFunc = CreateFunction(nFunc);
     
     return pFunc;
   } // GausFitCache::Get()

TF1 * hit::details::CompiledGausFitCacheBaseStruct::GetClone ( size_t nGaus )

virtual

Throws an exception (ROOT does not support cloning compiled functions)

Reimplemented from hit::GausFitCache.

Definition at line 95 of file GausFitCache.cxx.

References art::errors::LogicError.

                                                                    {
     throw art::Exception(art::errors::LogicError)
       << "CompiledGausFitCacheBaseStruct: compiled functions can't be cloned";
   } // CompiledGausFitCacheBaseStruct::GetClone()

std::string hit::GausFitCache::GetName ( ) const

inlineinherited

Return the name of this cache.

Definition at line 59 of file GausFitCache.h.

References hit::GausFitCache::FunctionName(), hit::GausFitCache::Get(), hit::GausFitCache::GetClone(), and hit::GausFitCache::name.

59 { return name; }

hit::GausFitCache::name

std::string name

name of the cache

Definition: GausFitCache.h:83

template<unsigned int NGaus>

void hit::details::CompiledGausFitCacheBaseStruct::InitializeCompiledGausFitVector ( )

protected

Returns a vector initialized with multigaussians.

Definition at line 387 of file GausFitCache.h.

                                                                          {
       if (NGaus > 0) InitializeCompiledGausFitVector<NGaus-1>();
       AppendFunction<NGaus>();
     } // CompiledGausFitCacheBaseStruct::InitializeCompiledGausFitVector()