#include <Chi2BoxSet.h>

Inheritance diagram for Chi2BoxSet:

Public Member Functions
	Chi2BoxSet ()

	Chi2BoxSet (const DalitzEventPattern &pat, const MINT::counted_ptr< IIntegrationCalculator > &integPtr=(MINT::counted_ptr< IIntegrationCalculator >) 0)

	Chi2BoxSet (const DalitzArea &area, const MINT::counted_ptr< IIntegrationCalculator > &integPtr=(MINT::counted_ptr< IIntegrationCalculator >) 0)

	Chi2BoxSet (const MINT::PolymorphVector< Chi2Box > &other)

	Chi2BoxSet (const Chi2BoxSet &other)

virtual	~Chi2BoxSet ()

void	setIIntegrationCalculator (MINT::counted_ptr< IIntegrationCalculator > fap)

void	add (const Chi2Box &box)

void	add (const Chi2BoxSet &boxSet)

void	resetEventCounts ()

bool	addData (const IDalitzEvent &evt)

bool	addData (const IDalitzEvent *evt)

bool	addMC (IDalitzEvent &evt, double weight)

bool	addMC (IDalitzEvent *evt, double weight)

void	printBoxInfo (std::ostream &os=std::cout) const

void	print (std::ostream &os=std::cout) const

int	nData () const

int	nMC () const

double	weightedMC () const

double	weightedMC2 () const

double	rmsMC (int Ntotal) const

void	setNormFactor (double nf)

double	normFactor () const

void	setHistoColour (Color_t fcolor)

void	setFillColour (Color_t fcolor)

void	setLineColour (Color_t fcolor)

DalitzHistoSet &	histoData ()

const DalitzHistoSet &	histoData () const

DalitzHistoSet &	histoMC ()

const DalitzHistoSet &	histoMC () const

double	chi2 (double normFactorPassed=-1) const

Public Member Functions inherited from MINT::PolymorphVector< Chi2Box >
	PolymorphVector ()

	PolymorphVector (unsigned int N)

	PolymorphVector (unsigned int N, const Chi2Box &c)

	PolymorphVector (const PolymorphVector &other)

	PolymorphVector (const typename std::vector< Chi2Box > &other)

virtual	~PolymorphVector ()

std::vector< Chi2Box > &	theVector ()

const std::vector< Chi2Box > &	theVector () const

Chi2Box &	operator[] (unsigned int i)

const Chi2Box &	operator[] (unsigned int i) const

Chi2Box &	at (unsigned int i)

const Chi2Box &	at (unsigned int i) const

std::vector< Chi2Box >::iterator	begin ()

std::vector< Chi2Box >::const_iterator	begin () const

std::vector< Chi2Box >::iterator	end ()

std::vector< Chi2Box >::const_iterator	end () const

std::vector< Chi2Box >::iterator	find (const Chi2Box &c)

std::vector< Chi2Box >::const_iterator	find (const Chi2Box &c) const

Chi2Box &	front ()

const Chi2Box &	front () const

Chi2Box &	back ()

const Chi2Box &	back () const

unsigned int	size () const

bool	empty () const

void	push_back (const Chi2Box &c)

void	pop_back ()

void	erase (typename std::vector< Chi2Box >::iterator pos)

void	erase (typename std::vector< Chi2Box >::iterator first, typename std::vector< Chi2Box >::iterator last)

PolymorphVector< Chi2Box > &	operator= (const PolymorphVector< Chi2Box > &other)

void	clear ()

void	resize (unsigned int N)

void	resize (unsigned int N, const Chi2Box &c)

	operator const typename std::vector< Chi2Box > & () const

	operator typename std::vector< Chi2Box > & ()

bool	operator== (const MINT::PolymorphVector< Chi2Box > &v2) const

bool	operator!= (const MINT::PolymorphVector< Chi2Box > &v2) const

bool	operator< (const MINT::PolymorphVector< Chi2Box > &v2) const

bool	operator> (const MINT::PolymorphVector< Chi2Box > &v2) const

Private Attributes
MINT::counted_ptr< IIntegrationCalculator >	_integCalc

DalitzHistoSet	_histoData

DalitzHistoSet	_histoMC

double	_normFactor

Additional Inherited Members
Protected Attributes inherited from MINT::PolymorphVector< Chi2Box >
std::vector< Chi2Box >	_vec

