AmpsPdfFlexiFastCPV Class Reference

#include <AmpsPdfFlexiFastCPV.h>

Inheritance diagram for AmpsPdfFlexiFastCPV:

Public Member Functions
void	parametersChanged ()
void	beginFit ()
void	endFit ()
void	printIntegralVals ()
double	un_normalised_noPs (IDalitzEvent &evt)
virtual double	getVal (IDalitzEvent &evt)
virtual double	getVal_withPs (IDalitzEvent &evt)
virtual double	getVal_noPs (IDalitzEvent &evt)
virtual double	getVal (IDalitzEvent *evt)
virtual double	getVal_withPs (IDalitzEvent *evt)
virtual double	getVal_noPs (IDalitzEvent *evt)
virtual DalitzHistoSet	histoSet ()
void	doFinalStatsAndSaveForAmp12 (MINT::Minimiser *min=0, const std::string &fname="FitAmpResults", const std::string &fnameROOT="fitFractions")
	AmpsPdfFlexiFastCPV (AmpsPdfFlexiFast *amps1, AmpsPdfFlexiFast *amps2, AmpsPdfFlexiFast *ampsSum, MINT::FitParameter &r, MINT::FitParameter &delta, MINT::FitParameter &gamma, MINT::FitParameter &tau, MINT::FitParameter &dGamma, MINT::FitParameter &dm, MINT::FitParameter &eff_tag, MINT::FitParameter &w)
void	parametersChanged ()
void	beginFit ()
void	endFit ()
void	printIntegralVals ()
double	un_normalised_noPs (IDalitzEvent &evt)
virtual double	getVal (IDalitzEvent &evt)
virtual double	getVal_withPs (IDalitzEvent &evt)
virtual double	getVal_noPs (IDalitzEvent &evt)
virtual double	getVal (IDalitzEvent *evt)
virtual double	getVal_withPs (IDalitzEvent *evt)
virtual double	getVal_noPs (IDalitzEvent *evt)
virtual DalitzHistoSet	histoSet ()
void	doFinalStatsAndSaveForAmp12 (MINT::Minimiser *min=0, const std::string &fname="FitAmpResults", const std::string &fnameROOT="fitFractions")
	AmpsPdfFlexiFastCPV (AmpsPdfFlexiFast *amps1, AmpsPdfFlexiFast *amps2, AmpsPdfFlexiFast *ampsSum, MINT::FitParameter &r, MINT::FitParameter &delta, MINT::FitParameter &gamma, MINT::FitParameter &tau, MINT::FitParameter &dGamma, MINT::FitParameter &dm, MINT::FitParameter &eff_tag, MINT::FitParameter &w)
Public Member Functions inherited from MINT::PdfBase< IDalitzEvent >
	PdfBase ()
	PdfBase (const PdfBase< IDalitzEvent > &)
virtual double	RealVal (IDalitzEvent &evt)
virtual double	getNewVal (IDalitzEvent &evt)
virtual	~PdfBase ()
Public Member Functions inherited from MINT::IPdf< IDalitzEvent >
virtual void	Gradient (IDalitzEvent &evt, std::vector< double > &grad, MINT::MinuitParameterSet *mps)
virtual bool	useAnalyticGradient ()
Public Member Functions inherited from MINT::IReturnRealForEvent< IDalitzEvent >
virtual	~IReturnRealForEvent ()

Protected Attributes
AmpsPdfFlexiFast *	_amps1
AmpsPdfFlexiFast *	_amps2
AmpsPdfFlexiFast *	_ampsSum
MINT::FitParameter &	_r
MINT::FitParameter &	_delta
MINT::FitParameter &	_gamma
MINT::FitParameter &	_tau
MINT::FitParameter &	_dGamma
MINT::FitParameter &	_dm
MINT::FitParameter &	_eff_tag
MINT::FitParameter &	_w
double	_intA
double	_intAbar
std::complex< double >	_intAAbar
FitParameter &	_r
FitParameter &	_delta
FitParameter &	_gamma
FitParameter &	_tau
FitParameter &	_dGamma
FitParameter &	_dm
FitParameter &	_eff_tag
FitParameter &	_w
complex< double >	_intAAbar

