#include <CoherenceFactorStoreAndEvaluate.h>
|
| CoherenceFactorStoreAndEvaluate (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 | varAngleDeg () const |
|
double | sigmaAngle () const |
|
double | sigmaAngleDeg () const |
|
void | printLong (std::ostream &os=std::cout) const |
|
void | print (std::ostream &os=std::cout, const std::string &pref="") const |
|
void | setPrecision (double prec) |
|
double | estimatedPrecision () const |
|
std::complex< double > | getCVal () |
|
double | absR () const |
|
double | phaseR () const |
|
double | phaseRdeg () const |
|
bool | addEvent (MINT::counted_ptr< IDalitzEvent > &evtPtr) |
|
const DalitzHistoSet & | getHistoA () const |
|
const DalitzHistoSet & | getHistoAbar () const |
|
const DalitzHistoSet & | getHistoBoth () const |
|
bool | saveHistos (const std::string &fname="CoherenceFactorHistos") const |
|
bool | drawHistos (const std::string &fname="CoherenceFactorHistos") const |
|
bool | drawHistosWith (const CoherenceFactorStoreAndEvaluate &other, const std::string &fname="CoherenceFactorHistosWith", double sf=1) const |
|
bool | drawHistoRatios (const CoherenceFactorStoreAndEvaluate &other, const std::string &fname="CoherenceFactorHistosWith", double sf=1) const |
|
void | scaleHistos (double sf=1) |
|
◆ CoherenceFactorStoreAndEvaluate()
Definition at line 18 of file CoherenceFactorStoreAndEvaluate.cpp.
MINT::IReturnComplexForEvent< IDalitzEvent > * A()
std::complex< double > _sumAxAbarStarSquared
DalitzHistoSet _histoBoth
MINT::IReturnRealForEvent< IDalitzEvent > * _eff
MINT::IReturnComplexForEvent< IDalitzEvent > * Abar()
std::complex< double > _sumAxAbarStar
DalitzHistoSet _histoAbar
◆ A()
◆ A_Value() [1/2]
complex< double > CoherenceFactorStoreAndEvaluate::A_Value |
( |
IDalitzEvent & |
evt | ) |
|
|
private |
◆ A_Value() [2/2]
std::complex<double> CoherenceFactorStoreAndEvaluate::A_Value |
( |
IDalitzEvent * |
evtPtr | ) |
|
|
inlineprivate |
◆ Abar()
◆ Abar_Value() [1/2]
complex< double > CoherenceFactorStoreAndEvaluate::Abar_Value |
( |
IDalitzEvent & |
evt | ) |
|
|
private |
Definition at line 234 of file CoherenceFactorStoreAndEvaluate.cpp.
237 if(dbThis) cout <<
"CoherenceFactorStoreAndEvaluate::Abar_Value" << endl;
239 if(dbThis) cout <<
"got value: " << val << endl;
virtual std::complex< double > ComplexVal(EVENT_TYPE &evt)=0
MINT::IReturnComplexForEvent< IDalitzEvent > * Abar()
◆ Abar_Value() [2/2]
std::complex<double> CoherenceFactorStoreAndEvaluate::Abar_Value |
( |
IDalitzEvent * |
evtPtr | ) |
|
|
inlineprivate |
◆ absR()
double CoherenceFactorStoreAndEvaluate::absR |
( |
| ) |
const |
◆ addEvent()
Definition at line 243 of file CoherenceFactorStoreAndEvaluate.cpp.
246 if(dbThis) cout <<
"Hello from CoherenceFactorStoreAndEvaluate::getCVal()" 252 bool printout =
false;
255 if(dbThis) cout <<
"WARNING in CoherenceFactorStoreAndEvaluate::getR()" 256 <<
" got evtPtr = " << evtPtr << endl;
261 cout <<
"=======================\n";
262 cout <<
"all A amplitudes\n";
264 cout <<
"all Abar ampltiudes\n";
280 complex<double> c_abar=
283 double d_a = norm(c_a);
284 double d_abar = norm(c_abar);
296 complex<double> c_prod = c_a * conj(c_abar);
298 complex<double> c_prodSquared(c_prod.real() * c_prod.real()
299 , c_prod.imag() * c_prod.imag()
306 cout <<
" DalitzEvent::EventCounter() " 309 cout <<
"INFO from CoherenceFactorStoreAndEvaluate::getR()" 310 <<
"\n\t _tries: " <<
_tries <<
", for " <<
_Nevents <<
" events." 318 cout <<
" --------------------------------- " << endl;
std::complex< double > Abar_Value(IDalitzEvent &evt)
virtual double getWeight() const =0
virtual void printAllAmps(std::ostream &os=std::cout) const
static long int eventCounter()
std::complex< double > _sumAxAbarStarSquared
DalitzHistoSet _histoBoth
double getEff(IDalitzEvent &evt)
void print(std::ostream &os=std::cout, const std::string &pref="") const
std::complex< double > A_Value(IDalitzEvent &evt)
virtual double getGeneratorPdfRelativeToPhaseSpace() const =0
void addEvent(const IDalitzEvent &evt, double weight=1)
std::complex< double > _sumAxAbarStar
DalitzHistoSet _histoAbar
◆ complexVar()
complex< double > CoherenceFactorStoreAndEvaluate::complexVar |
( |
std::complex< double > |
s, |
|
|
std::complex< double > |
s_sq |
|
) |
| const |
|
private |
◆ cvarAxAbarStar()
complex< double > CoherenceFactorStoreAndEvaluate::cvarAxAbarStar |
( |
| ) |
const |
|
private |
Definition at line 75 of file CoherenceFactorStoreAndEvaluate.cpp.
std::complex< double > complexVar(std::complex< double > s, std::complex< double > s_sq) const
std::complex< double > _sumAxAbarStarSquared
std::complex< double > _sumAxAbarStar
◆ drawHistoRatios()
bool CoherenceFactorStoreAndEvaluate::drawHistoRatios |
( |
const CoherenceFactorStoreAndEvaluate & |
other, |
|
|
const std::string & |
fname = "CoherenceFactorHistosWith" , |
|
|
double |
sf = 1 |
|
) |
| const |
◆ drawHistos()
bool CoherenceFactorStoreAndEvaluate::drawHistos |
( |
const std::string & |
fname = "CoherenceFactorHistos" | ) |
const |
Definition at line 341 of file CoherenceFactorStoreAndEvaluate.cpp.
DalitzHistoSet _histoBoth
bool draw(const std::string &baseName="", const std::string &drawOpt="", const std::string &format="eps") const
DalitzHistoSet _histoAbar
◆ drawHistosWith()
bool CoherenceFactorStoreAndEvaluate::drawHistosWith |
( |
const CoherenceFactorStoreAndEvaluate & |
other, |
|
|
const std::string & |
fname = "CoherenceFactorHistosWith" , |
|
|
double |
sf = 1 |
|
) |
| const |
Definition at line 349 of file CoherenceFactorStoreAndEvaluate.cpp.
bool drawWithFit(const DalitzHistoSet &fit, const std::string &baseName="", const std::string &format="eps", const std::string &fitDrawOpt="HIST C SAME") const
DalitzHistoSet _histoBoth
DalitzHistoSet _histoAbar
◆ estimatedPrecision()
double CoherenceFactorStoreAndEvaluate::estimatedPrecision |
( |
| ) |
const |
◆ getCVal()
std::complex<double> CoherenceFactorStoreAndEvaluate::getCVal |
( |
| ) |
|
◆ getEff() [1/2]
double CoherenceFactorStoreAndEvaluate::getEff |
( |
IDalitzEvent & |
evt | ) |
|
|
private |
◆ getEff() [2/2]
◆ getHistoA()
const DalitzHistoSet & CoherenceFactorStoreAndEvaluate::getHistoA |
( |
| ) |
const |
◆ getHistoAbar()
const DalitzHistoSet & CoherenceFactorStoreAndEvaluate::getHistoAbar |
( |
| ) |
const |
◆ getHistoBoth()
const DalitzHistoSet & CoherenceFactorStoreAndEvaluate::getHistoBoth |
( |
| ) |
const |
◆ integA()
double CoherenceFactorStoreAndEvaluate::integA |
( |
| ) |
const |
|
private |
◆ integAbar()
double CoherenceFactorStoreAndEvaluate::integAbar |
( |
| ) |
const |
|
private |
◆ integAbarSq()
double CoherenceFactorStoreAndEvaluate::integAbarSq |
( |
| ) |
const |
|
private |
◆ integASq()
double CoherenceFactorStoreAndEvaluate::integASq |
( |
| ) |
const |
|
private |
◆ integAxAbarStar()
complex< double > CoherenceFactorStoreAndEvaluate::integAxAbarStar |
( |
| ) |
const |
|
private |
◆ phaseR()
double CoherenceFactorStoreAndEvaluate::phaseR |
( |
| ) |
const |
◆ phaseRdeg()
double CoherenceFactorStoreAndEvaluate::phaseRdeg |
( |
| ) |
const |
◆ print()
void CoherenceFactorStoreAndEvaluate::print |
( |
std::ostream & |
os = std::cout , |
|
|
const std::string & |
pref = "" |
|
) |
| const |
Definition at line 212 of file CoherenceFactorStoreAndEvaluate.cpp.
214 os << pref <<
"\t int |A| = " << this->
integA()
215 <<
"\t int |Abar| = " << this->
integAbar()
217 <<
'\n' << pref <<
"\tR exp(i xi) = " << abs(
Rval())
218 <<
" * exp( i * " << arg(
Rval())*180.0/
pi <<
" deg )" 219 <<
"\n" << pref <<
"\t sig(R) = " <<
sigmaAbs()
std::complex< double > Rval() const
std::complex< double > integAxAbarStar() const
double sigmaAngle() const
◆ printLong()
void CoherenceFactorStoreAndEvaluate::printLong |
( |
std::ostream & |
os = std::cout | ) |
const |
Definition at line 182 of file CoherenceFactorStoreAndEvaluate.cpp.
183 os <<
"Integrals from CoherenceFactorStoreAndEvaluate\n" 188 os <<
"\n\t int |Abar|^2/rB^2 " <<
integAbarSq()/rBSq;
194 os <<
"\n\n Coherence factor result (same result, different formats):" 195 <<
"\n\tR exp(i xi) = " <<
Rval() <<
" +/- " <<
sigma()
196 <<
"\n\tR exp(i xi) = " << abs(
Rval())
197 <<
" * exp( i * " << arg(
Rval()) <<
" )" 200 <<
"\n\tR exp(i xi) = " << abs(
Rval())
201 <<
" * exp( i * " << arg(
Rval())*180.0/
pi <<
" deg )" 205 cout <<
"\n\tR exp(i xi) = " << abs(
Rval())
206 <<
" * exp( i * " << arg(
Rval())*180.0/
pi + 360. <<
" deg )" std::complex< double > sigma() const
std::complex< double > Rval() const
double sigmaAngle() const
double integAbarSq() const
◆ realAorAbarVar()
double CoherenceFactorStoreAndEvaluate::realAorAbarVar |
( |
double |
s, |
|
|
double |
sq |
|
) |
| const |
|
private |
◆ Rval()
complex< double > CoherenceFactorStoreAndEvaluate::Rval |
( |
| ) |
const |
|
private |
◆ saveHistos()
bool CoherenceFactorStoreAndEvaluate::saveHistos |
( |
const std::string & |
fname = "CoherenceFactorHistos" | ) |
const |
◆ scaleHistos()
void CoherenceFactorStoreAndEvaluate::scaleHistos |
( |
double |
sf = 1 | ) |
|
◆ setPrecision()
void CoherenceFactorStoreAndEvaluate::setPrecision |
( |
double |
prec | ) |
|
|
inline |
◆ sigma()
complex< double > CoherenceFactorStoreAndEvaluate::sigma |
( |
| ) |
const |
◆ sigmaAbs()
double CoherenceFactorStoreAndEvaluate::sigmaAbs |
( |
| ) |
const |
◆ sigmaAngle()
double CoherenceFactorStoreAndEvaluate::sigmaAngle |
( |
| ) |
const |
◆ sigmaAngleDeg()
double CoherenceFactorStoreAndEvaluate::sigmaAngleDeg |
( |
| ) |
const |
◆ var()
complex< double > CoherenceFactorStoreAndEvaluate::var |
( |
| ) |
const |
Definition at line 120 of file CoherenceFactorStoreAndEvaluate.cpp.
123 complex<double>
m =
Rval();
124 complex<double> som(
m.real()*
m.real(),
m.imag()*
m.imag());
127 if(den <= 0)
return complex<double>(-9999,-9999);
131 complex<double> a(ar, ai);
137 cout <<
"var() = " << a <<
" + " << b <<
" + " << c
138 <<
" = " << a + b + c << endl;
std::complex< double > cvarAxAbarStar() const
std::complex< double > Rval() const
◆ varA()
double CoherenceFactorStoreAndEvaluate::varA |
( |
| ) |
const |
|
private |
◆ varAbar()
double CoherenceFactorStoreAndEvaluate::varAbar |
( |
| ) |
const |
|
private |
◆ varAbarSq()
double CoherenceFactorStoreAndEvaluate::varAbarSq |
( |
| ) |
const |
|
private |
◆ varAbs()
double CoherenceFactorStoreAndEvaluate::varAbs |
( |
| ) |
const |
◆ varAngle()
double CoherenceFactorStoreAndEvaluate::varAngle |
( |
| ) |
const |
Definition at line 156 of file CoherenceFactorStoreAndEvaluate.cpp.
157 double a= fabs(
Rval().real());
158 double b= fabs(
Rval().imag());
159 double den = a*a + b*b;
160 if(den <= 0)
return -9999;
162 return var().real()*b/den +
var().imag()*a/den;
std::complex< double > Rval() const
std::complex< double > var() const
◆ varAngleDeg()
double CoherenceFactorStoreAndEvaluate::varAngleDeg |
( |
| ) |
const |
◆ varASq()
double CoherenceFactorStoreAndEvaluate::varASq |
( |
| ) |
const |
|
private |
◆ CoherenceFactorCalculator
◆ _A
FitAmpSum* CoherenceFactorStoreAndEvaluate::_A |
|
private |
◆ _A_plus_Abar
FitAmpSum CoherenceFactorStoreAndEvaluate::_A_plus_Abar |
|
private |
◆ _Abar
FitAmpSum * CoherenceFactorStoreAndEvaluate::_Abar |
|
private |
◆ _CSAbs
double CoherenceFactorStoreAndEvaluate::_CSAbs |
|
private |
◆ _CSPhase
double CoherenceFactorStoreAndEvaluate::_CSPhase |
|
private |
◆ _eff
◆ _histoA
◆ _histoAbar
◆ _histoBoth
◆ _Nevents
long int CoherenceFactorStoreAndEvaluate::_Nevents |
|
private |
◆ _precision
double CoherenceFactorStoreAndEvaluate::_precision |
|
private |
◆ _sumAbarSq
double CoherenceFactorStoreAndEvaluate::_sumAbarSq |
|
private |
◆ _sumAbarSqSquared
double CoherenceFactorStoreAndEvaluate::_sumAbarSqSquared |
|
private |
◆ _sumASq
double CoherenceFactorStoreAndEvaluate::_sumASq |
|
private |
◆ _sumASqSquared
double CoherenceFactorStoreAndEvaluate::_sumASqSquared |
|
private |
◆ _sumAxAbarStar
std::complex<double> CoherenceFactorStoreAndEvaluate::_sumAxAbarStar |
|
private |
◆ _sumAxAbarStarSquared
std::complex<double> CoherenceFactorStoreAndEvaluate::_sumAxAbarStarSquared |
|
private |
◆ _tries
long int CoherenceFactorStoreAndEvaluate::_tries |
|
private |
The documentation for this class was generated from the following files: