MINT2
Public Member Functions | Protected Attributes | List of all members
PhaseSpaceIntegral4bodyWith_s123s12 Class Reference

#include <phaseSpaceIntegrals.h>

Public Member Functions

 PhaseSpaceIntegral4bodyWith_s123s12 (IGenFct *s123f, IGenFct *s12f, const DalitzEventPattern &pat)
 
double getVal () const
 

Protected Attributes

IGenFct_s123Fct
 
IGenFct_s12Fct
 
DalitzEventPattern _pat
 

Detailed Description

Definition at line 78 of file phaseSpaceIntegrals.h.

Constructor & Destructor Documentation

◆ PhaseSpaceIntegral4bodyWith_s123s12()

PhaseSpaceIntegral4bodyWith_s123s12::PhaseSpaceIntegral4bodyWith_s123s12 ( IGenFct s123f,
IGenFct s12f,
const DalitzEventPattern pat 
)

Definition at line 352 of file phaseSpaceIntegrals.cpp.

356  : _s123Fct(s123f)
357  , _s12Fct(s12f)
358  , _pat(pat)
359 {
360 }

Member Function Documentation

◆ getVal()

double PhaseSpaceIntegral4bodyWith_s123s12::getVal ( ) const

Definition at line 424 of file phaseSpaceIntegrals.cpp.

424  {
425  bool dbThis=true;
426 
427  double xmi = _s123Fct->getRhoMi();
428  double xma = _s123Fct->getRhoMa();
429 
430  double ymi = _s12Fct->getRhoMi();
431  double yma = _s12Fct->getRhoMa();
432 
433  const int numParticles=5;
434 
437 
438  TF2 intFunc("d_ps4_by_ds13ds12_withFct"
440  , xmi, xma
441  , ymi, yma
442  , numParticles
443  );
444 
445  Double_t p[numParticles];
446  for(int i=0; i < numParticles; i++){
447  p[i] = _pat[i].mass();
448  }
449 
450  intFunc.SetParameters(p);
451 
452  Double_t a[2];
453  Double_t b[2];
454  a[0]=xmi; b[0]=xma;
455  a[1]=ymi; b[1]=yma;
456 
457  Int_t nDim = 2;
458  Int_t mipt = 4000000;
459  Int_t mapt = 400000000;
460  Double_t epsilon = 1.e-8;
461  Double_t relerr = -1;
462  Int_t nfnevl = 9999;
463  Int_t iFail = -9999;
464  double returnVal = intFunc.IntegralMultiple(nDim, a, b
465  , mipt, mapt
466  , epsilon, relerr
467  , nfnevl, iFail);
468 
469  if(dbThis){
470  cout << "123, 12: just called : "
471  << " intFunc.IntegralMultiple( " << nDim
472  << ", " << a[0] << "-" << b[0]
473  << ", " << a[1] << "-" << b[1]
474  << ", " << mipt <<", " << mapt
475  << ", " << epsilon << ", " << relerr
476  << ", " << nfnevl << ", " << iFail
477  << " );"
478  << " = " << returnVal << endl;
479  // cout << " ..... compare with: " << intFunc.Integral(xmi, xma, ymi, yma) << endl;
480  }
481 
482  return returnVal;
483 
484  // return intFunc.Integral(xmi, xma, ymi, yma);
485 
486 }
IGenFct * Global_s12Fct
virtual double getRhoMi() const
Definition: IGenFct.h:38
double d_ps4_by_ds123ds12_withFct(Double_t *x, Double_t *p)
IGenFct * Global_s123Fct
virtual double getRhoMa() const
Definition: IGenFct.h:39

Member Data Documentation

◆ _pat

DalitzEventPattern PhaseSpaceIntegral4bodyWith_s123s12::_pat
protected

Definition at line 81 of file phaseSpaceIntegrals.h.

◆ _s123Fct

IGenFct* PhaseSpaceIntegral4bodyWith_s123s12::_s123Fct
protected

Definition at line 80 of file phaseSpaceIntegrals.h.

◆ _s12Fct

IGenFct * PhaseSpaceIntegral4bodyWith_s123s12::_s12Fct
protected

Definition at line 80 of file phaseSpaceIntegrals.h.


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