#include <FitAmpSum.h>

Inheritance diagram for FitAmpSum:

Public Member Functions
	FitAmpSum (const DalitzEventPattern &pat, const char fname=0, MINT::MinuitParameterSet pset=0, const std::string &prefix="", const std::string &lineshapePrefix="", const std::string &opt="")

	FitAmpSum (const DalitzEventPattern &pat, MINT::MinuitParameterSet *pset, const std::string &prefix="", const std::string &lineshapePrefix="", const std::string &opt="")

	FitAmpSum (const DalitzEventPattern &pat, const std::string &prefix, const std::string &lineshapePrefix="", const std::string &opt="")

	FitAmpSum (const FitAmpSum &other)

	FitAmpSum (const FitAmpList &other)

virtual MINT::counted_ptr< FitAmpListBase >	GetCloneSameFitParameters () const

virtual MINT::counted_ptr< FitAmpListBase >	GetCloneOfSubsetSameFitParameters (std::string name) const

virtual DalitzBoxSet	makeBoxes (const DalitzEventPattern &pat, double nSigma=2)

virtual DalitzBWBoxSet	makeBWBoxes (const DalitzEventPattern &pat, TRandom *rnd=gRandom)

virtual std::complex< double >	getVal (IDalitzEvent &evt)

virtual std::complex< double >	getVal (IDalitzEvent *evt)

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

virtual bool	useAnalyticGradient ()

virtual MINT::counted_ptr< IIntegrationCalculator >	makeIntegrationCalculator ()

virtual MINT::counted_ptr< IntegCalculator >	makeIntegCalculator ()

virtual MINT::counted_ptr< FitAmpPairList >	makeFitAmpPairList ()

virtual double	Prob (IDalitzEvent &evt)

virtual double	RealVal (IDalitzEvent &evt)

virtual std::complex< double >	ComplexVal (IDalitzEvent &evt)

double	getAmpSqr (IDalitzEvent &evt, std::vector< std::string > ampNames, bool CC=false)

virtual MINT::counted_ptr< MINT::IUnweightedEventGenerator< IDalitzEvent > >	makeEventGenerator (const DalitzEventPattern &pat, TRandom *rnd=gRandom)

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

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

virtual void	printValues (IDalitzEvent &evt, std::ostream &os=std::cout)

virtual	~FitAmpSum ()

FitAmpSum &	operator *= (double r)

FitAmpSum &	operator *= (const std::complex< double > &z)

FitAmpSum &	operator *= (const MINT::counted_ptr< MINT::IReturnComplex > &irc)

FitAmpSum &	operator *= (const MINT::counted_ptr< MINT::IReturnComplexForEvent< IDalitzEvent > > &ircfe)

FitAmpSum	operator * (double r) const

FitAmpSum	operator * (const std::complex< double > &z) const

FitAmpSum	operator * (const MINT::counted_ptr< MINT::IReturnComplex > &irc) const

FitAmpSum	operator * (const MINT::counted_ptr< MINT::IReturnComplexForEvent< IDalitzEvent > > &ircfe) const

FitAmpSum &	operator= (const FitAmpSum &other)

FitAmpSum &	operator= (const FitAmpList &other)

FitAmpSum &	operator+= (const FitAmpSum &other)

FitAmpSum	operator+ (const FitAmpSum &other) const

Public Member Functions inherited from FitAmpList
	FitAmpList (const DalitzEventPattern &pat, const char fname=0, MINT::MinuitParameterSet pset=0, const std::string &prefix="", const std::string &lineshapePrefix="", const std::string &opt="")

	FitAmpList (const DalitzEventPattern &pat, MINT::MinuitParameterSet *pset, const std::string &prefix="", const std::string &lineshapePrefix="", const std::string &opt="")

	FitAmpList (const DalitzEventPattern &pat, const std::string &prefix, const std::string &lineshapePrefix="", const std::string &opt="")

	FitAmpList (const FitAmpList &other)

virtual MINT::MinuitParameterSet *	getMPS ()

virtual MINT::counted_ptr< FitAmpListBase >	GetCPConjugateSameFitParameters () const

virtual bool	createAllAmps (const DalitzEventPattern &thePattern, const std::string &prefix="", const std::string &lineshapePrefix="")

virtual bool	createAmpsFromFile ()

virtual	~FitAmpList ()

virtual void	printLargestAmp (IDalitzEvent &evt, std::ostream &os=std::cout)

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

virtual void	printAllAmps (IDalitzEvent &evt, std::ostream &os=std::cout)

virtual void	printNonZeroWithValue (IDalitzEvent &evt, std::ostream &os=std::cout)

FitAmpList &	operator *= (double r)

FitAmpList &	operator *= (const std::complex< double > &z)

