#include <iostream>
#include "T3piEtaCrossPart.h"
#include "TConstant.h"
#include <complex>
Go to the source code of this file.
Typedefs | |
typedef complex< double > | complex_t |
Functions | |
complex_t | operator+ (const double &x, const complex_t &y) |
complex_t | operator- (const double &x, const complex_t &y) |
complex_t | Amp (double s) |
complex_t | bw_rho (double s) |
complex_t | bw_om (double s) |
Definition at line 9 of file T3piEtaCrossPart.C.
complex_t Amp | ( | double | s | ) | [inline] |
Definition at line 18 of file T3piEtaCrossPart.C.
References alpha, exp(), rb::Geepsi, rb::Geepsip, rb::Geepsipp, rb::Gpsi, rb::Gpsip, rb::Gpsipp, I, M_PI, rb::Mpsi, rb::Mpsi2, rb::Mpsip, rb::Mpsip2, rb::Mpsipp, and rb::Mpsipp2.
00018 { 00019 const complex_t I(0,1); 00020 const double Apsi = alpha*Gpsi/(3*Geepsi); 00021 const double phipsi = -M_PI/2; 00022 const double Apsip = alpha*Gpsip/(3*Geepsip); 00023 const double phipsip = -M_PI/2; 00024 const double Apsipp = alpha*Gpsipp/(3*Geepsipp); 00025 const double phipsipp = -M_PI/2; 00026 complex_t res = 3*sqrt(s)/alpha* 00027 ( 00028 Geepsi*(1 + Apsi*exp(I*phipsi))/(s - Mpsi2 + I*Mpsi*Gpsi) + 00029 Geepsip*(1 + Apsip*exp(I*phipsip))/(s - Mpsip2 + I*Mpsip*Gpsip) + 00030 Geepsipp*(1 + Apsipp*exp(I*phipsipp))/(s - Mpsipp2 + I*Mpsipp*Gpsipp)); 00031 return res; 00032 }
complex_t bw_om | ( | double | s | ) | [inline] |
Definition at line 50 of file T3piEtaCrossPart.C.
References rb::Gomega, and rb::Momega.
Referenced by J5pi(), T4piCrossPart::SetJ(), and T3piEtaCrossPart::SetJ().
00050 { 00051 complex_t Domega = complex_t(Momega*Momega-s,-Momega*Gomega); 00052 complex_t bw = 1./Domega; 00053 return bw; 00054 }
complex_t bw_rho | ( | double | s | ) | [inline] |
Definition at line 34 of file T3piEtaCrossPart.C.
References rb::Grho, rb::Grhop, rb::Mrho, and rb::Mrhop.
00034 { 00035 const double sig = -0.1; 00036 complex_t Drho = complex_t(Mrho*Mrho-s,-Mrho*Grho); 00037 complex_t Drhop = complex_t(Mrhop*Mrhop-s,-Mrhop*Grhop); 00038 complex_t bw = (1./Drho + sig/Drhop); 00039 // 00040 // dumping factor added in omega part: to be substituted with a formfactor 00041 // 00042 const double Et2 = (2.2*1e3)*(2.2*1e3); 00043 if(s>Et2){ 00044 double r = Et2/s; 00045 bw *= r*r; 00046 } 00047 return bw; 00048 }