#include <PRunSingleHisto.h>

Inheritance diagram for PRunSingleHisto:

Collaboration diagram for PRunSingleHisto:

[legend]

Public Member Functions
	PRunSingleHisto ()

	PRunSingleHisto (PMsrHandler msrInfo, PRunDataHandler rawData, UInt_t runNo, EPMusrHandleTag tag, Bool_t theoAsData)

virtual	~PRunSingleHisto ()

virtual Double_t	CalcChiSquare (const std::vector< Double_t > &par)

virtual Double_t	CalcChiSquareExpected (const std::vector< Double_t > &par)

virtual Double_t	CalcMaxLikelihood (const std::vector< Double_t > &par)

virtual Double_t	CalcMaxLikelihoodExpected (const std::vector< Double_t > &par)

virtual void	CalcTheory ()

virtual UInt_t	GetNoOfFitBins ()

virtual void	SetFitRangeBin (const TString fitRange)

virtual Double_t	GetBackground ()

virtual Int_t	GetStartTimeBin ()

virtual Int_t	GetEndTimeBin ()

virtual Int_t	GetPacking ()

virtual Bool_t	GetScaleN0AndBkg ()

virtual void	CalcNoOfFitBins ()

Public Member Functions inherited from PRunBase
	PRunBase ()

	PRunBase (PMsrHandler msrInfo, PRunDataHandler rawData, UInt_t runNo, EPMusrHandleTag tag)

virtual	~PRunBase ()

virtual void	SetFitRange (PDoublePairVector fitRange)

virtual UInt_t	GetRunNo ()
	returns the number of runs of the msr-file More...

virtual PRunData *	GetData ()
	returns the data to be fitted More...

virtual void	CleanUp ()

virtual Bool_t	IsValid ()
	returns if the state is valid More...

Protected Member Functions
virtual Bool_t	PrepareData ()

virtual Bool_t	PrepareFitData (PRawRunData *runData, const UInt_t histoNo)

virtual Bool_t	PrepareRawViewData (PRawRunData *runData, const UInt_t histoNo)

virtual Bool_t	PrepareViewData (PRawRunData *runData, const UInt_t histoNo)

Protected Member Functions inherited from PRunBase
virtual void	CalculateKaiserFilterCoeff (Double_t wc, Double_t A, Double_t dw)

virtual void	FilterTheo ()

Private Member Functions
virtual Bool_t	GetProperT0 (PRawRunData runData, PMsrGlobalBlock globalBlock, PUIntVector &histoNo)

virtual Bool_t	GetProperDataRange ()

virtual void	GetProperFitRange (PMsrGlobalBlock *globalBlock)

virtual void	EstimateN0 ()

virtual Bool_t	EstimateBkg (UInt_t histoNo)

virtual Bool_t	IsScaleN0AndBkg ()

Private Attributes
Bool_t	fScaleN0AndBkg
	true=scale N0 and background to 1/ns, otherwise 1/bin More...

UInt_t	fNoOfFitBins
	number of bins to be fitted More...

Double_t	fBackground
	needed if background range is given (units: 1/bin) More...

Int_t	fPacking
	packing for this particular run. Either given in the RUN- or GLOBAL-block. More...

Bool_t	fTheoAsData
	true=only calculate the theory points at the data points, false=calculate more points for the theory as compared to data are calculated which lead to 'nicer' Fouriers More...

Int_t	fGoodBins [2]
	keep first/last good bins. 0=fgb, 1=lgb More...

PDoubleVector	fForward
	forward histo data More...

Int_t	fStartTimeBin
	bin at which the fit starts More...

Int_t	fEndTimeBin
	bin at which the fit ends More...

Additional Inherited Members
Protected Attributes inherited from PRunBase
Bool_t	fValid
	flag showing if the state of the class is valid More...

EPMusrHandleTag	fHandleTag
	tag telling whether this is used for fit, view, ... More...

Int_t	fRunNo
	number of the run within the msr-file More...

PMsrHandler *	fMsrInfo
	msr-file handler More...

PMsrRunBlock *	fRunInfo
	run info used to filter out needed infos of a run More...

PRunDataHandler *	fRawData
	holds the raw run data More...

PRunData	fData
	data to be fitted, viewed, i.e. binned data More...

Double_t	fTimeResolution
	time resolution in (us) More...

PMetaData	fMetaData
	keeps the meta data from the data file like field, temperature, energy, ... More...

