CoherenceFactorCalculator Class Reference

#include <CoherenceFactorCalculator.h>

Public Member Functions
	CoherenceFactorCalculator (FitAmpSum &A, FitAmpSum &Abar, double CSAbs=1, double CSPhase=0.0, MINT::IReturnRealForEvent< IDalitzEvent > *eff=0, double prec=1.e-3, const std::string &name="coherenceFactorC")
void	printResult (std::ostream &os=std::cout) const
std::complex< double >	evaluate ()
void	setPrecision (double prec)
double	estimatedPrecision () const
double	effA () const
double	effAbar () const
double	effAVar () const
double	effAbarVar () const
double	effASigma () const
double	effAbarSigma () const

Private Member Functions
MINT::counted_ptr< MINT::IEventGenerator< IDalitzEvent > >	getGenerator ()
MINT::counted_ptr< IDalitzEvent >	newEvent ()

Private Attributes
FitAmpSum	_A_plus_Abar
MINT::counted_ptr< MINT::IEventGenerator< IDalitzEvent > >	_myOwnGenAplusAbar
double	_precision
std::string	_name
std::vector< MINT::counted_ptr< CoherenceFactorStoreAndEvaluate > >	_cfList

Detailed Description

Definition at line 27 of file CoherenceFactorCalculator.h.

Constructor & Destructor Documentation

CoherenceFactorCalculator()

CoherenceFactorCalculator::CoherenceFactorCalculator	(	FitAmpSum &	A,
		FitAmpSum &	Abar,
		double	CSAbs = 1,
		double	CSPhase = 0.0,
		MINT::IReturnRealForEvent< IDalitzEvent > *	eff = 0,
		double	prec = 1.e-3,
		const std::string &	name = "coherenceFactorC"
	)

Definition at line 21 of file CoherenceFactorCalculator.cpp.

  : _A_plus_Abar(A+Abar)
  , _myOwnGenAplusAbar(0)
  , _precision(prec)
  , _name(name)
{
   counted_ptr< CoherenceFactorStoreAndEvaluate>
     ptr_flat(new CoherenceFactorStoreAndEvaluate(A, Abar
                          , CSAbs, CSPhase
                          , 0, prec));
   _cfList.push_back(ptr_flat);

   counted_ptr< CoherenceFactorStoreAndEvaluate>
     ptr_eff(new CoherenceFactorStoreAndEvaluate(A, Abar
                         , CSAbs, CSPhase
                         , eff, prec));
   _cfList.push_back(ptr_eff);
 
}

Member Function Documentation

effA()

double CoherenceFactorCalculator::effA

(

)

const

Definition at line 66 of file CoherenceFactorCalculator.cpp.

                                            {
  if(_cfList.size() < 2) return 1;
  double d = _cfList[0]->integASq();
  if(d <=0){
    return -1;
  }else{
    return _cfList[1]->integASq()/d;
  }
}

effAbar()

double CoherenceFactorCalculator::effAbar

(

)

const

Definition at line 75 of file CoherenceFactorCalculator.cpp.

                                               {
  if(_cfList.size() < 2) return 1;
  double d = _cfList[0]->integAbarSq();
  if(d <=0){
    return -1;
  }else{
    return _cfList[1]->integAbarSq()/d;
  }
}

effAbarSigma()

double CoherenceFactorCalculator::effAbarSigma

(

)

const

Definition at line 110 of file CoherenceFactorCalculator.cpp.

                                                    {
  double v = effAbarVar();
  if(v < 0) return -1;
  return sqrt(v);
}

effAbarVar()

double CoherenceFactorCalculator::effAbarVar

(

)

const

Definition at line 95 of file CoherenceFactorCalculator.cpp.

                                                  {
  double d = _cfList[0]->integAbarSq();
  if(d <=0){
    return -1;
  }else{
    double vn = _cfList[1]->varAbarSq();
    double vd = _cfList[0]->varAbarSq();
    return vn/(d*d) + effAbar()*effAbar()*vd/(d*d);
  }
}

effASigma()

double CoherenceFactorCalculator::effASigma

(

)

const

Definition at line 105 of file CoherenceFactorCalculator.cpp.

                                                 {
  double v = effAVar();
  if(v < 0) return -1;
  return sqrt(v);
}

effAVar()

double CoherenceFactorCalculator::effAVar

(

)

const

Definition at line 85 of file CoherenceFactorCalculator.cpp.

                                               {
  double d = _cfList[0]->integASq();
  if(d <=0){
    return -1;
  }else{
    double vn = _cfList[1]->varASq();
    double vd = _cfList[0]->varASq();
    return vn/(d*d) + effA()*effA()*vd/(d*d);
  }
}

estimatedPrecision()

double CoherenceFactorCalculator::estimatedPrecision

(

)

const

Definition at line 214 of file CoherenceFactorCalculator.cpp.

                                                          {
  double prec = -9999;
  for(unsigned int i=0; i < _cfList.size(); i++){
    double tp = _cfList[i]->estimatedPrecision();
    if(tp > prec) prec = tp;
  }
  return prec;
}

evaluate()

complex< double > CoherenceFactorCalculator::evaluate

(

)

Definition at line 150 of file CoherenceFactorCalculator.cpp.

                                                   {
  bool dbThis=false;
  if(dbThis) cout << "Hello from CoherenceFactorCalculator::evaluate()" << endl;

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

      for(unsigned int i=0; i < _cfList.size(); i++){
    _cfList[i]->addEvent(evtPtr);
      }
      
      if(printout){
    cout << " DalitzEvent::EventCounter() "
         << DalitzEvent::eventCounter()
         << endl;
    cout << "INFO from CoherenceFactor::getR()"
         << "\n\t tries: " << tries
         << 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;
      }
      if(makeHistos){
    /*
    for(unsigned int i=0; i < _cfList.size(); i++){
      char fname[1000];
      sprintf(fname, "%s_%d", _name.c_str(), i);
      _cfList[i]->saveHistos(fname);
      _cfList[i]->drawHistos(fname);
    }
    */
    if(_cfList.size() >=2){
      _cfList[1]->drawHistosWith(*_cfList[0], _name + "AmpCompare", effA());
      _cfList[1]->drawHistoRatios(*_cfList[0], _name + "Eff");
    }
      }
    }
  }
  return _cfList[1]->Rval();
}

getGenerator()

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

private

Definition at line 48 of file CoherenceFactorCalculator.cpp.

                                       {
  bool dbThis=false;
  if(0 == _myOwnGenAplusAbar){
    if(dbThis) cout << "CoherenceFactorCalculator::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;
}

newEvent()

counted_ptr< IDalitzEvent > CoherenceFactorCalculator::newEvent ( )

private

Definition at line 62 of file CoherenceFactorCalculator.cpp.

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

printResult()

void CoherenceFactorCalculator::printResult

(

std::ostream &

os = std::cout

)

const

Definition at line 115 of file CoherenceFactorCalculator.cpp.

                                                               {
  double R0 = -9999;
  if(_cfList.size() >=2){
    R0 = _cfList[0]->absR();
    os << "Efficiencies from CoherenceFactorCalculater\n";
    os << " A: " << effA() << " +/- " << effASigma();
    os << ", Abar: " << effAbar() << " +/- " << effAbarSigma();
  }
  os << "\n Coherence Factors";
  for(unsigned int i=0; i < _cfList.size(); i++){
    if(0 == i){
      os << "\n * No efficiency effects "
         << " - 100% efficient detector\n";
    }else{
      os << "\n * Including efficiency effects\n";
    }
    _cfList[i]->print(os, "\t" + anythingToString((int) i) + ")");
    if(i > 0){
      os << "\n ....difference..relative..to..no..efficiency..effects:....";
      os << "\n              diff R=  " << _cfList[i]->absR() - R0
     << " +/- " << sqrt(_cfList[i]->varAbs() + _cfList[0]->varAbs());
      if(R0 > 0){
    os << "\n     relative diff R=  " 
       << ((_cfList[i]->absR() - _cfList[0]->absR())/R0)
       << " +/- " << (sqrt(_cfList[i]->varAbs() + _cfList[0]->varAbs())/R0);
      os << "\n          diff phase=  " 
         << _cfList[i]->phaseRdeg() - _cfList[0]->phaseRdeg()
         << " +/- " << sqrt(_cfList[i]->varAngleDeg() + _cfList[0]->varAngleDeg())
         << "\n";
      }
    }
  }
  os << endl;
}

setPrecision()

void CoherenceFactorCalculator::setPrecision

(

double

prec

)

Definition at line 208 of file CoherenceFactorCalculator.cpp.

                                                       {
  _precision = prec;
  for(unsigned int i=0; i < _cfList.size(); i++){
    _cfList[i]->setPrecision(prec);
  }
}

Member Data Documentation

_A_plus_Abar

FitAmpSum CoherenceFactorCalculator::_A_plus_Abar

private

Definition at line 31 of file CoherenceFactorCalculator.h.

_cfList

std::vector< MINT::counted_ptr< CoherenceFactorStoreAndEvaluate> > CoherenceFactorCalculator::_cfList

private

Definition at line 40 of file CoherenceFactorCalculator.h.

_myOwnGenAplusAbar

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

private

Definition at line 32 of file CoherenceFactorCalculator.h.

_name

std::string CoherenceFactorCalculator::_name

private

Definition at line 34 of file CoherenceFactorCalculator.h.

_precision

double CoherenceFactorCalculator::_precision

private

Definition at line 33 of file CoherenceFactorCalculator.h.

---

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

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