#include <AmpsPdfFlexiFast.h>

Inheritance diagram for AmpsPdfFlexiFast:

Public Member Functions
double	un_normalised_noPs (IDalitzEvent &evt)

std::complex< double >	ComplexVal_un_normalised_noPs (IDalitzEvent &evt)

	AmpsPdfFlexiFast (const DalitzEventPattern &pat, IFastAmplitudeIntegrable amps, MINT::MinuitParameterSet mps, double precision=1.e-4, std::string method="efficient", std::string fname="SignalIntegrationEvents.root", bool genMoreEvents=false)

IFastAmplitudeIntegrable *	getAmpSum ()

	~AmpsPdfFlexiFast ()

double	un_normalised_noPs (IDalitzEvent &evt)

std::complex< double >	ComplexVal_un_normalised_noPs (IDalitzEvent &evt)

	AmpsPdfFlexiFast (const DalitzEventPattern &pat, IFastAmplitudeIntegrable amps, MinuitParameterSet mps, double precision=1.e-4, std::string method="efficient", std::string fname="SignalIntegrationEvents.root", bool genMoreEvents=false)

IFastAmplitudeIntegrable *	getAmpSum ()

	~AmpsPdfFlexiFast ()

Public Member Functions inherited from DalitzPdfBaseFlexiFastInteg
	DalitzPdfBaseFlexiFastInteg (const DalitzEventPattern &pat, MINT::IEventGenerator< IDalitzEvent > generator, IFastAmplitudeIntegrable amps, double precision=1.e-3, MINT::MinuitParameterSet *mps=0)

	DalitzPdfBaseFlexiFastInteg (const DalitzEventPattern &pat, MINT::IEventGenerator< IDalitzEvent > generator=0, MINT::MinuitParameterSet mps=0, double precision=1.e-3)

	DalitzPdfBaseFlexiFastInteg (const DalitzPdfBaseFlexiFastInteg &other)

virtual	~DalitzPdfBaseFlexiFastInteg ()

void	setIntegrationPrecision (double prec)

void	setEventGenerator (MINT::IEventGenerator< IDalitzEvent > *g)

MINT::IEventGenerator< IDalitzEvent > *	getEventGenerator ()

const MINT::IEventGenerator< IDalitzEvent > *	getEventGenerator () const

void	parametersChanged ()

virtual double	getVal (IDalitzEvent &evt)

virtual double	getVal_withPs (IDalitzEvent &evt)

virtual double	getVal_noPs (IDalitzEvent &evt)

virtual double	getNewVal (IDalitzEvent &evt)

virtual double	getVal (IDalitzEvent *evt)

virtual double	getVal_withPs (IDalitzEvent *evt)

virtual double	getVal_noPs (IDalitzEvent *evt)

virtual double	getNewVal (IDalitzEvent *evt)

virtual double	RealVal (IDalitzEvent &evt)

MINT::MinuitParameterSet *	getMPS ()

const MINT::MinuitParameterSet *	getMPS () const

bool	saveIntegrator (const std::string &fname) const

bool	makePlots (const std::string &filename) const

virtual DalitzHistoSet	histoSet () const

virtual DalitzHistoSet	histoSet ()

void	saveEachAmpsHistograms (const std::string &prefix) const

std::vector< DalitzHistoSet >	GetEachAmpsHistograms ()

virtual DalitzHistoSet	interferenceHistoSet () const

virtual DalitzHistoSet	interferenceHistoSet ()

void	saveInterferenceHistograms (const std::string &prefix) const

std::vector< DalitzHistoSet >	GetInterferenceHistograms ()

IFastAmplitudeIntegrable *	getFitAmpSum ()

virtual void	beginFit ()

virtual void	endFit ()

double	getIntegralValue () const

double	redoIntegrator ()

double	integralForMatchingPatterns (bool match, int pattern_sign)

std::complex< double >	ComplexSumForMatchingPatterns (bool match)

std::complex< double >	ComplexIntegralForTags (int tag1, int tag2)

double	sumOfFitFractions ()

double	absSumOfInterferenceFractions ()

double	absSumOfSqrtInterferenceFractions ()

double	sumOfSqrtFitFractions ()

int	numberOfFitFractionsLargerThanThreshold (double threshold)