Detailed Description

Definition at line 19 of file Chi2BoxSet.h.

Constructor & Destructor Documentation

◆ Chi2BoxSet() [1/5]

Chi2BoxSet::Chi2BoxSet ( )

Definition at line 10 of file Chi2BoxSet.cpp.

   : MINT::PolymorphVector<Chi2Box>()
   , _integCalc(0)
   , _normFactor(1)
 {
 }

◆ Chi2BoxSet() [2/5]

Chi2BoxSet::Chi2BoxSet	(	const DalitzEventPattern &	pat,
		const MINT::counted_ptr< IIntegrationCalculator > &	integPtr = `(MINT::counted_ptr<IIntegrationCalculator>)0`
	)

Definition at line 16 of file Chi2BoxSet.cpp.

   : MINT::PolymorphVector<Chi2Box>()
   , _integCalc(integPtr)
   , _normFactor(1)
 {
   Chi2Box b(pat);
   add(b);
 }

◆ Chi2BoxSet() [3/5]

Chi2BoxSet::Chi2BoxSet	(	const DalitzArea &	area,
		const MINT::counted_ptr< IIntegrationCalculator > &	integPtr = `(MINT::counted_ptr<IIntegrationCalculator>)0`
	)

Definition at line 25 of file Chi2BoxSet.cpp.

   : MINT::PolymorphVector<Chi2Box>()
   , _integCalc(integPtr)
   , _normFactor(1)
 {
   Chi2Box b(area);
   add(b);
 }

◆ Chi2BoxSet() [4/5]

Chi2BoxSet::Chi2BoxSet ( const MINT::PolymorphVector< Chi2Box > & other )

Definition at line 34 of file Chi2BoxSet.cpp.

   : MINT::PolymorphVector<Chi2Box>(other)
   , _integCalc(0)
   , _normFactor(1)
 {
 }

◆ Chi2BoxSet() [5/5]

Chi2BoxSet::Chi2BoxSet ( const Chi2BoxSet & other )

Definition at line 40 of file Chi2BoxSet.cpp.

   : MINT::PolymorphVector<Chi2Box>(other)
   , _integCalc(0)
   , _histoData(other._histoData)
   , _histoMC(other._histoMC)
   , _normFactor(other._normFactor)
 {
   if(0 != other._integCalc){
     counted_ptr<IIntegrationCalculator> newCP(other._integCalc->clone_IIntegrationCalculator());
     _integCalc = newCP;
   }
 }

◆ ~Chi2BoxSet()

virtual Chi2BoxSet::~Chi2BoxSet ( )

inlinevirtual

Definition at line 36 of file Chi2BoxSet.h.

36 {}

Member Function Documentation

◆ add() [1/2]

void Chi2BoxSet::add ( const Chi2Box & box )

Definition at line 53 of file Chi2BoxSet.cpp.

                                       {
   this->push_back(box);
 }

◆ add() [2/2]

void Chi2BoxSet::add ( const Chi2BoxSet & boxSet )

Definition at line 57 of file Chi2BoxSet.cpp.

                                             {
   for(unsigned int i=0; i<boxSet.size(); i++){
     this->add(boxSet[i]);
   }
 }

◆ addData() [1/2]

bool Chi2BoxSet::addData ( const IDalitzEvent & evt )

Definition at line 68 of file Chi2BoxSet.cpp.

                                                {
   for(unsigned int i=0; i < this->size(); i++){
     if ((*this)[i].addData(evt)){
       _histoData.addEvent(evt);
       return true;
     }
   }
   return false;
 }

◆ addData() [2/2]

bool Chi2BoxSet::addData ( const IDalitzEvent * evt )

Definition at line 77 of file Chi2BoxSet.cpp.

                                                {
   for(unsigned int i=0; i < this->size(); i++){
     if ((*this)[i].addData(evt)){
       _histoData.addEvent(*evt);
       return true;
     }
   }
   return false;
 }

◆ addMC() [1/2]

bool Chi2BoxSet::addMC	(	IDalitzEvent &	evt,
		double	weight
	)

Definition at line 86 of file Chi2BoxSet.cpp.

                                                       {
   for(unsigned int i=0; i < this->size(); i++){
     if((*this)[i].addMC(evt, weight)){
       _histoMC.addEvent(evt, weight);
       if(0 != _integCalc)_integCalc->addEvent(&evt);
       return true;
     }
   }
   return false;
 }