FitAmpList &	operator *= (const MINT::counted_ptr< MINT::IReturnComplex > &irc)

FitAmpList	operator * (double r) const

FitAmpList	operator * (const std::complex< double > &z) const

FitAmpList	operator * (const MINT::counted_ptr< MINT::IReturnComplex > &irc) const

FitAmpList &	operator= (const FitAmpList &other)

FitAmpList &	operator+= (const FitAmpList &other)

FitAmpList	operator+ (const FitAmpList &other) const

Public Member Functions inherited from FitAmpListBase
	FitAmpListBase ()

	FitAmpListBase (const FitAmpListBase &other)

	FitAmpListBase (const FitAmpListBase &other, std::string name)

virtual int	addCopyWithSameFitParameters (const FitAmpListBase &other, double factor=1)

virtual int	addCopyOfSubsetWithSameFitParameters (const FitAmpListBase &other, std::string name, double factor=1)

virtual bool	addAmplitude (FitAmplitude *fa)

virtual int	add (const FitAmpListBase &other, double factor=1)

virtual int	addAsList (const FitAmpListBase &other, double factor=1)

virtual unsigned int	size () const

virtual FitAmplitude *	getAmpPtr (unsigned int i)

const FitAmplitude *	getAmpPtr (unsigned int i) const

virtual bool	CPConjugateSameFitParameters ()

virtual bool	CConjugateFinalStateSameFitParameters ()

virtual MINT::counted_ptr< FitAmpListBase >	GetCConjugateFinalStateSameFitParameters () const

virtual bool	CConjugateInitialStateSameFitParameters ()

virtual MINT::counted_ptr< FitAmpListBase >	GetCConjugateInitialStateSameFitParameters () const

virtual bool	setLSameFitParameters (int L)

virtual MINT::counted_ptr< FitAmpListBase >	GetDifferentLSameFitParameters (int L) const

DalitzBoxSet	makeBoxes (const DalitzEventPattern &pat, MINT::IReturnRealForEvent< IDalitzEvent > *pdf, double nSigma=2)

DalitzBWBoxSet	makeBWBoxes (const DalitzEventPattern &pat, MINT::IReturnRealForEvent< IDalitzEvent > pdf, TRandom rnd=gRandom)

void	setAllAmpsTo (const std::complex< double > z)

virtual	~FitAmpListBase ()

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

double	efficiency (IDalitzEvent &evt)

virtual void	multiply (double r)

virtual void	multiply (const std::complex< double > &z)

virtual void	multiply (const MINT::counted_ptr< MINT::IReturnComplex > irc)

virtual void	multiply (MINT::counted_ptr< MINT::IReturnComplexForEvent< IDalitzEvent > > irc)

FitAmpListBase &	operator *= (double r)

FitAmpListBase &	operator *= (const std::complex< double > &z)

FitAmpListBase &	operator *= (const MINT::counted_ptr< MINT::IReturnComplex > &irc)

FitAmpListBase	operator * (double r) const

FitAmpListBase	operator * (const std::complex< double > &z) const

FitAmpListBase	operator * (const MINT::counted_ptr< MINT::IReturnComplex > &irc) const

void	normalizeAmps (DalitzEventList &evtList)

std::vector< double >	normFactors (DalitzEventList &evtList)

void	randomizeStartVals (int seed=0)

void	randomizePhaseStartVals (int seed=0)

void	setTag (int tag)

FitAmpListBase &	operator= (const FitAmpListBase &other)

FitAmpListBase &	operator+= (const FitAmpListBase &other)

FitAmpListBase	operator+ (const FitAmpListBase &other) const

Public Member Functions inherited from ILookLikeFitAmpSum
virtual	~ILookLikeFitAmpSum ()

Public Member Functions inherited from IFastAmplitudeIntegrable
virtual	~IFastAmplitudeIntegrable ()

Public Member Functions inherited from MINT::IReturnRealForEvent< IDalitzEvent >
virtual	~IReturnRealForEvent ()

Public Member Functions inherited from MINT::IReturnComplexForEvent< IDalitzEvent >
virtual	~IReturnComplexForEvent ()

Protected Attributes
MINT::NamedParameter< int >	_useAnalyticGradient

Protected Attributes inherited from FitAmpList
DalitzEventPattern	_pat

std::string	_paraFName

MINT::MinuitParameterSet *	_minuitParaSet

std::string	_opt

Protected Attributes inherited from FitAmpListBase
std::vector< FitAmplitude * >	_fitAmps

std::vector< MINT::counted_ptr< FitAmpListBase > >	_fitAmpLists

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

Friends
class	FitAmplitude

Additional Inherited Members
Protected Member Functions inherited from FitAmpListBase
virtual void	deleteAll ()

Protected Member Functions inherited from MINT::IReturnRealForEvent< IDalitzEvent >
	IReturnRealForEvent ()

