Logo ROOT  
Reference Guide
Loading...
Searching...
No Matches
TSpectrumFit Class Reference

Advanced 1-dimensional spectra fitting functions.

Author
Miroslav Morhac

Legacy Code

TSpectrumFit is a legacy interface: there will be no bug fixes nor new developments. Therefore it is not recommended to use it in new long-term production code. But, depending on the context, using TSpectrumFit might still be a valid solution. For modeling a spectrum fitting and estimating the background one can use RooFit while for deconvolution and unfolding one can use TUnfold.

Class for fitting 1D spectra using AWMI (algorithm without matrix inversion) and conjugate gradient algorithms for symmetrical matrices (Stiefel-Hestens method). AWMI method allows to fit simultaneously 100s up to 1000s peaks. Stiefel method is very stable, it converges faster, but is more time consuming

The algorithms in this class have been published in the following references:

  1. M. Morhac et al.: Efficient fitting algorithms applied to analysis of coincidence gamma-ray spectra. Computer Physics Communications, Vol 172/1 (2005) pp. 19-41.
  2. M. Morhac et al.: Study of fitting algorithms applied to simultaneous analysis of large number of peaks in gamma-ray spectra. Applied Spectroscopy, Vol. 57, No. 7, pp. 753-760, 2003.

Definition at line 18 of file TSpectrumFit.h.

Public Types

enum  {
  kFitOptimChiCounts =0 , kFitOptimChiFuncValues =1 , kFitOptimMaxLikelihood =2 , kFitAlphaHalving =0 ,
  kFitAlphaOptimal =1 , kFitPower2 =2 , kFitPower4 =4 , kFitPower6 =6 ,
  kFitPower8 =8 , kFitPower10 =10 , kFitPower12 =12 , kFitTaylorOrderFirst =0 ,
  kFitTaylorOrderSecond =1 , kFitNumRegulCycles =100
}
enum  { kSingleKey = (1ULL << (0)) , kOverwrite = (1ULL << (1)) , kWriteDelete = (1ULL << (2)) }
enum  {
  kIsOnHeap = 0x01000000 , kNotDeleted = 0x02000000 , kZombie = 0x04000000 , kInconsistent = 0x08000000 ,
  kBitMask = 0x00ffffff
}
enum  EDeprecatedStatusBits { kObjInCanvas = (1ULL << (3)) }
enum  EStatusBits {
  kCanDelete = (1ULL << (0)) , kMustCleanup = (1ULL << (3)) , kIsReferenced = (1ULL << (4)) , kHasUUID = (1ULL << (5)) ,
  kCannotPick = (1ULL << (6)) , kNoContextMenu = (1ULL << (8)) , kInvalidObject = (1ULL << (13))
}

Public Member Functions

 TSpectrumFit (Int_t numberPeaks)
 numberPeaks: number of fitted peaks (must be greater than zero)
 TSpectrumFit (void)
 Default constructor.
 ~TSpectrumFit () override
 Destructor.
void AbstractMethod (const char *method) const
 Call this function within a function that you don't want to define as purely virtual, in order not to force all users deriving from that class to implement that maybe (on their side) unused function; but at the same time, emit a run-time warning if they try to call it, telling that it is not implemented in the derived class: action must thus be taken on the user side to override it.
virtual void AppendPad (Option_t *option="")
 Append graphics object to current pad.
virtual void Browse (TBrowser *b)
 Browse object. May be overridden for another default action.
ULong_t CheckedHash ()
 Check and record whether this class has a consistent Hash/RecursiveRemove setup (*) and then return the regular Hash value for this object.
virtual const char * ClassName () const
 Returns name of class to which the object belongs.
void Clear (Option_t *option="") override
 Set name and title to empty strings ("").
TObjectClone (const char *newname="") const override
 Make a clone of an object using the Streamer facility.
Int_t Compare (const TObject *obj) const override
 Compare two TNamed objects.
void Copy (TObject &named) const override
 Copy this to obj.
virtual void Delete (Option_t *option="")
 Delete this object.
virtual Int_t DistancetoPrimitive (Int_t px, Int_t py)
 Computes distance from point (px,py) to the object.
virtual void Draw (Option_t *option="")
 Default Draw method for all objects.
virtual void DrawClass () const
 Draw class inheritance tree of the class to which this object belongs.
virtual TObjectDrawClone (Option_t *option="") const
 Draw a clone of this object in the current selected pad with: gROOT->SetSelectedPad(c1).
virtual void Dump () const
 Dump contents of object on stdout.
virtual void Error (const char *method, const char *msgfmt,...) const
 Issue error message.
virtual void Execute (const char *method, const char *params, Int_t *error=nullptr)
 Execute method on this object with the given parameter string, e.g.
virtual void Execute (TMethod *method, TObjArray *params, Int_t *error=nullptr)
 Execute method on this object with parameters stored in the TObjArray.
virtual void ExecuteEvent (Int_t event, Int_t px, Int_t py)
 Execute action corresponding to an event at (px,py).
virtual void Fatal (const char *method, const char *msgfmt,...) const
 Issue fatal error message.
virtual void FillBuffer (char *&buffer)
 Encode TNamed into output buffer.
virtual TObjectFindObject (const char *name) const
 Must be redefined in derived classes.
virtual TObjectFindObject (const TObject *obj) const
 Must be redefined in derived classes.
void FitAwmi (Double_t *source)
 This function fits the source spectrum.
void FitStiefel (Double_t *source)
 This function fits the source spectrum.
Double_tGetAmplitudes () const
Double_tGetAmplitudesErrors () const
Double_tGetAreas () const
Double_tGetAreasErrors () const
void GetBackgroundParameters (Double_t &a0, Double_t &a0Err, Double_t &a1, Double_t &a1Err, Double_t &a2, Double_t &a2Err)
 This function gets the background parameters and their errors.
Double_t GetChi () const
virtual Option_tGetDrawOption () const
 Get option used by the graphics system to draw this object.
virtual const char * GetIconName () const
 Returns mime type name of object.
const char * GetName () const override
 Returns name of object.
virtual char * GetObjectInfo (Int_t px, Int_t py) const
 Returns string containing info about the object at position (px,py).
virtual Option_tGetOption () const
Double_tGetPositions () const
Double_tGetPositionsErrors () const
void GetSigma (Double_t &sigma, Double_t &sigmaErr)
 This function gets the sigma parameter and its error.
void GetTailParameters (Double_t &t, Double_t &tErr, Double_t &b, Double_t &bErr, Double_t &s, Double_t &sErr)
 This function gets the tail parameters and their errors.
const char * GetTitle () const override
 Returns title of object.
virtual UInt_t GetUniqueID () const
 Return the unique object id.
virtual Bool_t HandleTimer (TTimer *timer)
 Execute action in response of a timer timing out.
ULong_t Hash () const override
 Return hash value for this object.
