#include <FitAmpPair.h>

Inheritance diagram for FitAmpPair:

Public Member Functions
	FitAmpPair ()

	FitAmpPair (FitAmplitude &a1, FitAmplitude &a2)

	FitAmpPair (const FitAmpPair &other)

bool	add (const FitAmpPair &other)

const std::string &	name () const

const std::string &	name ()

const std::string &	dirName () const

const std::string &	dirName ()

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

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

bool	isSingleAmp () const

bool	hasMatchingPattern () const

double	reAdd (IDalitzEvent &evt, double weight=1, double efficiency=1)

double	add (IDalitzEvent &evt, double weight=1, double efficiency=1)

double	add (IDalitzEvent *evt, double weight=1, double efficiency=1)

double	add (const MINT::counted_ptr< IDalitzEvent > &evt, double weight=1, double efficiency=1)

std::complex< double >	valNoFitPars () const

std::complex< double >	complexVal () const

double	integral () const

double	variance () const

double	weightSum () const

long int	N () const

void	reset ()

bool	needToReIntegrate () const

void	startReIntegration ()

void	startIntegration ()

void	endIntegration ()

bool	acceptEvents () const

std::complex< double >	lastEntry () const

DalitzHistoSet	histoSetRe () const

DalitzHistoSet	histoSetIm () const

DalitzHistoSet	histoSet () const

FitAmplitude &	fitAmp1 ()

FitAmplitude &	fitAmp2 ()

const FitAmplitude &	fitAmp1 () const

const FitAmplitude &	fitAmp2 () const

Amplitude &	rawAmp1 ()

Amplitude &	rawAmp2 ()

const Amplitude &	rawAmp1 () const

const Amplitude &	rawAmp2 () const

void	setSlow ()

void	setFast ()

bool	slow () const

bool	fast () const

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

virtual	~FitAmpPair ()

FitAmpPair &	operator+= (const FitAmpPair &other)

FitAmpPair	operator+ (const FitAmpPair &other) const

Public Member Functions inherited from MINT::FitParDependent
virtual unsigned int	size () const

virtual const FitParRef &	operator[] (unsigned int i) const

virtual FitParRef &	operator[] (unsigned int i)

virtual bool	changedSinceLastCall () const

virtual void	rememberFitParValues ()

virtual bool	registerFitParDependence (const IFitParDependent &fpd)

bool	registerFitParDependence (const FitParRef &fpr)

void	removeAllFitParDependencies ()

	FitParDependent (IFitParRegister *daddy=0)

	FitParDependent (const FitParDependent &other, IFitParRegister *newDaddy=0)

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

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

	PolymorphVector (unsigned int N)

	PolymorphVector (unsigned int N, const FitParRef &c)

	PolymorphVector (const PolymorphVector &other)

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

virtual	~PolymorphVector ()

std::vector< FitParRef > &	theVector ()

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

FitParRef &	operator[] (unsigned int i)

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

FitParRef &	at (unsigned int i)

const FitParRef &	at (unsigned int i) const

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

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

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

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

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

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

FitParRef &	front ()

const FitParRef &	front () const

FitParRef &	back ()

const FitParRef &	back () const

unsigned int	size () const

bool	empty () const

void	push_back (const FitParRef &c)

void	pop_back ()

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

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

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

void	clear ()

void	resize (unsigned int N)

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

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

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

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

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

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

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

Private Member Functions
std::complex< double >	ampValue (IDalitzEvent &evtPtr)

std::complex< double >	fitParValue () const

int	oneOrTwo () const

const DalitzHistoSet &	histosRe () const

const DalitzHistoSet &	histosIm () const

DalitzHistoSet &	histosRe ()

DalitzHistoSet &	histosIm ()

void	addToHistograms (IDalitzEvent *evtPtr, const std::complex< double > &c)

const std::string &	makeName ()

const std::string &	makeDirName ()

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

std::string	valueFileName (const std::string &asSubdirOf) const

std::string	histoReFileName (const std::string &asSubdirOf) const

std::string	histoImFileName (const std::string &asSubdirOf) const

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

bool	retrieveValues (const std::string &fromDirectory=".")

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

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

bool	isCompatibleWith (const FitAmpPair &other) const

void	rememberEventDependentParameters ()

Private Attributes
FitAmplitude *	_fitA1

FitAmplitude *	_fitA2

bool	_slow

bool	_beingIntegrated

MINT::FitParDependent	_eventDependentParameters

std::complex< double >	_sum

std::string	_sumName

std::complex< double >	_sumsq

std::string	_sumSqName

long int	_Nevents

std::string	_NName

double	_weightSum

std::string	_weightSumName

DalitzHistoSet	_hsRe

DalitzHistoSet	_hsIm

std::string	_name

std::string	_dirName

std::complex< double >	_lastEntry

Friends
class	FitAmpPairCovariance

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

Detailed Description

Definition at line 22 of file FitAmpPair.h.

Constructor & Destructor Documentation

◆ FitAmpPair() [1/3]

