PhaseSpaceIntegral3body Class Reference

#include <phaseSpaceIntegrals.h>

Public Member Functions
	PhaseSpaceIntegral3body ()
double	getVal (const DalitzEventPattern &_pat)
double	getVal (double mum, double d1, double d2, double d3)

Static Protected Attributes
static TF1 *	_f = 0

Detailed Description

Definition at line 40 of file phaseSpaceIntegrals.h.

Constructor & Destructor Documentation

PhaseSpaceIntegral3body()

PhaseSpaceIntegral3body::PhaseSpaceIntegral3body

(

)

Definition at line 187 of file phaseSpaceIntegrals.cpp.

                                                {
  if(0 == _f){
    _f = new TF1("d3body_by_ds12", d3body_by_ds12
         , 0, 10*GeV*GeV
         , 4
         );
  }
  if(0 == _f){
    cout << "phaseSpaceIntegral3body::phaseSpaceIntegral3body()"
     << " failed to get TF1!!!" << endl;
  }
}

Member Function Documentation

getVal() [1/2]

double PhaseSpaceIntegral3body::getVal

(

const DalitzEventPattern &

_pat

)

Definition at line 201 of file phaseSpaceIntegrals.cpp.

                                                                    {
  if(_pat.size() != 4){
    cout << "phaseSpaceIntegral3body: wrong pattern " << _pat << endl;
  }
  double mum = _pat[0].mass();
  double d1  = _pat[1].mass();
  double d2  = _pat[2].mass();
  double d3  = _pat[3].mass();

  return getVal(mum, d1, d2, d3);
}

getVal() [2/2]

double PhaseSpaceIntegral3body::getVal	(	double	mum,
		double	d1,
		double	d2,
		double	d3
	)

Definition at line 213 of file phaseSpaceIntegrals.cpp.

                                              {
  bool dbThis=false;
  double min_m12 = (d1 + d2);
  double max_m12 = (mum - d3);
  if(dbThis)cout << "PhaseSpaceIntegral3body::getVal : from, to " 
         << min_m12 << ", " << max_m12 << endl;
  if(min_m12 >= max_m12) return 0;

  double min_s12 = min_m12*min_m12;
  double max_s12 = max_m12*max_m12;

  Double_t para[4]={mum, d1, d2, d3};

  _f->SetParameters(para);


  /*
  Int_t np = 1000;
  double x[np];
  double w[np];
  _f->CalcGaussLegendreSamplingPoints(np,x,w,1e-15);
  double returnVal = _f->IntegralFast(np,x,w, min_s12, max_s12, para);
  */

  double returnVal =  _f->Integral(min_s12, max_s12);//, para);
  if(dbThis) cout << "PhaseSpaceIntegral3body::getVal() returning " 
          << returnVal << endl;
  return returnVal;
}

Member Data Documentation

_f

TF1 * PhaseSpaceIntegral3body::_f = 0

staticprotected

Definition at line 42 of file phaseSpaceIntegrals.h.

---

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

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