#include <FitAmpPairList.h>

Inheritance diagram for FitAmpPairList:

Public Member Functions
	FitAmpPairList ()

	FitAmpPairList (const FitAmpPairList &other)

virtual	~FitAmpPairList ()

MINT::counted_ptr< IIntegrationCalculator >	clone_IIntegrationCalculator () const

virtual void	addAmps (FitAmplitude a1, FitAmplitude a2)

virtual void	addAmps (FitAmplitude &a1, FitAmplitude &a2)

virtual void	addEvent (IDalitzEvent *evtPtr, double weight=1)

virtual void	addEvent (IDalitzEvent &evt, double weight=1)

virtual void	addEvent (MINT::counted_ptr< IDalitzEvent > evtPtr, double weight=1)

virtual void	reAddEvent (IDalitzEvent &evt, double weight=1)

bool	isCompatibleWith (const FitAmpPairList &other) const

virtual bool	add (const FitAmpPairList &otherList)

virtual bool	add (const FitAmpPairList *otherListPtr)

virtual bool	add (MINT::const_counted_ptr< FitAmpPairList > otherListPtr)

virtual bool	append (const FitAmpPairList &otherListPtr)

virtual bool	append (const FitAmpPairList *otherListPtr)

virtual bool	append (MINT::const_counted_ptr< FitAmpPairList > otherListPtr)

virtual int	numEvents () const

virtual double	integral () const

std::complex< double >	ComplexIntegralForTags (int tag1, int tag2) const

double	integralForMatchingPatterns (bool match, int pattern_sign) const

std::complex< double >	ComplexSumForMatchingPatterns (bool match) const

virtual double	variance () const

virtual double	sumOfVariances () const

virtual std::complex< double >	ComplexSum () const

virtual void	Gradient (MINT::MinuitParameterSet *mps, std::vector< double > &grad)

virtual void	GradientForLasso (MINT::MinuitParameterSet *mps, std::vector< double > &grad)

double	sumOfSqrtFitFractions ()

double	absSumOfInterferenceFractions ()

double	absSumOfSqrtInterferenceFractions ()

double	sumOfFitFractions ()

int	numberOfFitFractionsLargerThanThreshold (double threshold)

FitFractionList	getFractions () const

FitFractionList	getInterferenceTerms () const

bool	doFractions ()

void	setEfficiency (MINT::counted_ptr< MINT::IReturnRealForEvent< IDalitzEvent > > eff)

void	unsetEfficiency ()

double	efficiency (IDalitzEvent *evtPtr)

bool	haveEfficiency () const

virtual bool	makeAndStoreFractions (MINT::Minimiser *mini=0)

virtual bool	makeAndStoreFractions (const std::string &fname, const std::string &fnameROOT, MINT::Minimiser *min=0)

virtual double	getFractionChi2 () const

virtual DalitzHistoSet	histoSet () const

virtual DalitzHistoSet	un_normalised_histoSetRe () const

virtual DalitzHistoSet	un_normalised_histoSetIm () const

void	saveEachAmpsHistograms (const std::string &prefix) const

std::vector< DalitzHistoSet >	GetEachAmpsHistograms ()

DalitzHistoSet	interferenceHistoSet () const

void	saveInterferenceHistograms (const std::string &prefix) const

std::vector< DalitzHistoSet >	GetInterferenceHistograms ()

virtual void	doFinalStats (MINT::Minimiser *min=0)

virtual void	doFinalStatsAndSave (MINT::Minimiser *min=0, const std::string &fname="FitAmpResults.txt", const std::string &fnameROOT="fitFractions.root")

virtual bool	save (const std::string &asSubdirOf=".") const

virtual bool	retrieve (const std::string &asSubdirOf=".")

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

bool	needToReIntegrate () const

void	startIntegration ()

void	startReIntegration ()

void	startForcedReIntegration ()

void	endIntegration ()

void	setSlow ()

void	setFast ()

bool	slow () const

bool	fast () const

FitAmpPairList &	operator+= (const FitAmpPairList &other)

FitAmpPairList	operator+ (const FitAmpPairList &other) const

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

	PolymorphVector (unsigned int N)

	PolymorphVector (unsigned int N, const FitAmpPair &c)

	PolymorphVector (const PolymorphVector &other)

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

virtual	~PolymorphVector ()

std::vector< FitAmpPair > &	theVector ()

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

FitAmpPair &	operator[] (unsigned int i)

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

FitAmpPair &	at (unsigned int i)

const FitAmpPair &	at (unsigned int i) const

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

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

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

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

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

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

FitAmpPair &	front ()

const FitAmpPair &	front () const

FitAmpPair &	back ()

const FitAmpPair &	back () const

unsigned int	size () const

bool	empty () const

void	push_back (const FitAmpPair &c)

void	pop_back ()

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

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

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

void	clear ()

void	resize (unsigned int N)

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

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

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

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

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

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

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

Public Member Functions inherited from IIntegrationCalculator
virtual	~IIntegrationCalculator ()

Private Member Functions
void	applyHistoOption ()

double	phaseSpaceIntegral () const

std::string	dirName () const

bool	makeDirectory (const std::string &asSubdirOf=".") const

virtual double	oldVariance () const

bool	reset ()

Private Attributes
MINT::NamedParameter< std::string >	HistoOption

int	_Nevents

double	_sum

double	_sumsq

double	_psSum

double	_psSumSq

bool	_slow

FitAmpPairCovariance	_cov

MINT::counted_ptr< MINT::IReturnRealForEvent< IDalitzEvent > >	_efficiency

FitFractionList	_singleAmpFractions

FitFractionList	_interferenceFractions

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

Detailed Description

Definition at line 27 of file FitAmpPairList.h.

Constructor & Destructor Documentation

◆ FitAmpPairList() [1/2]

FitAmpPairList::FitAmpPairList ( )

Definition at line 50 of file FitAmpPairList.cpp.

   : MINT::PolymorphVector<FitAmpPair>()
   , _Nevents(0)
   , _sum(0)
   , _sumsq(0)
   , _psSum(0)
   , _psSumSq(0)
   , _slow(true)
   , _cov(this)
   , _efficiency(0)
   , HistoOption("FitAmpPairList::HistoOption", (std::string) "default")
 {
   applyHistoOption();
 }

◆ FitAmpPairList() [2/2]

FitAmpPairList::FitAmpPairList ( const FitAmpPairList & other )

Definition at line 64 of file FitAmpPairList.cpp.

   : IIntegrationCalculator()
   , MINT::PolymorphVector<FitAmpPair>(other)
   , _Nevents(other._Nevents)
   , _sum(other._sum)
   , _sumsq(other._sumsq)
   , _psSum(other._psSum)
   , _psSumSq(other._psSumSq)
   , _slow(other._slow)
   , _cov(other._cov, this)
   , _efficiency(other._efficiency)
   , HistoOption("FitAmpPairList::HistoOption", (std::string) "default")
 {
   applyHistoOption();
 }

◆ ~FitAmpPairList()

virtual FitAmpPairList::~FitAmpPairList ( )

inlinevirtual

Definition at line 60 of file FitAmpPairList.h.

60 {}

Member Function Documentation

◆ absSumOfInterferenceFractions()

double FitAmpPairList::absSumOfInterferenceFractions ( )

Definition at line 430 of file FitAmpPairList.cpp.

                                                      {
     double sum = 0;
     for(unsigned int i=0; i< this->size(); i++){
         if(!this->at(i).isSingleAmp())sum += abs(this->at(i).integral());
     }
     return sum/integral();
 }