FitAmpPair::FitAmpPair ( )

Definition at line 24 of file FitAmpPair.cpp.

   : _fitA1(0)
   , _fitA2(0)
   , _slow(true)
   , _beingIntegrated(true)
   , _eventDependentParameters(this)
   , _sum(0)
   , _sumName("FitAmpPair._sum")
   , _sumsq(0)
   , _sumSqName("FitAmpPair._sumSq")
   , _Nevents(0)
   , _NName("FitAmpPair._N")
   , _weightSum(0)
   , _weightSumName("FitAmpPair._weightSum")
   , _hsRe()
   , _hsIm()
   , _name("")
   , _dirName("")
   , _lastEntry(0)
 {
 }

◆ FitAmpPair() [2/3]

FitAmpPair::FitAmpPair	(	FitAmplitude &	a1,
		FitAmplitude &	a2
	)

Definition at line 46 of file FitAmpPair.cpp.

   : _fitA1(&a1)
   , _fitA2(&a2)
   , _slow(true)
   , _beingIntegrated(true)
   , _eventDependentParameters(this)
   , _sum(0)
   , _sumName("FitAmpPair._sum")
   , _sumsq(0)
   , _sumSqName("FitAmpPair._sumSq")
   , _Nevents(0)
   , _NName("FitAmpPair._N")
   , _weightSum(0)
   , _weightSumName("FitAmpPair._weightSum")
   , _hsRe()
   , _hsIm()
   , _name("")
   , _dirName("")
   , _lastEntry(0)
 {
 
   std::string name_1 = _fitA1->name();
   std::string name_2 = _fitA2->name();
   if(name_1 > name_2) std::swap(_fitA1, _fitA2);
   // this ensures that the amplitudes are always stored in 
   // a well-defined order. Determines sign of imagninary part.
   // Since we want Re(A1 A2*), this doesn't matter, however
   // we have to be consistent when we save and retrive and then add to this.
   // Also ensures consistent file and directory naming etc.
   makeName();
   makeDirName();
 
   _eventDependentParameters.registerFitParDependence(rawAmp1());
   _eventDependentParameters.registerFitParDependence(rawAmp2());
   _eventDependentParameters.registerFitParDependence(fitAmp1().eventDependentPrefactors());
   _eventDependentParameters.registerFitParDependence(fitAmp2().eventDependentPrefactors());
 }

◆ FitAmpPair() [3/3]

FitAmpPair::FitAmpPair ( const FitAmpPair & other )

Definition at line 83 of file FitAmpPair.cpp.

   : FitParDependent(other)
   , _fitA1(other._fitA1)
   , _fitA2(other._fitA2)
   , _slow(other._slow)
   , _beingIntegrated(other._beingIntegrated)
   , _eventDependentParameters(other._eventDependentParameters)
   , _sum(other._sum)
   , _sumName(other._sumName)
   , _sumsq(other._sumsq)
   , _sumSqName(other._sumSqName)
   , _Nevents(other._Nevents)
   , _NName(other._NName)
   , _weightSum(other._weightSum)
   , _weightSumName(other._weightSumName)
   , _hsRe(other._hsRe)
   , _hsIm(other._hsIm)
   , _name(other._name)
   , _dirName(other._dirName)
   , _lastEntry(other._lastEntry)
 {
 }

◆ ~FitAmpPair()

virtual FitAmpPair::~FitAmpPair ( )

inlinevirtual

Definition at line 184 of file FitAmpPair.h.

184 {}

Member Function Documentation

◆ acceptEvents()

bool FitAmpPair::acceptEvents ( ) const

Definition at line 520 of file FitAmpPair.cpp.

                                    {
   return _beingIntegrated;
 }

◆ add() [1/4]

bool FitAmpPair::add ( const FitAmpPair & other )

Definition at line 110 of file FitAmpPair.cpp.

                                            {
   if(! isCompatibleWith(other)){
     cout << "ERROR in FitAmpPair::add "
      << "trying to add two incompatible FitAmpPairs: "
      << "\n " << this->name() << " + " << other.name() 
      << endl;
     return false;
   }
   
   _sum        += other._sum;
   _sumsq      += other._sumsq;
   _Nevents    += other._Nevents;
   _weightSum  += other._weightSum;
 
   _hsRe       += other._hsRe;
   _hsIm       += other._hsIm;
 
   return true;
 }

◆ add() [2/4]

double FitAmpPair::add	(	IDalitzEvent &	evt,
		double	weight = `1`,
		double	efficiency = `1`
	)