FitFractionList	getFractions () const

FitFractionList	getInterferenceTerms () const

bool	doFractions ()

void	doFinalStats (MINT::Minimiser *mini=0)

void	doFinalStatsAndSave (MINT::Minimiser *min=0, const std::string &fname="FitAmpResults.txt", const std::string &fnameROOT="fitFractions.root")

Public Member Functions inherited from MINT::PdfBase< IDalitzEvent >
	PdfBase ()

	PdfBase (const PdfBase< IDalitzEvent > &)

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
TRandom *	_localRnd

SignalGenerator *	_sgGen

FromFileGenerator *	_fileGen

MINT::IEventGenerator< IDalitzEvent > *	_chosenGen

MINT::NamedParameter< std::string >	_integratorSource

std::string	_integratorFileName

IEventGenerator< IDalitzEvent > *	_chosenGen

NamedParameter< std::string >	_integratorSource

Protected Attributes inherited from DalitzPdfBaseFlexiFastInteg
MINT::MinuitParameterSet *	_mps

DalitzEventPattern	_pat

double	_norm

double	_precision

FlexiFastAmplitudeIntegrator	_faint

IFastAmplitudeIntegrable *	_amps

MINT::counted_ptr< IFastAmplitudeIntegrable >	_countedAmps

MINT::IEventGenerator< IDalitzEvent > *	_generator

bool	_integrating

MINT::counted_ptr< MINT::IEventGenerator< IDalitzEvent > >	_defaultGenerator

std::string	_commaSepList_of_SavedIntegrators

Additional Inherited Members
Protected Member Functions inherited from DalitzPdfBaseFlexiFastInteg
void	setup ()

bool	makeAmps ()

virtual double	phaseSpace (IDalitzEvent &evt)

virtual bool	getNorm ()

virtual bool	integrating ()

IFastAmplitudeIntegrable *	getAmps ()

MINT::IEventGenerator< IDalitzEvent > *	makeDefaultGenerator ()

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

Detailed Description

Definition at line 18 of file AmpsPdfFlexiFast.h.

Constructor & Destructor Documentation

◆ AmpsPdfFlexiFast() [1/2]

AmpsPdfFlexiFast::AmpsPdfFlexiFast	(	const DalitzEventPattern &	pat,
		IFastAmplitudeIntegrable *	amps,
		MINT::MinuitParameterSet *	mps,
		double	precision = `1.e-4`,
		std::string	method = `"efficient"`,
		std::string	fname = `"SignalIntegrationEvents.root"`,
		bool	genMoreEvents = `false`
	)

inline

Definition at line 38 of file AmpsPdfFlexiFast.h.

     : DalitzPdfBaseFlexiFastInteg(pat, 0, amps, precision, mps)
     , _localRnd(0)
     , _sgGen(0)
     , _fileGen(0)
     , _chosenGen(0)
     , _integratorSource("IntegratorSource", (std::string) "new", (char*) 0)
     , _integratorFileName(fname)
     {
         std::cout << " AmpsPdfFlexiFast with integ method " << method << std::endl;
         bool nonFlat = "efficient" == method;
         bool generateNew = ((std::string)_integratorSource == (std::string)"new");
         if(nonFlat){
             std::cout << "AmpsPdfFlexiFast uses nonFlat integration." << std::endl;
             if(generateNew){
                 _sgGen =  new SignalGenerator(pat);
                 _sgGen->setWeighted();
                 //_sgGen->setSaveEvents(_integratorFileName);
                 _chosenGen = _sgGen;
             }else{
                 // here, SignalGenerator is used by FromFileGenerator, to fill
                 // up missing events in case more are needed than found in the
                 // file.  Since we don't know with which random seed the
                 // events in the file were generated, we supply a random
                 // number generator with randomised seed.
                 _localRnd = new TRandom3(time(0));
                 _sgGen =  new SignalGenerator(pat, _localRnd);
                 _sgGen->setWeighted();
                 _sgGen->dontSaveEvents();// saving events is done by FromFileGenerator
                 if(genMoreEvents) _fileGen   = new FromFileGenerator(_integratorFileName, _sgGen);
                 else{
                     _fileGen = new FromFileGenerator(_integratorFileName, 0, "UPDATE");
                     std::cout << "not going to generate any more events" << std::endl;
                 }
                 _chosenGen = _fileGen;
             }
             this->setEventGenerator(_chosenGen);
         }else{
             std::cout << "AmpsPdfFlexiFast uses flat integration." << std::endl;
         }
     }