Protected Member Functions inherited from MINT::IReturnComplexForEvent< IDalitzEvent >
	IReturnComplexForEvent ()

Detailed Description

Definition at line 26 of file FitAmpSum.h.

Constructor & Destructor Documentation

◆ FitAmpSum() [1/5]

FitAmpSum::FitAmpSum	(	const DalitzEventPattern &	pat,
		const char *	fname = `0`,
		MINT::MinuitParameterSet *	pset = `0`,
		const std::string &	prefix = `""`,
		const std::string &	lineshapePrefix = `""`,
		const std::string &	opt = `""`
	)

Definition at line 20 of file FitAmpSum.cpp.

   : FitAmpList(pat, fname, pset, prefix, lineshapePrefix, opt)
   , _useAnalyticGradient("useAnalyticGradient",0)  
 {
   /*
     //Important! Ensures everything is initialised
     DalitzEventList eventTest;
     eventTest.generatePhaseSpaceEvents(1,pat);
     this->getVal(eventTest[0]);  
   */
 }

◆ FitAmpSum() [2/5]

FitAmpSum::FitAmpSum	(	const DalitzEventPattern &	pat,
		MINT::MinuitParameterSet *	pset,
		const std::string &	prefix = `""`,
		const std::string &	lineshapePrefix = `""`,
		const std::string &	opt = `""`
	)

Definition at line 38 of file FitAmpSum.cpp.

   : FitAmpList(pat, pset, prefix, lineshapePrefix, opt)
   , _useAnalyticGradient("useAnalyticGradient",0)  
 {
 
   // cout << "pset pointer in FitAmpSum::FitAmpSum " << pset << " = " << getMPS() << endl;
 
   /*
     //Important! Ensures everything is initialised
     DalitzEventList eventTest;
     eventTest.generatePhaseSpaceEvents(1,pat);
     this->getVal(eventTest[0]);  
   */
 }

◆ FitAmpSum() [3/5]

FitAmpSum::FitAmpSum	(	const DalitzEventPattern &	pat,
		const std::string &	prefix,
		const std::string &	lineshapePrefix = `""`,
		const std::string &	opt = `""`
	)

Definition at line 57 of file FitAmpSum.cpp.

   : FitAmpList(pat, prefix, lineshapePrefix, opt)
   , _useAnalyticGradient("useAnalyticGradient",0)
 {
   /*
     //Important! Ensures everything is initialised
     DalitzEventList eventTest;
     eventTest.generatePhaseSpaceEvents(1,pat);
     this->getVal(eventTest[0]); 
   */
 }

◆ FitAmpSum() [4/5]

FitAmpSum::FitAmpSum ( const FitAmpSum & other )

Definition at line 73 of file FitAmpSum.cpp.

   : IReturnRealForEvent<IDalitzEvent>()
   , IReturnComplexForEvent<IDalitzEvent>()
   , IFastAmplitudeIntegrable()
   , ILookLikeFitAmpSum()
   , FitAmpList(other)
   , _useAnalyticGradient("useAnalyticGradient",0)
 {
   /*
     //Important! Ensures everything is initialised
     DalitzEventList eventTest;
     eventTest.generatePhaseSpaceEvents(1,_pat);
     this->getVal(eventTest[0]); 
   */
 }

◆ FitAmpSum() [5/5]

FitAmpSum::FitAmpSum ( const FitAmpList & other )

Definition at line 89 of file FitAmpSum.cpp.

   : IReturnRealForEvent<IDalitzEvent>()
   , IReturnComplexForEvent<IDalitzEvent>()
   , IFastAmplitudeIntegrable()
   , ILookLikeFitAmpSum()
   , FitAmpList(other)
   , _useAnalyticGradient("useAnalyticGradient",0)
 {    
   /*
   //Important! Ensures everything is initialised
     DalitzEventList eventTest;
     eventTest.generatePhaseSpaceEvents(1,_pat);
     this->getVal(eventTest[0]); 
   */
 }

◆ ~FitAmpSum()

FitAmpSum::~FitAmpSum ( )

virtual

Definition at line 476 of file FitAmpSum.cpp.

                      {
   bool dbThis=true;
   if(dbThis) cout << "FitAmpSum::~FitAmpSum() " << endl;
   deleteAll();
   if(dbThis) cout << "FitAmpSum::~FitAmpSum() done" << endl;
 
 }

Member Function Documentation

◆ ComplexVal()

virtual std::complex<double> FitAmpSum::ComplexVal ( IDalitzEvent & evt )

inlinevirtual

Reimplemented from IFastAmplitudeIntegrable.

Definition at line 122 of file FitAmpSum.h.

122 {return getVal(evt);}

FitAmpSum::getVal

virtual std::complex< double > getVal(IDalitzEvent &evt)