Definition at line 343 of file FitAmpPair.cpp.

                 {
 //   bool dbThis=false;
   
   _Nevents++;
   
   double ps = evt.phaseSpace();
   if(ps <= 0.0000){
     if(!(_Nevents%100000)){
       cout << "WARNING in FitAmpPair::addToHistograms"
        << " event with phase space = " << ps << endl;
     }
     return 0; // should not happen.
   }
 
   double w = evt.getWeight()
     /evt.getGeneratorPdfRelativeToPhaseSpace();
   w *= weight;
   
   _weightSum += w;// / ps;
   
 //   if(dbThis){
 //     cout << " FitAmpPair::add, for pair "
 //   << fitAmp1().name() << " / " << fitAmp2().name()
 //   << endl;
 //   }
 
   complex<double> c=ampValue(evt) * efficiency * w;
   _lastEntry = c;
   _sum   += c;
 
   if(slow() && 0.0 != c) this->addToHistograms(&evt, c);
 
 //   if(dbThis){
 //     cout << "\t c = " << c
 //   << " _sum " << _sum
 //   << endl;
 //   }
   complex<double> csq(c.real()*c.real(), c.imag()*c.imag());
   _sumsq += csq;
 
   return (c * fitParValue()).real();
 }

◆ add() [3/4]

double FitAmpPair::add	(	IDalitzEvent *	evt,
		double	weight = `1`,
		double	efficiency = `1`
	)

Definition at line 336 of file FitAmpPair.cpp.

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

◆ add() [4/4]

double FitAmpPair::add	(	const MINT::counted_ptr< IDalitzEvent > &	evt,
		double	weight = `1`,
		double	efficiency = `1`
	)

Definition at line 437 of file FitAmpPair.cpp.

                 {
   return add(evtPtr.get(), weight, efficiency);
 }

◆ addToHistograms()

void FitAmpPair::addToHistograms	(	IDalitzEvent *	evtPtr,
		const std::complex< double > &	c
	)

private

Definition at line 129 of file FitAmpPair.cpp.

                                                {
   _hsRe.addEvent(*evtPtr, c.real());
   _hsIm.addEvent(*evtPtr, c.imag());
 }

◆ ampValue()

std::complex<double> FitAmpPair::ampValue ( IDalitzEvent & evtPtr )

inlineprivate

Definition at line 51 of file FitAmpPair.h.

                                                           {
     std::complex<double> c1 = rawAmp1().getVal(evtPtr);
     if(0.0 == c1) return 0;
     std::complex<double> c2 = rawAmp2().getVal(evtPtr);
     return (c1 * conj(c2));  // c1 x c2*
   }

◆ complexVal()

std::complex<double> FitAmpPair::complexVal ( ) const

inline

Definition at line 126 of file FitAmpPair.h.

                                               {
     return valNoFitPars() * fitParValue();
   }

◆ dirName() [1/2]

const std::string & FitAmpPair::dirName ( ) const

Definition at line 613 of file FitAmpPair.cpp.

                                           {
   return _dirName;
 }

◆ dirName() [2/2]

const std::string & FitAmpPair::dirName ( )

Definition at line 609 of file FitAmpPair.cpp.

                                     {
   if("" == _dirName) makeDirName();
   return _dirName;
 }

◆ endIntegration()

void FitAmpPair::endIntegration ( )

Definition at line 515 of file FitAmpPair.cpp.

                                {
   _beingIntegrated=false;
   rememberEventDependentParameters();
 }

◆ fast()

bool FitAmpPair::fast ( ) const

inline

Definition at line 180 of file FitAmpPair.h.

180 {return !_slow;}

FitAmpPair::_slow

bool _slow

Definition: FitAmpPair.h:27

◆ fitAmp1() [1/2]

FitAmplitude& FitAmpPair::fitAmp1 ( )

inline

Definition at line 153 of file FitAmpPair.h.

                          {
     if(0 == _fitA1){
       std::cout << "ERROR in FitAmpPair::fitAmp1: zero pointer " << std::endl;
       throw "zero ptr";
     }
     return *_fitA1;
   }

◆ fitAmp1() [2/2]

const FitAmplitude& FitAmpPair::fitAmp1 ( ) const

inline

Definition at line 167 of file FitAmpPair.h.

167 {return *_fitA1;}

FitAmpPair::_fitA1

FitAmplitude * _fitA1

Definition: FitAmpPair.h:24

◆ fitAmp2() [1/2]

FitAmplitude& FitAmpPair::fitAmp2 ( )

inline

Definition at line 160 of file FitAmpPair.h.

                          {
     if(0 == _fitA2){
       std::cout << "ERROR in FitAmpPair::fitAmp2: zero pointer " << std::endl;
       throw "zero ptr";
     }
     return *_fitA2;
   }

◆ fitAmp2() [2/2]

const FitAmplitude& FitAmpPair::fitAmp2 ( ) const

inline

Definition at line 168 of file FitAmpPair.h.

168 {return *_fitA2;}

FitAmpPair::_fitA2

FitAmplitude * _fitA2

Definition: FitAmpPair.h:25

◆ fitParValue()

std::complex<double> FitAmpPair::fitParValue ( ) const

inlineprivate

Definition at line 58 of file FitAmpPair.h.

                                               {
     std::complex<double> c1 = fitAmp1().AmpPhase();
     if(0.0 == c1) return 0;
     std::complex<double> c2 = fitAmp2().AmpPhase();
     return c1* conj(c2); // c1 x c2*
   }

