CoherenceFactor Class Reference

#include <CoherenceFactor.h>

Public Member Functions
	CoherenceFactor (FitAmpSum &A, FitAmpSum &Abar, double CSAbs=1, double CSPhase=0.0, MINT::IReturnRealForEvent< IDalitzEvent > *eff=0, double prec=1.e-3)
std::complex< double >	var () const
std::complex< double >	sigma () const
double	varAbs () const
double	sigmaAbs () const
double	varAngle () const
double	sigmaAngle () const
void	printResult (std::ostream &os=std::cout) const
void	setPrecision (double prec)
double	estimatedPrecision () const
std::complex< double >	getCVal ()

Private Member Functions
MINT::IReturnComplexForEvent< IDalitzEvent > *	A ()
MINT::IReturnComplexForEvent< IDalitzEvent > *	Abar ()
std::complex< double >	A_Value (IDalitzEvent &evt)
std::complex< double >	Abar_Value (IDalitzEvent &evt)
MINT::counted_ptr< MINT::IEventGenerator< IDalitzEvent > >	getGenerator ()
MINT::counted_ptr< IDalitzEvent >	newEvent ()
std::complex< double >	integAxAbarStar () const
double	integA () const
double	integAbar () const
double	integASq () const
double	integAbarSq () const
std::complex< double >	complexVar (std::complex< double > s, std::complex< double > s_sq) const
double	realAorAbarVar (double s, double sq) const
std::complex< double >	cvarAxAbarStar () const
double	varA () const
double	varAbar () const
std::complex< double >	Rval () const
double	getEff (IDalitzEvent &evt)
double	getEff (MINT::counted_ptr< IDalitzEvent > evtPtr)

Private Attributes
FitAmpSum *	_A
FitAmpSum *	_Abar
double	_CSAbs
double	_CSPhase
FitAmpSum	_A_plus_Abar
MINT::counted_ptr< MINT::IEventGenerator< IDalitzEvent > >	_myOwnGenAplusAbar
MINT::IReturnRealForEvent< IDalitzEvent > *	_eff
double	_precision
double	_sumASq
double	_sumASqSquared
double	_sumAbarSq
double	_sumAbarSqSquared
std::complex< double >	_sumAxAbarStar
std::complex< double >	_sumAxAbarStarSquared
long int	_Nevents

Detailed Description

Definition at line 25 of file CoherenceFactor.h.

Constructor & Destructor Documentation

CoherenceFactor()

CoherenceFactor::CoherenceFactor	(	FitAmpSum &	A,
		FitAmpSum &	Abar,
		double	CSAbs = 1,
		double	CSPhase = 0.0,
		MINT::IReturnRealForEvent< IDalitzEvent > *	eff = 0,
		double	prec = 1.e-3
	)

Definition at line 17 of file CoherenceFactor.cpp.

  : _A(&A)
  , _Abar(&Abar)
  , _CSAbs(CSAbs)
  , _CSPhase(CSPhase)
  , _A_plus_Abar(A+Abar)
  , _myOwnGenAplusAbar(0)
  , _eff(eff)
  , _precision(prec)
  , _sumASq(0), _sumASqSquared(0)
  , _sumAbarSq(0), _sumAbarSqSquared(0)
  , _sumAxAbarStar(0,0), _sumAxAbarStarSquared(0,0)
  , _Nevents(0)
{
}

Member Function Documentation

A()

MINT::IReturnComplexForEvent<IDalitzEvent>* CoherenceFactor::A ( )

inlineprivate

Definition at line 45 of file CoherenceFactor.h.

                                               {
    return (MINT::IReturnComplexForEvent<IDalitzEvent>*)_A;}

A_Value()

complex< double > CoherenceFactor::A_Value ( IDalitzEvent &  evt )

private

Definition at line 205 of file CoherenceFactor.cpp.

                                                         {
  //  bool dbThis=false;
  return A()->ComplexVal(evt);
}

