TimeBinning Class Reference

#include <TimeBinning.h>

Classes
class	Bin

Public Types
typedef std::deque< Bin >	Bins
typedef std::deque< Bins >	Bins2D

Public Member Functions
	TimeBinning (const std::vector< double > &, HadronicParameters::BinningPtr, double)
	TimeBinning (const std::string &, const std::string &fname="")
virtual	~TimeBinning ()
void	add (IDalitzEvent &, int, double, double weight=1.)
HadronicParameters::BinningPtr	phaseBinning () const
int	timeBin (double) const
void	Print (const std::string &, std::ostream &os=std::cout) const
void	write (const std::string &, const std::string &) const
double	chiSquared (unsigned, unsigned, double, double) const
unsigned	nBinsTime () const
unsigned	nBinsPhase () const
const Bin &	integratedBin (unsigned) const
const Bin &	bin (unsigned, unsigned) const
const Bin &	binBar (unsigned, unsigned) const
std::deque< TH1F >	plotVsTime (const std::string &, unsigned, int) const
std::deque< std::deque< TH1F > >	plotsVsTime (const std::string &) const
void	savePlotsVsTime (const std::string &, TFile &) const
double	meanUnmixedTime (unsigned) const
double	meanUnmixedTime2 (unsigned) const
void	setLifetime (double)
double	getLifetime () const
bool	isConsistent (const TimeBinning &) const
	Check that two TimeBinnings use the same binning scheme. More...
TimeBinning	operator+ (TimeBinning) const
TimeBinning &	operator+= (const TimeBinning &)

Private Member Functions
Bin &	_integratedBin (unsigned)
Bin &	_bin (unsigned, unsigned)
Bin &	_binBar (unsigned, unsigned)
virtual double	unmixedTimeMoment (unsigned, double, int) const

Private Attributes
std::vector< double >	m_timeBins
std::vector< double >	m_meant
std::vector< double >	m_meant2
double	m_lifetime
Bins2D	m_bins
Bins2D	m_binsBar
Bins	m_binsInt
HadronicParameters::BinningPtr	m_phaseBinning

Detailed Description

Definition at line 12 of file TimeBinning.h.

Member Typedef Documentation

Bins

typedef std::deque<Bin> TimeBinning::Bins

Definition at line 41 of file TimeBinning.h.

Bins2D

typedef std::deque<Bins> TimeBinning::Bins2D

Definition at line 42 of file TimeBinning.h.

Constructor & Destructor Documentation

TimeBinning() [1/2]

TimeBinning::TimeBinning	(	const std::vector< double > &	,
		HadronicParameters::BinningPtr	,
		double
	)

TimeBinning() [2/2]

TimeBinning::TimeBinning	(	const std::string &	,
		const std::string &	fname = ""
	)

~TimeBinning()

virtual TimeBinning::~TimeBinning ( )

inlinevirtual

Definition at line 46 of file TimeBinning.h.

{};

Member Function Documentation

_bin()

TimeBinning::Bin & TimeBinning::_bin ( unsigned  iTimeBin, unsigned  iPhaseBin  )

private

Definition at line 274 of file TimeBinning.cpp.

                                                                       {
  return m_bins.at(iTimeBin).at(iPhaseBin-1) ;
}

_binBar()

TimeBinning::Bin & TimeBinning::_binBar ( unsigned  iTimeBin, unsigned  iPhaseBin  )

private

Definition at line 278 of file TimeBinning.cpp.

                                                                          {
  return m_binsBar.at(iTimeBin).at(iPhaseBin-1) ;
}

_integratedBin()

TimeBinning::Bin & TimeBinning::_integratedBin ( unsigned  i )

private

Definition at line 270 of file TimeBinning.cpp.

                                                      {
  return m_binsInt.at(i) ;
}

add()

void TimeBinning::add	(	IDalitzEvent &	evt,
		int	tag,
		double	t,
		double	weight = 1.
	)