Definition: FitAmpSum.cpp:182

◆ getAmpSqr()

double FitAmpSum::getAmpSqr	(	IDalitzEvent &	evt,
		std::vector< std::string >	ampNames,
		bool	CC = `false`
	)

Definition at line 214 of file FitAmpSum.cpp.

                                                                                  {
     double dbThis=false;
     
     std::complex<double> sum(0);
     if(0 == size()) createAllAmps(evt.eventPattern());
     for(unsigned int i=0; i< this->size(); i++){
         for (unsigned int n=0; n <ampNames.size(); n++) {
             if(A_is_in_B(ampNames[n], this->getAmpPtr(i)->name())){
                 if(A_is_in_B("Cconj_FS", this->getAmpPtr(i)->name()) && !CC) continue;
                 if(dbThis) cout << "found amp " << this->getAmpPtr(i)->name() << endl;
                 sum += this->getAmpPtr(i)->getVal(evt);
             }
         }
     }
     std::complex<double> result(sqrt(fabs(efficiency(evt)))*sum);
 
     return norm(result);
 }

◆ GetCloneOfSubsetSameFitParameters()

counted_ptr< FitAmpListBase > FitAmpSum::GetCloneOfSubsetSameFitParameters ( std::string name ) const

virtual

Reimplemented from FitAmpListBase.

Definition at line 147 of file FitAmpSum.cpp.

                                                                                          { 
     // need to reform these one day...
     // ... it all relies on the copy-constructur/=operator in FitAmpitude
     // not copying the fit parameters, but just their pointers
     // which will need to be reviewed.
     //
     bool dbThis=true;
     if(dbThis) cout << "FitAmpSum::GetCloneOfSubsetSameFitParameters()" << endl;
     /* 
      There'll be 'physical' copies of all Amplitudes, but the
      FitParameters remain the same (pointers to the same
      FitParameter Object).  This is useful for the CP-con coding
      as it is now, but perhaps a bit counter-intuitive.  needs to
      be reviewed at some point. This behaviour is defined in the
      copy constructor of the FitAmplitude class.
      */
     
     /*
      counted_ptr<FitAmpSum> 
      newList(new FitAmpSum((IDalitzEventList*) this->getEventRecord()
      , _paraFName.c_str(), _minuitParaSet));
      */
     counted_ptr<FitAmpSum> newList(new FitAmpSum(*this));
     newList->deleteAll();
     
     if(dbThis) cout << "FitAmpSum::GetCloneOfSubsetSameFitParameters() made ptr" << endl;
     newList->addCopyOfSubsetWithSameFitParameters(*this,(string) name);
     if(dbThis) cout << "cloning subset of FitAmpSum " << newList->size() << endl;
     return newList;
 }

◆ GetCloneSameFitParameters()

counted_ptr< FitAmpListBase > FitAmpSum::GetCloneSameFitParameters ( ) const

virtual

Reimplemented from FitAmpList.

Definition at line 116 of file FitAmpSum.cpp.

                                                                       { 
 // need to reform these one day...
 // ... it all relies on the copy-constructur/=operator in FitAmpitude
 // not copying the fit parameters, but just their pointers
 // which will need to be reviewed.
 //
   bool dbThis=false;
   if(dbThis) cout << "FitAmpSum::GetCloneSameFitParameters()" << endl;
   /* 
      There'll be 'physical' copies of all Amplitudes, but the
      FitParameters remain the same (pointers to the same
      FitParameter Object).  This is useful for the CP-con coding
      as it is now, but perhaps a bit counter-intuitive.  needs to
      be reviewed at some point. This behaviour is defined in the
      copy constructor of the FitAmplitude class.
   */
 
   /*
   counted_ptr<FitAmpSum> 
     newList(new FitAmpSum((IDalitzEventList*) this->getEventRecord()
                , _paraFName.c_str(), _minuitParaSet));
   */
   counted_ptr<FitAmpSum> newList(new FitAmpSum(*this));
   newList->deleteAll();
 
   if(dbThis) cout << "FitAmpSum::GetCloneSameFitParameters() made ptr" << endl;
   newList->add(*this);
   if(dbThis) cout << "cloning FitAmpSum " << newList->size() << endl;
   return newList;
 }

◆ getVal() [1/2]

std::complex< double > FitAmpSum::getVal ( IDalitzEvent & evt )

virtual