◆ ~AmpsPdfFlexiFast() [1/2]

AmpsPdfFlexiFast::~AmpsPdfFlexiFast ( )

inline

Definition at line 88 of file AmpsPdfFlexiFast.h.

                        {
         if(0 != _fileGen)  delete _fileGen;
         if(0 != _sgGen)    delete _sgGen;
         if(0 != _localRnd) delete _localRnd;
     }

◆ AmpsPdfFlexiFast() [2/2]

AmpsPdfFlexiFast::AmpsPdfFlexiFast	(	const DalitzEventPattern &	pat,
		IFastAmplitudeIntegrable *	amps,
		MinuitParameterSet *	mps,
		double	precision = `1.e-4`,
		std::string	method = `"efficient"`,
		std::string	fname = `"SignalIntegrationEvents.root"`,
		bool	genMoreEvents = `false`
	)

inline

Definition at line 100 of file BsDsKpipi_TD_ampFit.cpp.

     : DalitzPdfBaseFlexiFastInteg(pat, 0, amps, precision, mps)
     , _localRnd(0)
     , _sgGen(0)
     , _fileGen(0)
     , _chosenGen(0)
     , _integratorSource("IntegratorSource", (std::string) "new", (char*) 0)
     , _integratorFileName(fname)
     {
         cout << " AmpsPdfFlexiFast with integ method " << method << endl;
         bool nonFlat = "efficient" == method;
         bool generateNew = ((string)_integratorSource == (string)"new");
         if(nonFlat){
             cout << "AmpsPdfFlexiFast uses nonFlat integration." << endl;
             if(generateNew){
                 _sgGen =  new SignalGenerator(pat);
                 _sgGen->setWeighted();
                 //_sgGen->setSaveEvents(_integratorFileName);
                 _chosenGen = _sgGen;
             }else{
                 // here, SignalGenerator is used by FromFileGenerator, to fill
                 // up missing events in case more are needed than found in the
                 // file.  Since we don't know with which random seed the
                 // events in the file were generated, we supply a random
                 // number generator with randomised seed.
                 _localRnd = new TRandom3(time(0));
                 _sgGen =  new SignalGenerator(pat, _localRnd);
                 _sgGen->setWeighted();
                 _sgGen->dontSaveEvents();// saving events is done by FromFileGenerator
                 if(genMoreEvents) _fileGen   = new FromFileGenerator(_integratorFileName, _sgGen);
                 else{
                     _fileGen = new FromFileGenerator(_integratorFileName, 0, "UPDATE");
                     cout << "not going to generate any more events" << endl;
                 }
                 _chosenGen = _fileGen;
             }
             this->setEventGenerator(_chosenGen);
         }else{
             cout << "AmpsPdfFlexiFast uses flat integration." << endl;
         }
     }

◆ ~AmpsPdfFlexiFast() [2/2]

AmpsPdfFlexiFast::~AmpsPdfFlexiFast ( )

inline

Definition at line 150 of file BsDsKpipi_TD_ampFit.cpp.

                        {
         if(0 != _fileGen)  delete _fileGen;
         if(0 != _sgGen)    delete _sgGen;
         if(0 != _localRnd) delete _localRnd;
     }

Member Function Documentation

◆ ComplexVal_un_normalised_noPs() [1/2]

std::complex<double> AmpsPdfFlexiFast::ComplexVal_un_normalised_noPs ( IDalitzEvent & evt )

inline

Definition at line 34 of file AmpsPdfFlexiFast.h.

                                                                        {
         return  _amps->ComplexVal(evt);
     }

◆ ComplexVal_un_normalised_noPs() [2/2]

std::complex<double> AmpsPdfFlexiFast::ComplexVal_un_normalised_noPs ( IDalitzEvent & evt )

inline

Definition at line 96 of file BsDsKpipi_TD_ampFit.cpp.

                                                                        {
         return  _amps->ComplexVal(evt);
     }