Definition at line 229 of file TimeBinning.cpp.

                                                                         {
  int timeBinNo = timeBin(t) ;
  if(timeBinNo < 0)
    return ;
  int phaseBin = m_phaseBinning->binNumber(evt) ;
  bool isPlus = phaseBin > 0 ;
  if(tag > 0){
    _bin(timeBinNo, abs(phaseBin)).add(t, isPlus, weight) ;
  }
  else{
    isPlus = !isPlus ;
    _binBar(timeBinNo, abs(phaseBin)).add(t, isPlus, weight) ;
  }
  _integratedBin(timeBinNo).add(t, isPlus, weight) ;
}

bin()

const TimeBinning::Bin & TimeBinning::bin	(	unsigned	iTimeBin,
		unsigned	iPhaseBin
	)		const

Definition at line 262 of file TimeBinning.cpp.

                                                                                  {
  return m_bins.at(iTimeBin).at(iPhaseBin-1) ;
}

binBar()

const TimeBinning::Bin & TimeBinning::binBar	(	unsigned	iTimeBin,
		unsigned	iPhaseBin
	)		const

Definition at line 266 of file TimeBinning.cpp.

                                                                                     {
  return m_binsBar.at(iTimeBin).at(iPhaseBin-1) ;
}

chiSquared()

double TimeBinning::chiSquared	(	unsigned	iTimeBin,
		unsigned	iPhaseBin,
		double	Rplus,
		double	Rminus
	)		const

Definition at line 245 of file TimeBinning.cpp.

                                                                                                       {
  return bin(iTimeBin, iPhaseBin).chiSquared(Rplus)
    + binBar(iTimeBin, iPhaseBin).chiSquared(Rminus) ;
}

getLifetime()

double TimeBinning::getLifetime

(

)

const

Definition at line 353 of file TimeBinning.cpp.

                                      {
  return m_lifetime ;
}

integratedBin()

const TimeBinning::Bin & TimeBinning::integratedBin

(

unsigned

i

)

const

Definition at line 258 of file TimeBinning.cpp.

                                                                 {
  return m_binsInt.at(i) ;
}

isConsistent()

bool TimeBinning::isConsistent

(

const TimeBinning &

other

)

const

Check that two TimeBinnings use the same binning scheme.

Definition at line 366 of file TimeBinning.cpp.

                                                             {
  return m_timeBins == other.m_timeBins && m_meant == other.m_meant && m_meant2 == other.m_meant2 \
    && m_lifetime == other.m_lifetime && *m_phaseBinning == *other.m_phaseBinning;
}

meanUnmixedTime()

double TimeBinning::meanUnmixedTime

(

unsigned

ibin

)

const

Definition at line 357 of file TimeBinning.cpp.

                                                       {
  return m_meant.at(ibin) ;
}

meanUnmixedTime2()

double TimeBinning::meanUnmixedTime2

(

unsigned

ibin

)

const

Definition at line 361 of file TimeBinning.cpp.

                                                        {
  return m_meant2.at(ibin) ;
}

nBinsPhase()

unsigned TimeBinning::nBinsPhase

(

)

const

Definition at line 254 of file TimeBinning.cpp.

                                       {
  return m_phaseBinning->nBins ;
}

nBinsTime()

unsigned TimeBinning::nBinsTime

(

)

const

Definition at line 250 of file TimeBinning.cpp.

                                      {
  return m_timeBins.size()-1 ;
}

operator+()

TimeBinning TimeBinning::operator+

(

TimeBinning

other

)

const

Definition at line 371 of file TimeBinning.cpp.

                                                          {
  return other += *this;
}

operator+=()

TimeBinning & TimeBinning::operator+=

(

const TimeBinning &

other

)