◆ absSumOfSqrtInterferenceFractions()

double FitAmpPairList::absSumOfSqrtInterferenceFractions ( )

Definition at line 422 of file FitAmpPairList.cpp.

                                                          {
     double sum = 0;
     for(unsigned int i=0; i< this->size(); i++){
         if(!this->at(i).isSingleAmp())sum += sqrt(abs(this->at(i).integral()));
     }
     return sum;
 }

◆ add() [1/3]

bool FitAmpPairList::add ( const FitAmpPairList & otherList )

virtual

Definition at line 177 of file FitAmpPairList.cpp.

                                                        {
   if(! isCompatibleWith(otherList)){
     cout << "ERROR in FitAmpPairList::add "
      << "trying to add incompatible lists"
      << endl;
     return false;
   }
 
   _Nevents  += otherList._Nevents;
   _sum      += otherList._sum;
   _sumsq    += otherList._sumsq;
   _psSum    += otherList._psSum;
   _psSumSq  += otherList._psSumSq;
 
   _cov      += otherList._cov;
 
   for(unsigned int i=0; i< this->size(); i++){
     this->at(i) += otherList[i];
   }
   
   return true;
 
 }

◆ add() [2/3]

bool FitAmpPairList::add ( const FitAmpPairList * otherListPtr )

virtual

Definition at line 165 of file FitAmpPairList.cpp.

                                                           {
   if(0 != otherListPtr){
     return add(*otherListPtr);
   }
   return true;
 }

◆ add() [3/3]

bool FitAmpPairList::add ( MINT::const_counted_ptr< FitAmpPairList > otherListPtr )

virtual

Definition at line 171 of file FitAmpPairList.cpp.

                                                                       {
   if(0 != otherListPtr){
     return add(*otherListPtr);
   }
   return true;
 }

◆ addAmps() [1/2]

void FitAmpPairList::addAmps	(	FitAmplitude *	a1,
		FitAmplitude *	a2
	)

virtual

Definition at line 85 of file FitAmpPairList.cpp.

                                                               {
   if(0 == a1 || 0 == a2){
     cout << "Error in FitAmpPairList::addAmps: " << a1 << ", " << a2 << endl;
     throw "zero pointer";
   }
   FitAmpPair amps(*a1, *a2);
   this->push_back(amps);
 }

◆ addAmps() [2/2]

void FitAmpPairList::addAmps	(	FitAmplitude &	a1,
		FitAmplitude &	a2
	)

virtual

Definition at line 93 of file FitAmpPairList.cpp.

                                                               {
   FitAmpPair amps(a1, a2);
   this->push_back(amps);
 }

◆ addEvent() [1/3]

void FitAmpPairList::addEvent	(	IDalitzEvent *	evtPtr,
		double	weight = `1`
	)

virtual

Implements IIntegrationCalculator.

Definition at line 98 of file FitAmpPairList.cpp.

                                                                 {
   if(0 == evtPtr) return;
   addEvent(*evtPtr, weight);
 }

◆ addEvent() [2/3]

void FitAmpPairList::addEvent	(	IDalitzEvent &	evt,
		double	weight = `1`
	)

virtual

Definition at line 103 of file FitAmpPairList.cpp.

                                                              {
   bool dbThis=false;
   static unsigned int nZeroPs=0;
   _Nevents++;
   double x=0;
   double ps = evt.phaseSpace();
   if(ps <= 0.0){
       if (!(nZeroPs%1000))
       {
         cout << "WARNING in FitAmpPairList::addToHistograms"
          << " " << ++nZeroPs << " th"
          << " event with phase space = " << ps << endl;
         return; // should not happen.
       }
   }
 
   for(unsigned int i=0; i< this->size(); i++){
     //    x += this->at(i).add(evtPtr, weight*efficiency(evtPtr));
     x += this->at(i).add(evt, weight, efficiency(&evt));
   }
   if(dbThis) cout << "FitAmpPairList::addEvent(): adding event" << endl;
   //if(_cov.isValid()) 
   _cov.addLastEventFromList();
 
   _sum   += x;
   _sumsq += x*x;
   
   //if(0 == evtPtr) return;
   //double ps = evtPtr->phaseSpace();
   _psSum   += ps;
   _psSumSq += ps*ps;
 
   return;
 }

◆ addEvent() [3/3]

void FitAmpPairList::addEvent	(	MINT::counted_ptr< IDalitzEvent >	evtPtr,
		double	weight = `1`
	)

virtual

Implements IIntegrationCalculator.

Definition at line 137 of file FitAmpPairList.cpp.

                                                                             {
   return addEvent(evtPtr.get(), weight);
 }

◆ append() [1/3]

bool FitAmpPairList::append ( const FitAmpPairList & otherListPtr )

virtual

Definition at line 207 of file FitAmpPairList.cpp.

                                                           {
   // WARNING: this will reset everything to zero.
 
   for(unsigned int i=0; i< otherList.size(); i++){
     this->push_back(otherList[i]);
   }
   this->reset();
 
   return true;
 
 }

◆ append() [2/3]

bool FitAmpPairList::append ( const FitAmpPairList * otherListPtr )

virtual

Definition at line 204 of file FitAmpPairList.cpp.

                                                              {
   return append(*otherListPtr);
 }

◆ append() [3/3]

bool FitAmpPairList::append ( MINT::const_counted_ptr< FitAmpPairList > otherListPtr )

virtual

Definition at line 201 of file FitAmpPairList.cpp.

                                                                              {
   return append(*otherListPtr);
 }

◆ applyHistoOption()

void FitAmpPairList::applyHistoOption ( )

private

Definition at line 45 of file FitAmpPairList.cpp.

                                      {
   if((string)HistoOption == (string)"fast") setFast();
   if((string)HistoOption == (string)"slow") setSlow();
 }

◆ clone_IIntegrationCalculator()

counted_ptr< IIntegrationCalculator > FitAmpPairList::clone_IIntegrationCalculator ( ) const

virtual

Implements IIntegrationCalculator.

Definition at line 80 of file FitAmpPairList.cpp.

                                                   {
   counted_ptr<IIntegrationCalculator> ptr(new FitAmpPairList(*this));
   return ptr;
 }

◆ ComplexIntegralForTags()

std::complex< double > FitAmpPairList::ComplexIntegralForTags	(	int	tag1,
		int	tag2
	)		const

Definition at line 242 of file FitAmpPairList.cpp.

                                                                                  {
     std::complex<double> sum = 0;
 
     for(unsigned int i=0; i< this->size(); i++){
         if((this->at(i).fitAmp1().getTag() == tag1) && ( this->at(i).fitAmp2().getTag() == tag2)){
             sum += conj(this->at(i).complexVal());       
         }
         else if((this->at(i).fitAmp1().getTag() == tag2) && ( this->at(i).fitAmp2().getTag() == tag1)){
             sum += this->at(i).complexVal();
         }
     }
     return sum;
 }

◆ ComplexSum()

std::complex< double > FitAmpPairList::ComplexSum ( ) const

virtual

Definition at line 281 of file FitAmpPairList.cpp.

                                                    { //laurenPsuedo
   bool dbThis = false;
   std::complex<double> sum = 0;
   for(unsigned int i=0; i< this->size(); i++){
     sum += this->at(i).complexVal();
   }
   return sum;
   (void)dbThis;
 }

