#include <DalitzPdfBaseFastInteg.h>

Inheritance diagram for DalitzPdfBaseFastInteg:

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

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

	DalitzPdfBaseFastInteg (const DalitzPdfBaseFastInteg &other)

virtual	~DalitzPdfBaseFastInteg ()

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 void	Gradient (IDalitzEvent &evt, std::vector< double > &grad, MINT::MinuitParameterSet *mps)

virtual void	GradientForLasso (std::vector< double > &grad)

bool	useAnalyticGradient ()

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	endFit ()

void	setIntegratorFileName (const std::string &commaSeparatedList)

double	getIntegralValue () const

double	getIntegralValue ()

double	redoIntegrator ()

double	sumOfFitFractions ()

double	absSumOfInterferenceFractions ()

double	absSumOfSqrtInterferenceFractions ()

double	sumOfSqrtFitFractions ()

int	numberOfFitFractionsLargerThanThreshold (double threshold)

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

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

	PdfBase (const PdfBase< IDalitzEvent > &)

virtual void	beginFit ()

virtual	~PdfBase ()

Public Member Functions inherited from MINT::IReturnRealForEvent< IDalitzEvent >
virtual	~IReturnRealForEvent ()

Protected Member Functions
void	setup ()

bool	makeAmps ()

virtual double	un_normalised_noPs (IDalitzEvent &evt)=0

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 ()

Protected Attributes
MINT::MinuitParameterSet *	_mps

DalitzEventPattern	_pat

double	_norm

double	_precision

FastAmplitudeIntegrator	_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

double	_val

std::vector< double >	_gradNorm

bool	_redoGradNorm

Detailed Description

Definition at line 52 of file DalitzPdfBaseFastInteg.h.

Constructor & Destructor Documentation

◆ DalitzPdfBaseFastInteg() [1/3]

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

Definition at line 124 of file DalitzPdfBaseFastInteg.cpp.

   : PdfBase<IDalitzEvent>()
   , _mps(mps)
   , _pat(pat)
   , _norm(-1)
   , _precision(prec)
   , _amps(amps)
   , _countedAmps(0)
     //, _efficiency(0)
   , _generator(generator)
   , _integrating(0)
   , _defaultGenerator(0)
   , _commaSepList_of_SavedIntegrators("")
   , _val(-1)
   , _gradNorm(1)
   , _redoGradNorm(true)
 {
   setup();
   _gradNorm.resize(this->getMPS()->size());
   //_gradNorm= new Double_t[this->getMPS()->size()];  
   for(unsigned int i=0; i < _gradNorm.size(); i++) _gradNorm[i]= 0.;
 }

◆ DalitzPdfBaseFastInteg() [2/3]

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

Definition at line 154 of file DalitzPdfBaseFastInteg.cpp.

   : PdfBase<IDalitzEvent>()
   , _mps(mps)
   , _pat(pat)
   , _norm(-1)
   , _precision(prec)
   , _amps(0)
   , _countedAmps(0)
     //  , _efficiency(eff)
   , _generator(generator)
   , _integrating(0)
   , _defaultGenerator(0)
   , _commaSepList_of_SavedIntegrators("")
   , _val(-1)
   , _gradNorm(1)
   , _redoGradNorm(true)
 {
   setup();
   _gradNorm.resize(this->getMPS()->size());//= new Double_t[this->getMPS()->size()];  
   for(unsigned int i=0; i < _gradNorm.size(); i++) _gradNorm[i]= 0.;
   //  cout << "pset pointer in DalitzPdfBaseFastInteg " << _mps << " = " << getMPS() << endl;
 }

◆ DalitzPdfBaseFastInteg() [3/3]

DalitzPdfBaseFastInteg::DalitzPdfBaseFastInteg ( const DalitzPdfBaseFastInteg & other )

Definition at line 183 of file DalitzPdfBaseFastInteg.cpp.

   : IReturnRealForEvent<IDalitzEvent>()
   , IPdf<IDalitzEvent>()
   , IDalitzPdf()
   , PdfBase<IDalitzEvent>(other)
   , _mps(other._mps)
   , _pat(other._pat)
   , _norm(other._norm)
   , _precision(other._precision)
   , _amps(other._amps)
   , _countedAmps(other._countedAmps)
     //, _efficiency(other._efficiency)
   , _generator(other._generator)
   , _integrating(other._integrating)
   , _defaultGenerator(other._defaultGenerator)
   , _commaSepList_of_SavedIntegrators(other._commaSepList_of_SavedIntegrators)
   , _val(-1)
   , _gradNorm(1)
   , _redoGradNorm(true)
 {
   setup();
   _gradNorm.resize(this->getMPS()->size());//= new Double_t[this->getMPS()->size()];  
   for(unsigned int i=0; i < _gradNorm.size(); i++) _gradNorm[i]= 0.;
 }