PDoubleVector	fT0s
	all t0 bins of a run! The derived classes will handle it. More...

std::vector< PDoubleVector >	fAddT0s
	all t0 bins of all addrun's of a run! The derived classes will handle it. More...

Double_t	fFitStartTime
	fit start time More...

Double_t	fFitEndTime
	fit end time More...

PDoubleVector	fFuncValues
	is keeping the values of the functions from the FUNCTIONS block More...

std::unique_ptr< PTheory >	fTheory
	theory needed to calculate chi-square More...

PDoubleVector	fKaiserFilter
	stores the Kaiser filter vector (needed for the RRF). More...

Detailed Description

Class handling single histogram fit type.

Definition at line 38 of file PRunSingleHisto.h.

Constructor & Destructor Documentation

◆ PRunSingleHisto() [1/2]

PRunSingleHisto::PRunSingleHisto ( )

Constructor

Definition at line 55 of file PRunSingleHisto.cpp.

References fBackground, fEndTimeBin, fGoodBins, fNoOfFitBins, fPacking, fScaleN0AndBkg, fStartTimeBin, and fTheoAsData.

◆ PRunSingleHisto() [2/2]

PRunSingleHisto::PRunSingleHisto	(	PMsrHandler *	msrInfo,
		PRunDataHandler *	rawData,
		UInt_t	runNo,
		EPMusrHandleTag	tag,
		Bool_t	theoAsData
	)

Constructor

Parameters

msrInfo	pointer to the msr-file handler
rawData	raw run data
runNo	number of the run within the msr-file
tag	tag showing what shall be done: kFit == fitting, kView == viewing

Definition at line 83 of file PRunSingleHisto.cpp.

References fBackground, fEndTimeBin, fGoodBins, PRunBase::fMsrInfo, fNoOfFitBins, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, fStartTimeBin, PRunBase::fValid, PMsrHandler::GetMsrGlobal(), PMsrGlobalBlock::GetPacking(), PMsrRunBlock::GetPacking(), IsScaleN0AndBkg(), and PrepareData().

◆ ~PRunSingleHisto()

PRunSingleHisto::~PRunSingleHisto ( )

virtual

Destructor

Definition at line 124 of file PRunSingleHisto.cpp.

References fForward.

Member Function Documentation

◆ CalcChiSquare()

Double_t PRunSingleHisto::CalcChiSquare ( const std::vector< Double_t > & par )

virtual

Calculate chi-square.

return:

chisq value

Parameters

par	parameter vector iterated by minuit2

Implements PRunBase.

Definition at line 140 of file PRunSingleHisto.cpp.

References PMsrHandler::EvalFunc(), fBackground, PRunBase::fData, fEndTimeBin, PRunBase::fFuncValues, PRunBase::fMetaData, PRunBase::fMsrInfo, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, fStartTimeBin, PRunBase::fTheory, PRunBase::fTimeResolution, PMsrRunBlock::GetBkgFitParamNo(), PMsrRunBlock::GetBkgFix(), PRunData::GetDataTimeStart(), PRunData::GetDataTimeStep(), PRunData::GetError(), PMsrHandler::GetFuncNo(), PMsrRunBlock::GetLifetimeParamNo(), PMsrRunBlock::GetMap(), PMsrHandler::GetNoOfFuncs(), PMsrRunBlock::GetNormParamNo(), PRunData::GetValue(), MSR_PARAM_FUN_OFFSET, PMUON_LIFETIME, and PMUSR_UNDEFINED.

◆ CalcChiSquareExpected()

Double_t PRunSingleHisto::CalcChiSquareExpected ( const std::vector< Double_t > & par )

virtual

Calculate expected chi-square.

return:

chisq value

Parameters

par	parameter vector iterated by minuit2

Definition at line 224 of file PRunSingleHisto.cpp.

References PMsrHandler::EvalFunc(), fBackground, PRunBase::fData, fEndTimeBin, PRunBase::fFuncValues, PRunBase::fMetaData, PRunBase::fMsrInfo, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, fStartTimeBin, PRunBase::fTheory, PRunBase::fTimeResolution, PMsrRunBlock::GetBkgFitParamNo(), PMsrRunBlock::GetBkgFix(), PRunData::GetDataTimeStart(), PRunData::GetDataTimeStep(), PMsrHandler::GetFuncNo(), PMsrRunBlock::GetLifetimeParamNo(), PMsrRunBlock::GetMap(), PMsrHandler::GetNoOfFuncs(), PMsrRunBlock::GetNormParamNo(), PRunData::GetValue(), MSR_PARAM_FUN_OFFSET, PMUON_LIFETIME, and PMUSR_UNDEFINED.