Abar()

MINT::IReturnComplexForEvent<IDalitzEvent>* CoherenceFactor::Abar ( )

inlineprivate

Definition at line 47 of file CoherenceFactor.h.

                                                  {
    return (MINT::IReturnComplexForEvent<IDalitzEvent>*) _Abar;}

Abar_Value()

complex< double > CoherenceFactor::Abar_Value ( IDalitzEvent &  evt )

private

Definition at line 209 of file CoherenceFactor.cpp.

                                                            {
  bool dbThis=false;
  if(dbThis) cout << "CoherenceFactor::Abar_Value" << endl;
  complex<double> val = Abar()->ComplexVal(evt);
  if(dbThis) cout << "got value: " << val << endl;
  return val;
}

complexVar()

complex< double > CoherenceFactor::complexVar ( std::complex< double >  s, std::complex< double >  s_sq  ) const

private

Definition at line 110 of file CoherenceFactor.cpp.

                              {
  if(_Nevents <=0) return 0;
  complex<double> m   = s        / ((double) _Nevents);
  complex<double> mos = s_sq / ((double) _Nevents);
  complex<double> som(m.real()*m.real(), m.imag()*m.imag());
  return (mos - som)/((double)_Nevents);
}

cvarAxAbarStar()

complex< double > CoherenceFactor::cvarAxAbarStar ( ) const

private

Definition at line 91 of file CoherenceFactor.cpp.

                                                    {
  return complexVar(_sumAxAbarStar, _sumAxAbarStarSquared);
}

estimatedPrecision()

double CoherenceFactor::estimatedPrecision

(

)

const

Definition at line 61 of file CoherenceFactor.cpp.

                                                {
  if(_Nevents <= 0) return 9999;
  double s = sigmaAbs();
  if(s < 0) s *=-1;
  return s;
}

getCVal()

complex< double > CoherenceFactor::getCVal

(

)

Definition at line 217 of file CoherenceFactor.cpp.

                                        {
  bool dbThis=false;
  if(dbThis) cout << "Hello from CoherenceFactor::getCVal()" << endl;
  long int rediculousNumber = 1000000000;
  long int tries=0;

  int printEvery = 10000;
  
  while(estimatedPrecision() > _precision && tries < rediculousNumber){
    for(int i=0; i < 1000; i++){
      tries++;
      bool printout = tries < 10 || (0 == tries%printEvery);
      
      counted_ptr<IDalitzEvent> evtPtr(newEvent());
      
      if(0 == evtPtr){
    if(dbThis) cout << "WARNING in CoherenceFactor::getR()"
            << " got evtPtr = " << evtPtr << endl;
    continue;
      }
      //cout << "evtPtr->phaseSpace() " << evtPtr->phaseSpace() << endl;

      double w = evtPtr->getWeight()
    /evtPtr->getGeneratorPdfRelativeToPhaseSpace();
      
      w*= getEff(*evtPtr);
      
      w = sqrt(w); // since we apply the weight to the amplitudes, 
      //              while the integrals (as one might expect) all 
      //              contain terms |A|^2, |Abar|^2, or |A*Abar|.

      complex<double> c_a=
    A_Value(*evtPtr) * w;

      complex<double> c_abar=
    Abar_Value(*evtPtr) *w*polar(_CSAbs, _CSPhase);
      
      double d_a    = norm(c_a);
      double d_abar = norm(c_abar);

      _sumASq += d_a;
      _sumASqSquared += d_a*d_a;

      _sumAbarSq += d_abar;
      _sumAbarSqSquared += d_abar*d_abar;
      
      complex<double> c_prod = c_a * conj(c_abar);
      _sumAxAbarStar += c_prod;
      complex<double> c_prodSquared(c_prod.real() * c_prod.real()
                    , c_prod.imag() * c_prod.imag()
                    );
      _sumAxAbarStarSquared += c_prodSquared;
      
      _Nevents++;
      
      if(printout){
    cout << " DalitzEvent::EventCounter() "
         << DalitzEvent::eventCounter()
         << endl;
    cout << "INFO from CoherenceFactor::getR()"
         << "\n\t tries: " << tries << ",  for " << _Nevents << " events."
         << endl;
      // << "\n\t c_a, c_abar: " << c_a << ", " << c_abar
      //  << "\n\t value: " << Rval()
      //  << " = " << abs(Rval()) << " * exp( i * " << arg(Rval()) << " )"
      //  << " +/- " << estimatedPrecision()
      //  << endl;
    printResult();
    cout << " --------------------------------- " << endl;
      }
    }
  }
  return Rval();
}

getEff() [1/2]

double CoherenceFactor::getEff ( IDalitzEvent &  evt )

private

Definition at line 57 of file CoherenceFactor.cpp.

                                               {
  if(0 == _eff) return 1.0;
  return _eff->RealVal(evt);
}

getEff() [2/2]

double CoherenceFactor::getEff ( MINT::counted_ptr< IDalitzEvent >  evtPtr )

inlineprivate

Definition at line 75 of file CoherenceFactor.h.

                                                     {
    if(0 == evtPtr) return 0;
    return getEff(*evtPtr);// just for backward compatibility.
  }

getGenerator()

counted_ptr< IEventGenerator< IDalitzEvent > > CoherenceFactor::getGenerator ( )

private

Definition at line 39 of file CoherenceFactor.cpp.

                             {
  bool dbThis=false;
  if(0 == _myOwnGenAplusAbar){
    if(dbThis) cout << "CoherenceFactor::getGenerator():  "
            << "making generator" << endl;
    if(_A_plus_Abar.size() <= 0) return 0;
    DalitzEventPattern pat = _A_plus_Abar.getAmpPtr(0)->getTreePattern();
    counted_ptr< IEventGenerator<IDalitzEvent> > 
      ptr( new DalitzBWBoxSet(_A_plus_Abar.makeBWBoxes(pat)));
    if(dbThis) cout << "made generator." << endl;
    _myOwnGenAplusAbar = ptr;
  }
  return _myOwnGenAplusAbar;
}

integA()

double CoherenceFactor::integA ( ) const

private

Definition at line 72 of file CoherenceFactor.cpp.

                                   {
  if(_Nevents <=0) return 0;
  if(_sumASq <=0) return -9999;
  return sqrt(_sumASq / ((double) _Nevents));
}

integAbar()

double CoherenceFactor::integAbar ( ) const

private

Definition at line 77 of file CoherenceFactor.cpp.

                                      {
  if(_Nevents <=0) return 0;
  if(_sumAbarSq <=0) return -9999;
  return sqrt(_sumAbarSq / ((double) _Nevents));
}

integAbarSq()

double CoherenceFactor::integAbarSq ( ) const

private

Definition at line 86 of file CoherenceFactor.cpp.

                                        {
  if(_Nevents <=0) return 0;
  return _sumAbarSq / ((double) _Nevents);
}

integASq()

double CoherenceFactor::integASq ( ) const

private

Definition at line 82 of file CoherenceFactor.cpp.

                                     {
  if(_Nevents <=0) return 0;
  return _sumASq / ((double) _Nevents);
}

integAxAbarStar()

complex< double > CoherenceFactor::integAxAbarStar ( ) const

private

Definition at line 68 of file CoherenceFactor.cpp.

                                                     {
  if(_Nevents <=0) return 0;
  return _sumAxAbarStar / ((double) _Nevents);
}

newEvent()

counted_ptr< IDalitzEvent > CoherenceFactor::newEvent ( )

private

Definition at line 53 of file CoherenceFactor.cpp.

                                                   {
  return getGenerator()->newEvent();
}

printResult()

void CoherenceFactor::printResult

(

std::ostream &

os = std::cout

)

const