◆ hasMatchingPattern()

bool FitAmpPair::hasMatchingPattern ( ) const

Definition at line 594 of file FitAmpPair.cpp.

                                          {
     if(_fitA1->theBareDecay().getVal() == _fitA2->theBareDecay().getVal()  )return true;
     return false;
 }

◆ histoImFileName()

std::string FitAmpPair::histoImFileName ( const std::string & asSubdirOf ) const

private

Definition at line 240 of file FitAmpPair.cpp.

                                                                      {
   return asSubdirOf + "/" + dirName() + "/histoIm";
 }

◆ histoReFileName()

std::string FitAmpPair::histoReFileName ( const std::string & asSubdirOf ) const

private

Definition at line 237 of file FitAmpPair.cpp.

                                                                      {
   return asSubdirOf + "/" + dirName() + "/histoRe";
 }

◆ histoSet()

DalitzHistoSet FitAmpPair::histoSet ( ) const

inline

Definition at line 151 of file FitAmpPair.h.

151 {return histoSetRe();} // ..this is.

FitAmpPair::histoSetRe

DalitzHistoSet histoSetRe() const

Definition: FitAmpPair.cpp:529

◆ histoSetIm()

DalitzHistoSet FitAmpPair::histoSetIm ( ) const

Definition at line 551 of file FitAmpPair.cpp.

                                            {
   // needed by Lauren's code, only
   bool dbThis=false;
   if(dbThis){
     cout << " FitAmpPair::histograms, for pair "
      << fitAmp1().name() << " / " << fitAmp2().name()
      << endl;
     cout << " r " << fitParValue().real() << ", " 
      << histosRe().begin()->second.histo()->Integral()
      << " i " << fitParValue().imag() << ", " 
      << histosIm().begin()->second.histo()->Integral()
      << endl;
   }
   
   return ( // this calculates a DalitzHistoSet:
       oneOrTwo() * 
       (fitParValue().real() * histosIm() 
        + fitParValue().imag() * histosRe())
        )/_weightSum;// / ((double)_Nevents);
 
   //  total /= (double)_Nevents;
 }

◆ histoSetRe()

DalitzHistoSet FitAmpPair::histoSetRe ( ) const

Definition at line 529 of file FitAmpPair.cpp.

                                            {
   bool dbThis=false;
   if(dbThis){
     cout << " FitAmpPair::histograms, for pair "
      << fitAmp1().name() << " / " << fitAmp2().name()
      << endl;
     cout << " r " << fitParValue().real() << ", " 
      << histosRe().begin()->second.histo()->Integral()
      << " i " << fitParValue().imag() << ", " 
      << histosIm().begin()->second.histo()->Integral()
      << endl;
   }
 
   
   return ( // this calculates a DalitzHistoSet:
       oneOrTwo() * 
       (fitParValue().real() * histosRe() 
        - fitParValue().imag() * histosIm())
        )/_weightSum;// / ((double)_Nevents);
 
   //  total /= (double)_Nevents;
 }

◆ histosIm() [1/2]

const DalitzHistoSet& FitAmpPair::histosIm ( ) const

inlineprivate

Definition at line 74 of file FitAmpPair.h.

74 {return _hsIm;}

FitAmpPair::_hsIm

DalitzHistoSet _hsIm

Definition: FitAmpPair.h:44

◆ histosIm() [2/2]

DalitzHistoSet& FitAmpPair::histosIm ( )

inlineprivate

Definition at line 76 of file FitAmpPair.h.

76 {return _hsIm;}

FitAmpPair::_hsIm

DalitzHistoSet _hsIm

Definition: FitAmpPair.h:44

◆ histosRe() [1/2]

const DalitzHistoSet& FitAmpPair::histosRe ( ) const

inlineprivate

Definition at line 73 of file FitAmpPair.h.

73 {return _hsRe;}

FitAmpPair::_hsRe

DalitzHistoSet _hsRe

Definition: FitAmpPair.h:43

◆ histosRe() [2/2]

DalitzHistoSet& FitAmpPair::histosRe ( )

inlineprivate

Definition at line 75 of file FitAmpPair.h.

75 {return _hsRe;}

FitAmpPair::_hsRe

DalitzHistoSet _hsRe

Definition: FitAmpPair.h:43

◆ integral()

double FitAmpPair::integral ( ) const

Definition at line 476 of file FitAmpPair.cpp.

                                  {
 
   return complexVal().real();
 }

◆ isCompatibleWith()

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

private

Definition at line 106 of file FitAmpPair.cpp.

                                                               {
   return (fitAmp1().name() == other.fitAmp1().name() && 
       fitAmp2().name() == other.fitAmp2().name());
 }

◆ isSingleAmp()

bool FitAmpPair::isSingleAmp ( ) const

Definition at line 590 of file FitAmpPair.cpp.

                                   {
   return _fitA1 == _fitA2;
 }

◆ lastEntry()