◆ ComplexSumForMatchingPatterns()

std::complex< double > FitAmpPairList::ComplexSumForMatchingPatterns ( bool match ) const

Definition at line 269 of file FitAmpPairList.cpp.

                                                                                 {
     std::complex<double> sum = 0;
     for(unsigned int i=0; i< this->size(); i++){
         if(this->at(i).hasMatchingPattern()==match){
             std::complex<double> val = this->at(i).complexVal();
             if(this->at(i).fitAmp1().theBareDecay().getVal() >  this->at(i).fitAmp2().theBareDecay().getVal() ) val = conj(val);
             sum += val;
         }
     }
     return sum;
 }

◆ dirName()

std::string FitAmpPairList::dirName ( ) const

private

Definition at line 1099 of file FitAmpPairList.cpp.

                                        {
   return "FitAmpPairList";
 }

◆ doFinalStats()

void FitAmpPairList::doFinalStats ( MINT::Minimiser * min = 0 )

virtual

Implements IIntegrationCalculator.

Definition at line 653 of file FitAmpPairList.cpp.

                                                 {
   bool dbThis=true;
   if(dbThis) cout << "FitAmpPairList::doFinalStats() called" << endl;
   makeAndStoreFractions(mini);
 }

◆ doFinalStatsAndSave()

virtual void FitAmpPairList::doFinalStatsAndSave	(	MINT::Minimiser *	min = `0`,
		const std::string &	fname = `"FitAmpResults.txt"`,
		const std::string &	fnameROOT = `"fitFractions.root"`
	)

inlinevirtual

Definition at line 129 of file FitAmpPairList.h.

                                                                                                                                                      {
       makeAndStoreFractions(fname,fnameROOT,min);
   }

◆ doFractions()

bool FitAmpPairList::doFractions ( )

Definition at line 660 of file FitAmpPairList.cpp.

                                 {
 
   if(this->empty()) return 0;
 
    _singleAmpFractions = FitFractionList() ;
    _interferenceFractions = FitFractionList();
 
   cout << "\n============================================"
        << "============================================"
        << "\n        Amplitude Fractions";
 
   double norm = integral();
 
   if(norm <= 0){
     cout << "ERROR in FitAmpPairList::makeAndStoreFractions()"
      << " integral = " << integral()
      << " won't do fractions."
      << endl;
     return false;
   }
 
   for(unsigned int i=0; i < this->size(); i++){    
     double frac = this->at(i).integral()/norm;
     string name;
     if(this->at(i).isSingleAmp()){
       name = this->at(i).fitAmp1().name();
     }else{
 //       name = this->at(i).name();
       name = this->at(i).fitAmp1().name() + "_x_" + this->at(i).fitAmp2().name();
     }
     FitFraction f(name, frac);
 
     if(this->at(i).isSingleAmp()){
       this->at(i).fitAmp1().setFraction(frac);
       _singleAmpFractions.add(f);
     }else{
       _interferenceFractions.add(f);
     }
   }
 
   cout << "filled FitFractionLists" << endl;
   FitFractionList interferenceFractionsSorted(_interferenceFractions);
   interferenceFractionsSorted.sortByMagnitudeDecending();
   
   cout <<   "================================================="
        << "\n================================================="
        << "\n FRACTIONS:"
        << "\n ^^^^^^^^^^" << endl;
   cout << _singleAmpFractions << endl;
   cout <<   "================================================="
        << "\n================================================="
        << "\n Interference terms (sorted by size)"
        << "\n ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^" << endl;
   cout << interferenceFractionsSorted << endl;
   cout <<   "================================================="
        << "\n=================================================" << endl;
 
   cout << "\n\n\t X-check: total sum of all terms"
        << " (amplitude fractions + interference terms): "
        << _singleAmpFractions.sum().frac() 
     + _interferenceFractions.sum().frac();
 
   return true;
 }

◆ efficiency()

double FitAmpPairList::efficiency ( IDalitzEvent * evtPtr )

Definition at line 529 of file FitAmpPairList.cpp.

                                                      {
   if(0 == _efficiency) return 1;
   if(0 == evtPtr) return 0;
   return _efficiency->RealVal(*evtPtr);
 }

◆ endIntegration()

void FitAmpPairList::endIntegration ( )

Definition at line 1218 of file FitAmpPairList.cpp.

                                    {
   for(unsigned int i=0; i < this->size(); i++){    
     this->at(i).endIntegration();
   }
 }

◆ fast()

bool FitAmpPairList::fast ( ) const

inline

Definition at line 150 of file FitAmpPairList.h.

150 {return !_slow;}

FitAmpPairList::_slow

bool _slow

Definition: FitAmpPairList.h:42

◆ GetEachAmpsHistograms()

std::vector< DalitzHistoSet > FitAmpPairList::GetEachAmpsHistograms ( )

Definition at line 576 of file FitAmpPairList.cpp.

                                                                {
     std::vector<DalitzHistoSet> HistoSet;
      DalitzHistoSet sum;
       int counter=0;
       for(unsigned int i=0; i< this->size(); i++){
         if(this->at(i).isSingleAmp()){
           counter++;
           std::string name =  anythingToString(counter) + ".root";
               double frac = this->at(i).integral()/integral();  
           DalitzHistoSet hs(this->at(i).histoSet());
           std::string title = this->at(i).fitAmp1().name();
           hs.setTitle(title);
           cout << "FitAmpPairList::saveEachAmpsHistograms: "
            << "saving " << title << " as " << name << endl;
               hs.scale(frac/hs.integral());  
           HistoSet.push_back(hs);
         }
       }
       return HistoSet;
 }

◆ getFractionChi2()

double FitAmpPairList::getFractionChi2 ( ) const

virtual

Implements IIntegrationCalculator.

Definition at line 1067 of file FitAmpPairList.cpp.

                                             {
   bool ErrorProportionalToSqrtTarget = true;
   double norm = integral();
   if(norm <= 0){
     cout << "ERROR in FitAmpPairList::getFractionChi2()"
      << " integral = " << integral()
      << " won't do fractions. Returning " << 9999.0e30
      << endl;
     return 9999.0e30;
   }
 
   double canonicalError=0.001;
 
   int counter=0;
   double sum=0;
   for(unsigned int i=0; i < this->size(); i++){
     if(this->at(i).isSingleAmp()){
       counter++;
       double frac = this->at(i).integral()/norm;
       double target_f = this->at(i).fitAmp1().getFraction();
       //if(i < 10){
     //cout << "(" << frac << " - " << target_f << ")" << endl;
       //}
       double dfs=(frac - target_f)/canonicalError;
       if(ErrorProportionalToSqrtTarget && target_f > 0) dfs /= sqrt(0.01*target_f);
       sum += dfs*dfs;
     }
   }
   return sum;
 }

◆ getFractions()

FitFractionList FitAmpPairList::getFractions ( ) const

inlinevirtual

Implements IIntegrationCalculator.

Definition at line 104 of file FitAmpPairList.h.

104 {return _singleAmpFractions;}

FitAmpPairList::_singleAmpFractions

FitFractionList _singleAmpFractions

Definition: FitAmpPairList.h:48

◆ GetInterferenceHistograms()

std::vector< DalitzHistoSet > FitAmpPairList::GetInterferenceHistograms ( )