Definition at line 175 of file CoherenceFactor.cpp.

                                                     {
  os << "Integrals from CoherenceFactor\n"
     << "\t int |A|^2: " << integASq() 
     << ", int |Abar|^2 " << integAbarSq();
  double rBSq = _CSAbs*_CSAbs;
  if(rBSq > 0){
    os << "\n\t int |Abar|^2/rB^2 " << integAbarSq()/rBSq;
  }
  if(integAbarSq() > 0){
    os << "\n ratio int |A|^2 / int |Abar|^2 = " << integASq()/integAbarSq()
       << ", ratio * rB^2 = " << rBSq*integASq()/integAbarSq();
  }
  os << "\n\n Coherence factor result (same result, different formats):"
     << "\n\tR exp(i xi) = " << Rval() << " +/- " << sigma()
     << "\n\tR exp(i xi) = " << abs(Rval()) 
     << " * exp( i * " << arg(Rval()) << " )"
     << "\n\t     sig(R) = " << sigmaAbs()
     << "  sig(xi) = " << sigmaAngle()
     << "\n\tR exp(i xi) = " << abs(Rval()) 
     << " * exp( i * " << arg(Rval())*180.0/pi << " deg )"
     << "\n\t     sig(R) = " << sigmaAbs()
     << "  sig(xi) = " << sigmaAngle()*180.0/pi << " deg";
  if(arg(Rval()) < 0){
    cout  << "\n\tR exp(i xi) = " << abs(Rval()) 
      << " * exp( i * " << arg(Rval())*180.0/pi + 360. << " deg )"
      << "\n\t     sig(R) = " << sigmaAbs()
      << "  sig(xi) = " << sigmaAngle()*180.0/pi << " deg";
  }
  cout << endl;
}

realAorAbarVar()

double CoherenceFactor::realAorAbarVar ( double  s, double  sq  ) const

private

Definition at line 94 of file CoherenceFactor.cpp.

                                                              {
  if(_Nevents <=0 ) return 0;
  if(s == 0) return -9999;
  double   m = s/((double) _Nevents);
  double msq = sq/((double) _Nevents); 
  double varAsq = (msq - m*m)/((double) _Nevents);
  
  return 0.5 * varAsq / fabs(s);
}

Rval()

complex< double > CoherenceFactor::Rval ( ) const

private

Definition at line 170 of file CoherenceFactor.cpp.

                                           {
  if(_Nevents <= 0 || integA() <= 0 || integAbar() <= 0) return -9999;
  return integAxAbarStar() / (integA() * integAbar());
}

setPrecision()

void CoherenceFactor::setPrecision ( double  prec )

inline

Definition at line 95 of file CoherenceFactor.h.

{_precision = prec;}

sigma()

complex< double > CoherenceFactor::sigma

(

)

const

Definition at line 144 of file CoherenceFactor.cpp.

                                            {
  return complex<double>(sqrt(fabs(var().real())), sqrt(fabs(var().imag())));
}

sigmaAbs()

double CoherenceFactor::sigmaAbs

(

)

const

Definition at line 152 of file CoherenceFactor.cpp.

                                     {
  double va = varAbs();
  if(va < 0) return -9999;
  return sqrt(va);
}

sigmaAngle()

double CoherenceFactor::sigmaAngle

(

)

const

Definition at line 165 of file CoherenceFactor.cpp.

                                       {
  double va = varAngle();
  if(va < 0) return -9999;
  return sqrt(va);
}

var()

complex< double > CoherenceFactor::var

(

)

const