complex< double > FitAmpPair::lastEntry ( ) const

Definition at line 443 of file FitAmpPair.cpp.

                                            {
   return _lastEntry;
 }

◆ makeDirectory()

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

private

Definition at line 171 of file FitAmpPair.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);
 }

◆ makeDirName()

const std::string & FitAmpPair::makeDirName ( )

private

Definition at line 140 of file FitAmpPair.cpp.

                                         {
   _dirName="";
   if("" == _name) makeName();
 
   std::string::const_iterator preEnd = _name.end();
   preEnd--;
   
   for(std::string::const_iterator it = _name.begin();
       it != _name.end(); it++){
     if(it != preEnd){
       std::string::const_iterator next = it; next++;
       if('-' == *it && '>' == *next){
     _dirName += "_to_";
     it++;
     continue;
       }
     }
     if('*' == *it) _dirName += "star";
     else if(' ' == *it) _dirName += "_";
     //    else if('(' == *it) _dirName +="_bra_";
     //    else if(')'==*it) _dirName +="_ket_";
     else if('['==*it) _dirName += "_sqbra_";
     else if(']'==*it) _dirName += "_sqket_";
     //    else if(','==*it) _dirName += "_comma_";
     //    else if('-' == *it) _dirName += "min";
     //    else if('+' == *it) _dirName += "pls";
     else _dirName += *it;
   }
   return _dirName;
 }

◆ makeName()

const std::string & FitAmpPair::makeName ( )

private

Definition at line 135 of file FitAmpPair.cpp.

                                      {
   _name = "A(" + fitAmp1().name() + ")_x_fitA(" + fitAmp2().name() + ")*";
   return _name;
 }

◆ N()

long int FitAmpPair::N ( ) const

Definition at line 485 of file FitAmpPair.cpp.

                             {
   return _Nevents;
 }

◆ name() [1/2]

const std::string & FitAmpPair::name ( ) const

Definition at line 605 of file FitAmpPair.cpp.

                                        {
   return _name;
 }

◆ name() [2/2]

const std::string & FitAmpPair::name ( )

Definition at line 599 of file FitAmpPair.cpp.

                                  {
   if("" == _name) makeName();
   return _name;
   //  if(isSingleAmp()) return "| " + _fitA1->name() + " |^2";
   //  else return "(" + _fitA1->name() +") * (" + _fitA2->name() + ")*";
 }

◆ needToReIntegrate()

bool FitAmpPair::needToReIntegrate ( ) const

Definition at line 500 of file FitAmpPair.cpp.

                                         {
   return _eventDependentParameters.changedSinceLastCall();
 }

◆ oneOrTwo()

int FitAmpPair::oneOrTwo ( ) const

inlineprivate

Definition at line 65 of file FitAmpPair.h.

                             {
     if(isSingleAmp()) return 1;
     else return 2;
   }  

◆ operator+()

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

Definition at line 666 of file FitAmpPair.cpp.

                                                              {
   FitAmpPair returnVal(*this);
   returnVal += other;
   return returnVal;
 }

◆ operator+=()

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

Definition at line 662 of file FitAmpPair.cpp.

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

◆ print()

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

virtual

Definition at line 634 of file FitAmpPair.cpp.

                                           {
   os << "FitAmpPair " << name()
      << ": N=" << N() << ", integral=" << integral();
 }

◆ rawAmp1() [1/2]

Amplitude& FitAmpPair::rawAmp1 ( )

inline

Definition at line 170 of file FitAmpPair.h.

170 {return fitAmp1().amp();}

FitAmpPair::fitAmp1

FitAmplitude & fitAmp1()

Definition: FitAmpPair.h:153

FitAmplitude::amp

Amplitude & amp()

Definition: FitAmplitude.h:163

◆ rawAmp1() [2/2]

const Amplitude& FitAmpPair::rawAmp1 ( ) const

inline

Definition at line 172 of file FitAmpPair.h.

172 {return fitAmp1().amp();}

FitAmpPair::fitAmp1

FitAmplitude & fitAmp1()

Definition: FitAmpPair.h:153

FitAmplitude::amp

Amplitude & amp()

Definition: FitAmplitude.h:163

◆ rawAmp2() [1/2]

Amplitude& FitAmpPair::rawAmp2 ( )

inline

Definition at line 171 of file FitAmpPair.h.

171 {return fitAmp2().amp();}

FitAmpPair::fitAmp2

FitAmplitude & fitAmp2()

Definition: FitAmpPair.h:160

FitAmplitude::amp

Amplitude & amp()

Definition: FitAmplitude.h:163

◆ rawAmp2() [2/2]

const Amplitude& FitAmpPair::rawAmp2 ( ) const

inline

Definition at line 173 of file FitAmpPair.h.

173 {return fitAmp2().amp();}

FitAmpPair::fitAmp2

FitAmplitude & fitAmp2()

Definition: FitAmpPair.h:160

FitAmplitude::amp

Amplitude & amp()

Definition: FitAmplitude.h:163

◆ reAdd()

double FitAmpPair::reAdd	(	IDalitzEvent &	evt,
		double	weight = `1`,
		double	efficiency = `1`
	)

Definition at line 389 of file FitAmpPair.cpp.

                 {
 //   bool dbThis=false;
   
   _Nevents++;
   
   double ps = evt.phaseSpace();
   if(ps <= 0.0000){
     if(!(_Nevents%100000)){
       cout << "WARNING in FitAmpPair::addToHistograms"
        << " event with phase space = " << ps << endl;
     }
     return 0; // should not happen.
   }
 
   double w = evt.getWeight()
     /evt.getGeneratorPdfRelativeToPhaseSpace();
   w *= weight;
   
   _weightSum += w;// / ps;
   
 //   if(dbThis){
 //     cout << " FitAmpPair::add, for pair "
 //   << fitAmp1().name() << " / " << fitAmp2().name()
 //   << endl;
 //   }
 
   complex<double> c=ampValue(evt) * efficiency * w;
   _lastEntry = c;
   _sum   += c;
 
   if(slow() && 0.0 != c) this->addToHistograms(&evt, c);
 
 //   if(dbThis){
 //     cout << "\t c = " << c
 //   << " _sum " << _sum
 //   << endl;
 //   }
 
   complex<double> csq(c.real()*c.real(), c.imag()*c.imag());
   _sumsq += csq;
 
   return (c * fitParValue()).real();
 }

◆ rememberEventDependentParameters()

void FitAmpPair::rememberEventDependentParameters ( )

private

Definition at line 503 of file FitAmpPair.cpp.

                                                  {
   _eventDependentParameters.rememberFitParValues();
 }

◆ reset()

void FitAmpPair::reset ( )

Definition at line 489 of file FitAmpPair.cpp.

                       {
   _beingIntegrated = true;
   _sum = 0;
   _sumsq = 0;
   _Nevents = 0;
   _weightSum = 0;
   if(slow()) _hsRe.clearAllHistos();
   if(slow()) _hsIm.clearAllHistos();
   _lastEntry = 0;
 }

◆ retrieve()

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

Definition at line 218 of file FitAmpPair.cpp.

                                                     {
   bool dbThis=false;
   bool success=true;
 
   success &= retrieveValues(asSubdirOf);
   if(slow()) success &= retrieveHistos(asSubdirOf);
 
   if(dbThis && slow()){
     cout << "after FitAmpPair::retrieve: pat = "
      << histosRe().begin()->second.pattern() << ", "
      << histosIm().begin()->second.pattern() << endl;
   }
   return success;
 }

◆ retrieveHistos()

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

private

Definition at line 250 of file FitAmpPair.cpp.

                                                           {
   bool sc = true;
   sc |= histosRe().retrieveFromDir(histoReFileName(asSubdirOf));
   sc |= histosIm().retrieveFromDir(histoImFileName(asSubdirOf));
   return sc;
 }

◆ retrieveValues()

bool FitAmpPair::retrieveValues ( const std::string & fromDirectory = "." )

private

Definition at line 299 of file FitAmpPair.cpp.

                                                              {
   bool dbThis=false;
   std::string fname = valueFileName(fromDirectory);
   if(dbThis)cout << "trying to retreive values from: " << fname << endl;
   if(dbThis)cout << "current sum " << _sum << endl;
 
   NamedParameter<double> n_sumRe(_sumName + ".real()", fname.c_str()
                  );
   n_sumRe.reloadFile(fname.c_str());
   NamedParameter<double> n_sumIm(_sumName + ".imag()", fname.c_str()
                  );//,NamedParameterBase::QUIET);
   NamedParameter<double> n_sumSqRe(_sumSqName + ".real()", fname.c_str()
                    );//, NamedParameterBase::QUIET);
   NamedParameter<double> n_sumSqIm(_sumSqName + ".imag()", fname.c_str()
                    );//, NamedParameterBase::QUIET);
   NamedParameter<long int> n_Nevents(_NName, fname.c_str()
                      );// , NamedParameterBase::QUIET);
   NamedParameter<double> n_weightSum(_weightSumName, fname.c_str()
                      );//, NamedParameterBase::QUIET);
 
   complex<double> n_sum(n_sumRe, n_sumIm);
 
   if(dbThis) cout << "n_sumRe " << n_sumRe << "n_sum " << n_sum << endl;
   _sum = n_sum;
   if(dbThis) cout << "_sum = " << _sum << endl;
   complex<double> n_sumsq(n_sumSqRe, n_sumSqIm);
   _sumsq = n_sumsq;
   _Nevents = n_Nevents;
   _weightSum = n_weightSum;
 
   return true;
 }

◆ save()

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

Definition at line 206 of file FitAmpPair.cpp.

                                                       {
   bool success=true;
   makeDirectory(asSubdirOf); // ignore error codes from mkdir, they could
               // be because directory already exists - if
               // not, we'll pick it up when we try to save
               // the files.
 
   success &= saveValues(asSubdirOf);
   success &= saveHistos(asSubdirOf);
 
   return success;
 }

◆ saveHistos()

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

private

Definition at line 244 of file FitAmpPair.cpp.

                                                             {
   bool sc = true;
   sc |= histosRe().saveAsDir(histoReFileName(asSubdirOf));
   sc |= histosIm().saveAsDir(histoImFileName(asSubdirOf));
   return sc;
 }

◆ saveValues()

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

private

Definition at line 257 of file FitAmpPair.cpp.

                                                             {
 
   NamedParameter<double> n_sumRe(_sumName + ".real()"
                  , NamedParameterBase::QUIET);
   NamedParameter<double> n_sumIm(_sumName + ".imag()"
                  ,NamedParameterBase::QUIET);
   NamedParameter<double> n_sumSqRe(_sumSqName + ".real()"
                    , NamedParameterBase::QUIET);
   NamedParameter<double> n_sumSqIm(_sumSqName + ".imag()"
                    , NamedParameterBase::QUIET);
   NamedParameter<long int> n_Nevents(_NName
                      , NamedParameterBase::QUIET);
   NamedParameter<double> n_weightSum(_weightSumName
                      , NamedParameterBase::QUIET);
 
   n_sumRe     = _sum.real();
   n_sumIm     = _sum.imag();
   n_sumSqRe   = _sumsq.real();
   n_sumSqIm   = _sumsq.imag();
   n_Nevents   = _Nevents;
   n_weightSum = _weightSum;
 
   std::string fname = valueFileName(asSubdirOf);
   ofstream os(fname.c_str());
   if(os.bad()){
     cout << "ERROR in FitAmpPair::saveValues of \n\t" << name()
      << "\n\t unable to create file: "
      << "\n\t" << fname << endl;
     return false;
   }
   os << setprecision(20)
      << name() 
      << '\n' << n_sumRe 
      << '\n' << n_sumIm
      << '\n' << n_sumSqRe 
      << '\n' << n_sumSqIm
      << '\n' << n_Nevents
      << '\n' << n_weightSum
      << endl;
   os.close();
   return true;
 }

◆ setFast()

void FitAmpPair::setFast ( )

inline

Definition at line 177 of file FitAmpPair.h.

177 {_slow = false;}

FitAmpPair::_slow

bool _slow

Definition: FitAmpPair.h:27

◆ setSlow()

void FitAmpPair::setSlow ( )

inline

Definition at line 176 of file FitAmpPair.h.

176 {_slow = true;}

FitAmpPair::_slow

bool _slow

Definition: FitAmpPair.h:27

◆ slow()

bool FitAmpPair::slow ( ) const

inline

Definition at line 179 of file FitAmpPair.h.

179 {return _slow;}

FitAmpPair::_slow

bool _slow

Definition: FitAmpPair.h:27

◆ startIntegration()

void FitAmpPair::startIntegration ( )

Definition at line 511 of file FitAmpPair.cpp.

                                  {
   _beingIntegrated=true;
 }

◆ startReIntegration()

void FitAmpPair::startReIntegration ( )

Definition at line 507 of file FitAmpPair.cpp.

                                    {
   reset();
   _beingIntegrated=true;
 }

◆ valNoFitPars()

std::complex< double > FitAmpPair::valNoFitPars ( ) const

Definition at line 447 of file FitAmpPair.cpp.

                                                  {
 /*  bool dbThis=false;
   if(dbThis){
     cout << " FitAmpPair::sumWithoutFitPars for pair "
      << fitAmp1().name() << " / " << fitAmp2().name()
      << endl;
   }*/
   double dN = (double) _Nevents;
   std::complex<double> total = ((std::complex<double>)oneOrTwo()) * _sum;
   std::complex<double> returnVal;
 
   if(0.0 == total) return 0.0;
 
   if(_weightSum > 0){
     returnVal = total/_weightSum;
   }else{
     returnVal = total/dN;
   }
 
 /*  if(dbThis){
     cout << "\t returning " << returnVal << endl;
     cout << "\t   = Real( " << oneOrTwo()
      << " * " << _sum
      << " / " << _Nevents
      << " )"
      << endl;
   }*/  
   return returnVal;
 }

◆ valueFileName()

std::string FitAmpPair::valueFileName ( const std::string & asSubdirOf ) const

private

Definition at line 234 of file FitAmpPair.cpp.

                                                                    {
   return asSubdirOf + "/" + dirName() + "/value.txt";
 }

◆ variance()

double FitAmpPair::variance ( ) const

Definition at line 574 of file FitAmpPair.cpp.

                                  {
   if(_Nevents <=0) return 0;
   double dN = (double) _Nevents;
   complex<double> mean   = _sum/dN;
   complex<double> meansq = _sumsq/dN;
   
   complex<double> var(meansq.real() - mean.real()*mean.real()
               , meansq.imag() - mean.imag()*mean.imag());
 
   var /= dN;
 
   if(_weightSum > 0) var *= dN*dN/(_weightSum*_weightSum);
 
   return oneOrTwo() * (fitParValue() * var).real();
 }