◆ getAmpSum() [1/2]

IFastAmplitudeIntegrable* AmpsPdfFlexiFast::getAmpSum ( )

inline

Definition at line 86 of file AmpsPdfFlexiFast.h.

86 { return _amps;}

DalitzPdfBaseFlexiFastInteg::_amps

IFastAmplitudeIntegrable * _amps

Definition: DalitzPdfBaseFlexiFastInteg.h:65

◆ getAmpSum() [2/2]

IFastAmplitudeIntegrable* AmpsPdfFlexiFast::getAmpSum ( )

inline

Definition at line 148 of file BsDsKpipi_TD_ampFit.cpp.

148 { return _amps;}

DalitzPdfBaseFlexiFastInteg::_amps

IFastAmplitudeIntegrable * _amps

Definition: DalitzPdfBaseFlexiFastInteg.h:65

◆ un_normalised_noPs() [1/2]

double AmpsPdfFlexiFast::un_normalised_noPs ( IDalitzEvent & evt )

inlinevirtual

Implements DalitzPdfBaseFlexiFastInteg.

Definition at line 29 of file AmpsPdfFlexiFast.h.

                                                 {
         double ampSq =  _amps->RealVal(evt);
         return ampSq;// * getEvent()->phaseSpace();
     }

◆ un_normalised_noPs() [2/2]

double AmpsPdfFlexiFast::un_normalised_noPs ( IDalitzEvent & evt )

inlinevirtual

Implements DalitzPdfBaseFlexiFastInteg.

Definition at line 91 of file BsDsKpipi_TD_ampFit.cpp.

                                                 {
         double ampSq =  _amps->RealVal(evt);
         return ampSq;// * getEvent()->phaseSpace();
     }

Member Data Documentation

◆ _chosenGen [1/2]

MINT::IEventGenerator<IDalitzEvent>* AmpsPdfFlexiFast::_chosenGen

protected

Definition at line 25 of file AmpsPdfFlexiFast.h.

◆ _chosenGen [2/2]

IEventGenerator<IDalitzEvent>* AmpsPdfFlexiFast::_chosenGen

protected

Definition at line 87 of file BsDsKpipi_TD_ampFit.cpp.

◆ _fileGen

FromFileGenerator * AmpsPdfFlexiFast::_fileGen

protected

Definition at line 24 of file AmpsPdfFlexiFast.h.

◆ _integratorFileName

std::string AmpsPdfFlexiFast::_integratorFileName

protected

Definition at line 27 of file AmpsPdfFlexiFast.h.

◆ _integratorSource [1/2]

MINT::NamedParameter<std::string> AmpsPdfFlexiFast::_integratorSource

protected

Definition at line 26 of file AmpsPdfFlexiFast.h.

◆ _integratorSource [2/2]

NamedParameter<std::string> AmpsPdfFlexiFast::_integratorSource

protected

Definition at line 88 of file BsDsKpipi_TD_ampFit.cpp.

◆ _localRnd

TRandom * AmpsPdfFlexiFast::_localRnd

protected

Definition at line 22 of file AmpsPdfFlexiFast.h.

◆ _sgGen

SignalGenerator * AmpsPdfFlexiFast::_sgGen

protected

Definition at line 23 of file AmpsPdfFlexiFast.h.

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

Mint/AmpsPdfFlexiFast.h
src/bin/BsDsKpipi_TD_ampFit.cpp

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ AmpsPdfFlexiFast() [1/2]

◆ ~AmpsPdfFlexiFast() [1/2]

◆ AmpsPdfFlexiFast() [2/2]

◆ ~AmpsPdfFlexiFast() [2/2]

Member Function Documentation

◆ ComplexVal_un_normalised_noPs() [1/2]

◆ ComplexVal_un_normalised_noPs() [2/2]

◆ getAmpSum() [1/2]

◆ getAmpSum() [2/2]

◆ un_normalised_noPs() [1/2]

◆ un_normalised_noPs() [2/2]

Member Data Documentation

◆ _chosenGen [1/2]

◆ _chosenGen [2/2]

◆ _fileGen

◆ _integratorFileName

◆ _integratorSource [1/2]

◆ _integratorSource [2/2]

◆ _localRnd

◆ _sgGen