Definition at line 182 of file FitAmpSum.cpp.

                                                      {
   //double dbThis=false;
 
   std::complex<double> sum(0);
   //if(0 == size()) createAllAmps(evt.eventPattern());
   for(unsigned int i=0; i< this->size(); i++){
     //if(dbThis) cout << "FitAmpSum::getVal: " << i << endl;
     sum = (! this->getAmpPtr(i)->isZero()) ? sum + this->getAmpPtr(i)->getVal(evt) : sum;
       //cout << sum << "," << endl;
   }
 //   if(dbThis){
 //     cout << "sum= " << sum << endl;
 //     cout << "efficiency(evt) " << efficiency(evt) << endl;
 //   }
   std::complex<double> result(sqrt(fabs(efficiency(evt)))*sum);
 //   if(dbThis) cout << "result " << result << endl;
   /* 
   double resultSq = norm(result);
   bool invalid = (! isfinite(resultSq)) || std::isnan(resultSq);
   if(dbThis || invalid){
     cout << "\n---------------------------" << endl;
     cout << "sqrt(efficiency())*sum = "
      << "sqrt(" << efficiency(evt) << ")*" << sum
      << endl;
     printValues(evt);
     cout << "\n---------------------------" << endl;
   }
    */
   return result;
 }

◆ getVal() [2/2]

std::complex< double > FitAmpSum::getVal ( IDalitzEvent * evt )

virtual

Definition at line 178 of file FitAmpSum.cpp.

                                                      {
   return getVal(*evt);
 }

◆ Gradient()

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

virtual

Reimplemented from IFastAmplitudeIntegrable.

Definition at line 233 of file FitAmpSum.cpp.

                                                                                       {
     
     std::complex<double> val = getVal(evt);
     std::complex<double> valConj= conj(val);
         
     for (unsigned int i=0; i<mps->size(); i++) {
         if(mps->getParPtr(i)->hidden())continue;
 
         string name_i= mps->getParPtr(i)->name();
         if(name_i.find("Inco")!=std::string::npos)continue;
 
         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,"");
             for(unsigned int j=0; j< this->size(); j++){
           if(A_is_in_B(name_i, this->getAmpPtr(j)->name())){
         if(i+1 >= grad.size()){
           cout << "WARNING in FitAmpSum::Gradient"
                << " have to increase size of grad to avoid memory issues" << endl;
           grad.resize(i+2);
         }
         
         std::complex<double> tmp = 2.*valConj* getAmpPtr(j)->getValWithoutFitParameters(evt);
                     grad[i]= tmp.real();
                     grad[i+1]= -tmp.imag();
                     i++;
                     break;
                 }
             }
         }
         // Doesn't work. Don't use! 
         /* 
         else if(name_i.find("_Amp")!=std::string::npos){
             name_i.replace(name_i.find("_Amp"),4,"");
             
             for(unsigned int j=0; j<_fitAmps.size(); j++){
                 if(A_is_in_B(name_i, _fitAmps[j]->name())){
                     std::complex<double> tmp = (_fitAmps[j])->getValWithoutFitParameters(evt);
                     grad[i]= 2.*(tmp*valConj).real()/std::abs((_fitAmps[j])->AmpPhase());
                     grad[i+1]= -2.*std::arg((_fitAmps[j])->AmpPhase())*(tmp*(valConj-conj((_fitAmps[j])->getVal(evt)))).imag();
                     i++;
                     break;
                 }
             }
         }
         */
         else if(mps->getParPtr(i)->iFixInit())continue;
         else {
             std::cout << "FitAmpSum::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";
         }
         
     }
     
 } 

◆ makeBoxes()

virtual DalitzBoxSet FitAmpSum::makeBoxes	(	const DalitzEventPattern &	pat,
		double	nSigma = `2`
	)

inlinevirtual

Definition at line 77 of file FitAmpSum.h.

78 {

79 return FitAmpList::makeBoxes(pat, this, nSigma);}

FitAmpListBase::makeBoxes

DalitzBoxSet makeBoxes(const DalitzEventPattern &pat, MINT::IReturnRealForEvent< IDalitzEvent > *pdf, double nSigma=2)

Definition: FitAmpListBase.cpp:374

◆ makeBWBoxes()

virtual DalitzBWBoxSet FitAmpSum::makeBWBoxes	(	const DalitzEventPattern &	pat,
		TRandom *	rnd = `gRandom`
	)

inlinevirtual

Implements ILookLikeFitAmpSum.

Definition at line 81 of file FitAmpSum.h.

82 {

83 return FitAmpList::makeBWBoxes(pat, this, rnd);}

FitAmpListBase::makeBWBoxes

DalitzBWBoxSet makeBWBoxes(const DalitzEventPattern &pat, MINT::IReturnRealForEvent< IDalitzEvent > *pdf, TRandom *rnd=gRandom)

Definition: FitAmpListBase.cpp:388

◆ makeEventGenerator()

virtual MINT::counted_ptr<MINT::IUnweightedEventGenerator<IDalitzEvent> > FitAmpSum::makeEventGenerator	(	const DalitzEventPattern &	pat,
		TRandom *	rnd = `gRandom`
	)

inlinevirtual