◆ addMC() [2/2]

bool Chi2BoxSet::addMC	(	IDalitzEvent *	evt,
		double	weight
	)

Definition at line 96 of file Chi2BoxSet.cpp.

                                                       {
   bool dbThis=false;
   if(dbThis) {
     cout << "Chi2BoxSet::addMC for pointers called with evt = " 
      << evt << endl;
   }
   for(unsigned int i=0; i < this->size(); i++){
     if((*this)[i].addMC(evt, weight)){
       _histoMC.addEvent(*evt, weight);
       if(0 != _integCalc)_integCalc->addEvent(evt);
       return true;
     }
   }
   return false;
 }

◆ chi2()

double Chi2BoxSet::chi2 ( double normFactorPassed = -1 ) const

Definition at line 209 of file Chi2BoxSet.cpp.

                                                     {
   double nf=normFactorPassed;
   if(nf < 0){
     nf = normFactor();
   }
 
   int n_data =  this->nData();
   double weight_mc = this->weightedMC() * nf;
 
   double var_mc = this->weightedMC2() * nf*nf;
   double varData_expected = weight_mc;
   
   double varData;
   if(0 != varData_expected) varData = varData_expected;
   else varData = n_data;
   
   double var = varData + var_mc;
 
   double delta_N = n_data - weight_mc;
   double dNSq = delta_N * delta_N;
   double chi2;
   if( dNSq < var * 1.e-20) chi2=0; // catches legit. 0-entry cases
   else chi2 = dNSq / var;
   
   return chi2;
 }

◆ histoData() [1/2]

DalitzHistoSet & Chi2BoxSet::histoData ( )

Definition at line 196 of file Chi2BoxSet.cpp.

                                      {
   return _histoData;
 }

◆ histoData() [2/2]

const DalitzHistoSet & Chi2BoxSet::histoData ( ) const

Definition at line 199 of file Chi2BoxSet.cpp.

                                                  {
   return _histoData;
 }

◆ histoMC() [1/2]

DalitzHistoSet & Chi2BoxSet::histoMC ( )

Definition at line 202 of file Chi2BoxSet.cpp.

                                    {
   return _histoMC;
 }

◆ histoMC() [2/2]

const DalitzHistoSet & Chi2BoxSet::histoMC ( ) const

Definition at line 205 of file Chi2BoxSet.cpp.

                                                {
   return _histoMC;
 }

◆ nData()

int Chi2BoxSet::nData ( ) const

Definition at line 127 of file Chi2BoxSet.cpp.

                            {
   int sum=0;
   for(unsigned int i=0; i < this->size(); i++){
     sum += (*this)[i].nData();
   }
   return sum;
 }

◆ nMC()

int Chi2BoxSet::nMC ( ) const

Definition at line 134 of file Chi2BoxSet.cpp.

                          {
   int sum=0;
   for(unsigned int i=0; i < this->size(); i++){
     sum += (*this)[i].nMC();
   }
   return sum;
 }

◆ normFactor()

double Chi2BoxSet::normFactor ( ) const

inline

Definition at line 64 of file Chi2BoxSet.h.

64 {return _normFactor;}

Chi2BoxSet::_normFactor

double _normFactor

Definition: Chi2BoxSet.h:24

◆ print()

void Chi2BoxSet::print ( std::ostream & os = std::cout ) const

inline

Definition at line 53 of file Chi2BoxSet.h.

53 { printBoxInfo(os); }

Chi2BoxSet::printBoxInfo

void printBoxInfo(std::ostream &os=std::cout) const

Definition: Chi2BoxSet.cpp:111

◆ printBoxInfo()

void Chi2BoxSet::printBoxInfo ( std::ostream & os = std::cout ) const

Definition at line 111 of file Chi2BoxSet.cpp.

                                                  {
   os << "Chi2BoxSet with " << this->size() << " sub-box";
   if(this->size() > 1) os << "es";
   os << ": ";
   for(unsigned int i=0; i < this->size(); i++){
     os << "\n " << i << ") " << (*this)[i];
   }
   if(0 != _integCalc){
     FitFractionList f = _integCalc->getFractions();
     f.sortByMagnitudeDecending();
     os << "\n sorted fractions in this box\n"
        << f << endl;
   }
   os << "\n--------------------------------------------------------\n" << endl;
 
 }

