phaseSpaceIntegrals.h File Reference

#include "Mint/DalitzEventPattern.h"
#include "Mint/IGenFct.h"
#include "Mint/lambda.h"
#include "Mint/CLHEPPhysicalConstants.h"
#include "TF1.h"
#include <iostream>
#include <cmath>

Go to the source code of this file.

Classes
class	PhaseSpaceIntegral3body
class	PhaseSpaceIntegral4body
class	PhaseSpaceIntegral4bodyWith_s12s34
class	PhaseSpaceIntegral4bodyWith_s123s12

Functions
double	phaseSpaceIntegral2body (double mum, double d1, double d2)
double	phaseSpaceIntegral2body (const DalitzEventPattern &_pat)
double	d3body_by_ds12 (Double_t *x, Double_t *p)
double	d4body_by_ds12 (Double_t *x, Double_t *p)
double	phaseSpaceIntegral_upTo4body (const DalitzEventPattern &pat)

Function Documentation

d3body_by_ds12()

double d3body_by_ds12	(	Double_t *	x,
		Double_t *	p
	)

Definition at line 52 of file phaseSpaceIntegrals.cpp.

                                               {
  bool dbThis=false;

  double s12 = x[0];
  if(s12 < 0) return 0;
  double m12 = sqrt(s12);
  double mum = p[0];
  double m1  = p[1];
  double m2  = p[2];
  double m3  = p[3];
  if(dbThis) cout << " d3body_by_ds12 called with " 
          << mum << ", "
          << m1 << ", "
          << m2 << ", "
          << m3 << endl;

  return phaseSpaceIntegral2body(mum, m12, m3)
    *    phaseSpaceIntegral2body(m12, m1, m2);
}

d4body_by_ds12()

double d4body_by_ds12	(	Double_t *	x,
		Double_t *	p
	)

Definition at line 22 of file phaseSpaceIntegrals.cpp.

                                               {
  bool dbThis=false;
  double s12 = x[0];
  if(s12 < 0) return 0;
  double m12 = sqrt(s12);
  double mum = p[0];
  double m1  = p[1];
  double m2  = p[2];
  double m3  = p[3];
  double m4  = p[4];

  if(dbThis) cout << " d4body_by_ds12 called with " 
          << mum << ", "
          << m1 << ", "
          << m2 << ", "
          << m3 << ", "
          << m4 << endl;

  PhaseSpaceIntegral3body p3;

  double p3val = p3.getVal(mum, m12, m3, m4);
  double p2val = phaseSpaceIntegral2body(m12, m1, m2);
  if(dbThis){
    cout << " d4body_by_ds12: returning " << p2val << " * " << p3val
     << " = " << p2val * p3val
     << endl;
  }
  return p2val * p3val;
}

phaseSpaceIntegral2body() [1/2]

double phaseSpaceIntegral2body ( double  mum, double  d1, double  d2  )

inline

Definition at line 20 of file phaseSpaceIntegrals.h.

                                                                       {
  if(mum <=0 ) return 0;
  if(mum < d1 + d2) return 0;

  double la = lambda(mum*mum, d1*d1, d2*d2);
  if(la <= 0) return 0;

  return pi * sqrt(la)/(2*mum*mum);
}

phaseSpaceIntegral2body() [2/2]

double phaseSpaceIntegral2body ( const DalitzEventPattern &  _pat )

inline

Definition at line 30 of file phaseSpaceIntegrals.h.

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

phaseSpaceIntegral_upTo4body()

double phaseSpaceIntegral_upTo4body

(

const DalitzEventPattern &

pat

)

Definition at line 488 of file phaseSpaceIntegrals.cpp.

                                                                  {

  if(pat.numDaughters() < 2 ){
    return 0;
  }else if(pat.numDaughters() == 2 ){
    return phaseSpaceIntegral2body(pat);
  }else if(pat.numDaughters() == 3){
    PhaseSpaceIntegral3body p3;
    return p3.getVal(pat);
  }else if(pat.numDaughters() == 4){
    PhaseSpaceIntegral4body p4;
    return p4.getVal(pat);
  }else{
    return -9999;
  }

}