Implements IFastAmplitudeIntegrable.

Definition at line 127 of file FitAmpSum.h.

                                                                            {
     MINT::counted_ptr<MINT::IUnweightedEventGenerator<IDalitzEvent> > 
       ptr(new DalitzBWBoxSet(makeBWBoxes(pat, rnd)));
     return ptr;
   }

◆ makeFitAmpPairList()

counted_ptr< FitAmpPairList > FitAmpSum::makeFitAmpPairList ( )

virtual

Implements IFastAmplitudeIntegrable.

Definition at line 404 of file FitAmpSum.cpp.

                                                          {
   counted_ptr<FitAmpPairList> list(new FitAmpPairList);
   for(unsigned int i=0; i < _fitAmps.size(); i++){
     if(_fitAmps[i]->canBeIgnored()) continue;
     for(unsigned int j=i; j < _fitAmps.size(); j++){
       if(_fitAmps[j]->canBeIgnored()) continue;
       // all terms within the list
       list->addAmps( (_fitAmps[i]),  (_fitAmps[j]));
     }
     // cross terms between list, and list-of-list
     for(unsigned int k=0; k < _fitAmpLists.size(); k++){
       for(unsigned int l=0; l < _fitAmpLists[k]->size(); l++){
     if(_fitAmpLists[k]->getAmpPtr(l)->canBeIgnored()) continue;
     list->addAmps(_fitAmps[i], _fitAmpLists[k]->getAmpPtr(l));
       }
     }
   }
 
 
   // Now it gets complicated. I want to leave all the 
   // terms within a list to whatever that list
   // wants to do. But we don't want to forget cross terms 
   // between lists.
 
   // cross terms between list, and list-of-list
   for(unsigned int i=1; i < _fitAmpLists.size(); i++){
     for(unsigned int j=0; j < i; j++){
       for(unsigned int ki = 0; ki < _fitAmpLists[i]->size(); ki++){
     if(_fitAmpLists[i]->getAmpPtr(ki)->canBeIgnored()) continue;
     for(unsigned int kj = 0; kj < _fitAmpLists[j]->size(); kj++){
       if(_fitAmpLists[j]->getAmpPtr(kj)->canBeIgnored()) continue;
       list->addAmps(_fitAmpLists[i]->getAmpPtr(ki), _fitAmpLists[j]->getAmpPtr(kj));
     }
       }
     }
   }
 
 
   // finally, the terms within each list
   // which is up to the list itself.
   // (this is where we save a lot of space and time
   // when it comes to the model-independent stuff
   // where there are no internal cross terms)
   for(unsigned int i=0; i < _fitAmpLists.size(); i++){
     list->append(_fitAmpLists[i]->makeFitAmpPairList());
   }
 
   cout << "FitAmpSum: setting efficiency POINTER "
        << " in integCalculator to " 
        << _efficiency.get();
   if(0 == _efficiency.get()){
     cout << " (0 means no pointer, 100% efficiency).";
   }
   cout << endl;
 
   list->setEfficiency(_efficiency);
 
   /*
   cout << "this->size()" << this->size() << endl;
   cout << "_fitAmpLists.size()" << _fitAmpLists.size() << endl;
   for(unsigned int i=0; i < _fitAmpLists.size(); i++){
     cout << "printing list number " << i << endl;
     _fitAmpLists[i]->print();
     cout << endl;
   }
   cout << " the pair list: " << endl;
   list->print();
   */
 
   return list;
 }

◆ makeIntegCalculator()

counted_ptr< IntegCalculator > FitAmpSum::makeIntegCalculator ( )

virtual

Implements IFastAmplitudeIntegrable.