◆ resetEventCounts()

void Chi2BoxSet::resetEventCounts ( )

Definition at line 63 of file Chi2BoxSet.cpp.

                                  {
   for(unsigned int i=0; i < this->size(); i++){
     (*this)[i].resetEventCounts();
   }
 }

◆ rmsMC()

double Chi2BoxSet::rmsMC ( int Ntotal ) const

Definition at line 157 of file Chi2BoxSet.cpp.

                                   {
   return weightedMC2();
 }

◆ setFillColour()

void Chi2BoxSet::setFillColour ( Color_t fcolor )

Definition at line 187 of file Chi2BoxSet.cpp.

                                             {
   histoData().setFillColour(fcolor);
   histoMC().setFillColour(fcolor);
 }

◆ setHistoColour()

void Chi2BoxSet::setHistoColour ( Color_t fcolor )

Definition at line 183 of file Chi2BoxSet.cpp.

                                              {
   setFillColour(fcolor);
   setLineColour(fcolor);
 }

◆ setIIntegrationCalculator()

void Chi2BoxSet::setIIntegrationCalculator ( MINT::counted_ptr< IIntegrationCalculator > fap )

inline

Definition at line 38 of file Chi2BoxSet.h.

                                                                              {
     _integCalc=fap;
   }

◆ setLineColour()

void Chi2BoxSet::setLineColour ( Color_t fcolor )

Definition at line 191 of file Chi2BoxSet.cpp.

                                             {
   histoData().setLineColour(fcolor);
   histoMC().setLineColour(fcolor);
 }

◆ setNormFactor()

void Chi2BoxSet::setNormFactor ( double nf )

inline

Definition at line 63 of file Chi2BoxSet.h.

63 {_normFactor = nf;}

Chi2BoxSet::_normFactor

double _normFactor

Definition: Chi2BoxSet.h:24

◆ weightedMC()

double Chi2BoxSet::weightedMC ( ) const

Definition at line 141 of file Chi2BoxSet.cpp.

                                    {
   double sum=0;
   for(unsigned int i=0; i < this->size(); i++){
     sum += (*this)[i].weightedMC();
   }
   return sum;
 }

◆ weightedMC2()

double Chi2BoxSet::weightedMC2 ( ) const

Definition at line 148 of file Chi2BoxSet.cpp.

                                     {
   double sum=0;
   for(unsigned int i=0; i < this->size(); i++){
     sum += (*this)[i].weightedMC2();
   }
   return sum;
 }

Member Data Documentation

◆ _histoData

DalitzHistoSet Chi2BoxSet::_histoData

private

Definition at line 22 of file Chi2BoxSet.h.

◆ _histoMC

DalitzHistoSet Chi2BoxSet::_histoMC

private

Definition at line 22 of file Chi2BoxSet.h.

◆ _integCalc

MINT::counted_ptr<IIntegrationCalculator> Chi2BoxSet::_integCalc

private

Definition at line 20 of file Chi2BoxSet.h.

◆ _normFactor

double Chi2BoxSet::_normFactor

private

Definition at line 24 of file Chi2BoxSet.h.

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

Mint/Chi2BoxSet.h
src/Mojito/Chi2Test/Chi2BoxSet.cpp

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Chi2BoxSet() [1/5]

◆ Chi2BoxSet() [2/5]

◆ Chi2BoxSet() [3/5]

◆ Chi2BoxSet() [4/5]

◆ Chi2BoxSet() [5/5]

◆ ~Chi2BoxSet()

Member Function Documentation

◆ add() [1/2]

◆ add() [2/2]

◆ addData() [1/2]

◆ addData() [2/2]

◆ addMC() [1/2]

◆ addMC() [2/2]

◆ chi2()

◆ histoData() [1/2]

◆ histoData() [2/2]

◆ histoMC() [1/2]

◆ histoMC() [2/2]

◆ nData()

◆ nMC()

◆ normFactor()

◆ print()

◆ printBoxInfo()

◆ resetEventCounts()

◆ rmsMC()

◆ setFillColour()

◆ setHistoColour()

◆ setIIntegrationCalculator()

◆ setLineColour()

◆ setNormFactor()

◆ weightedMC()

◆ weightedMC2()

Member Data Documentation

◆ _histoData

◆ _histoMC

◆ _integCalc

◆ _normFactor