#include <AmpPair.h>

Inheritance diagram for AmpPair:

Public Member Functions
	AmpPair (Amplitude a1=0, Amplitude a2=0)

	AmpPair (const AmpPair &other)

bool	add (const AmpPair &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

std::complex< double >	add (IDalitzEvent *evt, double weight=1, double efficiency=1)

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

std::complex< double >	integral () const

std::complex< double >	variance () const

double	weightSum () const

long int	N () const

std::complex< double >	lastEntry () const

DalitzHistoSet	histoSet () const

Amplitude *	amp1 ()

Amplitude *	amp2 ()

const Amplitude *	amp1 () const

const Amplitude *	amp2 () const

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

virtual	~AmpPair ()

AmpPair &	operator+= (const AmpPair &other)

AmpPair	operator+ (const AmpPair &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)

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 AmpPair &other) const

Private Attributes
Amplitude *	_a1

Amplitude *	_a2

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

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

Detailed Description

Definition at line 20 of file AmpPair.h.

Constructor & Destructor Documentation

◆ AmpPair() [1/2]

AmpPair::AmpPair	(	Amplitude *	a1 = `0`,
		Amplitude *	a2 = `0`
	)

Definition at line 24 of file AmpPair.cpp.

   : _a1(a1)
   , _a2(a2)
   , _sum(0)
   , _sumName("AmpPair._sum")
   , _sumsq(0)
   , _sumSqName("AmpPair._sumSq")
   , _Nevents(0)
   , _NName("AmpPair._N")
   , _weightSum(0)
   , _weightSumName("AmpPair._weightSum")
   , _hsRe()
   , _hsIm()
   , _name("")
   , _dirName("")
   , _lastEntry(0)
 {
   if(0 == _a1 || 0 == _a2){
     cout << "ERROR in AmpPair::AmpPair !"
      << "\n\t One of the amplitude pointers is zero:"
      << "\n\t _a1 = " << _a1 << ", _a2 = " << _a2
      << "\n\t ... will crash now" << endl;
     throw "zero amp pointer in AmpPair constructor";
   }
 
   std::string name_1 = _a1->name();
   std::string name_2 = _a2->name();
   if(name_1 > name_2) std::swap(_a1, _a2);
   // 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();
 
   registerFitParDependence(*_a1);
   if(_a1 != _a2)  registerFitParDependence(*_a2);
 }

◆ AmpPair() [2/2]

AmpPair::AmpPair ( const AmpPair & other )

Definition at line 63 of file AmpPair.cpp.

   : FitParDependent(other)
   , _a1(other._a1)
   , _a2(other._a2)
   , _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)
 {
 
   removeAllFitParDependencies();
   registerFitParDependence(*_a1);
   if(_a1 != _a2)  registerFitParDependence(*_a2);
 }

◆ ~AmpPair()

virtual AmpPair::~AmpPair ( )

inlinevirtual

Definition at line 109 of file AmpPair.h.

109 {}

Member Function Documentation

◆ add() [1/3]

bool AmpPair::add ( const AmpPair & other )