◆ ~DalitzPdfBaseFastInteg()

DalitzPdfBaseFastInteg::~DalitzPdfBaseFastInteg ( )

virtual

Definition at line 208 of file DalitzPdfBaseFastInteg.cpp.

                                                {
   bool dbThis=false;
   if(dbThis)cout << "DalitzPdfBaseFastInteg::~DalitzPdfBaseFastInteg() called" << endl;
   //delete[] _gradNorm;
   if(dbThis)cout << "DalitzPdfBaseFastInteg::~DalitzPdfBaseFastInteg() done" << endl;
 }

Member Function Documentation

◆ absSumOfInterferenceFractions()

double DalitzPdfBaseFastInteg::absSumOfInterferenceFractions ( )

inline

Definition at line 197 of file DalitzPdfBaseFastInteg.h.

                                          {
         if(!_faint.initialised())getNorm();
         return _faint.absSumOfInterferenceFractions();
   }  

◆ absSumOfSqrtInterferenceFractions()

double DalitzPdfBaseFastInteg::absSumOfSqrtInterferenceFractions ( )

inline

Definition at line 202 of file DalitzPdfBaseFastInteg.h.

                                              {
         if(!_faint.initialised())getNorm();
         return _faint.absSumOfSqrtInterferenceFractions();
   } 

◆ doFinalStats()

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

Definition at line 371 of file DalitzPdfBaseFastInteg.cpp.

                                                         {
   _faint.doFinalStats(mini);
 }

◆ endFit()

void DalitzPdfBaseFastInteg::endFit ( )

virtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 366 of file DalitzPdfBaseFastInteg.cpp.

                                    {
 
   //_faint.doFinalStats();
 }

◆ getAmps()

IFastAmplitudeIntegrable* DalitzPdfBaseFastInteg::getAmps ( )

inlineprotected

Definition at line 91 of file DalitzPdfBaseFastInteg.h.

91 { return _amps; }

DalitzPdfBaseFastInteg::_amps

IFastAmplitudeIntegrable * _amps

Definition: DalitzPdfBaseFastInteg.h:65

◆ GetEachAmpsHistograms()

std::vector< DalitzHistoSet > DalitzPdfBaseFastInteg::GetEachAmpsHistograms ( )

Definition at line 347 of file DalitzPdfBaseFastInteg.cpp.

                                                                        {
   return _faint.GetEachAmpsHistograms();
 }

◆ getEventGenerator() [1/2]

IEventGenerator< IDalitzEvent > * DalitzPdfBaseFastInteg::getEventGenerator ( )

Definition at line 98 of file DalitzPdfBaseFastInteg.cpp.

                                                                         {
   if(0 == _generator) makeDefaultGenerator();
   return _generator;
 }

◆ getEventGenerator() [2/2]

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

Definition at line 102 of file DalitzPdfBaseFastInteg.cpp.

                                                                                    {
   if(0 == _generator){
     cout << "WARNING in DalitzPdfBaseFastInteg::getEventGenerator() const"
      << " returning 0 pointer." << endl;
   }
   return _generator;
 }

◆ getFitAmpSum()

IFastAmplitudeIntegrable* DalitzPdfBaseFastInteg::getFitAmpSum ( )

inline

Definition at line 174 of file DalitzPdfBaseFastInteg.h.

174 { return _amps;}

DalitzPdfBaseFastInteg::_amps

IFastAmplitudeIntegrable * _amps

Definition: DalitzPdfBaseFastInteg.h:65

◆ getIntegralValue() [1/2]

double DalitzPdfBaseFastInteg::getIntegralValue ( ) const

inline

Definition at line 180 of file DalitzPdfBaseFastInteg.h.

                                  {
     return _norm;
   }

◆ getIntegralValue() [2/2]

double DalitzPdfBaseFastInteg::getIntegralValue ( )

inline

