BinFlipChi2 Class Reference

#include <BinFlipChi2.h>

Inheritance diagram for BinFlipChi2:

Public Member Functions
	BinFlipChi2 (BinFlipParSet *, const HadronicParameters &, const TimeBinning &)
	BinFlipChi2 (BinFlipParSet *, const std::string &, const std::string &, const std::string &fname="")
virtual double	getVal () override
std::deque< std::deque< TH1F > >	plotsVsTime (const std::string &) const
void	savePlotsVsTime (const std::string &, TFile &) const
Public Member Functions inherited from BinFlipChi2Base
	BinFlipChi2Base (BinFlipParSet *)
virtual void	parametersChanged () override
TH2D	scan2D (unsigned, unsigned, float nSigmaRange=5., unsigned nBins=300, bool zeroCentre=true, float scale=1.)
BinFlipParSet *	getBinFlipParSet ()
Public Member Functions inherited from MINT::Minimisable
	Minimisable (MinuitParameterSet *mps=0)
	Minimisable (const Minimisable &other)
virtual	~Minimisable ()
void	setPset (MinuitParameterSet *mps)
MinuitParameterSet *	getParSet ()
virtual void	beginFit ()
virtual void	endFit ()
double	getNewVal ()
Public Member Functions inherited from MINT::IMinimisable
virtual void	Gradient (std::vector< double > &grad)
virtual bool	useAnalyticGradient ()
virtual void	setUseAnalyticGradient (bool useAnalyticGradient)
virtual	~IMinimisable ()

Private Member Functions
void	checkPhaseBinning () const
void	init ()

Private Attributes
HadronicParameters	m_hadronicPars
TimeBinning	m_timeBinning
double	m_lifetime

Additional Inherited Members
Protected Attributes inherited from BinFlipChi2Base
BinFlipParSet *	m_fitPars
std::complex< double >	m_zcp
std::complex< double >	m_dz

Detailed Description

Definition at line 50 of file BinFlipChi2.h.

Constructor & Destructor Documentation

BinFlipChi2() [1/2]

BinFlipChi2::BinFlipChi2	(	BinFlipParSet *	fitPars,
		const HadronicParameters &	hadronicPars,
		const TimeBinning &	timeBins
	)

Definition at line 145 of file BinFlipChi2.cpp.

                                          :
  BinFlipChi2Base(fitPars),
  m_hadronicPars(hadronicPars),
  m_timeBinning(timeBins)
{
  init() ;
}

BinFlipChi2() [2/2]

BinFlipChi2::BinFlipChi2	(	BinFlipParSet *	,
		const std::string &	,
		const std::string &	,
		const std::string &	fname = ""
	)

Member Function Documentation

checkPhaseBinning()

void BinFlipChi2::checkPhaseBinning ( ) const

private

Definition at line 168 of file BinFlipChi2.cpp.

                                          {
  if(m_hadronicPars.binning() != *m_timeBinning.phaseBinning())
    throw invalid_argument("The phase binning of the hadronic parameters and binned data don't match!") ;
}

getVal()

double BinFlipChi2::getVal ( )

overridevirtual

Implements MINT::Minimisable.

Definition at line 173 of file BinFlipChi2.cpp.

                           {
  double chi2 = 0. ;
  for(unsigned iTimeBin = 0 ; iTimeBin < m_timeBinning.nBinsTime() ; ++iTimeBin){
    for(unsigned iPhaseBin = 1 ; iPhaseBin < m_timeBinning.nBinsPhase() + 1; ++iPhaseBin){
      // Needs the unmixed time.
      double Rplus = m_hadronicPars.bin(iPhaseBin).R(m_timeBinning.meanUnmixedTime(iTimeBin),
                             m_timeBinning.meanUnmixedTime2(iTimeBin),
                             m_timeBinning.getLifetime(), m_zcp, m_dz) ;
      double Rminus = m_hadronicPars.bin(iPhaseBin).Rbar(m_timeBinning.meanUnmixedTime(iTimeBin),
                             m_timeBinning.meanUnmixedTime2(iTimeBin),
                             m_timeBinning.getLifetime(), m_zcp, m_dz) ;
      chi2 += m_timeBinning.chiSquared(iTimeBin, iPhaseBin, Rplus, Rminus) ;
    }
  }
  return chi2 ;
}

init()

void BinFlipChi2::init ( )

private

Definition at line 164 of file BinFlipChi2.cpp.

                       {
  checkPhaseBinning() ;
}

plotsVsTime()

deque< deque< TH1F > > BinFlipChi2::plotsVsTime

(

const std::string &

)

const

Definition at line 190 of file BinFlipChi2.cpp.

                                                                     {
  deque<deque<TH1F> > histos = m_timeBinning.plotsVsTime(name) ;
  unsigned iPhaseBin = 1 ;
  for(auto& binhistos : histos) {
    TH1F& rminus = binhistos[2] ;
    TH1F& rplus = binhistos[5] ;
    TH1F& rdiff = binhistos[6] ;
    TH1F rminusfit = TH1F(rminus) ;
    rminusfit.SetName((string(rminus.GetName()) + "_fit").c_str()) ;
    TH1F rplusfit = TH1F(rplus) ;
    rplusfit.SetName((string(rplus.GetName()) + "_fit").c_str()) ;
    TH1F rdifffit = TH1F(rdiff) ;
    rdifffit.SetName((string(rdiff.GetName()) + "_fit").c_str()) ;
    for(unsigned iTimeBin = 0 ; iTimeBin < m_timeBinning.nBinsTime() ; ++iTimeBin){
      double t = m_timeBinning.meanUnmixedTime(iTimeBin) ;
      double t2 = m_timeBinning.meanUnmixedTime2(iTimeBin) ;
      double rplus = m_hadronicPars.bin(iPhaseBin).R(t, t2, m_timeBinning.getLifetime(), m_zcp, m_dz) ;
      double rminus = m_hadronicPars.bin(iPhaseBin).Rbar(t, t2, m_timeBinning.getLifetime(), m_zcp, m_dz) ;
      double rdiff = rplus-rminus ;
      rminusfit.SetBinContent(iTimeBin+1, rminus) ;
      rminusfit.SetBinError(iTimeBin+1, 0.) ;
      rplusfit.SetBinContent(iTimeBin+1, rplus) ;
      rplusfit.SetBinError(iTimeBin+1, 0.) ;
      rdifffit.SetBinContent(iTimeBin+1, rdiff) ;
      rdifffit.SetBinError(iTimeBin+1, 0.) ;
    }
    binhistos.push_back(rminusfit) ;
    binhistos.push_back(rplusfit) ;
    binhistos.push_back(rdifffit) ;
    ++iPhaseBin ;
  }
  return histos ;
}

savePlotsVsTime()

void BinFlipChi2::savePlotsVsTime	(	const std::string &	,
		TFile &
	)		const

Definition at line 224 of file BinFlipChi2.cpp.

                                                                             {
  outputfile.cd() ;
  deque<deque<TH1F> > histos = plotsVsTime(name) ;
  for(auto& binhistos : histos)
    for(auto& histo : binhistos)
      histo.Write() ;
}

Member Data Documentation

m_hadronicPars

HadronicParameters BinFlipChi2::m_hadronicPars

private

Definition at line 60 of file BinFlipChi2.h.

m_lifetime

double BinFlipChi2::m_lifetime

private

Definition at line 62 of file BinFlipChi2.h.

m_timeBinning

TimeBinning BinFlipChi2::m_timeBinning

private

Definition at line 61 of file BinFlipChi2.h.

---

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

• Mint/BinFlipChi2.h
• src/BinFlip/BinFlipChi2.cpp