Definition at line 91 of file AmpPair.cpp.

                                      {
   if(! isCompatibleWith(other)){
     cout << "ERROR in AmpPair::add "
      << "trying to add two incompatible AmpPairs: "
      << "\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/3]

std::complex< double > AmpPair::add	(	IDalitzEvent *	evt,
		double	weight = `1`,
		double	efficiency = `1`
	)

Definition at line 323 of file AmpPair.cpp.

                 {
   bool dbThis=false;
   
   _Nevents++;
   if(0 == evtPtr) return 0;
   
   double ps = evtPtr->phaseSpace();
   if(ps <= 0.0000){
       if(!(_Nevents%100000)){
     cout << "WARNING in AmpPair::addToHistograms"
      << " event with phase space = " << ps << endl;
       }
     return 0; // should not happen.
   }
 
   double w = evtPtr->getWeight()
     /evtPtr->getGeneratorPdfRelativeToPhaseSpace();
   w *= weight;
   
   _weightSum += w;// / ps;
   
   if(dbThis){
     cout << " AmpPair::add, for pair "
      << _a1->name() << " / " << _a2->name()
      << endl;
   }
 
   complex<double> c=ampValue(evtPtr) * efficiency * w;
   _lastEntry = c;
   _sum   += c;
   this->addToHistograms(evtPtr, 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;
 }

◆ add() [3/3]

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

Definition at line 369 of file AmpPair.cpp.

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

◆ addToHistograms()

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

private

Definition at line 110 of file AmpPair.cpp.

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

◆ amp1() [1/2]

Amplitude* AmpPair::amp1 ( )

inline

Definition at line 102 of file AmpPair.h.

102 {return _a1;}

AmpPair::_a1

Amplitude * _a1

Definition: AmpPair.h:22

◆ amp1() [2/2]

const Amplitude* AmpPair::amp1 ( ) const

inline

Definition at line 104 of file AmpPair.h.

104 {return _a1;}

AmpPair::_a1

Amplitude * _a1

Definition: AmpPair.h:22

◆ amp2() [1/2]

Amplitude* AmpPair::amp2 ( )

inline

Definition at line 103 of file AmpPair.h.

103 {return _a2;}

AmpPair::_a2

Amplitude * _a2

Definition: AmpPair.h:23

◆ amp2() [2/2]

const Amplitude* AmpPair::amp2 ( ) const

inline

Definition at line 105 of file AmpPair.h.

105 {return _a2;}

AmpPair::_a2

Amplitude * _a2

Definition: AmpPair.h:23

◆ ampValue()

complex< double > AmpPair::ampValue ( IDalitzEvent * evtPtr )

private

Definition at line 378 of file AmpPair.cpp.

                                                      {
   if(0 == _a1 || 0 == _a2 || 0 == evtPtr) return 0;
 
 
   complex<double> c1 = _a1->getVal(*evtPtr);
   complex<double> c2 = _a2->getVal(*evtPtr);
   complex<double> c2star = conj(c2);
   
   complex<double> val = (c1 * c2star);
 
   return val;
 }

◆ dirName() [1/2]

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

Definition at line 472 of file AmpPair.cpp.

                                        {
   return _dirName;
 }

◆ dirName() [2/2]

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

Definition at line 468 of file AmpPair.cpp.

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

◆ hasMatchingPattern()

bool AmpPair::hasMatchingPattern ( ) const

Definition at line 453 of file AmpPair.cpp.

                                       {
     if(_a1->theBareDecay().getVal() == _a2->theBareDecay().getVal()  )return true;
     return false;
 }

◆ histoImFileName()

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

private

Definition at line 228 of file AmpPair.cpp.

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

◆ histoReFileName()

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

private

Definition at line 225 of file AmpPair.cpp.

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

◆ histoSet()

DalitzHistoSet AmpPair::histoSet ( ) const

◆ histosIm() [1/2]

const DalitzHistoSet& AmpPair::histosIm ( ) const

inlineprivate

Definition at line 49 of file AmpPair.h.

49 {return _hsIm;}

AmpPair::_hsIm

DalitzHistoSet _hsIm

Definition: AmpPair.h:38

◆ histosIm() [2/2]

DalitzHistoSet& AmpPair::histosIm ( )

inlineprivate

Definition at line 51 of file AmpPair.h.

51 {return _hsIm;}

AmpPair::_hsIm

DalitzHistoSet _hsIm

Definition: AmpPair.h:38

◆ histosRe() [1/2]

const DalitzHistoSet& AmpPair::histosRe ( ) const

inlineprivate

Definition at line 48 of file AmpPair.h.

48 {return _hsRe;}

AmpPair::_hsRe

DalitzHistoSet _hsRe

Definition: AmpPair.h:37

◆ histosRe() [2/2]

DalitzHistoSet& AmpPair::histosRe ( )

inlineprivate

Definition at line 50 of file AmpPair.h.

50 {return _hsRe;}

AmpPair::_hsRe

DalitzHistoSet _hsRe

Definition: AmpPair.h:37

◆ integral()

std::complex< double > AmpPair::integral ( ) const

Definition at line 396 of file AmpPair.cpp.

                                           {
   bool dbThis=false;
   if(dbThis){
     cout << " AmpPair::sumWithoutFitPars for pair "
      << _a1->name() << " / " << _a2->name()
      << endl;
   }
   double dN = (double) _Nevents;
   std::complex<double> total = ((std::complex<double>)oneOrTwo()) * _sum;
   std::complex<double> returnVal;
 
   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;
 }

◆ isCompatibleWith()

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

private

Definition at line 87 of file AmpPair.cpp.

                                                         {
   return (amp1()->name() == other.amp1()->name() && 
       amp2()->name() == other.amp2()->name());
 }

◆ isSingleAmp()

bool AmpPair::isSingleAmp ( ) const

Definition at line 449 of file AmpPair.cpp.

                                {
   return _a1 == _a2;
 }

◆ lastEntry()

complex< double > AmpPair::lastEntry ( ) const

Definition at line 375 of file AmpPair.cpp.

                                         {
   return _lastEntry;
 }

◆ makeDirectory()

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

private

Definition at line 159 of file AmpPair.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 & AmpPair::makeDirName ( )

private

Definition at line 128 of file AmpPair.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 & AmpPair::makeName ( )

private

Definition at line 116 of file AmpPair.cpp.

                                   {
   if(0 == amp1() || 0 == amp2()){
     cout << "ERROR in AmpPair::makeName():"
      << " At least one of the amplitude pointers is 0:"
      << " _a1 = " << _a1 << ", _a2 = " << _a2 << endl;
     cout << "I'll crash. " << endl;
     throw "zero amplitude pointers";
   }
   _name = "A(" + amp1()->name() + ")_x_A(" + amp2()->name() + ")*";
   return _name;
 }

◆ N()

long int AmpPair::N ( ) const

Definition at line 429 of file AmpPair.cpp.

                          {
   return _Nevents;
 }

◆ name() [1/2]

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

Definition at line 464 of file AmpPair.cpp.

                                     {
   return _name;
 }

◆ name() [2/2]

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

Definition at line 458 of file AmpPair.cpp.

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

◆ oneOrTwo()

int AmpPair::oneOrTwo ( ) const

private

Definition at line 391 of file AmpPair.cpp.

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

◆ operator+()

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

Definition at line 525 of file AmpPair.cpp.

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

◆ operator+=()

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

Definition at line 521 of file AmpPair.cpp.

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

◆ print()

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

virtual

Definition at line 493 of file AmpPair.cpp.

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

◆ retrieve()

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

Definition at line 206 of file AmpPair.cpp.

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

◆ retrieveHistos()

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

private

Definition at line 238 of file AmpPair.cpp.

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

◆ retrieveValues()

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

private

Definition at line 287 of file AmpPair.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 AmpPair::save ( const std::string & asSubdirOf = "." ) const

Definition at line 194 of file AmpPair.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 AmpPair::saveHistos ( const std::string & asSubdirOf = "." ) const

private

Definition at line 232 of file AmpPair.cpp.

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

◆ saveValues()

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

private

Definition at line 245 of file AmpPair.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 AmpPair::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;
 }

◆ valueFileName()

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

private

Definition at line 222 of file AmpPair.cpp.

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

◆ variance()

std::complex< double > AmpPair::variance ( ) const

Definition at line 433 of file AmpPair.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 ((double) oneOrTwo()) * var;
 }

◆ weightSum()

double AmpPair::weightSum ( ) const

Definition at line 424 of file AmpPair.cpp.

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

Member Data Documentation

◆ _a1

Amplitude* AmpPair::_a1

private

Definition at line 22 of file AmpPair.h.

◆ _a2

Amplitude* AmpPair::_a2

private

Definition at line 23 of file AmpPair.h.

◆ _dirName

std::string AmpPair::_dirName

private

Definition at line 41 of file AmpPair.h.

◆ _hsIm

DalitzHistoSet AmpPair::_hsIm

private

Definition at line 38 of file AmpPair.h.

◆ _hsRe

DalitzHistoSet AmpPair::_hsRe

private

Definition at line 37 of file AmpPair.h.

◆ _lastEntry

std::complex<double> AmpPair::_lastEntry

private

Definition at line 43 of file AmpPair.h.

◆ _name

std::string AmpPair::_name

private

Definition at line 40 of file AmpPair.h.

◆ _Nevents

long int AmpPair::_Nevents

private

Definition at line 31 of file AmpPair.h.

◆ _NName

std::string AmpPair::_NName

private

Definition at line 32 of file AmpPair.h.

◆ _sum

std::complex<double> AmpPair::_sum

private

Definition at line 25 of file AmpPair.h.

◆ _sumName

std::string AmpPair::_sumName

private

Definition at line 26 of file AmpPair.h.

◆ _sumsq

std::complex<double> AmpPair::_sumsq

private

Definition at line 28 of file AmpPair.h.

◆ _sumSqName

std::string AmpPair::_sumSqName

private

Definition at line 29 of file AmpPair.h.

◆ _weightSum

double AmpPair::_weightSum

private

Definition at line 34 of file AmpPair.h.

◆ _weightSumName

std::string AmpPair::_weightSumName

private

Definition at line 35 of file AmpPair.h.

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

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

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ AmpPair() [1/2]

◆ AmpPair() [2/2]

◆ ~AmpPair()

Member Function Documentation

◆ add() [1/3]

◆ add() [2/3]

◆ add() [3/3]

◆ addToHistograms()

◆ amp1() [1/2]

◆ amp1() [2/2]

◆ amp2() [1/2]

◆ amp2() [2/2]

◆ ampValue()

◆ dirName() [1/2]

◆ dirName() [2/2]

◆ hasMatchingPattern()

◆ histoImFileName()

◆ histoReFileName()

◆ histoSet()

◆ 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]

◆ oneOrTwo()

◆ operator+()

◆ operator+=()

◆ print()

◆ retrieve()

◆ retrieveHistos()

◆ retrieveValues()

◆ save()

◆ saveHistos()

◆ saveValues()

◆ valueFileName()

◆ variance()

◆ weightSum()

Member Data Documentation

◆ _a1

◆ _a2

◆ _dirName

◆ _hsIm

◆ _hsRe

◆ _lastEntry

◆ _name

◆ _Nevents

◆ _NName

◆ _sum

◆ _sumName

◆ _sumsq

◆ _sumSqName

◆ _weightSum

◆ _weightSumName