Definition at line 121 of file CoherenceFactor.cpp.

                                         {
  bool dbThis=false;

  complex<double> m   = Rval();
  complex<double> som(m.real()*m.real(), m.imag()*m.imag());
  
  double den = integA()*integA() * integAbar()*integAbar();
  if(den <= 0) return complex<double>(-9999,-9999);

  double ar = cvarAxAbarStar().real() / den;
  double ai = cvarAxAbarStar().imag() / den;
  complex<double> a(ar, ai);
  complex<double> b = varA()/(integA()*integA()) * som;
  complex<double> c = varAbar()/(integAbar()*integAbar()) * som;

  if(dbThis){
    cout << " cvarAxAbarStar() << " << cvarAxAbarStar() << endl;
    cout << "var() = " << a << " + " << b << " + " << c
     << "  = " << a + b + c << endl;
  }
  return a + b + c;
}

varA()

double CoherenceFactor::varA ( ) const

private

Definition at line 103 of file CoherenceFactor.cpp.

                                 {
  return realAorAbarVar(_sumASq, _sumASqSquared);
}

varAbar()

double CoherenceFactor::varAbar ( ) const

private

Definition at line 106 of file CoherenceFactor.cpp.

                                    {
  return realAorAbarVar(_sumAbarSq, _sumAbarSqSquared);
}

varAbs()

double CoherenceFactor::varAbs

(

)

const

Definition at line 147 of file CoherenceFactor.cpp.

                                   {
  double a = var().real() * 2*fabs(Rval().real()/norm(Rval()));
  double b = var().imag() * 2*fabs(Rval().imag()/norm(Rval()));
  return a+b;
}

varAngle()

double CoherenceFactor::varAngle

(

)

const

Definition at line 157 of file CoherenceFactor.cpp.

                                     {
  double a= fabs(Rval().real());
  double b= fabs(Rval().imag());
  double den = a*a + b*b;
  if(den <= 0)return -9999;
  
  return var().real()*b/den + var().imag()*a/den;
}

Member Data Documentation

_A

FitAmpSum* CoherenceFactor::_A

private

Definition at line 29 of file CoherenceFactor.h.

_A_plus_Abar

FitAmpSum CoherenceFactor::_A_plus_Abar

private

Definition at line 33 of file CoherenceFactor.h.

_Abar

FitAmpSum * CoherenceFactor::_Abar

private

Definition at line 29 of file CoherenceFactor.h.

_CSAbs

double CoherenceFactor::_CSAbs

private

Definition at line 30 of file CoherenceFactor.h.

_CSPhase

double CoherenceFactor::_CSPhase

private

Definition at line 31 of file CoherenceFactor.h.

_eff

MINT::IReturnRealForEvent<IDalitzEvent>* CoherenceFactor::_eff

private

Definition at line 35 of file CoherenceFactor.h.

_myOwnGenAplusAbar

MINT::counted_ptr< MINT::IEventGenerator<IDalitzEvent> > CoherenceFactor::_myOwnGenAplusAbar

private

Definition at line 34 of file CoherenceFactor.h.

_Nevents

long int CoherenceFactor::_Nevents

private

Definition at line 42 of file CoherenceFactor.h.

_precision

double CoherenceFactor::_precision

private

Definition at line 36 of file CoherenceFactor.h.

_sumAbarSq

double CoherenceFactor::_sumAbarSq

private

Definition at line 39 of file CoherenceFactor.h.

_sumAbarSqSquared

double CoherenceFactor::_sumAbarSqSquared

private

Definition at line 39 of file CoherenceFactor.h.

_sumASq

double CoherenceFactor::_sumASq

private

Definition at line 38 of file CoherenceFactor.h.

_sumASqSquared

double CoherenceFactor::_sumASqSquared

private

Definition at line 38 of file CoherenceFactor.h.

_sumAxAbarStar

std::complex<double> CoherenceFactor::_sumAxAbarStar

private

Definition at line 40 of file CoherenceFactor.h.

_sumAxAbarStarSquared

std::complex<double> CoherenceFactor::_sumAxAbarStarSquared

private

Definition at line 40 of file CoherenceFactor.h.

---

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

• Mint/CoherenceFactor.h
• src/Mojito/CoherenceFactor/CoherenceFactor.cpp