Definition at line 627 of file FitAmpPairList.cpp.

                                                                    {
     std::vector<DalitzHistoSet> HistoSet;
     DalitzHistoSet sum;
     int counter=0;
     for(unsigned int i=0; i< this->size(); i++){
         if(!this->at(i).isSingleAmp()){
             counter++;
             std::string name =  anythingToString(counter) + ".root";
             double frac = this->at(i).integral()/integral();  
             DalitzHistoSet hs(this->at(i).histoSet());
             std::string title = this->at(i).name();
             hs.setTitle(title);
             cout << "FitAmpPairList::saveEachAmpsHistograms: "
             << "saving " << title << " as " << name << endl;
             hs.scale(frac/hs.integral());  
             HistoSet.push_back(hs);
         }
     }
     return HistoSet;
 }

◆ getInterferenceTerms()

FitFractionList FitAmpPairList::getInterferenceTerms ( ) const

inline

Definition at line 105 of file FitAmpPairList.h.

105 {return _interferenceFractions;}

FitAmpPairList::_interferenceFractions

FitFractionList _interferenceFractions

Definition: FitAmpPairList.h:48

◆ Gradient()

void FitAmpPairList::Gradient	(	MINT::MinuitParameterSet *	mps,
		std::vector< double > &	grad
	)

virtual

Definition at line 291 of file FitAmpPairList.cpp.

                                                                              {
     
   for (unsigned int i=0; i<mps->size(); i++) {
     if(mps->getParPtr(i)->hidden())continue;
     
     if(i+1 >= grad.size()){
       cout << "WARNING in FitAmpPairList::Gradient:"
        << " have to increase size of grad to avoid memory issues" << endl;
       grad.resize(i+2);
     }
  
     grad.at(i)=0; grad.at(i+1)=0;
 
     string name_i= mps->getParPtr(i)->name();
     if(name_i.find("_Re")!=std::string::npos){
       if(mps->getParPtr(i)->iFixInit() && mps->getParPtr(i+1)->iFixInit()){
     i++;
     continue;
       }
       name_i.replace(name_i.find("_Re"),3,"");
       complex<double> sum(0);
       for(unsigned int j=0; j< this->size(); j++){
     if(!A_is_in_B(name_i,this->at(j).name())) continue;
     if(A_is_in_B("Inco", name_i) != A_is_in_B("Inco",this->at(j).name()) ) continue;
     double singleAmpCorrection= 1.;
     if(this->at(j).isSingleAmp()) singleAmpCorrection = 2.;
     // 2 a_j^* A_i A_j^*
     if(A_is_in_B(name_i,this->at(j).fitAmp1().name())){
       sum += singleAmpCorrection* this->at(j).valNoFitPars()* conj(this->at(j).fitAmp2().AmpPhase());
     }
     else {
       sum += singleAmpCorrection* conj( this->at(j).valNoFitPars()* this->at(j).fitAmp1().AmpPhase() );
     }
         
       }
       grad.at(i)= sum.real();
       grad.at(i+1)= -sum.imag();
       i++;
     } 
     // Doesn't work. Don't use! 
     /* 
        else if(name_i.find("_Amp")!=std::string::npos){
        name_i.replace(name_i.find("_Amp"),4,""); 
        complex<double> sumAmp(0);
        complex<double> sumPhase(0);
        
        for(unsigned int j=0; j< this->size(); j++){
        if(!A_is_in_B(name_i,this->at(j).name())) continue;
        double singleAmpCorrection= 1.;
        if(this->at(j).isSingleAmp()) singleAmpCorrection = 2.;
        // 2 a_j^* A_i A_j^*
        if(A_is_in_B(name_i,this->at(j).fitAmp1().name())){
        sumAmp += singleAmpCorrection*this->at(j).complexVal()/std::abs(this->at(j).fitAmp1().AmpPhase());
        if(!this->at(j).isSingleAmp())sumPhase += this->at(j).complexVal()*std::arg(this->at(j).fitAmp1().AmpPhase());
        }
        else {
        sumAmp += singleAmpCorrection*conj(this->at(j).complexVal())/std::abs(this->at(j).fitAmp2().AmpPhase());
        if(!this->at(j).isSingleAmp())sumPhase += conj(this->at(j).complexVal())*std::arg(this->at(j).fitAmp2().AmpPhase());
        }
        
        }
        grad[i]= sumAmp.real();
        grad[i+1]= -sumPhase.imag();
        i++;
        }
     */
     else if(mps->getParPtr(i)->iFixInit())continue;
     else {
       std::cout << "FitAmpPairList::Gradient() called. Sorry, I don't know how to calculate the derivative with respect to the fit parameter " << mps->getParPtr(i)->name() << " ! Please implement me or set useAnalytic Gradient to 0 in your options file. I'll crash now. " << std::endl;
             throw "crash";
         }
     }
     
 }

◆ GradientForLasso()

void FitAmpPairList::GradientForLasso	(	MINT::MinuitParameterSet *	mps,
		std::vector< double > &	grad
	)

virtual

Definition at line 366 of file FitAmpPairList.cpp.

                                                 {
     
     for (unsigned int i=0; i<mps->size(); i++) {
         if(mps->getParPtr(i)->hidden())continue;
 
         grad[i]=0; grad[i+1]=0;
         string name_i= mps->getParPtr(i)->name();
         if(name_i.find("_Re")!=std::string::npos){
             if(mps->getParPtr(i)->iFixInit() && mps->getParPtr(i+1)->iFixInit())continue;
             name_i.replace(name_i.find("_Re"),3,"");
             for(unsigned int j=0; j< this->size(); j++){
                 if(!A_is_in_B(name_i,this->at(j).name())) continue;
                 if(A_is_in_B("Inco", name_i) != A_is_in_B("Inco",this->at(j).name()) ) continue;
                 if(!this->at(j).isSingleAmp()) continue;
 
         if(i+1 >= grad.size()){
           cout << "WARNING in FitAmpPairList::GradientForLasso:"
                << " have to increase size of grad to avoid memory issues" << endl;
           grad.resize(i+2);
     }
 
                 double integral = this->at(j).valNoFitPars().real()/sqrt(this->at(j).integral());
                 // Re(a_j) A_j A_j^* dphi
                 grad[i]= integral * this->at(j).fitAmp1().AmpPhase().real();
                 // Im(a_j) A_j A_j^* dphi
                 grad[i+1]= integral * this->at(j).fitAmp1().AmpPhase().imag();
                 i++;
                 break;
             }
         }
         else if(mps->getParPtr(i)->iFixInit())continue;
         else {
             std::cout << "Sorry, I don't know how to calculate the derivative with respect to the fit parameter " << mps->getParPtr(i)->name() << " ! Please implement me or set useAnalytic Gradient to 0 in your options file. I'll crash now. " << std::endl;
             throw "crash";
         }
     }
     
 }

◆ haveEfficiency()

bool FitAmpPairList::haveEfficiency ( ) const

inline

Definition at line 111 of file FitAmpPairList.h.

111 {return 0 != _efficiency;}

FitAmpPairList::_efficiency

MINT::counted_ptr< MINT::IReturnRealForEvent< IDalitzEvent > > _efficiency

Definition: FitAmpPairList.h:46

◆ histoSet()

DalitzHistoSet FitAmpPairList::histoSet ( ) const

virtual

Implements IIntegrationCalculator.

Definition at line 540 of file FitAmpPairList.cpp.

                                              {
   return un_normalised_histoSetRe()/integral();
 }

◆ integral()