Definition at line 183 of file DalitzPdfBaseFastInteg.h.

                            {
     if(_norm < 0) getNorm(); return _norm;
   }

◆ GetInterferenceHistograms()

std::vector< DalitzHistoSet > DalitzPdfBaseFastInteg::GetInterferenceHistograms ( )

Definition at line 362 of file DalitzPdfBaseFastInteg.cpp.

                                                                            {
     return _faint.GetInterferenceHistograms();
 }

◆ getMPS() [1/2]

MinuitParameterSet * DalitzPdfBaseFastInteg::getMPS ( )

Definition at line 13 of file DalitzPdfBaseFastInteg.cpp.

                                                   {
   return (0 == _mps ? MinuitParameterSet::getDefaultSet() : _mps);
 }

◆ getMPS() [2/2]

const MinuitParameterSet * DalitzPdfBaseFastInteg::getMPS ( ) const

Definition at line 16 of file DalitzPdfBaseFastInteg.cpp.

                                                               {
   return (0 == _mps ? MinuitParameterSet::getDefaultSet() : _mps);
 }

◆ getNewVal() [1/2]

virtual double DalitzPdfBaseFastInteg::getNewVal ( IDalitzEvent & evt )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 125 of file DalitzPdfBaseFastInteg.h.

                                              {
     this->parametersChanged();
     return getVal(evt);
   };

◆ getNewVal() [2/2]

virtual double DalitzPdfBaseFastInteg::getNewVal ( IDalitzEvent * evt )

inlinevirtual

Definition at line 143 of file DalitzPdfBaseFastInteg.h.

                                              {
     this->parametersChanged();
     if(0 == evt) return 0;
     return getNewVal(*evt);
   }

◆ getNorm()

bool DalitzPdfBaseFastInteg::getNorm ( )

protectedvirtual

Definition at line 24 of file DalitzPdfBaseFastInteg.cpp.

                                     {
   _integrating = true;
   if( ! _faint.initialised()){
     if(_pat.empty()){
       cout << "DalitzPdfBaseFastInteg::getNorm: don't know event pattern!"
        << endl;
       _integrating = false;
       return false;
     }
     if(0 == getAmps()){
       cout << "DalitzPdfBaseFastInteg::getNorm: don't have amps!"
        << endl;
       _integrating = false;
       return false;
     }
 
     /*
     if("" != _commaSepList_of_SavedIntegrators){
       _faint.initialiseFromFile(evtPtr->eventPattern()
                 , _amps
                 , _commaSepList_of_SavedIntegrators);
     }else{
     */
     _faint.initialise(_commaSepList_of_SavedIntegrators
               , _pat
               , getAmps()
               //, _efficiency
               , getEventGenerator()
               , gRandom
               , _precision
               );
       //}
     
 #ifdef CHECK_INTEGRATOR
     cout << "WARNING in DalitzPdfBaseFastInteg::getNorm!!!"
      << " using old, inefficient integrator."
      << " Won't work for efficiencies!!"
      << endl;
     _mcint.initialise(evtPtr->eventPattern()
               , getAmps()//this
               , getEventGenerator()
               , gRandom
               , _precision
               );
 #endif
   }
 
   _norm = _faint.getVal();
 
 
 #ifdef CHECK_INTEGRATOR
   double fastVal = _norm;
   double oldVal = _mcint.getVal();
   cout << "DalitzPdfBaseFastInteg::getNorm:"
        << " compare mcint integral with fast: "
        << "\n\t old: " << oldVal << " , fast: " << fastVal
        << ", ratio " << oldVal/fastVal
        << endl;
   _norm = oldVal;
 #endif
 
 
   _integrating = false;
   return _norm > 0;
 }

◆ getVal() [1/2]

double DalitzPdfBaseFastInteg::getVal ( IDalitzEvent & evt )

virtual

Implements IDalitzPdf.

Definition at line 238 of file DalitzPdfBaseFastInteg.cpp.

                                                       {
   if(_pat.empty()) _pat = evt.eventPattern();
   if(_integrating) return getVal_withPs(evt);
   else return getVal_noPs(evt);
 }

◆ getVal() [2/2]

