Dalitz.h

Go to the documentation of this file.

// -*- C++ -*-
 
//
// Author:      Dan Ambrose
// Created:     Mon Mar 5 2012
// based on Cleo package by Eric White and Warner Sun
 
#ifndef DALITZ_H
#define DALITZ_H
 
// #include "EvtGen/EvtComplex.h"
#include <complex>
using namespace std;
typedef complex<double> TComplex;
 
// For using "string"
/*
#include "DataHandler/Frame.h"
#include "DataHandler/Record.h"
#include "Experiment/Experiment.h"
#include "Experiment/report.h"
#include "Experiment/units.h"
#include "FrameAccess/FAItem.h"
#include "FrameAccess/FATable.h"
#include "FrameAccess/extract.h"
*/
using std::string;
 
// Class Dalitz
class Dalitz {
 
public:
  Dalitz();
  Dalitz( int binNum );
  TComplex Amplitude( double x, double y, double z );
  double   Phase( double x, double y, double z,
                  int Babar = 1 ); // 1 for Babar model, 2 for CLEO model, else for Belle model
  bool     Point_on_DP( double x, double y );
  bool     Point_on_DP2( double x, double y );
  TComplex CLEO_resAmp( double mAC, double mBC, double mAB, double mA, double mB, double mC,
                        double _ampl, double _theta, double _gamma, double _bwm, int _spin );
  TComplex resAmp( double mAC, double mBC, double mAB, double mA, double mB, double mC,
                   double _ampl, double _theta, double _gamma, double _bwm, int _spin );
  TComplex f_980( double mPP, double mR, double _ampl, double _theta );
  TComplex sakurai( double mkp, double mkm, double mpp, double _ampl, double _theta,
                    double gamma_r, double m_r );
  TComplex Babar_sakurai( double mkp, double mkm, double mpp, double _ampl, double _theta,
                          double gamma_r, double m_r );
  TComplex Babar_resAmp( double mAC, double mBC, double mAB, double mB, double mA, double mC,
                         double _ampl, double _theta, double _gamma, double _bwm, int _spin );
  TComplex Babar_Amplitude( double x, double y, double z );
  TComplex CLEO_Amplitude( double x, double y, double z );
  int      getBin( double mx, double my, double mz );
 
private:
  static double    PI;          // pi
  int              N;           // Bin number, set in constructor
  static const int Nbin = 1000; // number of dx for integral
};
 
#endif