Definition at line 333 of file FitAmpSum.cpp.

                                                            {
   counted_ptr<IntegCalculator> list(new IntegCalculator);
   for(unsigned int i=0; i < _fitAmps.size(); i++){
     if(_fitAmps[i]->canBeIgnored()) continue;
     for(unsigned int j=i; j < _fitAmps.size(); j++){
       if(_fitAmps[j]->canBeIgnored()) continue;
       // all terms within the list
       list->addAmps( (_fitAmps[i]),  (_fitAmps[j]));
     }
     // cross terms between list, and list-of-list
     for(unsigned int k=0; k < _fitAmpLists.size(); k++){
       for(unsigned int l=0; l < _fitAmpLists[k]->size(); l++){
     if(_fitAmpLists[k]->getAmpPtr(l)->canBeIgnored()) continue;
     list->addAmps(_fitAmps[i], _fitAmpLists[k]->getAmpPtr(l));
       }
     }
   }
 
 
   // Now it gets complicated. I want to leave all the 
   // terms within a list to whatever that list
   // wants to do. But we don't want to forget cross terms 
   // between lists.
 
   // cross terms between list, and list-of-list
   for(unsigned int i=1; i < _fitAmpLists.size(); i++){
     for(unsigned int j=0; j < i; j++){
       for(unsigned int ki = 0; ki < _fitAmpLists[i]->size(); ki++){
     if(_fitAmpLists[i]->getAmpPtr(ki)->canBeIgnored()) continue;
     for(unsigned int kj = 0; kj < _fitAmpLists[j]->size(); kj++){
       if(_fitAmpLists[j]->getAmpPtr(kj)->canBeIgnored()) continue;
       list->addAmps(_fitAmpLists[i]->getAmpPtr(ki), _fitAmpLists[j]->getAmpPtr(kj));
     }
       }
     }
   }
 
 
   // finally, the terms within each list
   // which is up to the list itself.
   // (this is where we save a lot of space and time
   // when it comes to the model-independent stuff
   // where there are no internal cross terms)
   for(unsigned int i=0; i < _fitAmpLists.size(); i++){
     list->append(_fitAmpLists[i]->makeIntegCalculator());
   }
 
   cout << "FitAmpSum: setting efficiency POINTER "
        << " in integCalculator to " 
        << _efficiency.get();
   if(0 == _efficiency.get()){
     cout << " (0 means no pointer, 100% efficiency).";
   }
   cout << endl;
 
   list->setEfficiency(_efficiency);
 
   
   cout << "this->size()" << this->size() << endl;
   cout << "_fitAmpLists.size()" << _fitAmpLists.size() << endl;
   for(unsigned int i=0; i < _fitAmpLists.size(); i++){
     cout << "printing list number " << i << endl;
     _fitAmpLists[i]->print();
     cout << endl;
   }
   cout << " the pair list: " << endl;
   list->print();
   
   return list;
 }

◆ makeIntegrationCalculator()

counted_ptr< IIntegrationCalculator > FitAmpSum::makeIntegrationCalculator ( )

virtual

Implements IFastAmplitudeIntegrable.

Definition at line 328 of file FitAmpSum.cpp.

                                                                         {
   return (counted_ptr<IIntegrationCalculator>) makeIntegCalculator();
 }

◆ operator *() [1/4]

FitAmpSum FitAmpSum::operator * ( double r ) const

Definition at line 511 of file FitAmpSum.cpp.

                                             {
   FitAmpSum fas(*this);
   fas.multiply(r);
   return fas;
 }

◆ operator *() [2/4]

FitAmpSum FitAmpSum::operator * ( const std::complex< double > & z ) const

◆ operator *() [3/4]

FitAmpSum FitAmpSum::operator * ( const MINT::counted_ptr< MINT::IReturnComplex > & irc ) const

Definition at line 521 of file FitAmpSum.cpp.

                                                                           {
   FitAmpSum fas(*this);
   fas.multiply(irc);
   return fas;
 }

◆ operator *() [4/4]

FitAmpSum FitAmpSum::operator * ( const MINT::counted_ptr< MINT::IReturnComplexForEvent< IDalitzEvent > > & ircfe ) const

Definition at line 526 of file FitAmpSum.cpp.

                                                                                                    {
   FitAmpSum fas(*this);
   fas.multiply(ircfe);
   return fas;
 }

◆ operator *=() [1/4]

FitAmpSum & FitAmpSum::operator *= ( double r )

Definition at line 494 of file FitAmpSum.cpp.

                                         {
   multiply(r);
   return *this;
 }

◆ operator *=() [2/4]

FitAmpSum& FitAmpSum::operator *= ( const std::complex< double > & z )

◆ operator *=() [3/4]

FitAmpSum & FitAmpSum::operator *= ( const MINT::counted_ptr< MINT::IReturnComplex > & irc )

Definition at line 502 of file FitAmpSum.cpp.

                                                                       {
   multiply(irc);
   return *this;
 }

◆ operator *=() [4/4]

FitAmpSum & FitAmpSum::operator *= ( const MINT::counted_ptr< MINT::IReturnComplexForEvent< IDalitzEvent > > & ircfe )

Definition at line 506 of file FitAmpSum.cpp.

                                                                                                {
   multiply(ircfe);
   return *this;
 }

◆ operator+()

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

Definition at line 488 of file FitAmpSum.cpp.

                                                         {
   FitAmpSum fas(*this);
   fas.add(rhs);
   return fas;
 }

◆ operator+=()

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

Definition at line 484 of file FitAmpSum.cpp.

                                                       {
   add(other);
   return *this;
 }

◆ operator=() [1/2]

FitAmpSum & FitAmpSum::operator= ( const FitAmpSum & other )

Definition at line 105 of file FitAmpSum.cpp.

                                                      {
   if(&other == this) return *this;
   (FitAmpList)(*this) = (FitAmpList) (other);
   return *this;
 }