◆ CalcMaxLikelihood()

Double_t PRunSingleHisto::CalcMaxLikelihood ( const std::vector< Double_t > & par )

virtual

Calculate log maximum-likelihood. See http://pdg.lbl.gov/index.html

return:

log maximum-likelihood value

Parameters

par	parameter vector iterated by minuit2

Implements PRunBase.

Definition at line 309 of file PRunSingleHisto.cpp.

References PMsrHandler::EvalFunc(), fBackground, PRunBase::fData, fEndTimeBin, PRunBase::fFuncValues, PRunBase::fMetaData, PRunBase::fMsrInfo, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, fStartTimeBin, PRunBase::fTheory, PRunBase::fTimeResolution, PMsrRunBlock::GetBkgFitParamNo(), PMsrRunBlock::GetBkgFix(), PRunData::GetDataTimeStart(), PRunData::GetDataTimeStep(), PMsrHandler::GetFuncNo(), PMsrRunBlock::GetLifetimeParamNo(), PMsrRunBlock::GetMap(), PMsrHandler::GetNoOfFuncs(), PMsrRunBlock::GetNormParamNo(), PRunData::GetValue(), MSR_PARAM_FUN_OFFSET, PMUON_LIFETIME, and PMUSR_UNDEFINED.

◆ CalcMaxLikelihoodExpected()

Double_t PRunSingleHisto::CalcMaxLikelihoodExpected ( const std::vector< Double_t > & par )

virtual

Calculate expected log maximum-likelihood.

return:

log maximum-likelihood value

Parameters

par	parameter vector iterated by minuit2

Definition at line 407 of file PRunSingleHisto.cpp.

References PMsrHandler::EvalFunc(), fBackground, PRunBase::fData, fEndTimeBin, PRunBase::fFuncValues, PRunBase::fMetaData, PRunBase::fMsrInfo, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, fStartTimeBin, PRunBase::fTheory, PRunBase::fTimeResolution, PMsrRunBlock::GetBkgFitParamNo(), PMsrRunBlock::GetBkgFix(), PRunData::GetDataTimeStart(), PRunData::GetDataTimeStep(), PMsrHandler::GetFuncNo(), PMsrRunBlock::GetLifetimeParamNo(), PMsrRunBlock::GetMap(), PMsrHandler::GetNoOfFuncs(), PMsrRunBlock::GetNormParamNo(), PRunData::GetValue(), MSR_PARAM_FUN_OFFSET, PMUON_LIFETIME, and PMUSR_UNDEFINED.

◆ CalcNoOfFitBins()

void PRunSingleHisto::CalcNoOfFitBins ( )

virtual

Calculate the number of fitted bins for the current fit range.

Definition at line 667 of file PRunSingleHisto.cpp.

References PRunBase::fData, fEndTimeBin, PRunBase::fFitEndTime, PRunBase::fFitStartTime, fNoOfFitBins, fStartTimeBin, PRunData::GetDataTimeStart(), PRunData::GetDataTimeStep(), and PRunData::GetValue().

Referenced by GetNoOfFitBins(), PrepareFitData(), PrepareRawViewData(), and PrepareViewData().

◆ CalcTheory()

void PRunSingleHisto::CalcTheory ( )

virtual

Calculate theory for a given set of fit-parameters.

Implements PRunBase.

Definition at line 498 of file PRunSingleHisto.cpp.

References PRunData::AppendTheoryValue(), PMsrHandler::EvalFunc(), fBackground, PRunBase::fData, PRunBase::fFuncValues, PRunBase::fMetaData, PRunBase::fMsrInfo, PRunBase::fRunInfo, PRunBase::fTheory, PMsrRunBlock::GetBkgFitParamNo(), PMsrRunBlock::GetBkgFix(), PRunData::GetDataTimeStart(), PRunData::GetDataTimeStep(), PMsrHandler::GetFuncNo(), PMsrRunBlock::GetLifetimeParamNo(), PMsrRunBlock::GetMap(), PMsrHandler::GetMsrParamList(), PMsrHandler::GetNoOfFuncs(), PMsrRunBlock::GetNormParamNo(), PRunData::GetValue(), MSR_PARAM_FUN_OFFSET, PMUON_LIFETIME, and PMUSR_UNDEFINED.