double FitAmpPairList::integral ( ) const

virtual

Implements IIntegrationCalculator.

Definition at line 234 of file FitAmpPairList.cpp.

                                      {
   double sum = 0;
   for(unsigned int i=0; i< this->size(); i++){
     sum += this->at(i).integral();
   }
   return sum;
 }

◆ integralForMatchingPatterns()

double FitAmpPairList::integralForMatchingPatterns	(	bool	match,
		int	pattern_sign
	)		const

Definition at line 256 of file FitAmpPairList.cpp.

                                                                                     {
     double sum = 0;
     for(unsigned int i=0; i< this->size(); i++){
         if(this->at(i).hasMatchingPattern()==match){
             if(match){
                 if(static_cast<double>(pattern_sign) * this->at(i).fitAmp1().theBareDecay().getVal() < 0 ) continue;
             }
             sum += this->at(i).integral();
         }
     }
     return sum;
 }

◆ interferenceHistoSet()

DalitzHistoSet FitAmpPairList::interferenceHistoSet ( ) const

Definition at line 597 of file FitAmpPairList.cpp.

                                                          {
     DalitzHistoSet sum;
     for(unsigned int i=0; i< this->size(); i++){
         if(!this->at(i).isSingleAmp())sum += this->at(i).histoSet();
     }
     sum /= integral();
     //  if(_Nevents > 0) sum /= (double) _Nevents;
     // above two lines normalise this to 1.
     
     return sum;
 }

◆ isCompatibleWith()

bool FitAmpPairList::isCompatibleWith ( const FitAmpPairList & other ) const

Definition at line 158 of file FitAmpPairList.cpp.

                                                                           {
   if(this->size() != otherList.size()) return false;
   for(unsigned int i=0; i < this->size(); i++){
     if(this->at(i).name() != otherList[i].name()) return false;
   }
   return true;
 }

◆ makeAndStoreFractions() [1/2]

virtual bool FitAmpPairList::makeAndStoreFractions ( MINT::Minimiser * mini = 0 )

inlinevirtual

Implements IIntegrationCalculator.

Definition at line 113 of file FitAmpPairList.h.

113 {

114 return makeAndStoreFractions("FitAmpResults.txt","fitFractions.root",mini);}

FitAmpPairList::makeAndStoreFractions

virtual bool makeAndStoreFractions(MINT::Minimiser *mini=0)

Definition: FitAmpPairList.h:113

◆ makeAndStoreFractions() [2/2]

bool FitAmpPairList::makeAndStoreFractions	(	const std::string &	fname,
		const std::string &	fnameROOT,
		MINT::Minimiser *	min = `0`
	)

virtual

latex table

Definition at line 726 of file FitAmpPairList.cpp.

                         {
   bool dbThis=true;
 
    _singleAmpFractions = FitFractionList() ;
    _interferenceFractions = FitFractionList();
 
   if(this->empty()) return 0;
   counted_ptr<FitAmpPairFitCovariance> fcov(0);
   if(0 != mini){
     if(mini->parSet() != this->at(0).fitAmp1().FitAmpPhase().p1().parSet()){
       cout << "ERROR in FitAmpPairList::makeAndStoreFractions"
        << "\n Minuit parameter set and fit parameters incompatible"
        << endl;
     }else{
       cout << "getting minuitCov" << endl;
       TMatrixTSym<double> minuitCov(mini->covMatrixFull());
       cout << "got minuitCov" << endl;
       counted_ptr<FitAmpPairFitCovariance> 
     fc(new FitAmpPairFitCovariance(this, minuitCov));
       fcov = fc;
       cout << "got fcov" << endl;
     }
   }
 
   bool printFitErrors = (fcov != 0 && fcov->isValid());
 
   cout << "\n============================================"
        << "============================================"
        << "\n        Amplitude Fractions";
   if(0 != fcov) cout << " +/- fit error ";
   cout << " +/- integration error "
        << endl;
   double norm = integral();
   if(norm <= 0){
     cout << "ERROR in FitAmpPairList::makeAndStoreFractions()"
      << " integral = " << integral()
      << " won't do fractions."
      << endl;
     return false;
   }
 
   for(unsigned int i=0; i < this->size(); i++){    
     double frac = this->at(i).integral()/norm;
     string name;
     if(this->at(i).isSingleAmp()){
       name = this->at(i).fitAmp1().name();
     }else{
       name = this->at(i).name();
     }
     FitFraction f(name, frac);
     if(printFitErrors){
       f.sigmaFit() = fcov->getFractionError(i);
     }
     if(_cov.isValid()){ 
       f.sigmaInteg() = _cov.getFractionError(i);
     }
 
     if(this->at(i).isSingleAmp()){
       this->at(i).fitAmp1().setFraction(frac);
       _singleAmpFractions.add(f);
     }else{
       _interferenceFractions.add(f);
     }
   }
 
   if(dbThis)cout << "filled FitFractionLists" << endl;
   if(printFitErrors){
     _singleAmpFractions.setSumFitError(fcov->getFractionSumError());
     _interferenceFractions.setSumFitError(fcov->getInterferenceFracSumError());
   }
   if(_cov.isValid()){
     _singleAmpFractions.setSumIntegError(_cov.getFractionSumError());
   }
 
   if(dbThis) cout << "... now including errors" << endl;
   if(dbThis)cout  << " sorting" << endl;
   _interferenceFractions.sortByMagnitudeDecending();
   if(dbThis)cout << "sorted" << endl;
 
   cout <<   "================================================="
        << "\n================================================="
        << "\n FRACTIONS:"
        << "\n ^^^^^^^^^^" << endl;
   cout << _singleAmpFractions << endl;
   cout <<   "================================================="
        << "\n================================================="
        << "\n Interference terms (sorted by size)"
        << "\n ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^" << endl;
   cout << _interferenceFractions << endl;
   cout <<   "================================================="
        << "\n=================================================" << endl;
 
   cout << "\n\n\t X-check: total sum of all terms"
        << " (amplitude fractions + interference terms): "
        << _singleAmpFractions.sum().frac() 
     + _interferenceFractions.sum().frac();
   cout << endl;
   ofstream os(fname.c_str(),std::ofstream::trunc);
 //   if(os){
 //     os << _singleAmpFractions << endl;
 //     os << _interferenceFractions << endl;
 //     os.close();
 //   }
   if(!os)return true;
   os << "\\begin{tabular}{l r}" << "\n";
   os << "\\hline" << "\n";
   os << "\\hline" << "\n";
   os << "Decay channel & Fraction [$\\%$] \\\\" << "\n";
   os << "\\hline" << "\n";
 
   for(unsigned int i=0; i < this->size(); i++){    
     double frac = this->at(i).integral()/norm;
     TString name;
     if(this->at(i).isSingleAmp()) name = this->at(i).fitAmp1().name();
     else continue;
     FitFraction f((string)name, frac);
     double frac_err = fcov->getFractionError(i);
     
     name.ReplaceAll("->"," \\to ");
     name.ReplaceAll("Bs0","B_s");
     name.ReplaceAll("Ds","D_s");
     name.ReplaceAll("pi","\\pi");
     name.ReplaceAll("rho","\\rho");
     name.ReplaceAll("sigma10","\\sigma");
     name.ReplaceAll("GS","");
     name.ReplaceAll("SBW","");
     name.ReplaceAll("RhoOmega","");
     name.ReplaceAll("GLass","");
     name.ReplaceAll("Lass","");
     name.ReplaceAll("Bugg","");
     name.ReplaceAll("Flatte","");
     name.ReplaceAll("+","^+");
     name.ReplaceAll("-","^-");
     name.ReplaceAll("*","^*");
     name.ReplaceAll(")0",")^0");
     name.ReplaceAll(","," \\, ");
     name.ReplaceAll("NonResS0( \\to D_s^- \\, \\pi^+)","( D_s^- \\, \\pi^+)_{S}");
     name.ReplaceAll("NonResV0( \\to D_s^- \\, \\pi^+)","( D_s^- \\, \\pi^+)_{P}");
     name.ReplaceAll("NonResS0( \\to D_s^- \\, K^+)","( D_s^- \\, K^+)_{S}");
     name.ReplaceAll("NonResV0( \\to D_s^- \\, K^+)","( D_s^- \\, K^+)_{P}");
 
     os << std::fixed << std::setprecision(2) << "$" << name << "$ & " << frac * 100. << " $\\pm$ "  << frac_err * 100. <<  " \\\\" << "\n";
   }
   os <<  " \\hline" << "\n";
   os << " Sum & " << std::fixed << std::setprecision(2) << _singleAmpFractions.sum().frac() * 100. << " $\\pm$ "  << fcov->getFractionSumError() * 100. <<  " \\\\" << "\n";
   os << "\\hline" << "\n";
   os << "\\hline" << "\n";
   os << "\\end{tabular}" << "\n";
   os.close();
 
   TFile* paraFile = new TFile((fnameROOT).c_str(), "RECREATE");
   TTree* tree = new TTree("fractions","fractions");
     
   for(unsigned int i=0; i < this->size(); i++){    
     double frac = this->at(i).integral()/norm;
     TString name;
     if(this->at(i).isSingleAmp())name = TString(this->at(i).fitAmp1().name());
     else name = TString(this->at(i).name());
     FitFraction f((string)name, frac);  
 
     name.ReplaceAll("A(","");
     name.ReplaceAll("fit","");
     name.ReplaceAll("[","_");
     name.ReplaceAll("]","_");
     name.ReplaceAll("(","_");
     name.ReplaceAll(")","_");
     name.ReplaceAll("->","_");
     name.ReplaceAll("+","p");
     name.ReplaceAll("-","m");
     name.ReplaceAll("*","s");
     name.ReplaceAll(",","");
     name.ReplaceAll("GS","");
     name.ReplaceAll("GLass","");
     name.ReplaceAll("Lass","");
     name.ReplaceAll("RhoOmega","");
     name.ReplaceAll("Bugg","");
     name.ReplaceAll("Flatte","");
 
     if(0 != tree){
         double val = frac;
         if(!this->at(i).isSingleAmp()) name = "int_" + name;
         TBranch* Bra_val = tree->Branch( ((string)name).c_str(), &val, ((string)(name+"/D")).c_str());
         if(0 != Bra_val)  Bra_val->Fill();
     }
   }
   double val = _singleAmpFractions.sum().frac();
   TBranch* Bra_val = tree->Branch( "Sum", &val, "Sum/D");
   if(0 != Bra_val) Bra_val->Fill();
 
   if(0 != tree) tree->Fill();
 
   tree->Write();  
   paraFile->Close();
   delete paraFile;
   
   return true;
 }