◆ operator=() [2/2]

FitAmpSum & FitAmpSum::operator= ( const FitAmpList & other )

Definition at line 110 of file FitAmpSum.cpp.

                                                       {
   if(&other == this) return *this;
   (FitAmpSum)(*this) = other;
   return *this;
 }

◆ print()

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

virtual

Implements IFastAmplitudeIntegrable.

Definition at line 292 of file FitAmpSum.cpp.

                                          {
   bool dbThis=false;
   os << "FitAmpSum::print " << this->size() << " amplitude components \n"
      << "======================================================" << endl;
   
   for(unsigned int i=0; i< this->size(); i++){
     if(dbThis){
       cout << "ampPtr( " << i << " ): " << endl;
       cout << getAmpPtr(i) << endl;
     }
     os << "\n\t" << getAmpPtr(i)->theBareDecay().oneLiner()
        << endl;
   }
 }

◆ printNonZero()

void FitAmpSum::printNonZero ( std::ostream & os = std::cout ) const

virtual

Reimplemented from FitAmpList.

Definition at line 306 of file FitAmpSum.cpp.

                                                 {
    os << "FitAmpSum::print\n====================";
 
   for(unsigned int i=0; i< this->size(); i++){
     if(getAmpPtr(i)->isZero()) continue;
     os << "\n\t" << getAmpPtr(i)->theBareDecay().oneLiner()
        << endl;
   }
 
 }

◆ printValues()

void FitAmpSum::printValues	(	IDalitzEvent &	evt,
		std::ostream &	os = `std::cout`
	)

virtual

Definition at line 316 of file FitAmpSum.cpp.

                                                             {
    os << "FitAmpSum::print\n====================";
 
   for(unsigned int i=0; i< this->size(); i++){
     complex<double> val = getAmpPtr(i)->getVal(evt);
     os << "\n\t" << getAmpPtr(i)->theBareDecay().oneLiner()
        << val;
     if(norm(val) > 1.e10) os << "  HUGE!!!"  << endl;
     os << endl;
   }
 }

◆ Prob()

virtual double FitAmpSum::Prob ( IDalitzEvent & evt )

inlinevirtual

Definition at line 100 of file FitAmpSum.h.

                                         {
     std::complex<double> res = getVal(evt);
     return res.real()*res.real() + res.imag()*res.imag();
   }

◆ RealVal()

virtual double FitAmpSum::RealVal ( IDalitzEvent & evt )

inlinevirtual

Implements ILookLikeFitAmpSum.

Definition at line 112 of file FitAmpSum.h.

                                            {
     return Prob(evt);
   }

◆ useAnalyticGradient()

virtual bool FitAmpSum::useAnalyticGradient ( )

inlinevirtual

Reimplemented from IFastAmplitudeIntegrable.

Definition at line 89 of file FitAmpSum.h.

89 {return _useAnalyticGradient;}

FitAmpSum::_useAnalyticGradient

MINT::NamedParameter< int > _useAnalyticGradient

Definition: FitAmpSum.h:34

Friends And Related Function Documentation

◆ FitAmplitude

friend class FitAmplitude

friend

Definition at line 139 of file FitAmpSum.h.

Member Data Documentation

◆ _useAnalyticGradient

MINT::NamedParameter<int> FitAmpSum::_useAnalyticGradient

protected

Definition at line 34 of file FitAmpSum.h.

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

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

Public Member Functions

Protected Attributes

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FitAmpSum() [1/5]

◆ FitAmpSum() [2/5]

◆ FitAmpSum() [3/5]

◆ FitAmpSum() [4/5]

◆ FitAmpSum() [5/5]

◆ ~FitAmpSum()

Member Function Documentation

◆ ComplexVal()

◆ getAmpSqr()

◆ GetCloneOfSubsetSameFitParameters()

◆ GetCloneSameFitParameters()

◆ getVal() [1/2]

◆ getVal() [2/2]

◆ Gradient()

◆ makeBoxes()

◆ makeBWBoxes()

◆ makeEventGenerator()

◆ makeFitAmpPairList()

◆ makeIntegCalculator()

◆ makeIntegrationCalculator()

◆ operator *() [1/4]

◆ operator *() [2/4]

◆ operator *() [3/4]

◆ operator *() [4/4]

◆ operator *=() [1/4]

◆ operator *=() [2/4]

◆ operator *=() [3/4]

◆ operator *=() [4/4]

◆ operator+()

◆ operator+=()

◆ operator=() [1/2]

◆ operator=() [2/2]

◆ print()

◆ printNonZero()

◆ printValues()

◆ Prob()

◆ RealVal()

◆ useAnalyticGradient()

Friends And Related Function Documentation

◆ FitAmplitude

Member Data Documentation

◆ _useAnalyticGradient