◆ EstimateBkg()

Bool_t PRunSingleHisto::EstimateBkg ( UInt_t histoNo )

privatevirtual

Estimate the background for a given interval.

return:

true, if everything went smooth
false, otherwise

Parameters

histoNo forward histogram number of the run

Definition at line 1750 of file PRunSingleHisto.cpp.

References ACCEL_PERIOD_PSI, ACCEL_PERIOD_RAL, ACCEL_PERIOD_TRIUMF, fBackground, fForward, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, PRunBase::fTimeResolution, PMsrRunBlock::GetBkgRange(), PMsrRunBlock::GetInstitute(), and PMsrRunBlock::SetBkgEstimated().

Referenced by PrepareFitData(), PrepareRawViewData(), and PrepareViewData().

◆ EstimateN0()

void PRunSingleHisto::EstimateN0 ( )

privatevirtual

Estimate the N0 for the given run.

Definition at line 1653 of file PRunSingleHisto.cpp.

References fForward, PRunBase::fMsrInfo, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, PRunBase::fT0s, PRunBase::fTimeResolution, PMsrRunBlock::GetBkgFitParamNo(), PMsrHandler::GetMsrParamList(), PMsrRunBlock::GetNormParamNo(), PMUON_LIFETIME, PMsrHandler::SetMsrParamStep(), and PMsrHandler::SetMsrParamValue().

Referenced by PrepareFitData().

◆ GetBackground()

virtual Double_t PRunSingleHisto::GetBackground ( )

inlinevirtual

Definition at line 55 of file PRunSingleHisto.h.

References fBackground.

◆ GetEndTimeBin()

virtual Int_t PRunSingleHisto::GetEndTimeBin ( )

inlinevirtual

Definition at line 58 of file PRunSingleHisto.h.

References fEndTimeBin.

◆ GetNoOfFitBins()

UInt_t PRunSingleHisto::GetNoOfFitBins ( )

virtual

Calculate the number of fitted bins for the current fit range.

return: number of fitted bins.

Definition at line 564 of file PRunSingleHisto.cpp.

References CalcNoOfFitBins(), and fNoOfFitBins.

◆ GetPacking()

virtual Int_t PRunSingleHisto::GetPacking ( )

inlinevirtual

Definition at line 59 of file PRunSingleHisto.h.

References fPacking.

◆ GetProperDataRange()

Bool_t PRunSingleHisto::GetProperDataRange ( )

privatevirtual

Get the proper data range, i.e. first/last good bin (fgb/lgb).

get fgb/lgb from the RUN block
if fgb/lgb still undefined, try to get it from the GLOBAL block
if fgb/lgb still undefined, try to estimate them.

return:

true if everthing went smooth
false, otherwise.

Definition at line 1533 of file PRunSingleHisto.cpp.

References fForward, fGoodBins, PRunBase::fMsrInfo, PRunBase::fRunInfo, PRunBase::fT0s, PRunBase::fTimeResolution, PMsrGlobalBlock::GetDataRange(), PMsrRunBlock::GetDataRange(), PMsrHandler::GetMsrGlobal(), and PMsrRunBlock::SetDataRange().

Referenced by PrepareData().

◆ GetProperFitRange()

void PRunSingleHisto::GetProperFitRange ( PMsrGlobalBlock * globalBlock )

privatevirtual

Get the proper fit range. There are two possible fit range commands: fit <start> <end> given in (usec), or fit fgb+offset_0 lgb-offset_1 given in (bins), therefore it works the following way:

get fit range assuming given in time from RUN block
if fit range in RUN block is given in bins, replace start/end
if fit range is NOT given yet, try fit range assuming given in time from GLOBAL block
if fit range in GLOBAL block is given in bins, replace start/end
if still no fit range is given, use fgb/lgb.

Parameters

globalBlock pointer to the GLOBAL block information form the msr-file.

Definition at line 1614 of file PRunSingleHisto.cpp.

References PRunBase::fFitEndTime, PRunBase::fFitStartTime, fGoodBins, PRunBase::fRunInfo, PRunBase::fT0s, PRunBase::fTimeResolution, PMsrGlobalBlock::GetFitRange(), PMsrRunBlock::GetFitRange(), PMsrGlobalBlock::GetFitRangeOffset(), PMsrRunBlock::GetFitRangeOffset(), PMsrGlobalBlock::IsFitRangeInBin(), PMsrRunBlock::IsFitRangeInBin(), PMUSR_UNDEFINED, PMsrGlobalBlock::SetFitRange(), and PMsrRunBlock::SetFitRange().