◆ weightSum()

double FitAmpPair::weightSum ( ) const

Definition at line 480 of file FitAmpPair.cpp.

                                   {
   if(_weightSum > 0) return _weightSum;
   else return (double) _Nevents;
 }

Friends And Related Function Documentation

◆ FitAmpPairCovariance

friend class FitAmpPairCovariance

friend

Definition at line 189 of file FitAmpPair.h.

Member Data Documentation

◆ _beingIntegrated

bool FitAmpPair::_beingIntegrated

private

Definition at line 28 of file FitAmpPair.h.

◆ _dirName

std::string FitAmpPair::_dirName

private

Definition at line 47 of file FitAmpPair.h.

◆ _eventDependentParameters

MINT::FitParDependent FitAmpPair::_eventDependentParameters

private

Definition at line 29 of file FitAmpPair.h.

◆ _fitA1

FitAmplitude* FitAmpPair::_fitA1

private

Definition at line 24 of file FitAmpPair.h.

◆ _fitA2

FitAmplitude* FitAmpPair::_fitA2

private

Definition at line 25 of file FitAmpPair.h.

◆ _hsIm

DalitzHistoSet FitAmpPair::_hsIm

private

Definition at line 44 of file FitAmpPair.h.

◆ _hsRe

DalitzHistoSet FitAmpPair::_hsRe

private

Definition at line 43 of file FitAmpPair.h.

◆ _lastEntry

std::complex<double> FitAmpPair::_lastEntry

private

Definition at line 49 of file FitAmpPair.h.

◆ _name

std::string FitAmpPair::_name

private

Definition at line 46 of file FitAmpPair.h.

◆ _Nevents

long int FitAmpPair::_Nevents

private

Definition at line 37 of file FitAmpPair.h.

◆ _NName

std::string FitAmpPair::_NName

private

Definition at line 38 of file FitAmpPair.h.

◆ _slow

bool FitAmpPair::_slow

private

Definition at line 27 of file FitAmpPair.h.

◆ _sum

std::complex<double> FitAmpPair::_sum

private

Definition at line 31 of file FitAmpPair.h.

◆ _sumName

std::string FitAmpPair::_sumName

private

Definition at line 32 of file FitAmpPair.h.

◆ _sumsq

std::complex<double> FitAmpPair::_sumsq

private

Definition at line 34 of file FitAmpPair.h.

◆ _sumSqName

std::string FitAmpPair::_sumSqName

private

Definition at line 35 of file FitAmpPair.h.

◆ _weightSum

double FitAmpPair::_weightSum

private

Definition at line 40 of file FitAmpPair.h.

◆ _weightSumName

std::string FitAmpPair::_weightSumName

private

Definition at line 41 of file FitAmpPair.h.

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

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

Public Member Functions

Private Member Functions

Private Attributes

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FitAmpPair() [1/3]

◆ FitAmpPair() [2/3]

◆ FitAmpPair() [3/3]

◆ ~FitAmpPair()

Member Function Documentation

◆ acceptEvents()

◆ add() [1/4]

◆ add() [2/4]

◆ add() [3/4]

◆ add() [4/4]

◆ addToHistograms()

◆ ampValue()

◆ complexVal()

◆ dirName() [1/2]

◆ dirName() [2/2]

◆ endIntegration()

◆ fast()

◆ fitAmp1() [1/2]

◆ fitAmp1() [2/2]

◆ fitAmp2() [1/2]

◆ fitAmp2() [2/2]

◆ fitParValue()

◆ hasMatchingPattern()

◆ histoImFileName()

◆ histoReFileName()

◆ histoSet()

◆ histoSetIm()

◆ histoSetRe()

◆ histosIm() [1/2]

◆ histosIm() [2/2]

◆ histosRe() [1/2]

◆ histosRe() [2/2]

◆ integral()

◆ isCompatibleWith()

◆ isSingleAmp()

◆ lastEntry()

◆ makeDirectory()

◆ makeDirName()

◆ makeName()

◆ N()

◆ name() [1/2]

◆ name() [2/2]

◆ needToReIntegrate()

◆ oneOrTwo()

◆ operator+()

◆ operator+=()

◆ print()

◆ rawAmp1() [1/2]

◆ rawAmp1() [2/2]

◆ rawAmp2() [1/2]

◆ rawAmp2() [2/2]

◆ reAdd()

◆ rememberEventDependentParameters()

◆ reset()

◆ retrieve()

◆ retrieveHistos()

◆ retrieveValues()

◆ save()

◆ saveHistos()

◆ saveValues()

◆ setFast()

◆ setSlow()

◆ slow()

◆ startIntegration()

◆ startReIntegration()

◆ valNoFitPars()

◆ valueFileName()

◆ variance()

◆ weightSum()

Friends And Related Function Documentation

◆ FitAmpPairCovariance

Member Data Documentation

◆ _beingIntegrated