PhaseSpaceIntegral4bodyWith_s12s34 Class Reference

#include <phaseSpaceIntegrals.h>

Public Member Functions
	PhaseSpaceIntegral4bodyWith_s12s34 (IGenFct *s12f, IGenFct *s34f, const DalitzEventPattern &pat)
double	getVal () const

Protected Attributes
IGenFct *	_s12Fct
IGenFct *	_s34Fct
DalitzEventPattern	_pat

Detailed Description

Definition at line 66 of file phaseSpaceIntegrals.h.

Constructor & Destructor Documentation

PhaseSpaceIntegral4bodyWith_s12s34()

PhaseSpaceIntegral4bodyWith_s12s34::PhaseSpaceIntegral4bodyWith_s12s34	(	IGenFct *	s12f,
		IGenFct *	s34f,
		const DalitzEventPattern &	pat
	)

Definition at line 341 of file phaseSpaceIntegrals.cpp.

  : _s12Fct(s12f)
  , _s34Fct(s34f)
  , _pat(pat)
{
}

Member Function Documentation

getVal()

double PhaseSpaceIntegral4bodyWith_s12s34::getVal

(

)

const

Definition at line 362 of file phaseSpaceIntegrals.cpp.

                                                      {
  bool dbThis=true;

  double xmi = _s12Fct->getRhoMi();
  double xma = _s12Fct->getRhoMa();

  double ymi = _s34Fct->getRhoMi();
  double yma = _s34Fct->getRhoMa();

  const int numParticles=5;

  Global_s12Fct = _s12Fct;
  Global_s34Fct = _s34Fct;

  Double_t p[numParticles];
  for(int i=0; i < numParticles; i++){
    p[i] = _pat[i].mass();
  }

  TF2 intFunc("d_ps4_by_ds12ds34_withFct"
          , d_ps4_by_ds12ds34_withFct
          , xmi, xma
          , ymi, yma
          , numParticles
          );
  

  intFunc.SetParameters(p);

  Double_t a[2];
  Double_t b[2];
  a[0]=xmi; b[0]=xma;
  a[1]=ymi; b[1]=yma;

  Int_t nDim    = 2;
  Int_t mipt    =    4000000;
  Int_t mapt    =  400000000;
  Double_t epsilon =  1.e-8;
  Double_t relerr  = -1;
  Int_t nfnevl = 9999;
  Int_t iFail = -9999;
  double returnVal = intFunc.IntegralMultiple(nDim, a, b
                          , mipt, mapt
                          , epsilon, relerr
                          , nfnevl, iFail);
  if(dbThis){
    cout << "12, 34:  just called : "
     << " intFunc.IntegralMultiple( " << nDim 
     << ", " << a[0] << "-" << b[0]
     << ", " << a[1] << "-" << b[1] 
     << ", " << mipt <<", " << mapt
     << ", " << epsilon << ", " << relerr 
     << ", " << nfnevl << ", " << iFail
     << " );" 
     << " = " << returnVal << endl;
    //    cout << " ..... compare with: " << intFunc.Integral(xmi, xma, ymi, yma) << endl;

  }
  return returnVal;

}

Member Data Documentation

_pat

DalitzEventPattern PhaseSpaceIntegral4bodyWith_s12s34::_pat

protected

Definition at line 69 of file phaseSpaceIntegrals.h.

_s12Fct

IGenFct* PhaseSpaceIntegral4bodyWith_s12s34::_s12Fct

protected

Definition at line 68 of file phaseSpaceIntegrals.h.

_s34Fct

IGenFct * PhaseSpaceIntegral4bodyWith_s12s34::_s34Fct

protected

Definition at line 68 of file phaseSpaceIntegrals.h.

---

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

• Mint/phaseSpaceIntegrals.h
• src/Mojito/BycklingKajantje/phaseSpaceIntegrals.cpp