Referenced by PrepareData().

◆ GetProperT0()

Bool_t PRunSingleHisto::GetProperT0	(	PRawRunData *	runData,
		PMsrGlobalBlock *	globalBlock,
		PUIntVector &	histoNo
	)

privatevirtual

Get the proper t0 for the single histogram run.

the t0 vector size = number of detectors (grouping) for forward.
initialize t0's with -1
fill t0's from RUN block
if t0's are missing (i.e. t0 == -1), try to fill from the GLOBAL block.
if t0's are missing, try t0's from the data file
if t0's are missing, try to estimate them

Parameters

runData	pointer to the current RUN block entry from the msr-file
globalBlock	pointer to the GLOBLA block entry from the msr-file
histoNo	histogram number vector of forward; histoNo = msr-file forward + redGreen_offset - 1

return:

true if everthing went smooth
false, otherwise.

Definition at line 1402 of file PRunSingleHisto.cpp.

References PRunBase::fAddT0s, PRunBase::fRawData, PRunBase::fRunInfo, PRunBase::fT0s, PMsrRunBlock::GetAddT0Bin(), PRawRunData::GetDataBin(), PMsrRunBlock::GetForwardHistoNoSize(), PRunDataHandler::GetRunData(), PMsrRunBlock::GetRunName(), PMsrRunBlock::GetRunNameSize(), PRawRunData::GetT0Bin(), PMsrGlobalBlock::GetT0Bin(), PMsrRunBlock::GetT0Bin(), PRawRunData::GetT0BinEstimated(), PMsrGlobalBlock::GetT0BinSize(), PMsrRunBlock::GetT0BinSize(), PMsrRunBlock::SetAddT0Bin(), and PMsrRunBlock::SetT0Bin().

Referenced by PrepareData().

◆ GetScaleN0AndBkg()

virtual Bool_t PRunSingleHisto::GetScaleN0AndBkg ( )

inlinevirtual

Definition at line 60 of file PRunSingleHisto.h.

References fScaleN0AndBkg.

◆ GetStartTimeBin()

virtual Int_t PRunSingleHisto::GetStartTimeBin ( )

inlinevirtual

Definition at line 57 of file PRunSingleHisto.h.

References fStartTimeBin.

◆ IsScaleN0AndBkg()

Bool_t PRunSingleHisto::IsScaleN0AndBkg ( )

privatevirtual

Checks if N0/Bkg normalization to 1/ns is whished. The default is yes, since most of the users want to have it that way. To overwrite this, one should add the line 'SCALE_N0_BKG FALSE' to the command block of the msr-file.

return:

true, if scaling of N0 and Bkg to 1/ns is whished
false, otherwise

Parameters

histoNo forward histogram number of the run

Definition at line 1834 of file PRunSingleHisto.cpp.

References PRunBase::fMsrInfo, and PMsrHandler::GetMsrCommands().

Referenced by PRunSingleHisto().

◆ PrepareData()

Bool_t PRunSingleHisto::PrepareData ( )

protectedvirtual

Prepare data for fitting or viewing. What is already processed at this stage:

get proper raw run data
get all needed forward histograms
get time resolution
get t0's and perform necessary cross checks (e.g. if t0 of msr-file (if present) are consistent with t0 of the data files, etc.)
add runs (if addruns are present)
group histograms (if grouping is present)

return:

true if everthing went smooth
false, otherwise.

Implements PRunBase.

Definition at line 699 of file PRunSingleHisto.cpp.

References PRunBase::fAddT0s, PMetaData::fEnergy, PMetaData::fField, fForward, PRunBase::fHandleTag, PRunBase::fMetaData, PRunBase::fMsrInfo, PRunBase::fRawData, PRunBase::fRunInfo, PRunBase::fT0s, PMetaData::fTemp, PRunBase::fTimeResolution, PRunBase::fValid, PRawRunData::GetDataBin(), PRawRunData::GetEnergy(), PRawRunData::GetField(), PMsrRunBlock::GetForwardHistoNo(), PMsrRunBlock::GetForwardHistoNoSize(), PMsrHandler::GetMsrGlobal(), PMsrHandler::GetMsrPlotList(), PRawRunData::GetNoOfTemperatures(), GetProperDataRange(), GetProperFitRange(), GetProperT0(), PRunDataHandler::GetRunData(), PMsrRunBlock::GetRunName(), PMsrRunBlock::GetRunNameSize(), PRawRunData::GetTemperature(), PRawRunData::GetTimeResolution(), PRawRunData::IsPresent(), kFit, kView, PrepareFitData(), PrepareRawViewData(), and PrepareViewData().