virtual double DalitzPdfBaseFastInteg::getVal ( IDalitzEvent * evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 131 of file DalitzPdfBaseFastInteg.h.

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

◆ getVal_noPs() [1/2]

double DalitzPdfBaseFastInteg::getVal_noPs ( IDalitzEvent & evt )

virtual

Implements IDalitzPdf.

Definition at line 243 of file DalitzPdfBaseFastInteg.cpp.

                                                            {
   bool dbThis=false;
   static double maxVal=0;
   if(_pat.empty()) _pat = evt.eventPattern();
 
   if(_integrating){
     // shouldn't really do that - use getVal or getVal_withPs
     // when you integrate (automatically done in getVal()):
     return un_normalised_noPs(evt);
   }else{
     if(_norm == -1) getNorm();
     double num = un_normalised_noPs(evt);
     if(dbThis){
       double val=num/_norm;
       if(fabs(val) > maxVal){
     maxVal=fabs(val);
     if(dbThis){
       cout  << "biggest un_normalised / norm, yet: " 
         << num << " / " <<_norm 
         << " = " << num/_norm
         << ". ps = " << phaseSpace(evt)
         << endl;
       cout << "for event: "; evt.print(); cout << endl;
     }
       }
     }
 
     /*
     cout << "\n DalitzPdfBaseFastInteg::getVal_noPs():"
      << " for event ptr " << getEvent()
      << " returning "
      << num << " / " << _norm << endl;
     */
       
     _val= num/_norm;
     if(dbThis) cout << "DalitzPdfBaseFastInteg::getVal_noPs: returning "
             << _val << endl;
     return num/_norm;
   }
 }

◆ getVal_noPs() [2/2]

virtual double DalitzPdfBaseFastInteg::getVal_noPs ( IDalitzEvent * evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 139 of file DalitzPdfBaseFastInteg.h.

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

◆ getVal_withPs() [1/2]

double DalitzPdfBaseFastInteg::getVal_withPs ( IDalitzEvent & evt )

virtual

Implements IDalitzPdf.

Definition at line 283 of file DalitzPdfBaseFastInteg.cpp.

                                                              {
   if(_pat.empty()) _pat = evt.eventPattern();
 
   if(_integrating){
     return un_normalised_noPs(evt)*phaseSpace(evt);
   }else{
     if(_norm == -1) getNorm();
     double num = un_normalised_noPs(evt)*phaseSpace(evt);
     return num/_norm;
   }
 }

◆ getVal_withPs() [2/2]

virtual double DalitzPdfBaseFastInteg::getVal_withPs ( IDalitzEvent * evt )

inlinevirtual

Implements IDalitzPdf.

Definition at line 135 of file DalitzPdfBaseFastInteg.h.

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

◆ Gradient()

void DalitzPdfBaseFastInteg::Gradient	(	IDalitzEvent &	evt,
		std::vector< double > &	grad,
		MINT::MinuitParameterSet *	mps
	)

virtual

Reimplemented from MINT::IPdf< IDalitzEvent >.

Definition at line 295 of file DalitzPdfBaseFastInteg.cpp.

                                                                                                          {
     
     if(_val == -1)un_normalised_noPs(evt);
     if(_norm == -1) getNorm();
     
     unsigned int N= mps->size();
     
     std::vector<double> gradAmpSum(N);
     for(unsigned int j=0; j < gradAmpSum.size(); j++) gradAmpSum[j]= 0.;   
     
     _amps->Gradient(evt,gradAmpSum,mps);
     
     if(_redoGradNorm){
       for(unsigned int j=0; j < _gradNorm.size(); j++) _gradNorm[j]= 0.;   
       _faint.Gradient(mps, _gradNorm);
       _redoGradNorm=false;
     }
     
     if(_gradNorm.size() != grad.size() || gradAmpSum.size() !=grad.size()
        ||  gradAmpSum.size() != _gradNorm.size()){
       cout << "array size mismatch in DalitzPdfBaseFastInteg::Gradient" << endl;
       throw "out";
     }
        
     for (unsigned int i=0; i<N; i++) {
       if(mps->getParPtr(i)->hidden() || mps->getParPtr(i)->iFixInit())continue;
       grad.at(i)= gradAmpSum.at(i)/_norm- _val/_norm*_gradNorm.at(i);
     }
     
     
 }

◆ GradientForLasso()

void DalitzPdfBaseFastInteg::GradientForLasso ( std::vector< double > & grad )

virtual

Definition at line 327 of file DalitzPdfBaseFastInteg.cpp.

                                                                  {
     _faint.GradientForLasso(this->getMPS(),grad);
 }

◆ histoSet() [1/2]

DalitzHistoSet DalitzPdfBaseFastInteg::histoSet ( ) const

virtual

Definition at line 334 of file DalitzPdfBaseFastInteg.cpp.

                                                      {
   return (_faint.histoSet());
 }

◆ histoSet() [2/2]

DalitzHistoSet DalitzPdfBaseFastInteg::histoSet ( )

virtual

Implements IDalitzPdf.

Definition at line 337 of file DalitzPdfBaseFastInteg.cpp.

                                                {
     // non-const version to satisfy IDalitzPdf
   return (_faint.histoSet());
 }

◆ integrating()

bool DalitzPdfBaseFastInteg::integrating ( )

protectedvirtual

Definition at line 231 of file DalitzPdfBaseFastInteg.cpp.

                                         {
   return _integrating;
 }

◆ interferenceHistoSet() [1/2]

DalitzHistoSet DalitzPdfBaseFastInteg::interferenceHistoSet ( ) const

virtual

Definition at line 351 of file DalitzPdfBaseFastInteg.cpp.

                                                                  {
     return (_faint.interferenceHistoSet());
 }

◆ interferenceHistoSet() [2/2]

DalitzHistoSet DalitzPdfBaseFastInteg::interferenceHistoSet ( )

virtual

Definition at line 354 of file DalitzPdfBaseFastInteg.cpp.

                                                            {
     // non-const version to satisfy IDalitzPdf
     return (_faint.interferenceHistoSet());
 }

◆ makeAmps()

bool DalitzPdfBaseFastInteg::makeAmps ( )

protected

Definition at line 218 of file DalitzPdfBaseFastInteg.cpp.

                                      {
   cout << "DalitzPdfBaseFastInteg::makeAmps() with MPS pointer " << getMPS() 
        << " and _amps pointer " << _amps << endl;
   if(0 == _amps && (! _pat.empty())){
       cout << "DalitzPdfBaseFastInteg::getAmps() making _amps with pointer " << getMPS() 
        << endl;
     _amps = new FitAmpSum(_pat, getMPS());
     counted_ptr<IFastAmplitudeIntegrable> cp(_amps);
     _countedAmps = cp;
   }
   return (bool) _amps;
 }

◆ makeDefaultGenerator()

IEventGenerator< IDalitzEvent > * DalitzPdfBaseFastInteg::makeDefaultGenerator ( )

protected

Definition at line 89 of file DalitzPdfBaseFastInteg.cpp.

                                                                            {
   if(_pat.empty()) return 0;
   if(0 == getAmps()) return 0;
   counted_ptr<IEventGenerator<IDalitzEvent> > 
     ptr(getAmps()->makeEventGenerator(_pat));
   _defaultGenerator = ptr;
   return _defaultGenerator.get();
 }

◆ makePlots()

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

Definition at line 331 of file DalitzPdfBaseFastInteg.cpp.

                                                                      {
   return histoSet().save(filename);
 }

◆ numberOfFitFractionsLargerThanThreshold()

int DalitzPdfBaseFastInteg::numberOfFitFractionsLargerThanThreshold ( double threshold )

inline

Definition at line 212 of file DalitzPdfBaseFastInteg.h.

                                                                {
         return _faint.numberOfFitFractionsLargerThanThreshold(threshold);
   }

◆ parametersChanged()

void DalitzPdfBaseFastInteg::parametersChanged ( )

virtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 118 of file DalitzPdfBaseFastInteg.cpp.

                                               {
   getNorm();
   _val= -1;
   _redoGradNorm=true;
 }

◆ phaseSpace()

double DalitzPdfBaseFastInteg::phaseSpace ( IDalitzEvent & evt )

protectedvirtual

Definition at line 235 of file DalitzPdfBaseFastInteg.cpp.

                                                           {
   return evt.phaseSpace();
 }

◆ RealVal()

virtual double DalitzPdfBaseFastInteg::RealVal ( IDalitzEvent & evt )

inlinevirtual

Reimplemented from MINT::PdfBase< IDalitzEvent >.

Definition at line 153 of file DalitzPdfBaseFastInteg.h.

153 {return getVal(evt);}

DalitzPdfBaseFastInteg::getVal

virtual double getVal(IDalitzEvent &evt)

Definition: DalitzPdfBaseFastInteg.cpp:238

◆ redoIntegrator()

double DalitzPdfBaseFastInteg::redoIntegrator ( )

inline

Definition at line 186 of file DalitzPdfBaseFastInteg.h.

                          {
     _norm = -1;
     getNorm();
     return _norm;
   }

◆ saveEachAmpsHistograms()

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

Definition at line 342 of file DalitzPdfBaseFastInteg.cpp.

                                                                                 {
   _faint.saveEachAmpsHistograms(prefix);
   return;
 }

◆ saveIntegrator()

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

Definition at line 20 of file DalitzPdfBaseFastInteg.cpp.

                                                                        {
   return _faint.save(fname);
 }

◆ saveInterferenceHistograms()

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

Definition at line 358 of file DalitzPdfBaseFastInteg.cpp.

                                                                                     {
     _faint.saveInterferenceHistograms(prefix);
     return;
 }

◆ setEventGenerator()

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

inline

Definition at line 116 of file DalitzPdfBaseFastInteg.h.

116 {_generator=g;}

DalitzPdfBaseFastInteg::_generator

MINT::IEventGenerator< IDalitzEvent > * _generator

Definition: DalitzPdfBaseFastInteg.h:69

g

static const double g

Definition: CLHEPSystemOfUnits.h:165

◆ setIntegrationPrecision()

void DalitzPdfBaseFastInteg::setIntegrationPrecision ( double prec )

Definition at line 110 of file DalitzPdfBaseFastInteg.cpp.

                                                                {
   _precision = prec;
 #ifdef CHECK_INTEGRATOR
   _mcint.setPrecision(_precision);
 #endif
   _faint.setPrecision(_precision);
 }

◆ setIntegratorFileName()

void DalitzPdfBaseFastInteg::setIntegratorFileName ( const std::string & commaSeparatedList )

Definition at line 375 of file DalitzPdfBaseFastInteg.cpp.

                                                                                      {
   _commaSepList_of_SavedIntegrators=commaSeparatedList;
 }

◆ setup()

void DalitzPdfBaseFastInteg::setup ( )

protected

Definition at line 215 of file DalitzPdfBaseFastInteg.cpp.

                                   {
   makeAmps(); // it's important to do this at construction time, otherwise the MinuitParameterSet will be empty.
 }

◆ sumOfFitFractions()

double DalitzPdfBaseFastInteg::sumOfFitFractions ( )

inline

Definition at line 192 of file DalitzPdfBaseFastInteg.h.

                              {
         if(!_faint.initialised())getNorm();
         return _faint.sumOfFitFractions();
   }  

◆ sumOfSqrtFitFractions()

double DalitzPdfBaseFastInteg::sumOfSqrtFitFractions ( )

inline

Definition at line 207 of file DalitzPdfBaseFastInteg.h.

                                  {
       if(!_faint.initialised())getNorm();
       return _faint.sumOfSqrtFitFractions();
   }   

◆ un_normalised_noPs()

virtual double DalitzPdfBaseFastInteg::un_normalised_noPs ( IDalitzEvent & evt )

protectedpure virtual

Implemented in AmpsPdf, and DalitzPdfSaveInteg.

◆ useAnalyticGradient()

bool DalitzPdfBaseFastInteg::useAnalyticGradient ( )

inlinevirtual

Reimplemented from MINT::IPdf< IDalitzEvent >.

Definition at line 151 of file DalitzPdfBaseFastInteg.h.

151 {return getAmps()->useAnalyticGradient();}

DalitzPdfBaseFastInteg::getAmps

IFastAmplitudeIntegrable * getAmps()

Definition: DalitzPdfBaseFastInteg.h:91

IFastAmplitudeIntegrable::useAnalyticGradient

virtual bool useAnalyticGradient()

Definition: IFastAmplitudeIntegrable.h:45

Member Data Documentation

◆ _amps

IFastAmplitudeIntegrable* DalitzPdfBaseFastInteg::_amps

protected

Definition at line 65 of file DalitzPdfBaseFastInteg.h.

◆ _commaSepList_of_SavedIntegrators

std::string DalitzPdfBaseFastInteg::_commaSepList_of_SavedIntegrators

protected

Definition at line 73 of file DalitzPdfBaseFastInteg.h.

◆ _countedAmps

MINT::counted_ptr<IFastAmplitudeIntegrable> DalitzPdfBaseFastInteg::_countedAmps

protected

Definition at line 66 of file DalitzPdfBaseFastInteg.h.

◆ _defaultGenerator

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

protected

Definition at line 71 of file DalitzPdfBaseFastInteg.h.

◆ _faint

FastAmplitudeIntegrator DalitzPdfBaseFastInteg::_faint

protected

Definition at line 64 of file DalitzPdfBaseFastInteg.h.

◆ _generator

MINT::IEventGenerator<IDalitzEvent>* DalitzPdfBaseFastInteg::_generator

protected

Definition at line 69 of file DalitzPdfBaseFastInteg.h.

◆ _gradNorm

std::vector<double> DalitzPdfBaseFastInteg::_gradNorm

protected

Definition at line 75 of file DalitzPdfBaseFastInteg.h.

◆ _integrating

bool DalitzPdfBaseFastInteg::_integrating

protected

Definition at line 70 of file DalitzPdfBaseFastInteg.h.

◆ _mps

MINT::MinuitParameterSet* DalitzPdfBaseFastInteg::_mps

protected

Definition at line 57 of file DalitzPdfBaseFastInteg.h.

◆ _norm

double DalitzPdfBaseFastInteg::_norm

protected

Definition at line 59 of file DalitzPdfBaseFastInteg.h.

◆ _pat

DalitzEventPattern DalitzPdfBaseFastInteg::_pat

protected

Definition at line 58 of file DalitzPdfBaseFastInteg.h.

◆ _precision

double DalitzPdfBaseFastInteg::_precision

protected

Definition at line 60 of file DalitzPdfBaseFastInteg.h.

◆ _redoGradNorm

bool DalitzPdfBaseFastInteg::_redoGradNorm

protected

Definition at line 76 of file DalitzPdfBaseFastInteg.h.

◆ _val

double DalitzPdfBaseFastInteg::_val

protected

Definition at line 74 of file DalitzPdfBaseFastInteg.h.

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

Mint/DalitzPdfBaseFastInteg.h
src/Mojito/DalitzFits/DalitzPdfBaseFastInteg.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ DalitzPdfBaseFastInteg() [1/3]

◆ DalitzPdfBaseFastInteg() [2/3]

◆ DalitzPdfBaseFastInteg() [3/3]

◆ ~DalitzPdfBaseFastInteg()

Member Function Documentation

◆ absSumOfInterferenceFractions()

◆ absSumOfSqrtInterferenceFractions()

◆ doFinalStats()

◆ endFit()

◆ getAmps()

◆ GetEachAmpsHistograms()

◆ getEventGenerator() [1/2]

◆ getEventGenerator() [2/2]

◆ getFitAmpSum()

◆ getIntegralValue() [1/2]

◆ getIntegralValue() [2/2]

◆ GetInterferenceHistograms()

◆ getMPS() [1/2]

◆ getMPS() [2/2]

◆ getNewVal() [1/2]

◆ getNewVal() [2/2]

◆ getNorm()

◆ getVal() [1/2]

◆ getVal() [2/2]

◆ getVal_noPs() [1/2]

◆ getVal_noPs() [2/2]

◆ getVal_withPs() [1/2]

◆ getVal_withPs() [2/2]

◆ Gradient()

◆ GradientForLasso()

◆ histoSet() [1/2]

◆ histoSet() [2/2]

◆ integrating()

◆ interferenceHistoSet() [1/2]

◆ interferenceHistoSet() [2/2]

◆ makeAmps()

◆ makeDefaultGenerator()

◆ makePlots()

◆ numberOfFitFractionsLargerThanThreshold()

◆ parametersChanged()

◆ phaseSpace()

◆ RealVal()

◆ redoIntegrator()

◆ saveEachAmpsHistograms()

◆ saveIntegrator()

◆ saveInterferenceHistograms()

◆ setEventGenerator()

◆ setIntegrationPrecision()

◆ setIntegratorFileName()

◆ setup()

◆ sumOfFitFractions()

◆ sumOfSqrtFitFractions()

◆ un_normalised_noPs()

◆ useAnalyticGradient()

Member Data Documentation

◆ _amps

◆ _commaSepList_of_SavedIntegrators

◆ _countedAmps

◆ _defaultGenerator

◆ _faint

◆ _generator

◆ _gradNorm

◆ _integrating

◆ _mps

◆ _norm

◆ _pat

◆ _precision

◆ _redoGradNorm

◆ _val