Additional Inherited Members
Protected Member Functions inherited from MINT::IReturnRealForEvent< IDalitzEvent >
	IReturnRealForEvent ()

Detailed Description

Definition at line 13 of file AmpsPdfFlexiFastCPV.h.

Constructor & Destructor Documentation

AmpsPdfFlexiFastCPV() [1/2]

AmpsPdfFlexiFastCPV::AmpsPdfFlexiFastCPV ( AmpsPdfFlexiFast *  amps1, AmpsPdfFlexiFast *  amps2, AmpsPdfFlexiFast *  ampsSum, MINT::FitParameter &  r, MINT::FitParameter &  delta, MINT::FitParameter &  gamma, MINT::FitParameter &  tau, MINT::FitParameter &  dGamma, MINT::FitParameter &  dm, MINT::FitParameter &  eff_tag, MINT::FitParameter &  w  )

inline

Definition at line 128 of file AmpsPdfFlexiFastCPV.h.

                                                            :
    _amps1(amps1),
      _amps2(amps2),
      _ampsSum(ampsSum),
      _r(r),
      _delta(delta),
      _gamma(gamma),
      _tau(tau),
      _dGamma(dGamma),
      _dm(dm),
      _eff_tag(eff_tag),
      _w(w),
      _intA(-1),
      _intAbar(-1),
      _intAAbar(-1)
    {}

AmpsPdfFlexiFastCPV() [2/2]

AmpsPdfFlexiFastCPV::AmpsPdfFlexiFastCPV ( AmpsPdfFlexiFast *  amps1, AmpsPdfFlexiFast *  amps2, AmpsPdfFlexiFast *  ampsSum, MINT::FitParameter &  r, MINT::FitParameter &  delta, MINT::FitParameter &  gamma, MINT::FitParameter &  tau, MINT::FitParameter &  dGamma, MINT::FitParameter &  dm, MINT::FitParameter &  eff_tag, MINT::FitParameter &  w  )

inline

Definition at line 985 of file BsDsKpipi_TD_ampFit.cpp.

                                                                                                                                               :
    _amps1(amps1),_amps2(amps2),_ampsSum(ampsSum),_r(r),_delta(delta),_gamma(gamma),_tau(tau),_dGamma(dGamma),_dm(dm),_eff_tag(eff_tag),_w(w),
    _intA(-1),_intAbar(-1),_intAAbar(-1)
    {;}

Member Function Documentation

beginFit() [1/2]

void AmpsPdfFlexiFastCPV::beginFit ( )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 42 of file AmpsPdfFlexiFastCPV.h.

                   {
        _ampsSum->beginFit();
        printIntegralVals();
    }

beginFit() [2/2]

void AmpsPdfFlexiFastCPV::beginFit ( )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 900 of file BsDsKpipi_TD_ampFit.cpp.

                   {
        _ampsSum->beginFit();
        printIntegralVals();
    }

doFinalStatsAndSaveForAmp12() [1/2]

void AmpsPdfFlexiFastCPV::doFinalStatsAndSaveForAmp12 ( MINT::Minimiser *  min = 0, const std::string &  fname = "FitAmpResults", const std::string &  fnameROOT = "fitFractions"  )

inline

Definition at line 120 of file AmpsPdfFlexiFastCPV.h.

                                                               {
        _amps1->redoIntegrator();
        _amps2->redoIntegrator();
        _amps1->doFinalStatsAndSave(min,((std::string)fname+".txt").c_str(),((std::string)fnameROOT+".root").c_str());
        _amps2->doFinalStatsAndSave(min,((std::string)fname+"_CC.txt").c_str(),((std::string)fnameROOT+"_CC.root").c_str());        
    }

doFinalStatsAndSaveForAmp12() [2/2]