Referenced by PRunSingleHisto().

◆ PrepareFitData()

Bool_t PRunSingleHisto::PrepareFitData	(	PRawRunData *	runData,
		const UInt_t	histoNo
	)

protectedvirtual

Take the pre-processed data (i.e. grouping and addrun are preformed) and form the histogram for fitting. The following steps are preformed:

get fit start/stop time
check that 'first good data bin', 'last good data bin', and 't0' make any sense
check how the background shall be handled, i.e. fitted, subtracted from background estimate data range, or subtacted from a given fixed background.
packing (i.e rebinning)

return:

true, if everything went smooth
false, otherwise

Parameters

runData	raw run data handler
histoNo	forward histogram number

Definition at line 851 of file PRunSingleHisto.cpp.

References PRunData::AppendErrorValue(), PRunData::AppendValue(), CalcNoOfFitBins(), EstimateBkg(), EstimateN0(), PMsrHandler::EstimateN0(), PRunBase::fData, fForward, fGoodBins, PRunBase::fMsrInfo, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, PRunBase::fT0s, PRunBase::fTimeResolution, PMsrRunBlock::GetBkgFitParamNo(), PMsrRunBlock::GetBkgFix(), PMsrRunBlock::GetBkgRange(), PMUSR_UNDEFINED, PMsrRunBlock::SetBkgRange(), PRunData::SetDataTimeStart(), and PRunData::SetDataTimeStep().

Referenced by PrepareData().

◆ PrepareRawViewData()

Bool_t PRunSingleHisto::PrepareRawViewData	(	PRawRunData *	runData,
		const UInt_t	histoNo
	)

protectedvirtual

Take the pre-processed data (i.e. grouping and addrun are preformed) and form the histogram for viewing without any life time correction.

The following steps are preformed:

check if view packing is whished.
check that 'first good data bin', 'last good data bin', and 't0' makes any sense
packing (i.e. rebinnig)
calculate theory

return:

true, if everything went smooth
false, otherwise.

Parameters

runData	raw run data handler
histoNo	forward histogram number

Definition at line 944 of file PRunSingleHisto.cpp.

References PRunData::AppendErrorValue(), PRunData::AppendTheoryValue(), PRunData::AppendValue(), CalcNoOfFitBins(), EstimateBkg(), PMsrHandler::EvalFunc(), fBackground, PRunBase::fData, fForward, PRunBase::fFuncValues, fGoodBins, PRunBase::fMetaData, PRunBase::fMsrInfo, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, PRunBase::fT0s, fTheoAsData, PRunBase::fTheory, PRunBase::fTimeResolution, PMsrRunBlock::GetBkgFitParamNo(), PMsrRunBlock::GetBkgFix(), PMsrRunBlock::GetBkgRange(), PRunData::GetDataTimeStart(), PRunData::GetDataTimeStep(), PMsrHandler::GetFuncNo(), PMsrRunBlock::GetLifetimeParamNo(), PMsrRunBlock::GetMap(), PMsrHandler::GetMsrParamList(), PMsrHandler::GetMsrPlotList(), PMsrHandler::GetNoOfFuncs(), PMsrRunBlock::GetNormParamNo(), PRunData::GetTheoryTimeStart(), PRunData::GetTheoryTimeStep(), MSR_PARAM_FUN_OFFSET, PMUON_LIFETIME, PMUSR_UNDEFINED, PMsrRunBlock::SetBkgRange(), PRunData::SetDataTimeStart(), PRunData::SetDataTimeStep(), PRunData::SetTheoryTimeStart(), and PRunData::SetTheoryTimeStep().

Referenced by PrepareData().

◆ PrepareViewData()

Bool_t PRunSingleHisto::PrepareViewData	(	PRawRunData *	runData,
		const UInt_t	histoNo
	)

protectedvirtual

Take the pre-processed data (i.e. grouping and addrun are preformed) and form the histogram for viewing with life time correction, i.e. the exponential decay is removed.

The following steps are preformed:

check if view packing is whished.
check that 'first good data bin', 'last good data bin', and 't0' makes any sense
transform data sets (see below).
calculate theory

Muon life time corrected data: Starting from

$N(t) = N_0 e^{-t/\tau} [ 1 + A(t) ] + \mathrm{Bkg}$

it follows that

$A(t) = (-1) + e^{+t/\tau}\, \frac{N(t)-\mathrm{Bkg}}{N_0}.$

For the error estimate only the statistical error of $ N(t) $ is used, and hence

$\Delta A(t) = \frac{e^{t/\tau}}{N_0}\,\sqrt{\frac{N(t)}{p}}$

where $ p $ is the packing, and $ N(t) $ are the packed data, i.e.

$N(t_i) = \frac{1}{p}\, \sum_{j=i}^{i+p} n_j$

with $ n_j $ the raw histogram data bins.

return:

true, if everything went smooth
false, otherwise

Parameters

runData	raw run data handler
histoNo	forward histogram number

Definition at line 1133 of file PRunSingleHisto.cpp.

References PRunData::AppendErrorValue(), PRunData::AppendTheoryValue(), PRunData::AppendValue(), CalcNoOfFitBins(), PRunBase::CalculateKaiserFilterCoeff(), EstimateBkg(), PMsrHandler::EvalFunc(), fBackground, PRunBase::fData, fForward, PRunBase::fFuncValues, fGoodBins, PRunBase::FilterTheo(), PRunBase::fMetaData, PRunBase::fMsrInfo, fPacking, PRunBase::fRunInfo, fScaleN0AndBkg, PRunBase::fT0s, fTheoAsData, PRunBase::fTheory, PRunBase::fTimeResolution, GAMMA_BAR_MUON, PMsrRunBlock::GetBkgFitParamNo(), PMsrRunBlock::GetBkgFix(), PMsrRunBlock::GetBkgRange(), PRunData::GetDataTimeStart(), PRunData::GetDataTimeStep(), PMsrHandler::GetFuncNo(), PMsrRunBlock::GetLifetimeParamNo(), PMsrRunBlock::GetMap(), PMsrHandler::GetMsrParamList(), PMsrHandler::GetMsrPlotList(), PMsrHandler::GetNoOfFuncs(), PMsrRunBlock::GetNormParamNo(), PRunData::GetTheory(), PRunData::GetTheoryTimeStart(), PRunData::GetTheoryTimeStep(), MSR_PARAM_FUN_OFFSET, PMUON_LIFETIME, PMUSR_UNDEFINED, PRunData::ReplaceTheory(), RRF_UNIT_G, RRF_UNIT_kHz, RRF_UNIT_Mcs, RRF_UNIT_MHz, RRF_UNIT_T, PMsrRunBlock::SetBkgRange(), PRunData::SetDataTimeStart(), PRunData::SetDataTimeStep(), PRunData::SetTheoryTimeStart(), and PRunData::SetTheoryTimeStep().

Referenced by PrepareData().

◆ SetFitRangeBin()

void PRunSingleHisto::SetFitRangeBin ( const TString fitRange )

virtual

Allows to change the fit range on the fly. Used in the COMMAND block. The syntax of the string is: FIT_RANGE fgb[+n00] lgb[-n01] [fgb[+n10] lgb[-n11] ... fgb[+nN0] lgb[-nN1]]. If only one pair of fgb/lgb is given, it is used for all runs in the RUN block section. If multiple fgb/lgb's are given, the number N has to be the number of RUN blocks in the msr-file.

nXY are offsets which can be used to shift, limit the fit range.

Parameters

fitRange string containing the necessary information.

Definition at line 585 of file PRunSingleHisto.cpp.

References PRunBase::fFitEndTime, PRunBase::fFitStartTime, fGoodBins, PRunBase::fRunNo, PRunBase::fT0s, and PRunBase::fTimeResolution.

Member Data Documentation

◆ fBackground

Double_t PRunSingleHisto::fBackground

private

needed if background range is given (units: 1/bin)

Definition at line 73 of file PRunSingleHisto.h.

Referenced by CalcChiSquare(), CalcChiSquareExpected(), CalcMaxLikelihood(), CalcMaxLikelihoodExpected(), CalcTheory(), EstimateBkg(), GetBackground(), PrepareRawViewData(), PrepareViewData(), and PRunSingleHisto().

◆ fEndTimeBin

Int_t PRunSingleHisto::fEndTimeBin

private

bin at which the fit ends

Definition at line 82 of file PRunSingleHisto.h.