◆ makeDirectory()

bool FitAmpPairList::makeDirectory ( const std::string & asSubdirOf = "." ) const

private

Definition at line 1155 of file FitAmpPairList.cpp.

                                                                   {
   /*
     A mode is created from or'd permission bit masks defined
      in <sys/stat.h>:
            #define S_IRWXU 0000700     RWX mask for owner 
            #define S_IRUSR 0000400     R for owner 
            #define S_IWUSR 0000200     W for owner 
            #define S_IXUSR 0000100     X for owner 
 
            #define S_IRWXG 0000070     RWX mask for group 
            #define S_IRGRP 0000040     R for group 
            #define S_IWGRP 0000020     W for group 
            #define S_IXGRP 0000010     X for group 
 
            #define S_IRWXO 0000007     RWX mask for other 
            #define S_IROTH 0000004     R for other 
            #define S_IWOTH 0000002     W for other 
            #define S_IXOTH 0000001     X for other 
 
            #define S_ISUID 0004000     set user id on execution 
            #define S_ISGID 0002000     set group id on execution 
            #define S_ISVTX 0001000     save swapped text even after use
    */
 
   mode_t mode = S_IRWXU | S_IRWXG | S_IRWXO 
               | S_ISUID | S_ISGID; 
   // see above for meaning. I want everybody to be allowed to read/write/exec.
   // Not sure about the last two bits.
 
   int zeroForSuccess = 0;
   zeroForSuccess |= mkdir( (asSubdirOf ).c_str(), mode );
   zeroForSuccess |= mkdir( (asSubdirOf + "/" + dirName() ).c_str(), mode );
   return (0 == zeroForSuccess);
 }

◆ needToReIntegrate()

bool FitAmpPairList::needToReIntegrate ( ) const

Definition at line 1197 of file FitAmpPairList.cpp.

                                             {
   for(unsigned int i=0; i < this->size(); i++){    
     if (this->at(i).needToReIntegrate() ) return true;
   }
   return false;
 }

◆ numberOfFitFractionsLargerThanThreshold()

int FitAmpPairList::numberOfFitFractionsLargerThanThreshold ( double threshold )

Definition at line 438 of file FitAmpPairList.cpp.

                                                                            {
     int n=0;
     double norm = integral();
     for(unsigned int i=0; i< this->size(); i++){
         if(this->at(i).isSingleAmp()){
             if(this->at(i).integral()/norm > threshold)n++ ;
         }
     }
     return n;
 }

◆ numEvents()

int FitAmpPairList::numEvents ( ) const

virtual

Implements IIntegrationCalculator.

Definition at line 230 of file FitAmpPairList.cpp.

                                    {
   return _Nevents;
 }

◆ oldVariance()

double FitAmpPairList::oldVariance ( ) const

privatevirtual

Definition at line 500 of file FitAmpPairList.cpp.

                                         {
   // this does not (yet) take into account
   // that the variance changes with the
   // fit paramters.
   bool dbThis=false;
   if(0 == _Nevents) return -9999.0;
 
   double dN = (double) _Nevents;
 
   double mean   = _sum/dN;
   double meansq = _sumsq/dN;
   double v = (meansq - mean*mean)/dN;
   double sf=1.0;
   if(0 != integral()){
     sf = integral()/mean;
     v        *= sf*sf;
     // this construct allows for a scaling factor (such as _weightSum)
     // in the integral.
   }
   if(dbThis){
     cout << "FitAmpPairList::oldVariance() "
      << " N = " << dN
      << " sum = " << _sum
      << " sumsq = " << _sumsq
      << " sf " << sf
      << " v " << v << endl;
   }
   return v;
 }

◆ operator+()

FitAmpPairList FitAmpPairList::operator+ ( const FitAmpPairList & other ) const

Definition at line 1241 of file FitAmpPairList.cpp.

                                                                          {
   FitAmpPairList returnVal(*this);
   returnVal.add(other);
   return returnVal;
 }

◆ operator+=()

FitAmpPairList & FitAmpPairList::operator+= ( const FitAmpPairList & other )

Definition at line 1237 of file FitAmpPairList.cpp.

                                                                      {
   this->add(other);
   return *this;
 }