void AmpsPdfFlexiFastCPV::doFinalStatsAndSaveForAmp12 ( MINT::Minimiser *  min = 0, const std::string &  fname = "FitAmpResults", const std::string &  fnameROOT = "fitFractions"  )

inline

Definition at line 978 of file BsDsKpipi_TD_ampFit.cpp.

                                                                                                                                              {
        _amps1->redoIntegrator();
        _amps2->redoIntegrator();
        _amps1->doFinalStatsAndSave(min,((string)fname+".txt").c_str(),((string)fnameROOT+".root").c_str());
        _amps2->doFinalStatsAndSave(min,((string)fname+"_CC.txt").c_str(),((string)fnameROOT+"_CC.root").c_str());        
    }

endFit() [1/2]

void AmpsPdfFlexiFastCPV::endFit ( )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 46 of file AmpsPdfFlexiFastCPV.h.

                 {
        printIntegralVals();
        _ampsSum->endFit();
    }

endFit() [2/2]

void AmpsPdfFlexiFastCPV::endFit ( )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 904 of file BsDsKpipi_TD_ampFit.cpp.

                 {
        printIntegralVals();
        _ampsSum->endFit();
    }

getVal() [1/4]

virtual double AmpsPdfFlexiFastCPV::getVal ( IDalitzEvent &  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 84 of file AmpsPdfFlexiFastCPV.h.

                                            {
        
        const double f = static_cast<double>((int)evt.getValueFromVector(4));
        const double val = un_normalised_noPs(evt);
        
        double r = (double)_r; // * sqrt(_intA/_intAbar);
        const double Gamma = 1./((double) _tau);
        const std::complex<double> phase_diff = std::polar((double)r,((double) _delta -(double)_gamma*f)/360.*2*pi);
        const double int_interference = (phase_diff*_intAAbar).real();
        
        if(_intA == -1 ){
            std::cout << "AmpsPdfFlexiFastCPV:: _norm = -1, should not have happened." << std::endl;
            throw "can't deal with that";
        }
        
        return val/(2.* ((_intA + r* r * _intAbar) * 4.*Gamma - int_interference * 2. * _dGamma )/ (4.*Gamma*Gamma-_dGamma*_dGamma));
    }

getVal() [2/4]

virtual double AmpsPdfFlexiFastCPV::getVal ( IDalitzEvent *  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 105 of file AmpsPdfFlexiFastCPV.h.

                                            {
        if(0 == evt) return 0;
        return getVal(*evt);
    }

getVal() [3/4]

virtual double AmpsPdfFlexiFastCPV::getVal ( IDalitzEvent &  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 942 of file BsDsKpipi_TD_ampFit.cpp.

                                            {
        
        const double f = static_cast<double>((int)evt.getValueFromVector(4));
        const double val = un_normalised_noPs(evt);
        
        double r = (double)_r; // * sqrt(_intA/_intAbar);
        const double Gamma = 1./((double) _tau);
        const complex<double> phase_diff = polar((double)r,((double) _delta -(double)_gamma*f)/360.*2*pi);
        const double int_interference = (phase_diff*_intAAbar).real();
        
        if(_intA == -1 ){
            cout << "AmpsPdfFlexiFastCPV:: _norm = -1, should not have happened." << endl;
            throw "can't deal with that";
        }
        
        return val/(2.* ((_intA + r* r * _intAbar) * 4.*Gamma - int_interference * 2. * _dGamma )/ (4.*Gamma*Gamma-_dGamma*_dGamma));
    }

getVal() [4/4]

virtual double AmpsPdfFlexiFastCPV::getVal ( IDalitzEvent *  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 963 of file BsDsKpipi_TD_ampFit.cpp.

                                            {
        if(0 == evt) return 0;
        return getVal(*evt);
    }

getVal_noPs() [1/4]

virtual double AmpsPdfFlexiFastCPV::getVal_noPs ( IDalitzEvent &  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 103 of file AmpsPdfFlexiFastCPV.h.

{return getVal(evt);}

getVal_noPs() [2/4]

virtual double AmpsPdfFlexiFastCPV::getVal_noPs ( IDalitzEvent *  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 113 of file AmpsPdfFlexiFastCPV.h.

                                                 {
        if(0 == evt) return 0;
        return getVal_noPs(*evt);
    }

getVal_noPs() [3/4]

virtual double AmpsPdfFlexiFastCPV::getVal_noPs ( IDalitzEvent &  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 961 of file BsDsKpipi_TD_ampFit.cpp.

{return getVal(evt);}

getVal_noPs() [4/4]

virtual double AmpsPdfFlexiFastCPV::getVal_noPs ( IDalitzEvent *  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 971 of file BsDsKpipi_TD_ampFit.cpp.

                                                 {
        if(0 == evt) return 0;
        return getVal_noPs(*evt);
    }

getVal_withPs() [1/4]

virtual double AmpsPdfFlexiFastCPV::getVal_withPs ( IDalitzEvent &  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 102 of file AmpsPdfFlexiFastCPV.h.

{return getVal(evt);}

getVal_withPs() [2/4]

virtual double AmpsPdfFlexiFastCPV::getVal_withPs ( IDalitzEvent *  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 109 of file AmpsPdfFlexiFastCPV.h.

                                                   {
        if(0 == evt) return 0;
        return getVal_withPs(*evt);
    }

getVal_withPs() [3/4]

virtual double AmpsPdfFlexiFastCPV::getVal_withPs ( IDalitzEvent &  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 960 of file BsDsKpipi_TD_ampFit.cpp.

{return getVal(evt);}

getVal_withPs() [4/4]

virtual double AmpsPdfFlexiFastCPV::getVal_withPs ( IDalitzEvent *  evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 967 of file BsDsKpipi_TD_ampFit.cpp.

                                                   {
        if(0 == evt) return 0;
        return getVal_withPs(*evt);
    }

histoSet() [1/2]

virtual DalitzHistoSet AmpsPdfFlexiFastCPV::histoSet ( )

inlinevirtual

Implements IDalitzPdf.

Definition at line 118 of file AmpsPdfFlexiFastCPV.h.

{return _ampsSum->histoSet();}

histoSet() [2/2]

virtual DalitzHistoSet AmpsPdfFlexiFastCPV::histoSet ( )

inlinevirtual

Implements IDalitzPdf.

Definition at line 976 of file BsDsKpipi_TD_ampFit.cpp.

{return _ampsSum->histoSet();}

parametersChanged() [1/2]

void AmpsPdfFlexiFastCPV::parametersChanged ( )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 36 of file AmpsPdfFlexiFastCPV.h.

                            {
        _ampsSum->parametersChanged();
        _intA = _ampsSum->integralForMatchingPatterns(true,1);
        _intAbar = _ampsSum->integralForMatchingPatterns(true,-1);        
        _intAAbar = _ampsSum->ComplexSumForMatchingPatterns(false);
    }

parametersChanged() [2/2]

void AmpsPdfFlexiFastCPV::parametersChanged ( )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 894 of file BsDsKpipi_TD_ampFit.cpp.

                            {
        _ampsSum->parametersChanged();
        _intA = _ampsSum->integralForMatchingPatterns(true,1);
        _intAbar = _ampsSum->integralForMatchingPatterns(true,-1);        
        _intAAbar = _ampsSum->ComplexSumForMatchingPatterns(false);
    }

printIntegralVals() [1/2]

void AmpsPdfFlexiFastCPV::printIntegralVals ( )

inline

Definition at line 51 of file AmpsPdfFlexiFastCPV.h.

                            {
        std::cout << "intSum = " << _ampsSum->getIntegralValue() << std::endl;
        std::cout << "intA = " << _intA << std::endl;
        std::cout << "intAbar = " << _intAbar << std::endl;
        std::cout << "intAAbar = " << _intAAbar.real() << std::endl;
    }

printIntegralVals() [2/2]

void AmpsPdfFlexiFastCPV::printIntegralVals ( )

inline

Definition at line 909 of file BsDsKpipi_TD_ampFit.cpp.

                            {
        cout << "intSum = " << _ampsSum->getIntegralValue() << endl;
        cout << "intA = " << _intA << endl;
        cout << "intAbar = " << _intAbar << endl;
        cout << "intAAbar = " << _intAAbar.real() << endl;
    }

un_normalised_noPs() [1/2]

double AmpsPdfFlexiFastCPV::un_normalised_noPs ( IDalitzEvent &  evt )

inline

Definition at line 58 of file AmpsPdfFlexiFastCPV.h.

                                                       {
        const double t = (double) evt.getValueFromVector(0);
        const double dt = (double) evt.getValueFromVector(1);
        const double q = static_cast<double>((int)evt.getValueFromVector(2)); 
        const double w = (double) evt.getValueFromVector(3);
        const double f = static_cast<double>((int)evt.getValueFromVector(4));
        
        const double e_eff = fabs(q)*_eff_tag + (1.-fabs(q))*(1.-_eff_tag);
        
        double r = (double)_r; // * sqrt(_intA/_intAbar);
        const std::complex<double> phase_diff = std::polar((double)r,((double) _delta -(double)_gamma*f)/360.*2*pi);
        
        const std::complex<double> amp = _amps1->ComplexVal_un_normalised_noPs(evt) ;
        const std::complex<double> amp_bar = _amps2->ComplexVal_un_normalised_noPs(evt) * phase_diff;
        
        const double val =  exp(-fabs(t)/(double)_tau) *
        (
         (2.-fabs(q))*(norm(amp) + norm(amp_bar))*cosh((double)_dGamma/2.*t)
         +f*q*(1.-2.*w)*(norm(amp) - norm(amp_bar)) *cos((double)_dm*t)
         -(2.-fabs(q))*2.0*real(amp_bar*conj(amp))*sinh((double)_dGamma/2.*t)
         -f*2.0*q*(1.-2.*w)*imag(amp_bar*conj(amp))*sin((double)_dm*t)
         )*e_eff;
        
        return val;
    }

un_normalised_noPs() [2/2]

double AmpsPdfFlexiFastCPV::un_normalised_noPs ( IDalitzEvent &  evt )

inline

Definition at line 916 of file BsDsKpipi_TD_ampFit.cpp.

                                                       {
        const double t = (double) evt.getValueFromVector(0);
        const double dt = (double) evt.getValueFromVector(1);
        const double q = static_cast<double>((int)evt.getValueFromVector(2)); 
        const double w = (double) evt.getValueFromVector(3);
        const double f = static_cast<double>((int)evt.getValueFromVector(4));
        
        const double e_eff = fabs(q)*_eff_tag + (1.-fabs(q))*(1.-_eff_tag);
        
        double r = (double)_r; // * sqrt(_intA/_intAbar);
        const complex<double> phase_diff = polar((double)r,((double) _delta -(double)_gamma*f)/360.*2*pi);
        
        const std::complex<double> amp = _amps1->ComplexVal_un_normalised_noPs(evt) ;
        const std::complex<double> amp_bar = _amps2->ComplexVal_un_normalised_noPs(evt) * phase_diff;
        
        const double val =  exp(-fabs(t)/(double)_tau) *
        (
         (2.-fabs(q))*(norm(amp) + norm(amp_bar))*cosh((double)_dGamma/2.*t)
         +f*q*(1.-2.*w)*(norm(amp) - norm(amp_bar)) *cos((double)_dm*t)
         -(2.-fabs(q))*2.0*real(amp_bar*conj(amp))*sinh((double)_dGamma/2.*t)
         -f*2.0*q*(1.-2.*w)*imag(amp_bar*conj(amp))*sin((double)_dm*t)
         )*e_eff;
        
        return val;
    }

Member Data Documentation

_amps1

AmpsPdfFlexiFast * AmpsPdfFlexiFastCPV::_amps1

protected

Definition at line 17 of file AmpsPdfFlexiFastCPV.h.

_amps2

AmpsPdfFlexiFast * AmpsPdfFlexiFastCPV::_amps2

protected

Definition at line 18 of file AmpsPdfFlexiFastCPV.h.

_ampsSum

AmpsPdfFlexiFast * AmpsPdfFlexiFastCPV::_ampsSum

protected

Definition at line 19 of file AmpsPdfFlexiFastCPV.h.

_delta [1/2]

MINT::FitParameter& AmpsPdfFlexiFastCPV::_delta

protected

Definition at line 22 of file AmpsPdfFlexiFastCPV.h.

_delta [2/2]

FitParameter& AmpsPdfFlexiFastCPV::_delta

protected

Definition at line 880 of file BsDsKpipi_TD_ampFit.cpp.

_dGamma [1/2]

MINT::FitParameter& AmpsPdfFlexiFastCPV::_dGamma

protected

Definition at line 26 of file AmpsPdfFlexiFastCPV.h.

_dGamma [2/2]

FitParameter& AmpsPdfFlexiFastCPV::_dGamma

protected

Definition at line 884 of file BsDsKpipi_TD_ampFit.cpp.

_dm [1/2]

MINT::FitParameter& AmpsPdfFlexiFastCPV::_dm

protected

Definition at line 27 of file AmpsPdfFlexiFastCPV.h.

_dm [2/2]

FitParameter& AmpsPdfFlexiFastCPV::_dm

protected

Definition at line 885 of file BsDsKpipi_TD_ampFit.cpp.

_eff_tag [1/2]

MINT::FitParameter& AmpsPdfFlexiFastCPV::_eff_tag

protected

Definition at line 28 of file AmpsPdfFlexiFastCPV.h.

_eff_tag [2/2]

FitParameter& AmpsPdfFlexiFastCPV::_eff_tag

protected

Definition at line 886 of file BsDsKpipi_TD_ampFit.cpp.

_gamma [1/2]

MINT::FitParameter& AmpsPdfFlexiFastCPV::_gamma

protected

Definition at line 23 of file AmpsPdfFlexiFastCPV.h.

_gamma [2/2]

FitParameter& AmpsPdfFlexiFastCPV::_gamma

protected

Definition at line 881 of file BsDsKpipi_TD_ampFit.cpp.

_intA

double AmpsPdfFlexiFastCPV::_intA

protected

Definition at line 31 of file AmpsPdfFlexiFastCPV.h.

_intAAbar [1/2]

std::complex<double> AmpsPdfFlexiFastCPV::_intAAbar

protected

Definition at line 33 of file AmpsPdfFlexiFastCPV.h.

_intAAbar [2/2]

complex<double> AmpsPdfFlexiFastCPV::_intAAbar

protected

Definition at line 891 of file BsDsKpipi_TD_ampFit.cpp.

_intAbar

double AmpsPdfFlexiFastCPV::_intAbar

protected

Definition at line 32 of file AmpsPdfFlexiFastCPV.h.

_r [1/2]

MINT::FitParameter& AmpsPdfFlexiFastCPV::_r

protected

Definition at line 21 of file AmpsPdfFlexiFastCPV.h.

_r [2/2]

FitParameter& AmpsPdfFlexiFastCPV::_r

protected

Definition at line 879 of file BsDsKpipi_TD_ampFit.cpp.

_tau [1/2]

MINT::FitParameter& AmpsPdfFlexiFastCPV::_tau

protected

Definition at line 25 of file AmpsPdfFlexiFastCPV.h.

_tau [2/2]

FitParameter& AmpsPdfFlexiFastCPV::_tau

protected

Definition at line 883 of file BsDsKpipi_TD_ampFit.cpp.

_w [1/2]

MINT::FitParameter& AmpsPdfFlexiFastCPV::_w

protected

Definition at line 29 of file AmpsPdfFlexiFastCPV.h.

_w [2/2]

FitParameter& AmpsPdfFlexiFastCPV::_w

protected

Definition at line 887 of file BsDsKpipi_TD_ampFit.cpp.

---

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

• Mint/AmpsPdfFlexiFastCPV.h
• src/bin/BsDsKpipi_TD_ampFit.cpp