Referenced by CalcChiSquare(), CalcChiSquareExpected(), CalcMaxLikelihood(), CalcMaxLikelihoodExpected(), CalcNoOfFitBins(), GetEndTimeBin(), and PRunSingleHisto().

◆ fForward

PDoubleVector PRunSingleHisto::fForward

private

forward histo data

Definition at line 79 of file PRunSingleHisto.h.

Referenced by EstimateBkg(), EstimateN0(), GetProperDataRange(), PrepareData(), PrepareFitData(), PrepareRawViewData(), PrepareViewData(), and ~PRunSingleHisto().

◆ fGoodBins

Int_t PRunSingleHisto::fGoodBins[2]

private

keep first/last good bins. 0=fgb, 1=lgb

Definition at line 77 of file PRunSingleHisto.h.

Referenced by GetProperDataRange(), GetProperFitRange(), PrepareFitData(), PrepareRawViewData(), PrepareViewData(), PRunSingleHisto(), and SetFitRangeBin().

◆ fNoOfFitBins

UInt_t PRunSingleHisto::fNoOfFitBins

private

number of bins to be fitted

Definition at line 72 of file PRunSingleHisto.h.

Referenced by CalcNoOfFitBins(), GetNoOfFitBins(), and PRunSingleHisto().

◆ fPacking

Int_t PRunSingleHisto::fPacking

private

packing for this particular run. Either given in the RUN- or GLOBAL-block.

Definition at line 74 of file PRunSingleHisto.h.

Referenced by CalcChiSquare(), CalcChiSquareExpected(), CalcMaxLikelihood(), CalcMaxLikelihoodExpected(), EstimateBkg(), EstimateN0(), GetPacking(), PrepareFitData(), PrepareRawViewData(), PrepareViewData(), and PRunSingleHisto().

◆ fScaleN0AndBkg

Bool_t PRunSingleHisto::fScaleN0AndBkg

private

true=scale N0 and background to 1/ns, otherwise 1/bin

Definition at line 71 of file PRunSingleHisto.h.

Referenced by CalcChiSquare(), CalcChiSquareExpected(), CalcMaxLikelihood(), CalcMaxLikelihoodExpected(), EstimateBkg(), EstimateN0(), GetScaleN0AndBkg(), PrepareFitData(), PrepareRawViewData(), PrepareViewData(), and PRunSingleHisto().

◆ fStartTimeBin

Int_t PRunSingleHisto::fStartTimeBin

private

bin at which the fit starts

Definition at line 81 of file PRunSingleHisto.h.

Referenced by CalcChiSquare(), CalcChiSquareExpected(), CalcMaxLikelihood(), CalcMaxLikelihoodExpected(), CalcNoOfFitBins(), GetStartTimeBin(), and PRunSingleHisto().

◆ fTheoAsData

Bool_t PRunSingleHisto::fTheoAsData

private

true=only calculate the theory points at the data points, false=calculate more points for the theory as compared to data are calculated which lead to 'nicer' Fouriers

Definition at line 75 of file PRunSingleHisto.h.

Referenced by PrepareRawViewData(), PrepareViewData(), and PRunSingleHisto().

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

/home/nemu/Apps/musrfit/src/include/PRunSingleHisto.h
/home/nemu/Apps/musrfit/src/classes/PRunSingleHisto.cpp

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PRunSingleHisto() [1/2]

◆ PRunSingleHisto() [2/2]

◆ ~PRunSingleHisto()

Member Function Documentation

◆ CalcChiSquare()

◆ CalcChiSquareExpected()

◆ CalcMaxLikelihood()

◆ CalcMaxLikelihoodExpected()

◆ CalcNoOfFitBins()

◆ CalcTheory()

◆ EstimateBkg()

◆ EstimateN0()

◆ GetBackground()

◆ GetEndTimeBin()

◆ GetNoOfFitBins()

◆ GetPacking()

◆ GetProperDataRange()

◆ GetProperFitRange()

◆ GetProperT0()

◆ GetScaleN0AndBkg()

◆ GetStartTimeBin()

◆ IsScaleN0AndBkg()

◆ PrepareData()

◆ PrepareFitData()

◆ PrepareRawViewData()

◆ PrepareViewData()

◆ SetFitRangeBin()

Member Data Documentation

◆ fBackground

◆ fEndTimeBin

◆ fForward

◆ fGoodBins

◆ fNoOfFitBins

◆ fPacking

◆ fScaleN0AndBkg

◆ fStartTimeBin

◆ fTheoAsData