Definition at line 390 of file TimeBinning.cpp.

                                                             {
  if(!isConsistent(other))
    throw std::invalid_argument("Attempt to combine incompatible TimeBinning instances!");
  //m_bins += other.m_bins;
  //m_binsBar += other.m_binsBar;
  //m_binsInt += other.m_binsInt;
  for(unsigned iTimeBin = 0; iTimeBin < nBinsTime(); ++iTimeBin){
    m_binsInt.at(iTimeBin) += other.m_binsInt.at(iTimeBin);
    for(unsigned iPhaseBin = 0; iPhaseBin < nBinsPhase(); ++iPhaseBin){
      m_bins.at(iTimeBin).at(iPhaseBin) += other.m_bins.at(iTimeBin).at(iPhaseBin);
      m_binsBar.at(iTimeBin).at(iPhaseBin) += other.m_binsBar.at(iTimeBin).at(iPhaseBin);
    }
  }
  return *this;
}

phaseBinning()

HadronicParameters::BinningPtr TimeBinning::phaseBinning

(

)

const

Definition at line 218 of file TimeBinning.cpp.

                                                             {
  return m_phaseBinning ;
}

plotsVsTime()

deque< deque< TH1F > > TimeBinning::plotsVsTime

(

const std::string &

)

const

Definition at line 308 of file TimeBinning.cpp.

                                                                     {
  deque<deque<TH1F> > histos ;
  for(unsigned iphase = 1 ; iphase < m_phaseBinning->nBins + 1 ; ++iphase){
    histos.push_back(deque<TH1F>()) ;
    deque<TH1F>& binhistos = histos.back() ;
    for(int tag = -1 ; tag < 2 ; tag += 2){
      deque<TH1F> taghistos = plotVsTime(name, iphase, tag) ;
      binhistos.insert(binhistos.end(), taghistos.begin(), taghistos.end()) ;
    }
    TH1F rdiff(binhistos.back()) ;
    string rdiffname(Bin::getName(name + "_phase", iphase) + "ratiodiff") ;
    rdiff.SetName(rdiffname.c_str()) ;
    rdiff.Add(&binhistos[2], -1.) ;
    binhistos.push_back(rdiff) ;
  }
  return histos ;
}

plotVsTime()

deque< TH1F > TimeBinning::plotVsTime	(	const std::string &	,
		unsigned	,
		int
	)		const

Definition at line 282 of file TimeBinning.cpp.

                                                                                         {
  ostringstream stag ;
  stag << tag ;
  string name = Bin::getName(_name + "_" + stag.str() + "_phase", phaseBin) ;
  
  TH1F hplus((name + "nPlus").c_str(), "", nBinsTime(), &m_timeBins[0]) ;
  TH1F hminus((name + "nMinus").c_str(), "", nBinsTime(), &m_timeBins[0]) ;
  TH1F hr((name + "ratio").c_str(), "", nBinsTime(), &m_timeBins[0]) ;
  for(unsigned i = 1 ; i < nBinsTime()+1 ; ++i) {
    const Bin& pbin = tag > 0 ? bin(i-1, phaseBin) : binBar(i-1, phaseBin) ;
    hplus.SetBinContent(i, pbin.nPlus()) ;
    hplus.SetBinError(i, pbin.nPlusErr()) ;
    hminus.SetBinContent(i, pbin.nMinus()) ;
    hminus.SetBinError(i, pbin.nMinusErr()) ;
    double r = pbin.nPlus() > 0. ? pbin.nMinus()/pbin.nPlus() : 0. ;
    double rerr = r > 0 && pbin.nMinus() > 0 ? r * sqrt(pow(pbin.nPlusErr()/pbin.nPlus(), 2) + pow(pbin.nMinusErr()/pbin.nMinus(), 2)) : 0. ;
    hr.SetBinContent(i, r) ;
    hr.SetBinError(i, rerr) ;
  }
  deque<TH1F> histos ;
  histos.push_back(hplus) ;
  histos.push_back(hminus) ;
  histos.push_back(hr) ;
  return histos ;
}

Print()

void TimeBinning::Print	(	const std::string &	,
		std::ostream &	os = std::cout
	)		const