◆ phaseSpaceIntegral()

double FitAmpPairList::phaseSpaceIntegral ( ) const

private

Definition at line 450 of file FitAmpPairList.cpp.

                                               {
   // only for debug
   if(_Nevents <= 0) return -9999;
   return _psSum/((double) _Nevents);
 }

◆ print()

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

virtual

Implements IIntegrationCalculator.

Definition at line 1190 of file FitAmpPairList.cpp.

                                               {
   os << "FitAmpPairList with " << this->size() << " FitAmpPairs:";
   for(unsigned int i=0; i < this->size(); i++){
     os << "\n\t" << i << ") " << this->at(i);
   }
 }

◆ reAddEvent()

void FitAmpPairList::reAddEvent	(	IDalitzEvent &	evt,
		double	weight = `1`
	)

virtual

Definition at line 141 of file FitAmpPairList.cpp.

                                                                {
   if(evt.phaseSpace() < 0) return;
 
   for(unsigned int i=0; i< this->size(); i++){
     if( this->at(i).acceptEvents()){
       this->at(i).reAdd(evt, weight, efficiency(&evt));
     }
   }
 
   //  _cov.addLastEventFromList(); (need to do something about it)
 
   return;
 }

◆ reset()

bool FitAmpPairList::reset ( )

private

Definition at line 218 of file FitAmpPairList.cpp.

                           {
   _Nevents=0;
   _sum=0;
   _sumsq=0;
   _psSum=0;
   _psSumSq=0;
   _cov.reset();
   _singleAmpFractions.clear();
   _interferenceFractions.clear();
 
   return true;
 }

◆ retrieve()

bool FitAmpPairList::retrieve ( const std::string & asSubdirOf = "." )

virtual

Implements IIntegrationCalculator.

Definition at line 1125 of file FitAmpPairList.cpp.

                                                         {
   bool sc=true;
   bool verbose=true;
   if(verbose){
     cout << "FitAmpPairList::retrieve: reading from "
      << asSubdirOf << endl; 
   }
 
   struct stat buf;
   if(stat(asSubdirOf.c_str(), &buf) != 0){
     if(verbose) cout << "FitAmpPairList::retrieve: " << asSubdirOf
              << " does not exist" << endl;
     return false;
   }
 
   ifstream is((asSubdirOf +"/"+ dirName() + "values.txt").c_str());
   std::string dummy;
   is >> dummy >> _Nevents;
   is >> dummy >> _sum;
   is >> dummy >> _sumsq;
   is >> dummy >> _psSum;
   is >> dummy >> _psSumSq;
   is.close();
   for(unsigned int i=0; i < this->size(); i++){
     sc &= this->at(i).retrieve(asSubdirOf +"/"+ dirName());
   }
   _cov.retrieve(asSubdirOf + "/" + dirName());
   return sc;
 }

◆ save()

bool FitAmpPairList::save ( const std::string & asSubdirOf = "." ) const

virtual

Implements IIntegrationCalculator.

Definition at line 1102 of file FitAmpPairList.cpp.

                                                           {
   bool sc=true;
   makeDirectory(asSubdirOf);
   ofstream os((asSubdirOf +"/"+ dirName() + "values.txt").c_str());
   os << "N " << _Nevents << endl;
   os << "sum " << setprecision(20) << _sum << endl;
   os << "sumsq " << setprecision(20) << _sumsq << endl;
   os << "psSum " << setprecision(20) << _psSum << endl;
   os << "psSumSq " << setprecision(20) << _psSumSq << endl;
   os.close();
 
   for(unsigned int i=0; i < this->size(); i++){
     sc &= this->at(i).save(asSubdirOf +"/"+ dirName());
   }
   if(_cov.isValid()){
     _cov.save(asSubdirOf + "/" + dirName());
   }else{
     cout << "WARNING IN FitAmpPairList::save(" << asSubdirOf
      << "): invalid _cov" << endl;
   }
 
   return sc;
 }

◆ saveEachAmpsHistograms()

void FitAmpPairList::saveEachAmpsHistograms ( const std::string & prefix ) const

virtual

Implements IIntegrationCalculator.

Definition at line 558 of file FitAmpPairList.cpp.

                                                                         {
   DalitzHistoSet sum;
   int counter=0;
   for(unsigned int i=0; i< this->size(); i++){
     if(this->at(i).isSingleAmp()){
       counter++;
       std::string name = prefix + "_" + anythingToString(counter) + ".root";
       DalitzHistoSet hs(this->at(i).histoSet());
       std::string title = this->at(i).fitAmp1().name();
       hs.setTitle(title);
       cout << "FitAmpPairList::saveEachAmpsHistograms: "
        << "saving " << title << " as " << name << endl;
       hs.save(name);
     }
   }
   return;
 }

◆ saveInterferenceHistograms()

void FitAmpPairList::saveInterferenceHistograms ( const std::string & prefix ) const

Definition at line 609 of file FitAmpPairList.cpp.

                                                                             {
     DalitzHistoSet sum;
     int counter=0;
     for(unsigned int i=0; i< this->size(); i++){
         if(!this->at(i).isSingleAmp()){
             counter++;
             std::string name = prefix + "_" + anythingToString(counter) + ".root";
             DalitzHistoSet hs(this->at(i).histoSet());
             std::string title = this->at(i).name();
             hs.setTitle(title);
             cout << "FitAmpPairList::saveEachAmpsHistograms: "
             << "saving " << title << " as " << name << endl;
             hs.save(name);
         }
     }
     return;
 }

◆ setEfficiency()

void FitAmpPairList::setEfficiency ( MINT::counted_ptr< MINT::IReturnRealForEvent< IDalitzEvent > > eff )

Definition at line 534 of file FitAmpPairList.cpp.

                                                                                      {
   _efficiency=eff;
 }

◆ setFast()

void FitAmpPairList::setFast ( )

Definition at line 1230 of file FitAmpPairList.cpp.

                             {
   _slow=false;
   for(unsigned int i=0; i < this->size(); i++){    
     this->at(i).setFast();
   }
 }

◆ setSlow()

void FitAmpPairList::setSlow ( )

Definition at line 1224 of file FitAmpPairList.cpp.

                             {
   _slow=true;
   for(unsigned int i=0; i < this->size(); i++){    
     this->at(i).setSlow();
   }
 }

◆ slow()

bool FitAmpPairList::slow ( ) const

inline

Definition at line 149 of file FitAmpPairList.h.

149 {return _slow;}

FitAmpPairList::_slow

bool _slow

Definition: FitAmpPairList.h:42

◆ startForcedReIntegration()

void FitAmpPairList::startForcedReIntegration ( )

Definition at line 1213 of file FitAmpPairList.cpp.

                                              {
   for(unsigned int i=0; i < this->size(); i++){    
     this->at(i).startReIntegration();
   }
 }

◆ startIntegration()

void FitAmpPairList::startIntegration ( )

Definition at line 1203 of file FitAmpPairList.cpp.

                                      {
   for(unsigned int i=0; i < this->size(); i++){    
     this->at(i).startIntegration();
   }
 }

◆ startReIntegration()

void FitAmpPairList::startReIntegration ( )

Definition at line 1208 of file FitAmpPairList.cpp.

                                        {
   for(unsigned int i=0; i < this->size(); i++){    
     if( this->at(i).needToReIntegrate() ) this->at(i).startReIntegration();
   }
 }

◆ sumOfFitFractions()

double FitAmpPairList::sumOfFitFractions ( )

Definition at line 414 of file FitAmpPairList.cpp.

                                          {
     double sum = 0;
     for(unsigned int i=0; i< this->size(); i++){
         if(this->at(i).isSingleAmp())sum += this->at(i).integral();
     }
     return sum/integral();
 }

◆ sumOfSqrtFitFractions()

double FitAmpPairList::sumOfSqrtFitFractions ( )

Definition at line 406 of file FitAmpPairList.cpp.

                                              {
     double sum = 0;
     for(unsigned int i=0; i< this->size(); i++){
         if(this->at(i).isSingleAmp())sum += sqrt(this->at(i).integral());
     }
     return sum;
 }

◆ sumOfVariances()

double FitAmpPairList::sumOfVariances ( ) const

virtual

Definition at line 493 of file FitAmpPairList.cpp.

                                            {
   double sum=0;
   for(unsigned int i=0; i < this->size(); i++){    
     sum += this->at(i).variance();
   }
   return sum;
 }

◆ un_normalised_histoSetIm()

DalitzHistoSet FitAmpPairList::un_normalised_histoSetIm ( ) const

virtual

Definition at line 550 of file FitAmpPairList.cpp.

                                                              {
   DalitzHistoSet sum;
   for(unsigned int i=0; i< this->size(); i++){
     sum += this->at(i).histoSetIm();
   }
   return sum;
 }

◆ un_normalised_histoSetRe()

DalitzHistoSet FitAmpPairList::un_normalised_histoSetRe ( ) const

virtual

Definition at line 543 of file FitAmpPairList.cpp.

                                                              {
   DalitzHistoSet sum;
   for(unsigned int i=0; i< this->size(); i++){
     sum += this->at(i).histoSetRe();
   }
   return sum;
 }

◆ unsetEfficiency()

void FitAmpPairList::unsetEfficiency ( )

Definition at line 537 of file FitAmpPairList.cpp.

                                     {
   _efficiency=counted_ptr<IReturnRealForEvent<IDalitzEvent> >(0);
 }

◆ variance()

double FitAmpPairList::variance ( ) const

virtual

Implements IIntegrationCalculator.

Definition at line 456 of file FitAmpPairList.cpp.

                                      {
   bool dbThis=false;
   static int printouts=0;
   static bool printedWarning=false;
   if(_Nevents <=0)return 0;
 
   if(_cov.isValid()){
     if(dbThis && printouts < 100){
       printouts++;
       cout << "FitAmpPairList::variance() with "
        << _Nevents << "(" << _cov.Nevents() << ") events. Compare: "
        << " old: " << oldVariance()
        << " sum: " << sumOfVariances()
        << ", new: " << _cov.getIntegralVariance()
        << endl;
     }
     if(_Nevents < 10000) return oldVariance();
     double var = _cov.getIntegralVariance();
     if(var < 1.e-6 * oldVariance()){
       if(dbThis) cout << "returning OLD variance " << var << endl;
       // (if var is tiny, it usually means something went wrong)
       return oldVariance(); 
     }
     if(dbThis) cout << "returning new variance " << var << endl;
     return var;
   }else{
     if(! printedWarning){
       cout << "WARNING in FitAmpPairList::variance()"
        << " ignoring correlations because _cov is invalid."
        << " (I'll print this only once)"
        << endl;
       printedWarning=true;
     }
     return sumOfVariances();
   }
 }

Member Data Documentation

◆ _cov

FitAmpPairCovariance FitAmpPairList::_cov

mutableprivate

Definition at line 44 of file FitAmpPairList.h.

◆ _efficiency

MINT::counted_ptr<MINT::IReturnRealForEvent<IDalitzEvent> > FitAmpPairList::_efficiency

private

Definition at line 46 of file FitAmpPairList.h.

◆ _interferenceFractions

FitFractionList FitAmpPairList::_interferenceFractions

private

Definition at line 48 of file FitAmpPairList.h.

◆ _Nevents

int FitAmpPairList::_Nevents

private

Definition at line 35 of file FitAmpPairList.h.

◆ _psSum

double FitAmpPairList::_psSum

private

Definition at line 39 of file FitAmpPairList.h.

◆ _psSumSq

double FitAmpPairList::_psSumSq

private

Definition at line 40 of file FitAmpPairList.h.

◆ _singleAmpFractions

FitFractionList FitAmpPairList::_singleAmpFractions

private

Definition at line 48 of file FitAmpPairList.h.

◆ _slow

bool FitAmpPairList::_slow

private

Definition at line 42 of file FitAmpPairList.h.

◆ _sum

double FitAmpPairList::_sum

private

Definition at line 36 of file FitAmpPairList.h.

◆ _sumsq

double FitAmpPairList::_sumsq

private

Definition at line 37 of file FitAmpPairList.h.

◆ HistoOption

MINT::NamedParameter<std::string> FitAmpPairList::HistoOption

private

Definition at line 32 of file FitAmpPairList.h.

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

Mint/FitAmpPairList.h
src/Mojito/FitAmplitude/FitAmpPairList.cpp

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FitAmpPairList() [1/2]

◆ FitAmpPairList() [2/2]

◆ ~FitAmpPairList()

Member Function Documentation

◆ absSumOfInterferenceFractions()

◆ absSumOfSqrtInterferenceFractions()

◆ add() [1/3]

◆ add() [2/3]

◆ add() [3/3]

◆ addAmps() [1/2]

◆ addAmps() [2/2]

◆ addEvent() [1/3]

◆ addEvent() [2/3]

◆ addEvent() [3/3]

◆ append() [1/3]

◆ append() [2/3]

◆ append() [3/3]

◆ applyHistoOption()

◆ clone_IIntegrationCalculator()

◆ ComplexIntegralForTags()

◆ ComplexSum()

◆ ComplexSumForMatchingPatterns()

◆ dirName()

◆ doFinalStats()

◆ doFinalStatsAndSave()

◆ doFractions()

◆ efficiency()

◆ endIntegration()

◆ fast()

◆ GetEachAmpsHistograms()

◆ getFractionChi2()

◆ getFractions()

◆ GetInterferenceHistograms()

◆ getInterferenceTerms()

◆ Gradient()

◆ GradientForLasso()

◆ haveEfficiency()

◆ histoSet()

◆ integral()

◆ integralForMatchingPatterns()

◆ interferenceHistoSet()

◆ isCompatibleWith()

◆ makeAndStoreFractions() [1/2]

◆ makeAndStoreFractions() [2/2]

◆ makeDirectory()

◆ needToReIntegrate()

◆ numberOfFitFractionsLargerThanThreshold()

◆ numEvents()

◆ oldVariance()

◆ operator+()

◆ operator+=()

◆ phaseSpaceIntegral()

◆ print()

◆ reAddEvent()

◆ reset()

◆ retrieve()

◆ save()

◆ saveEachAmpsHistograms()

◆ saveInterferenceHistograms()

◆ setEfficiency()

◆ setFast()

◆ setSlow()

◆ slow()

◆ startForcedReIntegration()

◆ startIntegration()

◆ startReIntegration()

◆ sumOfFitFractions()

◆ sumOfSqrtFitFractions()

◆ sumOfVariances()

◆ un_normalised_histoSetIm()

◆ un_normalised_histoSetRe()

◆ unsetEfficiency()

◆ variance()

Member Data Documentation