Bool_t HasInconsistentHash () const
 Return true is the type of this object is known to have an inconsistent setup for Hash and RecursiveRemove (i.e.
virtual void Info (const char *method, const char *msgfmt,...) const
 Issue info message.
virtual Bool_t InheritsFrom (const char *classname) const
 Returns kTRUE if object inherits from class "classname".
virtual Bool_t InheritsFrom (const TClass *cl) const
 Returns kTRUE if object inherits from TClass cl.
virtual void Inspect () const
 Dump contents of this object in a graphics canvas.
void InvertBit (UInt_t f)
TClassIsA () const override
Bool_t IsDestructed () const
 IsDestructed.
virtual Bool_t IsEqual (const TObject *obj) const
 Default equal comparison (objects are equal if they have the same address in memory).
virtual Bool_t IsFolder () const
 Returns kTRUE in case object contains browsable objects (like containers or lists of other objects).
Bool_t IsOnHeap () const
Bool_t IsSortable () const override
Bool_t IsZombie () const
void ls (Option_t *option="") const override
 List TNamed name and title.
void MayNotUse (const char *method) const
 Use this method to signal that a method (defined in a base class) may not be called in a derived class (in principle against good design since a child class should not provide less functionality than its parent, however, sometimes it is necessary).
virtual Bool_t Notify ()
 This method must be overridden to handle object notification (the base implementation is no-op).
void Obsolete (const char *method, const char *asOfVers, const char *removedFromVers) const
 Use this method to declare a method obsolete.
void operator delete (void *, size_t)
 Operator delete for sized deallocation.
void operator delete (void *ptr)
 Operator delete.
void operator delete (void *ptr, void *vp)
 Only called by placement new when throwing an exception.
void operator delete[] (void *, size_t)
 Operator delete [] for sized deallocation.
void operator delete[] (void *ptr)
 Operator delete [].
void operator delete[] (void *ptr, void *vp)
 Only called by placement new[] when throwing an exception.
void * operator new (size_t sz)
void * operator new (size_t sz, void *vp)
void * operator new[] (size_t sz)
void * operator new[] (size_t sz, void *vp)
virtual void Paint (Option_t *option="")
 This method must be overridden if a class wants to paint itself.
virtual void Pop ()
 Pop on object drawn in a pad to the top of the display list.
void Print (Option_t *option="") const override
 Print TNamed name and title.
virtual Int_t Read (const char *name)
 Read contents of object with specified name from the current directory.
virtual void RecursiveRemove (TObject *obj)
 Recursively remove this object from a list.
void ResetBit (UInt_t f)
virtual void SaveAs (const char *filename="", Option_t *option="") const
 Save this object in the file specified by filename.
virtual void SavePrimitive (std::ostream &out, Option_t *option="")
 Save a primitive as a C++ statement(s) on output stream "out".
void SetBackgroundParameters (Double_t a0Init, Bool_t fixA0, Double_t a1Init, Bool_t fixA1, Double_t a2Init, Bool_t fixA2)
 This function sets the following fitting parameters of background:
void SetBit (UInt_t f)
void SetBit (UInt_t f, Bool_t set)
 Set or unset the user status bits as specified in f.
virtual void SetDrawOption (Option_t *option="")
 Set drawing option for object.
void SetFitParameters (Int_t xmin, Int_t xmax, Int_t numberIterations, Double_t alpha, Int_t statisticType, Int_t alphaOptim, Int_t power, Int_t fitTaylor)
 This function sets the following fitting parameters:
virtual void SetName (const char *name)
 Set the name of the TNamed.
virtual void SetNameTitle (const char *name, const char *title)
 Set all the TNamed parameters (name and title).
void SetPeakParameters (Double_t sigma, Bool_t fixSigma, const Double_t *positionInit, const Bool_t *fixPosition, const Double_t *ampInit, const Bool_t *fixAmp)
 This function sets the following fitting parameters of peaks:
void SetTailParameters (Double_t tInit, Bool_t fixT, Double_t bInit, Bool_t fixB, Double_t sInit, Bool_t fixS)
 This function sets the following fitting parameters of tails of peaks.
virtual void SetTitle (const char *title="")
 Set the title of the TNamed.
virtual void SetUniqueID (UInt_t uid)
 Set the unique object id.
virtual Int_t Sizeof () const
 Return size of the TNamed part of the TObject.
void Streamer (TBuffer &) override
 Stream an object of class TObject.
void StreamerNVirtual (TBuffer &ClassDef_StreamerNVirtual_b)
virtual void SysError (const char *method, const char *msgfmt,...) const
 Issue system error message.
Bool_t TestBit (UInt_t f) const
Int_t TestBits (UInt_t f) const
virtual void UseCurrentStyle ()
 Set current style settings in this object This function is called when either TCanvas::UseCurrentStyle or TROOT::ForceStyle have been invoked.
virtual void Warning (const char *method, const char *msgfmt,...) const
 Issue warning message.
virtual Int_t Write (const char *name=nullptr, Int_t option=0, Int_t bufsize=0)
 Write this object to the current directory.
virtual Int_t Write (const char *name=nullptr, Int_t option=0, Int_t bufsize=0) const
 Write this object to the current directory.

Static Public Member Functions

static TClassClass ()
static const char * Class_Name ()
static constexpr Version_t Class_Version ()
static const char * DeclFileName ()
static Longptr_t GetDtorOnly ()
 Return destructor only flag.
static Bool_t GetObjectStat ()
 Get status of object stat flag.
static void SetDtorOnly (void *obj)
 Set destructor only flag.
static void SetObjectStat (Bool_t stat)
 Turn on/off tracking of objects in the TObjectTable.

Protected Types

enum  { kOnlyPrepStep = (1ULL << (3)) }

Protected Member Functions

Double_t Area (Double_t a, Double_t sigma, Double_t t, Double_t b)
 This function calculates area of a peak Function parameters:
Double_t Dera1 (Double_t i)
 Derivative of background according to a1.
Double_t Dera2 (Double_t i)
 Derivative of background according to a2.
Double_t Deramp (Double_t i, Double_t i0, Double_t sigma, Double_t t, Double_t s, Double_t b)
 This function calculates derivative of peak shape function (see manual) according to amplitude of peak.
Double_t Derb (Int_t num_of_fitted_peaks, Double_t i, const Double_t *parameter, Double_t sigma, Double_t t, Double_t b)
 This function calculates derivative of peaks shape function (see manual) according to slope b.
Double_t Derderi0 (Double_t i, Double_t amp, Double_t i0, Double_t sigma)
 This function calculates second derivative of peak shape function (see manual) according to peak position.
Double_t Derdersigma (Int_t num_of_fitted_peaks, Double_t i, const Double_t *parameter, Double_t sigma)
 This function calculates second derivative of peaks shape function (see manual) according to sigma of peaks.
Double_t Derfc (Double_t x)
 This function calculates derivative of error function of x.
Double_t Deri0 (Double_t i, Double_t amp, Double_t i0, Double_t sigma, Double_t t, Double_t s, Double_t b)
 This function calculates derivative of peak shape function (see manual) according to peak position.
Double_t Derpa (Double_t sigma, Double_t t, Double_t b)
 This function calculates derivative of the area of peak according to its amplitude.
Double_t Derpb (Double_t a, Double_t sigma, Double_t t, Double_t b)
 This function calculates derivative of the area of peak according to b parameter.
Double_t Derpsigma (Double_t a, Double_t t, Double_t b)
 This function calculates derivative of the area of peak according to sigma of peaks.
Double_t Derpt (Double_t a, Double_t sigma, Double_t b)
 This function calculates derivative of the area of peak according to t parameter.
Double_t Ders (Int_t num_of_fitted_peaks, Double_t i, const Double_t *parameter, Double_t sigma)
 This function calculates derivative of peaks shape function (see manual) according to relative amplitude s.
Double_t Dersigma (Int_t num_of_fitted_peaks, Double_t i, const Double_t *parameter, Double_t sigma, Double_t t, Double_t s, Double_t b)
 This function calculates derivative of peaks shape function (see manual) according to sigma of peaks.
Double_t Dert (Int_t num_of_fitted_peaks, Double_t i, const Double_t *parameter, Double_t sigma, Double_t b)
 This function calculates derivative of peaks shape function (see manual) according to relative amplitude t.
virtual void DoError (int level, const char *location, const char *fmt, va_list va) const
 Interface to ErrorHandler (protected).
Double_t Erfc (Double_t x)
void MakeZombie ()
Double_t Ourpowl (Double_t a, Int_t pw)
 Power function.
void SavePrimitiveNameTitle (std::ostream &out, const char *variable_name)
 Save object name and title into the output stream "out".
Double_t Shape (Int_t num_of_fitted_peaks, Double_t i, const Double_t *parameter, Double_t sigma, Double_t t, Double_t s, Double_t b, Double_t a0, Double_t a1, Double_t a2)
 This function calculates peaks shape function (see manual) Function parameters:
void StiefelInversion (Double_t **a, Int_t rozmer)
 This function calculates solution of the system of linear equations.

Static Protected Member Functions

static void SavePrimitiveConstructor (std::ostream &out, TClass *cl, const char *variable_name, const char *constructor_agrs="", Bool_t empty_line=kTRUE)
 Save object constructor in the output stream "out".
static void SavePrimitiveDraw (std::ostream &out, const char *variable_name, Option_t *option=nullptr)
 Save invocation of primitive Draw() method Skipped if option contains "nodraw" string.
static TString SavePrimitiveVector (std::ostream &out, const char *prefix, Int_t len, Double_t *arr, Int_t flag=0)
 Save array in the output stream "out" as vector.

Protected Attributes

Double_t fA0Calc
 calculated value of background a0 parameter
Double_t fA0Err
 error value of background a0 parameter
Double_t fA0Init
 initial value of background a0 parameter(backgroud is estimated as a0+a1*x+a2*x*x)
Double_t fA1Calc
 calculated value of background a1 parameter
Double_t fA1Err
 error value of background a1 parameter
Double_t fA1Init
 initial value of background a1 parameter(backgroud is estimated as a0+a1*x+a2*x*x)
Double_t fA2Calc
 calculated value of background a2 parameter
Double_t fA2Err
 error value of background a2 parameter
Double_t fA2Init
 initial value of background a2 parameter(backgroud is estimated as a0+a1*x+a2*x*x)
Double_t fAlpha
 convergence coefficient, input parameter, it should be positive number and <=1, for details see references
Int_t fAlphaOptim
 optimization of convergence algorithm, possible values kFitAlphaHalving, kFitAlphaOptimal
Double_tfAmpCalc
 [fNPeaks] array of calculated values of fitted amplitudes, output parameters
Double_tfAmpErr
 [fNPeaks] array of amplitude errors
Double_tfAmpInit
 [fNPeaks] array of initial values of peaks amplitudes, input parameters
Double_tfArea
 [fNPeaks] array of calculated areas of peaks
Double_tfAreaErr
 [fNPeaks] array of errors of peak areas
Double_t fBCalc
 calculated value of b parameter
Double_t fBErr
 error value of b parameter
Double_t fBInit
 initial value of b parameter (slope), for details see html manual and references
Double_t fChi
 here the fitting functions return resulting chi square
Int_t fFitTaylor
 order of Taylor expansion, possible values kFitTaylorOrderFirst, kFitTaylorOrderSecond. It applies only for Awmi fitting function.
Bool_t fFixA0
 logical value of a0 parameter, which allows to fix the parameter (not to fit).
Bool_t fFixA1
 logical value of a1 parameter, which allows to fix the parameter (not to fit).
Bool_t fFixA2
 logical value of a2 parameter, which allows to fix the parameter (not to fit).
Bool_tfFixAmp
 [fNPeaks] array of logical values which allow to fix appropriate amplitudes (not fit). However they are present in the estimated functional
Bool_t fFixB
 logical value of b parameter, which allows to fix the parameter (not to fit).
Bool_tfFixPosition
 [fNPeaks] array of logical values which allow to fix appropriate positions (not fit). However they are present in the estimated functional
Bool_t fFixS
 logical value of s parameter, which allows to fix the parameter (not to fit).
Bool_t fFixSigma
 logical value of sigma parameter, which allows to fix the parameter (not to fit).
Bool_t fFixT
 logical value of t parameter, which allows to fix the parameter (not to fit).
TString fName
Int_t fNPeaks
 number of peaks present in fit, input parameter, it should be > 0
Int_t fNumberIterations
 number of iterations in fitting procedure, input parameter, it should be > 0
Double_tfPositionCalc
 [fNPeaks] array of calculated values of fitted positions, output parameters
Double_tfPositionErr
 [fNPeaks] array of position errors
Double_tfPositionInit
 [fNPeaks] array of initial values of peaks positions, input parameters
Int_t fPower
 possible values kFitPower2,4,6,8,10,12, for details see references. It applies only for Awmi fitting function.
Double_t fSCalc
 calculated value of s parameter
Double_t fSErr
 error value of s parameter
Double_t fSigmaCalc
 calculated value of sigma parameter
Double_t fSigmaErr
 error value of sigma parameter
Double_t fSigmaInit
 initial value of sigma parameter
Double_t fSInit
 initial value of s parameter (relative amplitude of step), for details see html manual and references
Int_t fStatisticType
 type of statistics, possible values kFitOptimChiCounts (chi square statistics with counts as weighting coefficients), kFitOptimChiFuncValues (chi square statistics with function values as weighting coefficients),kFitOptimMaxLikelihood
Double_t fTCalc
 calculated value of t parameter
Double_t fTErr
 error value of t parameter
Double_t fTInit
 initial value of t parameter (relative amplitude of tail), for details see html manual and references
TString fTitle
Int_t fXmax
 last fitted channel
Int_t fXmin
 first fitted channel

Static Private Member Functions

static void AddToTObjectTable (TObject *)
 Private helper function which will dispatch to TObjectTable::AddObj.

Private Attributes

UInt_t fBits
 bit field status word
UInt_t fUniqueID
 object unique identifier

Static Private Attributes

static Longptr_t fgDtorOnly = 0
 object for which to call dtor only (i.e. no delete)
static Bool_t fgObjectStat = kTRUE
 if true keep track of objects in TObjectTable

#include <TSpectrumFit.h>

Inheritance diagram for TSpectrumFit:
TNamed TObject

Member Enumeration Documentation

◆ anonymous enum

anonymous enum
Enumerator
kFitOptimChiCounts 
kFitOptimChiFuncValues 
kFitOptimMaxLikelihood 
kFitAlphaHalving 
kFitAlphaOptimal 
kFitPower2 
kFitPower4 
kFitPower6 
kFitPower8 
kFitPower10 
kFitPower12 
kFitTaylorOrderFirst 
kFitTaylorOrderSecond 
kFitNumRegulCycles 

Definition at line 70 of file TSpectrumFit.h.

◆ anonymous enum

anonymous enum
protectedinherited
Enumerator
kOnlyPrepStep 

Used to request that the class specific implementation of TObject::Write just prepare the objects to be ready to be written but do not actually write them into the TBuffer.

This is just for example by TBufferMerger to request that the TTree inside the file calls TTree::FlushBaskets (outside of the merging lock) and TBufferMerger will later ask for the write (inside the merging lock). To take advantage of this feature the class needs to overload TObject::Write and use this enum value accordingly. (See TTree::Write and TObject::Write) Do not use, this feature will be migrate to the Merge function (See TClass and TTree::Merge)

Definition at line 106 of file TObject.h.

◆ anonymous enum

anonymous enum
inherited
Enumerator
kSingleKey 

write collection with single key

kOverwrite 

overwrite existing object with same name

kWriteDelete 

write object, then delete previous key with same name

Definition at line 99 of file TObject.h.

◆ anonymous enum

anonymous enum
inherited
Enumerator
kIsOnHeap 

object is on heap

kNotDeleted 

object has not been deleted

kZombie 

object ctor failed

kInconsistent 

class overload Hash but does call RecursiveRemove in destructor

kBitMask 

Definition at line 89 of file TObject.h.

◆ EDeprecatedStatusBits

Enumerator
kObjInCanvas 

for backward compatibility only, use kMustCleanup

Definition at line 84 of file TObject.h.

◆ EStatusBits

enum TObject::EStatusBits
inherited
Enumerator
kCanDelete 

if object in a list can be deleted

kMustCleanup 

if object destructor must call RecursiveRemove()

kIsReferenced 

if object is referenced by a TRef or TRefArray

kHasUUID 

if object has a TUUID (its fUniqueID=UUIDNumber)

kCannotPick 

if object in a pad cannot be picked

kNoContextMenu 

if object does not want context menu

kInvalidObject 

if object ctor succeeded but object should not be used

Definition at line 70 of file TObject.h.

Constructor & Destructor Documentation

◆ TSpectrumFit() [1/2]

TSpectrumFit::TSpectrumFit ( void )

Default constructor.

Definition at line 36 of file TSpectrumFit.cxx.

◆ TSpectrumFit() [2/2]

TSpectrumFit::TSpectrumFit ( Int_t numberPeaks)

numberPeaks: number of fitted peaks (must be greater than zero)

the constructor allocates arrays for all fitted parameters (peak positions, amplitudes etc) and sets the member variables to their default values. One can change these variables by member functions (setters) of TSpectrumFit class.

Shape function of the fitted peaks is

where a represents vector of fitted parameters (positions p(j), amplitudes A(j), sigma, relative amplitudes T, S and slope B).

Definition at line 102 of file TSpectrumFit.cxx.

◆ ~TSpectrumFit()

TSpectrumFit::~TSpectrumFit ( )
override

Destructor.

Definition at line 161 of file TSpectrumFit.cxx.

Member Function Documentation

◆ AbstractMethod()

void TObject::AbstractMethod ( const char * method) const
inherited

Call this function within a function that you don't want to define as purely virtual, in order not to force all users deriving from that class to implement that maybe (on their side) unused function; but at the same time, emit a run-time warning if they try to call it, telling that it is not implemented in the derived class: action must thus be taken on the user side to override it.

In other word, this method acts as a "runtime purely virtual" warning instead of a "compiler purely virtual" error.

Warning
This interface is a legacy function that is no longer recommended to be used by new development code.
Note
The name "AbstractMethod" does not imply that it's an abstract method in the strict C++ sense.

Definition at line 1149 of file TObject.cxx.

◆ AddToTObjectTable()

void TObject::AddToTObjectTable ( TObject * op)
staticprivateinherited

Private helper function which will dispatch to TObjectTable::AddObj.

Included here to avoid circular dependency between header files.

Definition at line 195 of file TObject.cxx.

◆ AppendPad()

void TObject::AppendPad ( Option_t * option = "")
virtualinherited

Append graphics object to current pad.

In case no current pad is set yet, create a default canvas with the name "c1".

Definition at line 204 of file TObject.cxx.

◆ Area()

Double_t TSpectrumFit::Area ( Double_t a,
Double_t sigma,
Double_t t,
Double_t b )
protected

This function calculates area of a peak Function parameters:

  • a-amplitude of the peak
  • sigma-sigma of peak
  • t-relative amplitude
  • b-slope

Definition at line 570 of file TSpectrumFit.cxx.

◆ Browse()

◆ CheckedHash()

ULong_t TObject::CheckedHash ( )
inlineinherited

Check and record whether this class has a consistent Hash/RecursiveRemove setup (*) and then return the regular Hash value for this object.

The intent is for this routine to be called instead of directly calling the function Hash during "insert" operations. See TObject::HasInconsistenTObjectHash();

(*) The setup is consistent when all classes in the class hierarchy that overload TObject::Hash do call ROOT::CallRecursiveRemoveIfNeeded in their destructor. i.e. it is safe to call the Hash virtual function during the RecursiveRemove operation.

Definition at line 332 of file TObject.h.

◆ Class()

TClass * TSpectrumFit::Class ( )
static
Returns
TClass describing this class

◆ Class_Name()

const char * TSpectrumFit::Class_Name ( )
static
Returns
Name of this class

◆ Class_Version()

constexpr Version_t TSpectrumFit::Class_Version ( )
inlinestaticconstexpr
Returns
Version of this class

Definition at line 131 of file TSpectrumFit.h.

◆ ClassName()

const char * TObject::ClassName ( ) const
virtualinherited

Returns name of class to which the object belongs.

Definition at line 227 of file TObject.cxx.

◆ Clear()

void TNamed::Clear ( Option_t * option = "")
overridevirtualinherited

Set name and title to empty strings ("").

Reimplemented from TObject.

Reimplemented in TPrincipal, TProcessID, TStreamerInfo, TTask, TVirtualFitter, and TVirtualStreamerInfo.

Definition at line 63 of file TNamed.cxx.

◆ Clone()

TObject * TNamed::Clone ( const char * newname = "") const
overridevirtualinherited

Make a clone of an object using the Streamer facility.

If newname is specified, this will be the name of the new object.

Reimplemented from TObject.

Reimplemented in TStreamerInfo, and TTreeIndex.

Definition at line 73 of file TNamed.cxx.

◆ Compare()

Int_t TNamed::Compare ( const TObject * obj) const
overridevirtualinherited

Compare two TNamed objects.

Returns 0 when equal, -1 when this is smaller and +1 when bigger (like strcmp).

Reimplemented from TObject.

Reimplemented in TStructNodeProperty.

Definition at line 84 of file TNamed.cxx.

◆ Copy()

void TNamed::Copy ( TObject & named) const
overridevirtualinherited

Copy this to obj.

Reimplemented from TObject.

Reimplemented in TPieSlice, TProfile2D, TProfile3D, TProfile, TStyle, TSystemDirectory, TSystemFile, TText, and TXTRU.

Definition at line 93 of file TNamed.cxx.

◆ DeclFileName()

const char * TSpectrumFit::DeclFileName ( )
inlinestatic
Returns
Name of the file containing the class declaration

Definition at line 131 of file TSpectrumFit.h.

◆ Delete()

void TObject::Delete ( Option_t * option = "")
virtualinherited

◆ Dera1()

Double_t TSpectrumFit::Dera1 ( Double_t i)
protected

Derivative of background according to a1.

Definition at line 493 of file TSpectrumFit.cxx.

◆ Dera2()

Double_t TSpectrumFit::Dera2 ( Double_t i)
protected

Derivative of background according to a2.

Definition at line 501 of file TSpectrumFit.cxx.

◆ Deramp()

Double_t TSpectrumFit::Deramp ( Double_t i,
Double_t i0,
Double_t sigma,
Double_t t,
Double_t s,
Double_t b )
protected

This function calculates derivative of peak shape function (see manual) according to amplitude of peak.

Function parameters:

  • i-channel
  • i0-position of peak
  • sigma-sigma of peak
  • t, s-relative amplitudes
  • b-slope

Definition at line 230 of file TSpectrumFit.cxx.

◆ Derb()

Double_t TSpectrumFit::Derb ( Int_t num_of_fitted_peaks,
Double_t i,
const Double_t * parameter,
Double_t sigma,
Double_t t,
Double_t b )
protected

This function calculates derivative of peaks shape function (see manual) according to slope b.

Function parameters:

  • num_of_fitted_peaks-number of fitted peaks
  • i-channel
  • parameter-array of peaks parameters (amplitudes and positions)
  • sigma-sigma of peak
  • t-relative amplitude
  • b-slope

Definition at line 464 of file TSpectrumFit.cxx.

◆ Derderi0()

Double_t TSpectrumFit::Derderi0 ( Double_t i,
Double_t amp,
Double_t i0,
Double_t sigma )
protected

This function calculates second derivative of peak shape function (see manual) according to peak position.

Function parameters:

  • i-channel
  • amp-amplitude of peak
  • i0-position of peak
  • sigma-width of peak

Definition at line 304 of file TSpectrumFit.cxx.

◆ Derdersigma()

Double_t TSpectrumFit::Derdersigma ( Int_t num_of_fitted_peaks,
Double_t i,
const Double_t * parameter,
Double_t sigma )
protected

This function calculates second derivative of peaks shape function (see manual) according to sigma of peaks.

Function parameters:

  • num_of_fitted_peaks-number of fitted peaks
  • i-channel
  • parameter-array of peaks parameters (amplitudes and positions)
  • sigma-sigma of peak

Definition at line 377 of file TSpectrumFit.cxx.

◆ Derfc()

Double_t TSpectrumFit::Derfc ( Double_t x)
protected

This function calculates derivative of error function of x.

Definition at line 201 of file TSpectrumFit.cxx.

◆ Deri0()

Double_t TSpectrumFit::Deri0 ( Double_t i,
Double_t amp,
Double_t i0,
Double_t sigma,
Double_t t,
Double_t s,
Double_t b )
protected

This function calculates derivative of peak shape function (see manual) according to peak position.

Function parameters:

  • i-channel
  • amp-amplitude of peak
  • i0-position of peak
  • sigma-sigma of peak
  • t, s-relative amplitudes
  • b-slope

Definition at line 267 of file TSpectrumFit.cxx.

◆ Derpa()

Double_t TSpectrumFit::Derpa ( Double_t sigma,
Double_t t,
Double_t b )
protected

This function calculates derivative of the area of peak according to its amplitude.

Function parameters:

  • sigma-sigma of peak
  • t-relative amplitudes
  • b-slope

Definition at line 593 of file TSpectrumFit.cxx.

◆ Derpb()

Double_t TSpectrumFit::Derpb ( Double_t a,
Double_t sigma,
Double_t t,
Double_t b )
protected

This function calculates derivative of the area of peak according to b parameter.

Function parameters:

  • sigma-sigma of peak
  • t-relative amplitudes
  • b-slope

Definition at line 659 of file TSpectrumFit.cxx.

◆ Derpsigma()

Double_t TSpectrumFit::Derpsigma ( Double_t a,
Double_t t,
Double_t b )
protected

This function calculates derivative of the area of peak according to sigma of peaks.

Function parameters:

  • a-amplitude of peak
  • t-relative amplitudes
  • b-slope

Definition at line 615 of file TSpectrumFit.cxx.

◆ Derpt()

Double_t TSpectrumFit::Derpt ( Double_t a,
Double_t sigma,
Double_t b )
protected

This function calculates derivative of the area of peak according to t parameter.

Function parameters:

  • sigma-sigma of peak
  • t-relative amplitudes
  • b-slope

Definition at line 637 of file TSpectrumFit.cxx.

◆ Ders()

Double_t TSpectrumFit::Ders ( Int_t num_of_fitted_peaks,
Double_t i,
const Double_t * parameter,
Double_t sigma )
protected

This function calculates derivative of peaks shape function (see manual) according to relative amplitude s.

Function parameters:

  • num_of_fitted_peaks-number of fitted peaks
  • i-channel
  • parameter-array of peaks parameters (amplitudes and positions)
  • sigma-sigma of peak

Definition at line 438 of file TSpectrumFit.cxx.

◆ Dersigma()

Double_t TSpectrumFit::Dersigma ( Int_t num_of_fitted_peaks,
Double_t i,
const Double_t * parameter,
Double_t sigma,
Double_t t,
Double_t s,
Double_t b )
protected

This function calculates derivative of peaks shape function (see manual) according to sigma of peaks.

Function parameters:

  • num_of_fitted_peaks-number of fitted peaks
  • i-channel
  • parameter-array of peaks parameters (amplitudes and positions)
  • sigma-sigma of peak
  • t, s-relative amplitudes
  • b-slope

Definition at line 332 of file TSpectrumFit.cxx.

◆ Dert()

Double_t TSpectrumFit::Dert ( Int_t num_of_fitted_peaks,
Double_t i,
const Double_t * parameter,
Double_t sigma,
Double_t b )
protected

This function calculates derivative of peaks shape function (see manual) according to relative amplitude t.

Function parameters:

  • num_of_fitted_peaks-number of fitted peaks
  • i-channel
  • parameter-array of peaks parameters (amplitudes and positions)
  • sigma-sigma of peak
  • b-slope

Definition at line 410 of file TSpectrumFit.cxx.

◆ DistancetoPrimitive()

◆ DoError()

void TObject::DoError ( int level,
const char * location,
const char * fmt,
va_list va ) const
protectedvirtualinherited

Interface to ErrorHandler (protected).

Reimplemented in TThread, and TTreeViewer.

Definition at line 1059 of file TObject.cxx.

◆ Draw()

◆ DrawClass()

void TObject::DrawClass ( ) const
virtualinherited

Draw class inheritance tree of the class to which this object belongs.

If a class B inherits from a class A, description of B is drawn on the right side of description of A. Member functions overridden by B are shown in class A with a blue line crossing-out the corresponding member function. The following picture is the class inheritance tree of class TPaveLabel:

Reimplemented in TGFrame, TSystemDirectory, and TSystemFile.

Definition at line 308 of file TObject.cxx.

◆ DrawClone()

TObject * TObject::DrawClone ( Option_t * option = "") const
virtualinherited

Draw a clone of this object in the current selected pad with: gROOT->SetSelectedPad(c1).

If pad was not selected - gPad will be used.

Note
For histograms, use the more specialised TH1::DrawCopy().

Reimplemented in TAxis, TCanvas, TGFrame, TSystemDirectory, and TSystemFile.

Definition at line 319 of file TObject.cxx.

◆ Dump()

void TObject::Dump ( ) const
virtualinherited

Dump contents of object on stdout.

Using the information in the object dictionary (class TClass) each data member is interpreted. If a data member is a pointer, the pointer value is printed

The following output is the Dump of a TArrow object:

fAngle 0 Arrow opening angle (degrees)
fArrowSize 0.2 Arrow Size
fOption.*fData
fX1 0.1 X of 1st point
fY1 0.15 Y of 1st point
fX2 0.67 X of 2nd point
fY2 0.83 Y of 2nd point
fUniqueID 0 object unique identifier
fBits 50331648 bit field status word
fLineColor 1 line color
fLineStyle 1 line style
fLineWidth 1 line width
fFillColor 19 fill area color
fFillStyle 1001 fill area style
#define X(type, name)
UInt_t fUniqueID
object unique identifier
Definition TObject.h:46
UInt_t fBits
bit field status word
Definition TObject.h:47
TLine * line
TCanvas * style()
Definition style.C:1

Reimplemented in TClass, TCollection, TGFrame, TGPack, and TSystemFile.

Definition at line 367 of file TObject.cxx.

◆ Erfc()

Double_t TSpectrumFit::Erfc ( Double_t x)
protected

Definition at line 178 of file TSpectrumFit.cxx.

◆ Error()

void TObject::Error ( const char * location,
const char * fmt,
... ) const
virtualinherited

Issue error message.

Use "location" to specify the method where the error occurred. Accepts standard printf formatting arguments.

Reimplemented in TFitResult.

Definition at line 1098 of file TObject.cxx.

◆ Execute() [1/2]

void TObject::Execute ( const char * method,
const char * params,
Int_t * error = nullptr )
virtualinherited

Execute method on this object with the given parameter string, e.g.

"3.14,1,\"text\"".

Reimplemented in ROOT::R::TRInterface, TCling, TContextMenu, TInterpreter, and TMethodCall.

Definition at line 378 of file TObject.cxx.

◆ Execute() [2/2]

void TObject::Execute ( TMethod * method,
TObjArray * params,
Int_t * error = nullptr )
virtualinherited

Execute method on this object with parameters stored in the TObjArray.

The TObjArray should contain an argv vector like:

argv[0] ... argv[n] = the list of TObjString parameters
Collectable string class.
Definition TObjString.h:28
const Int_t n
Definition legend1.C:16

Reimplemented in ROOT::R::TRInterface, TCling, TContextMenu, TInterpreter, and TMethodCall.

Definition at line 398 of file TObject.cxx.

◆ ExecuteEvent()

◆ Fatal()

void TObject::Fatal ( const char * location,
const char * fmt,
... ) const
virtualinherited

Issue fatal error message.

Use "location" to specify the method where the fatal error occurred. Accepts standard printf formatting arguments.

Definition at line 1126 of file TObject.cxx.

◆ FillBuffer()

void TNamed::FillBuffer ( char *& buffer)
virtualinherited

Encode TNamed into output buffer.

Reimplemented in TDirectoryFile, TFile, TKey, TKeySQL, TKeyXML, TSQLFile, and TXMLFile.

Definition at line 103 of file TNamed.cxx.

◆ FindObject() [1/2]

TObject * TObject::FindObject ( const char * name) const
virtualinherited

Must be redefined in derived classes.

This function is typically used with TCollections, but can also be used to find an object by name inside this object.

Reimplemented in RooAbsCollection, RooLinkedList, TBtree, TCollection, TDirectory, TFolder, TGeometry, TGraph2D, TGraph, TH1, THashList, THashTable, THbookFile, TList, TListOfDataMembers, TListOfEnums, TListOfEnumsWithLock, TListOfFunctions, TListOfFunctionTemplates, TListOfTypes, TMap, TObjArray, TPad, TROOT, TViewPubDataMembers, and TViewPubFunctions.

Definition at line 425 of file TObject.cxx.

◆ FindObject() [2/2]

TObject * TObject::FindObject ( const TObject * obj) const
virtualinherited

Must be redefined in derived classes.

This function is typically used with TCollections, but can also be used to find an object inside this object.

Reimplemented in RooAbsCollection, RooLinkedList, TBtree, TCollection, TDirectory, TFolder, TGeometry, TGraph2D, TGraph, TH1, THashList, THashTable, THbookFile, TList, TListOfDataMembers, TListOfEnums, TListOfEnumsWithLock, TListOfFunctions, TListOfFunctionTemplates, TListOfTypes, TMap, TObjArray, TPad, TROOT, TViewPubDataMembers, and TViewPubFunctions.

Definition at line 435 of file TObject.cxx.

◆ FitAwmi()

void TSpectrumFit::FitAwmi ( Double_t * source)

This function fits the source spectrum.

The calling program should fill in input parameters of the TSpectrumFit class The fitted parameters are written into TSpectrumFit class output parameters and fitted data are written into source spectrum.

Function parameters:

  • source-pointer to the vector of source spectrum

Fitting

Goal: to estimate simultaneously peak shape parameters in spectra with large number of peaks

  • peaks can be fitted separately, each peak (or multiplets) in a region or together all peaks in a spectrum. To fit separately each peak one needs to determine the fitted region. However it can happen that the regions of neighbouring peaks are overlapping. Then the results of fitting are very poor. On the other hand, when fitting together all peaks found in a spectrum, one needs to have a method that is stable (converges) and fast enough to carry out fitting in reasonable time
  • we have implemented the non-symmetrical semi-empirical peak shape function [1]
  • it contains the symmetrical Gaussian as well as non-symmetrical terms.

where T and S are relative amplitudes and B is slope.

  • algorithm without matrix inversion (AWMI) allows fitting tens, hundreds of peaks simultaneously that represent sometimes thousands of parameters [2], [5].

References:

[1] Phillps G.W., Marlow K.W., NIM 137 (1976) 525.

[2] I. A. Slavic: Nonlinear least-squares fitting without matrix inversion applied to complex Gaussian spectra analysis. NIM 134 (1976) 285-289.

[3] T. Awaya: A new method for curve fitting to the data with low statistics not using chi-square method. NIM 165 (1979) 317-323.

[4] T. Hauschild, M. Jentschel: Comparison of maximum likelihood estimation and chi-square statistics applied to counting experiments. NIM A 457 (2001) 384-401.

[5] M. Morhac, J. Kliman, M. Jandel, L. Krupa, V. Matouoek: Study of fitting algorithms applied to simultaneous analysis of large number of peaks in -ray spectra. Applied Spectroscopy, Vol. 57, No. 7, pp. 753-760, 2003

Example - script FitAwmi.c:

Fig. 1 Original spectrum (black line) and fitted spectrum using AWMI algorithm (red line) and number of iteration steps = 1000. Positions of fitted peaks are denoted by markers

Script:

Example to illustrate fitting function using AWMI algorithm. To execute this example, do:

root > .x FitAwmi.C

void FitAwmi() {
Int_t i,nfound=0,bin;
Int_t nbins = 256;
Int_t xmin = 0;
Int_t xmax = nbins;
Double_t * source = new Double_t[nbins];
Double_t * dest = new Double_t[nbins];
TH1F *h = new TH1F("h","Fitting using AWMI algorithm",nbins,xmin,xmax);
TH1F *d = new TH1F("d","",nbins,xmin,xmax);
TFile *f = new TFile("TSpectrum.root");
h=(TH1F*) f->Get("fit;1");
for (i = 0; i < nbins; i++) source[i]=h->GetBinContent(i + 1);
TCanvas *Fit1 = gROOT->GetListOfCanvases()->FindObject("Fit1");
if (!Fit1) Fit1 = new TCanvas("Fit1","Fit1",10,10,1000,700);
h->Draw("L");
TSpectrum *s = new TSpectrum();
//searching for candidate peaks positions
nfound = s->SearchHighRes(source, dest, nbins, 2, 0.1, kFALSE, 10000, kFALSE, 0);
Bool_t *FixPos =new Bool_t[nfound];
Bool_t *FixAmp = new Bool_t[nfound];
for(i = 0; i< nfound ; i++){
FixPos[i] = kFALSE;
FixAmp[i] = kFALSE;
}
//filling in the
initial estimates of the input parameters
Double_t *PosX = new Double_t[nfound];
Double_t *PosY = new Double_t[nfound];
PosX = s->GetPositionX();
for (i = 0; i < nfound; i++) {
a=PosX[i];
bin = 1 + Int_t(a + 0.5);
PosY[i] = h->GetBinContent(bin);
}
TSpectrumFit *pfit=new TSpectrumFit(nfound);
pfit->SetFitParameters(xmin, xmax-1, 1000, 0.1, pfit->kFitOptimChiCounts,
pfit->SetPeakParameters(2, kFALSE, PosX, (Bool_t *) FixPos, PosY, (Bool_t *) FixAmp);
pfit->FitAwmi(source);
Double_t *CalcPositions = new Double_t[nfound];
Double_t *CalcAmplitudes = new Double_t[nfound];
CalcPositions=pfit->GetPositions();
CalcAmplitudes=pfit->GetAmplitudes();
for (i = 0; i < nbins; i++) d->SetBinContent(i + 1,source[i]);
d->SetLineColor(kRed);
d->Draw("SAME L");
for (i = 0; i < nfound; i++) {
a=CalcPositions[i];
bin = 1 + Int_t(a + 0.5);
PosX[i] = d->GetBinCenter(bin);
PosY[i] = d->GetBinContent(bin);
}
TPolyMarker * pm = (TPolyMarker*)h->GetListOfFunctions()->FindObject("TPolyMarker");
if (pm) {
h->GetListOfFunctions()->Remove(pm);
delete pm;
}
pm = new TPolyMarker(nfound, PosX, PosY);
h->GetListOfFunctions()->Add(pm);
pm->SetMarkerStyle(23);
pm->SetMarkerSize(1);
}
#define d(i)
Definition RSha256.hxx:102
#define f(i)
Definition RSha256.hxx:104
#define a(i)
Definition RSha256.hxx:99
#define h(i)
Definition RSha256.hxx:106
int Int_t
Signed integer 4 bytes (int).
Definition RtypesCore.h:59
bool Bool_t
Boolean (0=false, 1=true) (bool).
Definition RtypesCore.h:77
constexpr Bool_t kFALSE
Definition RtypesCore.h:108
double Double_t
Double 8 bytes.
Definition RtypesCore.h:73
@ kRed
Definition Rtypes.h:67
float xmin
float xmax
#define gROOT
Definition TROOT.h:417
virtual void SetMarkerColor(Color_t mcolor=1)
Set the marker color.
Definition TAttMarker.h:41
virtual void SetMarkerStyle(Style_t mstyle=1)
Set the marker style.
Definition TAttMarker.h:43
virtual void SetMarkerSize(Size_t msize=1)
Set the marker size.
Definition TAttMarker.h:48
The Canvas class.
Definition TCanvas.h:23
A file, usually with extension .root, that stores data and code in the form of serialized objects in ...
Definition TFile.h:130
1-D histogram with a float per channel (see TH1 documentation)
Definition TH1.h:878
void SetPeakParameters(Double_t sigma, Bool_t fixSigma, const Double_t *positionInit, const Bool_t *fixPosition, const Double_t *ampInit, const Bool_t *fixAmp)
This function sets the following fitting parameters of peaks:
void FitAwmi(Double_t *source)
This function fits the source spectrum.
TSpectrumFit(void)
Default constructor.
Double_t * GetAmplitudes() const
void SetFitParameters(Int_t xmin, Int_t xmax, Int_t numberIterations, Double_t alpha, Int_t statisticType, Int_t alphaOptim, Int_t power, Int_t fitTaylor)
This function sets the following fitting parameters:
Double_t * GetPositions() const
Advanced Spectra Processing.
Definition TSpectrum.h:18
Int_t SearchHighRes(Double_t *source, Double_t *destVector, Int_t ssize, Double_t sigma, Double_t threshold, bool backgroundRemove, Int_t deconIterations, bool markov, Int_t averWindow)
Double_t * GetPositionX() const
Definition TSpectrum.h:58

Definition at line 821 of file TSpectrumFit.cxx.

◆ FitStiefel()

void TSpectrumFit::FitStiefel ( Double_t * source)

This function fits the source spectrum.

The calling program should fill in input parameters The fitted parameters are written into output parameters and fitted data are written into source spectrum.

Function parameters:

  • source-pointer to the vector of source spectrum

Example - script FitStiefel.c:

Fig. 2 Original spectrum (black line) and fitted spectrum using Stiefel-Hestens method (red line) and number of iteration steps = 100. Positions of fitted peaks are denoted by markers

Script:

Example to illustrate fitting function using Stiefel-Hestens method. To execute this example, do:

root > .x FitStiefel.C

void FitStiefel() {
Int_t i,nfound=0,bin;
Int_t nbins = 256;
Int_t xmin = 0;
Int_t xmax = nbins;
Double_t * source = new Double_t[nbins];
Double_t * dest = new Double_t[nbins];
TH1F *h = new TH1F("h","Fitting using AWMI algorithm",nbins,xmin,xmax);
TH1F *d = new TH1F("d","",nbins,xmin,xmax);
TFile *f = new TFile("TSpectrum.root");
h=(TH1F*) f->Get("fit;1");
for (i = 0; i < nbins;i++) source[i]=h->GetBinContent(i + 1);
TCanvas *Fit1 = gROOT->GetListOfCanvases()->FindObject("Fit1");
if (!Fit1) Fit1 = new TCanvas("Fit1","Fit1",10,10,1000,700);
h->Draw("L");
TSpectrum *s = new TSpectrum();
//searching for candidate peaks positions
nfound = s->SearchHighRes(source, dest, nbins, 2, 0.1, kFALSE, 10000, kFALSE, 0);
Bool_t *FixPos = new Bool_t[nfound];
Bool_t *FixAmp = new Bool_t[nfound];
for(i = 0; i< nfound ; i++){
FixPos[i] = kFALSE;
FixAmp[i] = kFALSE;
}
//filling in the initial estimates of the input parameters
Double_t *PosX = new Double_t[nfound];
Double_t *PosY = new Double_t[nfound];
PosX = s->GetPositionX();
for (i = 0; i < nfound; i++) {
a=PosX[i];
bin = 1 + Int_t(a + 0.5);
PosY[i] = h->GetBinContent(bin);
}
TSpectrumFit *pfit = new TSpectrumFit(nfound);
pfit->SetFitParameters(xmin, xmax-1, 1000, 0.1, pfit->kFitOptimChiCounts,
pfit->SetPeakParameters(2, kFALSE, PosX, (Bool_t *) FixPos, PosY, (Bool_t *) FixAmp);
pfit->FitStiefel(source);
Double_t *CalcPositions = new Double_t[nfound];
Double_t *CalcAmplitudes = new Double_t[nfound];
CalcPositions=pfit->GetPositions();
CalcAmplitudes=pfit->GetAmplitudes();
for (i = 0; i < nbins; i++) d->SetBinContent(i + 1,source[i]);
d->SetLineColor(kRed);
d->Draw("SAMEL");
for (i = 0; i < nfound; i++) {
a=CalcPositions[i];
bin = 1 + Int_t(a + 0.5);
PosX[i] = d->GetBinCenter(bin);
PosY[i] = d->GetBinContent(bin);
}
TPolyMarker * pm = (TPolyMarker*)h->GetListOfFunctions()->FindObject("TPolyMarker");
if (pm) {
h->GetListOfFunctions()->Remove(pm);
delete pm;
}
pm = new TPolyMarker(nfound, PosX, PosY);
h->GetListOfFunctions()->Add(pm);
pm->SetMarkerStyle(23);
pm->SetMarkerSize(1);
}
void FitStiefel(Double_t *source)
This function fits the source spectrum.

Definition at line 1858 of file TSpectrumFit.cxx.

◆ GetAmplitudes()

Double_t * TSpectrumFit::GetAmplitudes ( ) const
inline

Definition at line 116 of file TSpectrumFit.h.

◆ GetAmplitudesErrors()

Double_t * TSpectrumFit::GetAmplitudesErrors ( ) const
inline

Definition at line 117 of file TSpectrumFit.h.

◆ GetAreas()

Double_t * TSpectrumFit::GetAreas ( ) const
inline

Definition at line 118 of file TSpectrumFit.h.

◆ GetAreasErrors()

Double_t * TSpectrumFit::GetAreasErrors ( ) const
inline

Definition at line 119 of file TSpectrumFit.h.

◆ GetBackgroundParameters()

void TSpectrumFit::GetBackgroundParameters ( Double_t & a0,
Double_t & a0Err,
Double_t & a1,
Double_t & a1Err,
Double_t & a2,
Double_t & a2Err )

This function gets the background parameters and their errors.

  • a0 - gets the fitted value of a0 parameter
  • a0Err - gets error value of a0 parameter
  • a1 - gets the fitted value of a1 parameter
  • a1Err - gets error value of a1 parameter
  • a2 - gets the fitted value of a2 parameter
  • a2Err - gets error value of a2 parameter

Definition at line 2741 of file TSpectrumFit.cxx.

◆ GetChi()

Double_t TSpectrumFit::GetChi ( ) const
inline

Definition at line 121 of file TSpectrumFit.h.

◆ GetDrawOption()

Option_t * TObject::GetDrawOption ( ) const
virtualinherited

Get option used by the graphics system to draw this object.

Note that before calling object.GetDrawOption(), you must have called object.Draw(..) before in the current pad.

Reimplemented in TBrowser, TFitEditor, TGedFrame, TGFileBrowser, TRootBrowser, and TRootBrowserLite.

Definition at line 445 of file TObject.cxx.

◆ GetDtorOnly()

Longptr_t TObject::GetDtorOnly ( )
staticinherited

Return destructor only flag.

Definition at line 1196 of file TObject.cxx.

◆ GetIconName()

const char * TObject::GetIconName ( ) const
virtualinherited

Returns mime type name of object.

Used by the TBrowser (via TGMimeTypes class). Override for class of which you would like to have different icons for objects of the same class.

Reimplemented in ROOT::Experimental::XRooFit::xRooNode, TASImage, TBranch, TBranchElement, TGeoVolume, TGMainFrame, TKey, TMethodBrowsable, TSystemFile, and TVirtualBranchBrowsable.

Definition at line 472 of file TObject.cxx.

◆ GetName()

const char * TNamed::GetName ( ) const
inlineoverridevirtualinherited

Returns name of object.

This default method returns the class name. Classes that give objects a name should override this method.

Reimplemented from TObject.

Definition at line 49 of file TNamed.h.

◆ GetObjectInfo()

char * TObject::GetObjectInfo ( Int_t px,
Int_t py ) const
virtualinherited

Returns string containing info about the object at position (px,py).

This method is typically overridden by classes of which the objects can report peculiarities for different positions. Returned string will be re-used (lock in MT environment).

Reimplemented in TASImage, TAxis3D, TColorWheel, TF1, TF2, TFileDrawMap, TGeoNode, TGeoTrack, TGeoVolume, TGL5DDataSet, TGLHistPainter, TGLParametricEquation, TGLTH3Composition, TGraph, TH1, THistPainter, TNode, TPaletteAxis, TParallelCoordVar, and TVirtualHistPainter.

Definition at line 491 of file TObject.cxx.

◆ GetObjectStat()

Bool_t TObject::GetObjectStat ( )
staticinherited

Get status of object stat flag.

Definition at line 1181 of file TObject.cxx.

◆ GetOption()

virtual Option_t * TObject::GetOption ( ) const
inlinevirtualinherited

◆ GetPositions()

Double_t * TSpectrumFit::GetPositions ( ) const
inline

Definition at line 122 of file TSpectrumFit.h.

◆ GetPositionsErrors()

Double_t * TSpectrumFit::GetPositionsErrors ( ) const
inline

Definition at line 123 of file TSpectrumFit.h.

◆ GetSigma()

void TSpectrumFit::GetSigma ( Double_t & sigma,
Double_t & sigmaErr )

This function gets the sigma parameter and its error.

  • sigma - gets the fitted value of sigma parameter
  • sigmaErr - gets error value of sigma parameter

Definition at line 2726 of file TSpectrumFit.cxx.

◆ GetTailParameters()

void TSpectrumFit::GetTailParameters ( Double_t & t,
Double_t & tErr,
Double_t & b,
Double_t & bErr,
Double_t & s,
Double_t & sErr )

This function gets the tail parameters and their errors.

  • t - gets the fitted value of t parameter
  • tErr - gets error value of t parameter
  • b - gets the fitted value of b parameter
  • bErr - gets error value of b parameter
  • s - gets the fitted value of s parameter
  • sErr - gets error value of s parameter

Definition at line 2761 of file TSpectrumFit.cxx.

◆ GetTitle()

const char * TNamed::GetTitle ( ) const
inlineoverridevirtualinherited

Returns title of object.

This default method returns the class title (i.e. description). Classes that give objects a title should override this method.

Reimplemented from TObject.

Definition at line 50 of file TNamed.h.

◆ GetUniqueID()

UInt_t TObject::GetUniqueID ( ) const
virtualinherited

Return the unique object id.

Definition at line 480 of file TObject.cxx.

◆ HandleTimer()

Bool_t TObject::HandleTimer ( TTimer * timer)
virtualinherited

Execute action in response of a timer timing out.

This method must be overridden if an object has to react to timers.

Reimplemented in TGCommandPlugin, TGDNDManager, TGFileContainer, TGHtml, TGLEventHandler, TGPopupMenu, TGraphTime, TGScrollBar, TGShutter, TGTextEdit, TGTextEditor, TGTextEntry, TGTextView, TGToolTip, TGuiBldDragManager, TGWindow, and TTreeViewer.

Definition at line 516 of file TObject.cxx.

◆ Hash()

ULong_t TNamed::Hash ( ) const
inlineoverridevirtualinherited

Return hash value for this object.

Note: If this routine is overloaded in a derived class, this derived class should also add

void CallRecursiveRemoveIfNeeded(TObject &obj)
call RecursiveRemove for obj if gROOT is valid and obj.TestBit(kMustCleanup) is true.
Definition TROOT.h:406

Otherwise, when RecursiveRemove is called (by ~TObject or example) for this type of object, the transversal of THashList and THashTable containers will will have to be done without call Hash (and hence be linear rather than logarithmic complexity). You will also see warnings like

Error in <ROOT::Internal::TCheckHashRecursiveRemoveConsistency::CheckRecursiveRemove>: The class SomeName overrides
TObject::Hash but does not call TROOT::RecursiveRemove in its destructor.
ULong_t Hash() const override
Return hash value for this object.
Definition TNamed.h:51
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
Definition TObject.cxx:1098
TObject()
TObject constructor.
Definition TObject.h:259
void RecursiveRemove(TObject *obj) override
Recursively remove this object from the list of Cleanups.
Definition TROOT.cxx:2651

Reimplemented from TObject.

Definition at line 51 of file TNamed.h.

◆ HasInconsistentHash()

Bool_t TObject::HasInconsistentHash ( ) const
inlineinherited

Return true is the type of this object is known to have an inconsistent setup for Hash and RecursiveRemove (i.e.

missing call to RecursiveRemove in destructor).

Note: Since the consistency is only tested for during inserts, this routine will return true for object that have never been inserted whether or not they have a consistent setup. This has no negative side-effect as searching for the object with the right or wrong Hash will always yield a not-found answer (Since anyway no hash can be guaranteed unique, there is always a check)

Definition at line 366 of file TObject.h.

◆ Info()

void TObject::Info ( const char * location,
const char * fmt,
... ) const
virtualinherited

Issue info message.

Use "location" to specify the method where the warning occurred. Accepts standard printf formatting arguments.

Definition at line 1072 of file TObject.cxx.

◆ InheritsFrom() [1/2]

Bool_t TObject::InheritsFrom ( const char * classname) const
virtualinherited

Returns kTRUE if object inherits from class "classname".

Reimplemented in TClass.

Definition at line 549 of file TObject.cxx.

◆ InheritsFrom() [2/2]

Bool_t TObject::InheritsFrom ( const TClass * cl) const
virtualinherited

Returns kTRUE if object inherits from TClass cl.

Reimplemented in TClass.

Definition at line 557 of file TObject.cxx.

◆ Inspect()

void TObject::Inspect ( ) const
virtualinherited

Dump contents of this object in a graphics canvas.

Same action as Dump but in a graphical form. In addition pointers to other objects can be followed.

The following picture is the Inspect of a histogram object:

Reimplemented in ROOT::Experimental::XRooFit::xRooNode, TGFrame, TInspectorObject, and TSystemFile.

Definition at line 570 of file TObject.cxx.

◆ InvertBit()

void TObject::InvertBit ( UInt_t f)
inlineinherited

Definition at line 206 of file TObject.h.

◆ IsA()

TClass * TSpectrumFit::IsA ( ) const
inlineoverridevirtual
Returns
TClass describing current object

Reimplemented from TNamed.

Definition at line 131 of file TSpectrumFit.h.

◆ IsDestructed()

Bool_t TObject::IsDestructed ( ) const
inlineinherited

IsDestructed.

Note
This function must be non-virtual as it can be used on destructed (but not yet modified) memory. This is used for example in TClonesArray to record the element that have been destructed but not deleted and thus are ready for re-use (by operator new with placement).
Returns
true if this object's destructor has been run.

Definition at line 186 of file TObject.h.

◆ IsEqual()

Bool_t TObject::IsEqual ( const TObject * obj) const
virtualinherited

Default equal comparison (objects are equal if they have the same address in memory).

More complicated classes might want to override this function.

Reimplemented in TGObject, TObjString, TPair, and TQCommand.

Definition at line 589 of file TObject.cxx.

◆ IsFolder()

◆ IsOnHeap()

Bool_t TObject::IsOnHeap ( ) const
inlineinherited

Definition at line 160 of file TObject.h.

◆ IsSortable()

Bool_t TNamed::IsSortable ( ) const
inlineoverridevirtualinherited

Reimplemented from TObject.

Reimplemented in TStructNodeProperty.

Definition at line 52 of file TNamed.h.

◆ IsZombie()

Bool_t TObject::IsZombie ( ) const
inlineinherited

Definition at line 161 of file TObject.h.

◆ ls()

void TNamed::ls ( Option_t * option = "") const
overridevirtualinherited

List TNamed name and title.

Reimplemented from TObject.

Reimplemented in ROOT::Experimental::XRooFit::xRooBrowser, TNode, TROOT, TStreamerBase, TStreamerElement, TStreamerInfo, TStreamerSTL, TTask, TText, and TVirtualStreamerInfo.

Definition at line 112 of file TNamed.cxx.

◆ MakeZombie()

void TObject::MakeZombie ( )
inlineprotectedinherited

Definition at line 55 of file TObject.h.

◆ MayNotUse()

void TObject::MayNotUse ( const char * method) const
inherited

Use this method to signal that a method (defined in a base class) may not be called in a derived class (in principle against good design since a child class should not provide less functionality than its parent, however, sometimes it is necessary).

Definition at line 1160 of file TObject.cxx.

◆ Notify()

Bool_t TObject::Notify ( )
virtualinherited

This method must be overridden to handle object notification (the base implementation is no-op).

Different objects in ROOT use the Notify method for different purposes, in coordination with other objects that call this method at the appropriate time.

For example, TLeaf uses it to load class information; TBranchRef to load contents of referenced branches TBranchRef; most notably, based on Notify, TChain implements a callback mechanism to inform interested parties when it switches to a new sub-tree.

Reimplemented in h1analysis, h1analysisTreeReader, TARInterruptHandler, TASInputHandler, TASInterruptHandler, TASLogHandler, TASSigPipeHandler, TBlinkTimer, TBranchElement, TBranchRef, TBreakLineCom, TBrowserTimer, TCollection, TDelCharCom, TDelTextCom, TFileHandler, TGContainerKeyboardTimer, TGContainerScrollTimer, TGInputHandler, TGLRedrawTimer, TGTextEditHist, TGuiBldDragManagerRepeatTimer, TIdleTimer, TInsCharCom, TInsTextCom, TInterruptHandler, TLeafObject, TMessageHandler, TNotifyLink< Type >, TNotifyLink< RNoCleanupNotifierHelper >, TNotifyLink< ROOT::Detail::TBranchProxy >, TNotifyLink< TTreeReader >, TPopupDelayTimer, TProcessEventTimer, TRefTable, TRepeatTimer, TSBRepeatTimer, TSelector, TSelectorDraw, TSelectorEntries, TSignalHandler, TSingleShotCleaner, TSocketHandler, TStdExceptionHandler, TSysEvtHandler, TTermInputHandler, TThreadTimer, TTimeOutTimer, TTimer, TTipDelayTimer, TTree, TTreeFormula, TTreeFormulaManager, TTreeReader, TViewTimer, and TViewUpdateTimer.

Definition at line 618 of file TObject.cxx.

◆ Obsolete()

void TObject::Obsolete ( const char * method,
const char * asOfVers,
const char * removedFromVers ) const
inherited

Use this method to declare a method obsolete.

Specify as of which version the method is obsolete and as from which version it will be removed.

Definition at line 1169 of file TObject.cxx.

◆ operator delete() [1/3]

void TObject::operator delete ( void * ptr,
size_t size )
inherited

Operator delete for sized deallocation.

Definition at line 1234 of file TObject.cxx.

◆ operator delete() [2/3]

void TObject::operator delete ( void * ptr)
inherited

Operator delete.

Definition at line 1212 of file TObject.cxx.

◆ operator delete() [3/3]

void TObject::operator delete ( void * ptr,
void * vp )
inherited

Only called by placement new when throwing an exception.

Definition at line 1266 of file TObject.cxx.

◆ operator delete[]() [1/3]

void TObject::operator delete[] ( void * ptr,
size_t size )
inherited

Operator delete [] for sized deallocation.

Definition at line 1245 of file TObject.cxx.

◆ operator delete[]() [2/3]

void TObject::operator delete[] ( void * ptr)
inherited

Operator delete [].

Definition at line 1223 of file TObject.cxx.

◆ operator delete[]() [3/3]

void TObject::operator delete[] ( void * ptr,
void * vp )
inherited

Only called by placement new[] when throwing an exception.

Definition at line 1274 of file TObject.cxx.

◆ operator new() [1/2]

void * TObject::operator new ( size_t sz)
inlineinherited

Definition at line 189 of file TObject.h.

◆ operator new() [2/2]

void * TObject::operator new ( size_t sz,
void * vp )
inlineinherited

Definition at line 191 of file TObject.h.

◆ operator new[]() [1/2]

void * TObject::operator new[] ( size_t sz)
inlineinherited

Definition at line 190 of file TObject.h.

◆ operator new[]() [2/2]

void * TObject::operator new[] ( size_t sz,
void * vp )
inlineinherited

Definition at line 192 of file TObject.h.

◆ Ourpowl()

Double_t TSpectrumFit::Ourpowl ( Double_t a,
Int_t pw )
protected

Power function.

Definition at line 675 of file TSpectrumFit.cxx.

◆ Paint()

void TObject::Paint ( Option_t * option = "")
virtualinherited

This method must be overridden if a class wants to paint itself.

The difference between Paint() and Draw() is that when a object draws itself it is added to the display list of the pad in which it is drawn (and automatically redrawn whenever the pad is redrawn). While paint just draws the object without adding it to the pad display list.

Reimplemented in ROOT::Experimental::RTreeMapPainter, ROOT::RGeoPainter, TAnnotation, TArrow, TASImage, TASPaletteEditor::LimitLine, TASPaletteEditor::PaintPalette, TAxis3D, TBits, TBox, TButton, TCanvas, TClassTree, TCollection, TColorWheel, TCrown, TDiamond, TDirectory, TEfficiency, TEllipse, TEveArrow, TEveCaloViz, TEveDigitSet, TEveGeoShape, TEveGeoTopNode, TEvePlot3D, TEvePointSet, TEveProjectionAxes, TEveScene, TEveShape, TEveStraightLineSet, TEveText, TEveTriangleSet, TExec, TF1, TF2, TF3, TFile, TFileDrawMap, TFrame, TGaxis, TGenerator, TGeoBoolNode, TGeoIntersection, TGeoNode, TGeoOverlap, TGeoPainter, TGeoPhysicalNode, TGeoShape, TGeoSubtraction, TGeoTrack, TGeoUnion, TGeoVGShape, TGeoVolume, TGL5DDataSet, TGLHistPainter, TGLParametricEquation, TGLTH3Composition, TGraph2D, TGraph2DPainter, TGraph, TGraphEdge, TGraphNode, TGraphPolargram, TGraphTime, TH1, THistPainter, THStack, TLatex, TLegend, TLine, TMacro, TMarker3DBox, TMarker, TMathText, TMultiGraph, TNode, TNodeDiv, TPad, TPaletteAxis, TParallelCoord, TParallelCoordRange, TParallelCoordVar, TParticle, TPave, TPaveLabel, TPaveStats, TPavesText, TPaveText, TPie, TPolyLine3D, TPolyLine, TPolyMarker3D, TPolyMarker, TPrimary, TRatioPlot, TScatter2D, TScatter, TShape, TSpectrum2Painter, TSpider, TSpline, TSQLFile, TStyle, TText, TTreePerfStats, TVirtualGeoPainter, TVirtualGeoTrack, TVirtualHistPainter, TVirtualPad, TWbox, and TXMLFile.

Definition at line 631 of file TObject.cxx.

◆ Pop()

void TObject::Pop ( )
virtualinherited

Pop on object drawn in a pad to the top of the display list.

I.e. it will be drawn last and on top of all other primitives.

Reimplemented in TFrame, TPad, and TVirtualPad.

Definition at line 640 of file TObject.cxx.

◆ Print()

◆ Read()

Int_t TObject::Read ( const char * name)
virtualinherited

Read contents of object with specified name from the current directory.

First the key with the given name is searched in the current directory, next the key buffer is deserialized into the object. The object must have been created before via the default constructor. See TObject::Write().

Reimplemented in TBuffer, TKey, TKeySQL, and TKeyXML.

Definition at line 673 of file TObject.cxx.

◆ RecursiveRemove()

◆ ResetBit()

void TObject::ResetBit ( UInt_t f)
inlineinherited

Definition at line 203 of file TObject.h.

◆ SaveAs()

void TObject::SaveAs ( const char * filename = "",
Option_t * option = "" ) const
virtualinherited

Save this object in the file specified by filename.

  • if "filename" contains ".root" the object is saved in filename as root binary file.
  • if "filename" contains ".xml" the object is saved in filename as a xml ascii file.
  • if "filename" contains ".cc" the object is saved in filename as C code independent from ROOT. The code is generated via SavePrimitive(). Specific code should be implemented in each object to handle this option. Like in TF1::SavePrimitive().
  • otherwise the object is written to filename as a CINT/C++ script. The C++ code to rebuild this object is generated via SavePrimitive(). The "option" parameter is passed to SavePrimitive. By default it is an empty string. It can be used to specify the Draw option in the code generated by SavePrimitive.

    The function is available via the object context menu.

Reimplemented in ROOT::Experimental::XRooFit::xRooNode, TClassTree, TFolder, TGeoVolume, TGObject, TGraph, TH1, TPad, TPaveClass, TSpline3, TSpline5, TSpline, TTreePerfStats, and TVirtualPad.

Definition at line 708 of file TObject.cxx.

◆ SavePrimitive()

void TObject::SavePrimitive ( std::ostream & out,
Option_t * option = "" )
virtualinherited

Save a primitive as a C++ statement(s) on output stream "out".

Reimplemented in TAnnotation, TArc, TArrow, TASImage, TAxis3D, TBox, TButton, TCanvas, TChain, TCrown, TCurlyArc, TCurlyLine, TCutG, TDiamond, TEfficiency, TEllipse, TExec, TF12, TF1, TF2, TF3, TFrame, TGaxis, TGButton, TGButtonGroup, TGCanvas, TGCheckButton, TGColorSelect, TGColumnLayout, TGComboBox, TGCompositeFrame, TGContainer, TGDockableFrame, TGDoubleHSlider, TGDoubleVSlider, TGedMarkerSelect, TGedPatternSelect, TGeoArb8, TGeoBBox, TGeoBoolNode, TGeoCombiTrans, TGeoCompositeShape, TGeoCone, TGeoConeSeg, TGeoCtub, TGeoDecayChannel, TGeoElementRN, TGeoEltu, TGeoGtra, TGeoHalfSpace, TGeoHMatrix, TGeoHype, TGeoIdentity, TGeoIntersection, TGeoMaterial, TGeoMedium, TGeoMixture, TGeoPara, TGeoParaboloid, TGeoPatternCylPhi, TGeoPatternCylR, TGeoPatternParaX, TGeoPatternParaY, TGeoPatternParaZ, TGeoPatternSphPhi, TGeoPatternSphR, TGeoPatternSphTheta, TGeoPatternTrapZ, TGeoPatternX, TGeoPatternY, TGeoPatternZ, TGeoPcon, TGeoPgon, TGeoRotation, TGeoScaledShape, TGeoShapeAssembly, TGeoSphere, TGeoSubtraction, TGeoTessellated, TGeoTorus, TGeoTranslation, TGeoTrap, TGeoTrd1, TGeoTrd2, TGeoTube, TGeoTubeSeg, TGeoUnion, TGeoVolume, TGeoXtru, TGFileContainer, TGFont, TGFrame, TGFSComboBox, TGGC, TGGroupFrame, TGHButtonGroup, TGHorizontal3DLine, TGHorizontalFrame, TGHorizontalLayout, TGHProgressBar, TGHScrollBar, TGHSlider, TGHSplitter, TGHtml, TGIcon, TGLabel, TGLayoutHints, TGLineStyleComboBox, TGLineWidthComboBox, TGListBox, TGListDetailsLayout, TGListLayout, TGListTree, TGListView, TGLVContainer, TGMainFrame, TGMatrixLayout, TGMdiFrame, TGMdiMainFrame, TGMdiMenuBar, TGMenuBar, TGMenuTitle, TGNumberEntry, TGNumberEntryField, TGPictureButton, TGPopupMenu, TGProgressBar, TGRadioButton, TGraph2D, TGraph2DAsymmErrors, TGraph2DErrors, TGraph, TGraphAsymmErrors, TGraphBentErrors, TGraphEdge, TGraphErrors, TGraphMultiErrors, TGraphNode, TGraphPolar, TGraphPolargram, TGraphStruct, TGroupButton, TGRowLayout, TGShapedFrame, TGShutter, TGShutterItem, TGSplitFrame, TGStatusBar, TGTab, TGTabLayout, TGTableLayout, TGTableLayoutHints, TGTextButton, TGTextEdit, TGTextEntry, TGTextLBEntry, TGTextView, TGTileLayout, TGToolBar, TGTransientFrame, TGTripleHSlider, TGTripleVSlider, TGVButtonGroup, TGVertical3DLine, TGVerticalFrame, TGVerticalLayout, TGVFileSplitter, TGVProgressBar, TGVScrollBar, TGVSlider, TGVSplitter, TGXYLayout, TGXYLayoutHints, TH1, TH2Poly, THelix, THStack, TLatex, TLegend, TLine, TMacro, TMarker3DBox, TMarker, TMathText, TMultiGraph, TPad, TPaletteAxis, TParallelCoord, TParallelCoordVar, TPave, TPaveClass, TPaveLabel, TPaveStats, TPavesText, TPaveText, TPie, TPieSlice, TPolyLine3D, TPolyLine, TPolyMarker3D, TPolyMarker, TProfile2D, TProfile3D, TProfile, TRootContainer, TRootEmbeddedCanvas, TScatter2D, TScatter, TSlider, TSliderBox, TSpline3, TSpline5, TStyle, TText, TTreePerfStats, and TWbox.

Definition at line 858 of file TObject.cxx.

◆ SavePrimitiveConstructor()

void TObject::SavePrimitiveConstructor ( std::ostream & out,
TClass * cl,
const char * variable_name,
const char * constructor_agrs = "",
Bool_t empty_line = kTRUE )
staticprotectedinherited

Save object constructor in the output stream "out".

Can be used as first statement when implementing SavePrimitive() method for the object

Definition at line 777 of file TObject.cxx.

◆ SavePrimitiveDraw()

void TObject::SavePrimitiveDraw ( std::ostream & out,
const char * variable_name,
Option_t * option = nullptr )
staticprotectedinherited

Save invocation of primitive Draw() method Skipped if option contains "nodraw" string.

Definition at line 845 of file TObject.cxx.

◆ SavePrimitiveNameTitle()

void TNamed::SavePrimitiveNameTitle ( std::ostream & out,
const char * variable_name )
protectedinherited

Save object name and title into the output stream "out".

Definition at line 135 of file TNamed.cxx.

◆ SavePrimitiveVector()

TString TObject::SavePrimitiveVector ( std::ostream & out,
const char * prefix,
Int_t len,
Double_t * arr,
Int_t flag = 0 )
staticprotectedinherited

Save array in the output stream "out" as vector.

Create unique variable name based on prefix value Returns name of vector which can be used in constructor or in other places of C++ code If flag === kTRUE, just add empty line If flag === 111, check if array is empty and return nullptr or <vectorname>.data()

Definition at line 796 of file TObject.cxx.

◆ SetBackgroundParameters()

void TSpectrumFit::SetBackgroundParameters ( Double_t a0Init,
Bool_t fixA0,
Double_t a1Init,
Bool_t fixA1,
Double_t a2Init,
Bool_t fixA2 )

This function sets the following fitting parameters of background:

  • a0Init - initial value of a0 parameter (background is estimated as a0+a1*x+a2*x*x)
  • fixA0 - logical value of a0 parameter, which allows to fix the parameter (not to fit)
  • a1Init - initial value of a1 parameter
  • fixA1 - logical value of a1 parameter, which allows to fix the parameter (not to fit)
  • a2Init - initial value of a2 parameter
  • fixA2 - logical value of a2 parameter, which allows to fix the parameter (not to fit)

Definition at line 2692 of file TSpectrumFit.cxx.

◆ SetBit() [1/2]

void TObject::SetBit ( UInt_t f)
inlineinherited

Definition at line 202 of file TObject.h.

◆ SetBit() [2/2]

void TObject::SetBit ( UInt_t f,
Bool_t set )
inherited

Set or unset the user status bits as specified in f.

Definition at line 888 of file TObject.cxx.

◆ SetDrawOption()

void TObject::SetDrawOption ( Option_t * option = "")
virtualinherited

Set drawing option for object.

This option only affects the drawing style and is stored in the option field of the TObjOptLink supporting a TPad's primitive list (TList). Note that it does not make sense to call object.SetDrawOption(option) before having called object.Draw().

Reimplemented in RooPlot, TAxis, TBrowser, TGedFrame, TGFrame, TPad, TPaveStats, TRootBrowserLite, TSystemDirectory, and TSystemFile.

Definition at line 871 of file TObject.cxx.

◆ SetDtorOnly()

void TObject::SetDtorOnly ( void * obj)
staticinherited

Set destructor only flag.

Definition at line 1204 of file TObject.cxx.

◆ SetFitParameters()

void TSpectrumFit::SetFitParameters ( Int_t xmin,
Int_t xmax,
Int_t numberIterations,
Double_t alpha,
Int_t statisticType,
Int_t alphaOptim,
Int_t power,
Int_t fitTaylor )

This function sets the following fitting parameters:

  • xmin, xmax - fitting region
  • numberIterations - # of desired iterations in the fit
  • alpha - convergence coefficient, it should be positive number and <=1, for details see references
  • statisticType - type of statistics, possible values kFitOptimChiCounts (chi square statistics with counts as weighting coefficients), kFitOptimChiFuncValues (chi square statistics with function values as weighting coefficients),kFitOptimMaxLikelihood
  • alphaOptim - optimization of convergence algorithm, possible values kFitAlphaHalving, kFitAlphaOptimal
  • power - possible values kFitPower2,4,6,8,10,12, for details see references. It applies only for Awmi fitting function.
  • fitTaylor - order of Taylor expansion, possible values kFitTaylorOrderFirst, kFitTaylorOrderSecond. It applies only for Awmi fitting function.

Definition at line 2609 of file TSpectrumFit.cxx.

◆ SetName()

void TNamed::SetName ( const char * name)
virtualinherited

Set the name of the TNamed.

WARNING: if the object is a member of a THashTable or THashList container the container must be Rehash()'ed after SetName(). For example the list of objects in the current directory is a THashList.

Reimplemented in RooAbsArg, RooAbsData, RooDataHist, RooDataSet, RooFitResult, RooPlot, ROOT::Experimental::XRooFit::xRooNode, TChain, TColor, TDirectory, TEfficiency, TEventList, TEveScene, TFormula, TGraph2D, TGraph, TH1, TNode, TRotMatrix, TShape, TSystemDirectory, TSystemFile, and TTree.

Definition at line 149 of file TNamed.cxx.

◆ SetNameTitle()

void TNamed::SetNameTitle ( const char * name,
const char * title )
virtualinherited

Set all the TNamed parameters (name and title).

WARNING: if the name is changed and the object is a member of a THashTable or THashList container the container must be Rehash()'ed after SetName(). For example the list of objects in the current directory is a THashList.

Reimplemented in RooAbsArg, RooAbsData, RooDataHist, RooDataSet, RooFitResult, RooPlot, TContextMenu, TGraph2D, TGraph, TH1, and TNode.

Definition at line 163 of file TNamed.cxx.

◆ SetObjectStat()

void TObject::SetObjectStat ( Bool_t stat)
staticinherited

Turn on/off tracking of objects in the TObjectTable.

Definition at line 1188 of file TObject.cxx.

◆ SetPeakParameters()

void TSpectrumFit::SetPeakParameters ( Double_t sigma,
Bool_t fixSigma,
const Double_t * positionInit,
const Bool_t * fixPosition,
const Double_t * ampInit,
const Bool_t * fixAmp )

This function sets the following fitting parameters of peaks:

  • sigma - initial value of sigma parameter
  • fixSigma - logical value of sigma parameter, which allows to fix the parameter (not to fit)
  • positionInit - array of initial values of peaks positions
  • fixPosition - array of logical values which allow to fix appropriate positions (not fit). However they are present in the estimated functional.
  • ampInit - array of initial values of peaks amplitudes
  • fixAmp - array of logical values which allow to fix appropriate amplitudes (not fit). However they are present in the estimated functional

Definition at line 2657 of file TSpectrumFit.cxx.

◆ SetTailParameters()

void TSpectrumFit::SetTailParameters ( Double_t tInit,
Bool_t fixT,
Double_t bInit,
Bool_t fixB,
Double_t sInit,
Bool_t fixS )

This function sets the following fitting parameters of tails of peaks.

  • tInit - initial value of t parameter
  • fixT - logical value of t parameter, which allows to fix the parameter (not to fit)
  • bInit - initial value of b parameter
  • fixB - logical value of b parameter, which allows to fix the parameter (not to fit)
  • sInit - initial value of s parameter
  • fixS - logical value of s parameter, which allows to fix the parameter (not to fit)

Definition at line 2711 of file TSpectrumFit.cxx.

◆ SetTitle()

void TNamed::SetTitle ( const char * title = "")
virtualinherited

◆ SetUniqueID()

void TObject::SetUniqueID ( UInt_t uid)
virtualinherited

Set the unique object id.

Definition at line 899 of file TObject.cxx.

◆ Shape()

Double_t TSpectrumFit::Shape ( Int_t num_of_fitted_peaks,
Double_t i,
const Double_t * parameter,
Double_t sigma,
Double_t t,
Double_t s,
Double_t b,
Double_t a0,
Double_t a1,
Double_t a2 )
protected

This function calculates peaks shape function (see manual) Function parameters:

  • num_of_fitted_peaks-number of fitted peaks
  • i-channel
  • parameter-array of peaks parameters (amplitudes and positions)
  • sigma-sigma of peak
  • t, s-relative amplitudes
  • b-slope
  • a0, a1, a2- background coefficients

Definition at line 517 of file TSpectrumFit.cxx.

◆ Sizeof()

Int_t TNamed::Sizeof ( ) const
virtualinherited

Return size of the TNamed part of the TObject.

Reimplemented in TDirectory, TDirectoryFile, TFile, TKey, TSQLFile, and TXMLFile.

Definition at line 182 of file TNamed.cxx.

◆ StiefelInversion()

void TSpectrumFit::StiefelInversion ( Double_t ** a,
Int_t size )
protected

This function calculates solution of the system of linear equations.

The matrix a should have a dimension size*(size+4) The calling function should fill in the matrix, the column size should contain vector y (right side of the system of equations). The result is placed into size+1 column of the matrix. according to sigma of peaks.

Function parameters:

  • a-matrix with dimension size*(size+4)
  • size-number of rows of the matrix

Definition at line 1722 of file TSpectrumFit.cxx.

◆ Streamer()

void TSpectrumFit::Streamer ( TBuffer & R__b)
overridevirtual

Stream an object of class TObject.

Reimplemented from TNamed.

◆ StreamerNVirtual()

void TSpectrumFit::StreamerNVirtual ( TBuffer & ClassDef_StreamerNVirtual_b)
inline

Definition at line 131 of file TSpectrumFit.h.

◆ SysError()

void TObject::SysError ( const char * location,
const char * fmt,
... ) const
virtualinherited

Issue system error message.

Use "location" to specify the method where the system error occurred. Accepts standard printf formatting arguments.

Definition at line 1112 of file TObject.cxx.

◆ TestBit()

Bool_t TObject::TestBit ( UInt_t f) const
inlineinherited

Definition at line 204 of file TObject.h.

◆ TestBits()

Int_t TObject::TestBits ( UInt_t f) const
inlineinherited

Definition at line 205 of file TObject.h.

◆ UseCurrentStyle()

void TObject::UseCurrentStyle ( )
virtualinherited

Set current style settings in this object This function is called when either TCanvas::UseCurrentStyle or TROOT::ForceStyle have been invoked.

Reimplemented in TAxis3D, TCanvas, TFrame, TGraph, TH1, TPad, TPaveStats, TPaveText, and TTree.

Definition at line 909 of file TObject.cxx.

◆ Warning()

void TObject::Warning ( const char * location,
const char * fmt,
... ) const
virtualinherited

Issue warning message.

Use "location" to specify the method where the warning occurred. Accepts standard printf formatting arguments.

Definition at line 1084 of file TObject.cxx.

◆ Write() [1/2]

Int_t TObject::Write ( const char * name = nullptr,
Int_t option = 0,
Int_t bufsize = 0 )
virtualinherited

Write this object to the current directory.

For more see the const version of this method.

Reimplemented in ROOT::TBufferMergerFile, TBuffer, TCollection, TDirectory, TDirectoryFile, TFile, TMap, TParallelMergingFile, TSQLFile, TTree, and TXMLFile.

Definition at line 989 of file TObject.cxx.

◆ Write() [2/2]

Int_t TObject::Write ( const char * name = nullptr,
Int_t option = 0,
Int_t bufsize = 0 ) const
virtualinherited

Write this object to the current directory.

The data structure corresponding to this object is serialized. The corresponding buffer is written to the current directory with an associated key with name "name".

Writing an object to a file involves the following steps:

  • Creation of a support TKey object in the current directory. The TKey object creates a TBuffer object.
  • The TBuffer object is filled via the class::Streamer function.
  • If the file is compressed (default) a second buffer is created to hold the compressed buffer.
  • Reservation of the corresponding space in the file by looking in the TFree list of free blocks of the file.
  • The buffer is written to the file.

Bufsize can be given to force a given buffer size to write this object. By default, the buffersize will be taken from the average buffer size of all objects written to the current file so far.

If a name is specified, it will be the name of the key. If name is not given, the name of the key will be the name as returned by GetName().

The option can be a combination of: kSingleKey, kOverwrite or kWriteDelete Using the kOverwrite option a previous key with the same name is overwritten. The previous key is deleted before writing the new object. Using the kWriteDelete option a previous key with the same name is deleted only after the new object has been written. This option is safer than kOverwrite but it is slower. NOTE: Neither kOverwrite nor kWriteDelete reduces the size of a TFile– the space is simply freed up to be overwritten; in the case of a TTree, it is more complicated. If one opens a TTree, appends some entries, then writes it out, the behaviour is effectively the same. If, however, one creates a new TTree and writes it out in this way, only the metadata is replaced, effectively making the old data invisible without deleting it. TTree::Delete() can be used to mark all disk space occupied by a TTree as free before overwriting its metadata this way. The kSingleKey option is only used by TCollection::Write() to write a container with a single key instead of each object in the container with its own key.

An object is read from the file into memory via TKey::Read() or via TObject::Read().

The function returns the total number of bytes written to the file. It returns 0 if the object cannot be written.

Reimplemented in TBuffer, TCollection, TDirectory, TDirectoryFile, TFile, TMap, TParallelMergingFile, TSQLFile, TTree, and TXMLFile.

Definition at line 964 of file TObject.cxx.

Member Data Documentation

◆ fA0Calc

Double_t TSpectrumFit::fA0Calc
protected

calculated value of background a0 parameter

Definition at line 51 of file TSpectrumFit.h.

◆ fA0Err

Double_t TSpectrumFit::fA0Err
protected

error value of background a0 parameter

Definition at line 52 of file TSpectrumFit.h.

◆ fA0Init

Double_t TSpectrumFit::fA0Init
protected

initial value of background a0 parameter(backgroud is estimated as a0+a1*x+a2*x*x)

Definition at line 50 of file TSpectrumFit.h.

◆ fA1Calc

Double_t TSpectrumFit::fA1Calc
protected

calculated value of background a1 parameter

Definition at line 54 of file TSpectrumFit.h.

◆ fA1Err

Double_t TSpectrumFit::fA1Err
protected

error value of background a1 parameter

Definition at line 55 of file TSpectrumFit.h.

◆ fA1Init

Double_t TSpectrumFit::fA1Init
protected

initial value of background a1 parameter(backgroud is estimated as a0+a1*x+a2*x*x)

Definition at line 53 of file TSpectrumFit.h.

◆ fA2Calc

Double_t TSpectrumFit::fA2Calc
protected

calculated value of background a2 parameter

Definition at line 57 of file TSpectrumFit.h.

◆ fA2Err

Double_t TSpectrumFit::fA2Err
protected

error value of background a2 parameter

Definition at line 58 of file TSpectrumFit.h.

◆ fA2Init

Double_t TSpectrumFit::fA2Init
protected

initial value of background a2 parameter(backgroud is estimated as a0+a1*x+a2*x*x)

Definition at line 56 of file TSpectrumFit.h.

◆ fAlpha

Double_t TSpectrumFit::fAlpha
protected

convergence coefficient, input parameter, it should be positive number and <=1, for details see references

Definition at line 28 of file TSpectrumFit.h.

◆ fAlphaOptim

Int_t TSpectrumFit::fAlphaOptim
protected

optimization of convergence algorithm, possible values kFitAlphaHalving, kFitAlphaOptimal

Definition at line 25 of file TSpectrumFit.h.

◆ fAmpCalc

Double_t* TSpectrumFit::fAmpCalc
protected

[fNPeaks] array of calculated values of fitted amplitudes, output parameters

Definition at line 34 of file TSpectrumFit.h.

◆ fAmpErr

Double_t* TSpectrumFit::fAmpErr
protected

[fNPeaks] array of amplitude errors

Definition at line 35 of file TSpectrumFit.h.

◆ fAmpInit

Double_t* TSpectrumFit::fAmpInit
protected

[fNPeaks] array of initial values of peaks amplitudes, input parameters

Definition at line 33 of file TSpectrumFit.h.

◆ fArea

Double_t* TSpectrumFit::fArea
protected

[fNPeaks] array of calculated areas of peaks

Definition at line 36 of file TSpectrumFit.h.

◆ fAreaErr

Double_t* TSpectrumFit::fAreaErr
protected

[fNPeaks] array of errors of peak areas

Definition at line 37 of file TSpectrumFit.h.

◆ fBCalc

Double_t TSpectrumFit::fBCalc
protected

calculated value of b parameter

Definition at line 45 of file TSpectrumFit.h.

◆ fBErr

Double_t TSpectrumFit::fBErr
protected

error value of b parameter

Definition at line 46 of file TSpectrumFit.h.

◆ fBInit

Double_t TSpectrumFit::fBInit
protected

initial value of b parameter (slope), for details see html manual and references

Definition at line 44 of file TSpectrumFit.h.

◆ fBits

UInt_t TObject::fBits
privateinherited

bit field status word

Definition at line 47 of file TObject.h.

◆ fChi

Double_t TSpectrumFit::fChi
protected

here the fitting functions return resulting chi square

Definition at line 29 of file TSpectrumFit.h.

◆ fFitTaylor

Int_t TSpectrumFit::fFitTaylor
protected

order of Taylor expansion, possible values kFitTaylorOrderFirst, kFitTaylorOrderSecond. It applies only for Awmi fitting function.

Definition at line 27 of file TSpectrumFit.h.

◆ fFixA0

Bool_t TSpectrumFit::fFixA0
protected

logical value of a0 parameter, which allows to fix the parameter (not to fit).

Definition at line 65 of file TSpectrumFit.h.

◆ fFixA1

Bool_t TSpectrumFit::fFixA1
protected

logical value of a1 parameter, which allows to fix the parameter (not to fit).

Definition at line 66 of file TSpectrumFit.h.

◆ fFixA2

Bool_t TSpectrumFit::fFixA2
protected

logical value of a2 parameter, which allows to fix the parameter (not to fit).

Definition at line 67 of file TSpectrumFit.h.

◆ fFixAmp

Bool_t* TSpectrumFit::fFixAmp
protected

[fNPeaks] array of logical values which allow to fix appropriate amplitudes (not fit). However they are present in the estimated functional

Definition at line 60 of file TSpectrumFit.h.

◆ fFixB

Bool_t TSpectrumFit::fFixB
protected

logical value of b parameter, which allows to fix the parameter (not to fit).

Definition at line 63 of file TSpectrumFit.h.

◆ fFixPosition

Bool_t* TSpectrumFit::fFixPosition
protected

[fNPeaks] array of logical values which allow to fix appropriate positions (not fit). However they are present in the estimated functional

Definition at line 59 of file TSpectrumFit.h.

◆ fFixS

Bool_t TSpectrumFit::fFixS
protected

logical value of s parameter, which allows to fix the parameter (not to fit).

Definition at line 64 of file TSpectrumFit.h.

◆ fFixSigma

Bool_t TSpectrumFit::fFixSigma
protected

logical value of sigma parameter, which allows to fix the parameter (not to fit).

Definition at line 61 of file TSpectrumFit.h.

◆ fFixT

Bool_t TSpectrumFit::fFixT
protected

logical value of t parameter, which allows to fix the parameter (not to fit).

Definition at line 62 of file TSpectrumFit.h.

◆ fgDtorOnly

Longptr_t TObject::fgDtorOnly = 0
staticprivateinherited

object for which to call dtor only (i.e. no delete)

Definition at line 49 of file TObject.h.

◆ fgObjectStat

Bool_t TObject::fgObjectStat = kTRUE
staticprivateinherited

if true keep track of objects in TObjectTable

Definition at line 50 of file TObject.h.

◆ fName

TString TNamed::fName
protectedinherited

Definition at line 32 of file TNamed.h.

◆ fNPeaks

Int_t TSpectrumFit::fNPeaks
protected

number of peaks present in fit, input parameter, it should be > 0

Definition at line 20 of file TSpectrumFit.h.

◆ fNumberIterations

Int_t TSpectrumFit::fNumberIterations
protected

number of iterations in fitting procedure, input parameter, it should be > 0

Definition at line 21 of file TSpectrumFit.h.

◆ fPositionCalc

Double_t* TSpectrumFit::fPositionCalc
protected

[fNPeaks] array of calculated values of fitted positions, output parameters

Definition at line 31 of file TSpectrumFit.h.

◆ fPositionErr

Double_t* TSpectrumFit::fPositionErr
protected

[fNPeaks] array of position errors

Definition at line 32 of file TSpectrumFit.h.

◆ fPositionInit

Double_t* TSpectrumFit::fPositionInit
protected

[fNPeaks] array of initial values of peaks positions, input parameters

Definition at line 30 of file TSpectrumFit.h.

◆ fPower

Int_t TSpectrumFit::fPower
protected

possible values kFitPower2,4,6,8,10,12, for details see references. It applies only for Awmi fitting function.

Definition at line 26 of file TSpectrumFit.h.

◆ fSCalc

Double_t TSpectrumFit::fSCalc
protected

calculated value of s parameter

Definition at line 48 of file TSpectrumFit.h.

◆ fSErr

Double_t TSpectrumFit::fSErr
protected

error value of s parameter

Definition at line 49 of file TSpectrumFit.h.

◆ fSigmaCalc

Double_t TSpectrumFit::fSigmaCalc
protected

calculated value of sigma parameter

Definition at line 39 of file TSpectrumFit.h.

◆ fSigmaErr

Double_t TSpectrumFit::fSigmaErr
protected

error value of sigma parameter

Definition at line 40 of file TSpectrumFit.h.

◆ fSigmaInit

Double_t TSpectrumFit::fSigmaInit
protected

initial value of sigma parameter

Definition at line 38 of file TSpectrumFit.h.

◆ fSInit

Double_t TSpectrumFit::fSInit
protected

initial value of s parameter (relative amplitude of step), for details see html manual and references

Definition at line 47 of file TSpectrumFit.h.

◆ fStatisticType

Int_t TSpectrumFit::fStatisticType
protected

type of statistics, possible values kFitOptimChiCounts (chi square statistics with counts as weighting coefficients), kFitOptimChiFuncValues (chi square statistics with function values as weighting coefficients),kFitOptimMaxLikelihood

Definition at line 24 of file TSpectrumFit.h.

◆ fTCalc

Double_t TSpectrumFit::fTCalc
protected

calculated value of t parameter

Definition at line 42 of file TSpectrumFit.h.

◆ fTErr

Double_t TSpectrumFit::fTErr
protected

error value of t parameter

Definition at line 43 of file TSpectrumFit.h.

◆ fTInit

Double_t TSpectrumFit::fTInit
protected

initial value of t parameter (relative amplitude of tail), for details see html manual and references

Definition at line 41 of file TSpectrumFit.h.

◆ fTitle

TString TNamed::fTitle
protectedinherited

Definition at line 33 of file TNamed.h.

◆ fUniqueID

UInt_t TObject::fUniqueID
privateinherited

object unique identifier

Definition at line 46 of file TObject.h.

◆ fXmax

Int_t TSpectrumFit::fXmax
protected

last fitted channel

Definition at line 23 of file TSpectrumFit.h.

◆ fXmin

Int_t TSpectrumFit::fXmin
protected

first fitted channel

Definition at line 22 of file TSpectrumFit.h.


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