Definition at line 189 of file TimeBinning.cpp.

                                                             {
  os.precision(16) ;
  os << name << "_timeBins" ;
  for(const auto& bin : m_timeBins)
    os << "\t" << bin ;
  os << endl ;
  os << name << "_lifetime\t" << m_lifetime << endl ;
  string nameint = name + "_int" ;
  string nameplus = name + "_plus" ;
  string nameminus = name + "_minus" ;
  for(unsigned iTimeBin = 0 ; iTimeBin < nBinsTime() ; ++iTimeBin){
    m_binsInt[iTimeBin].Print(nameint, iTimeBin, os) ;
    string nameplusbin = Bin::getName(nameplus + "_time", iTimeBin) + "phase" ;
    string nameminusbin = Bin::getName(nameminus + "_time", iTimeBin) + "phase" ;
    for(unsigned iPhaseBin = 1 ; iPhaseBin < m_phaseBinning->nBins + 1; ++iPhaseBin){
      bin(iTimeBin, iPhaseBin).Print(nameplusbin, iPhaseBin, os) ;
      binBar(iTimeBin, iPhaseBin).Print(nameminusbin, iPhaseBin, os) ;
    }
  }
  m_phaseBinning->Print(name, os) ;
}

savePlotsVsTime()

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

Definition at line 326 of file TimeBinning.cpp.

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

setLifetime()

void TimeBinning::setLifetime

(

double

lifetime

)

Definition at line 345 of file TimeBinning.cpp.

                                             {
  m_lifetime = lifetime ;
  for(unsigned ibin = 0 ; ibin < nBinsTime() ; ++ibin){
    m_meant.at(ibin) = unmixedTimeMoment(ibin, lifetime, 1) ;
    m_meant2.at(ibin) = unmixedTimeMoment(ibin, lifetime, 2) ;
  }
}

timeBin()

int TimeBinning::timeBin

(

double

t

)

const

Definition at line 222 of file TimeBinning.cpp.

                                       {
  for(unsigned i = 0 ; i < nBinsTime() ; ++i)
    if(m_timeBins[i] <= t && t < m_timeBins[i+1])
      return i ;
  return -1 ;
}

unmixedTimeMoment()

double TimeBinning::unmixedTimeMoment ( unsigned  ibin, double  lifetime, int  moment  ) const

privatevirtual

Definition at line 334 of file TimeBinning.cpp.

                                                                                      {
  const double& tmin = m_timeBins.at(ibin) ;
  const double& tmax = m_timeBins.at(ibin+1) ;
  double norm = exp(-tmin/lifetime) - exp(-tmax/lifetime) ;
  double timeMoment = (pow(tmin, moment) * exp(-tmin/lifetime) - pow(tmax, moment) * exp(-tmax/lifetime))/norm ;
  if(moment == 0)
    return timeMoment ;
  timeMoment += moment * lifetime * unmixedTimeMoment(ibin, lifetime, moment-1) ;
  return timeMoment ;
}

write()

void TimeBinning::write	(	const std::string &	,
		const std::string &
	)		const

Definition at line 211 of file TimeBinning.cpp.

                                                                     {
  ofstream outfile ;
  outfile.open(fname) ;
  Print(name, outfile) ;
  outfile.close() ;
}

Member Data Documentation

m_bins

Bins2D TimeBinning::m_bins

private

Definition at line 76 of file TimeBinning.h.

m_binsBar

Bins2D TimeBinning::m_binsBar

private

Definition at line 77 of file TimeBinning.h.

m_binsInt

Bins TimeBinning::m_binsInt

private

Definition at line 78 of file TimeBinning.h.

m_lifetime

double TimeBinning::m_lifetime

private

Definition at line 75 of file TimeBinning.h.

m_meant

std::vector<double> TimeBinning::m_meant

private

Definition at line 73 of file TimeBinning.h.

m_meant2

std::vector<double> TimeBinning::m_meant2

private

Definition at line 74 of file TimeBinning.h.

m_phaseBinning

HadronicParameters::BinningPtr TimeBinning::m_phaseBinning

private

Definition at line 79 of file TimeBinning.h.

m_timeBins

std::vector<double> TimeBinning::m_timeBins

private

Definition at line 72 of file TimeBinning.h.

---

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

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