#include <EvtVVS.hh>
Public Member Functions | |
VVS (EvtVector4R pd1, EvtVector4R pd2, EvtVector4R pd3, double ResonanceMass, double ResonanceWidth, double r1, double r2, double phase1, double phase2) | |
double | Fij (int i, int j, double r) |
double | R00 (double r) |
EvtComplex | amps1 (int m, int i, int j) |
double | amps () |
Private Attributes | |
EvtVector4R | _pd [5] |
double | _res [5] |
double | _par [5] |
Definition at line 35 of file EvtVVS.hh.
VVS::VVS | ( | EvtVector4R | pd1, | |
EvtVector4R | pd2, | |||
EvtVector4R | pd3, | |||
double | ResonanceMass, | |||
double | ResonanceWidth, | |||
double | r1, | |||
double | r2, | |||
double | phase1, | |||
double | phase2 | |||
) | [inline] |
Definition at line 37 of file EvtVVS.hh.
References _par, _pd, and _res.
00037 { 00039 // For the decay 1-- -->R(1--) +Pseudoscalar-->3 Pseudoscalars 00040 // pd1,pd2,pd3: four momentum vector for three daughters 00041 // ResonanceMass, ResonanceWidth : resonance parameter for intermediate state R(1--) 00042 // r1,r2 : the relative coupling strength for two resonaces R R.P.T another R 00043 // phase1,phase2: the relatev phase of two phases R.P.T another R 00045 _pd[0]=pd1; 00046 _pd[1]=pd2; 00047 _pd[2]=pd3; 00048 _res[0]=ResonanceMass; 00049 _res[1]=ResonanceWidth; 00050 _par[0]=r1; 00051 _par[1]=phase1; 00052 _par[2]=r2; 00053 _par[3]=phase2; 00054 }
double VVS::amps | ( | ) |
Definition at line 71 of file EvtVVS.cc.
References _par, abs, amps1(), exp(), and subSeperate::temp.
Referenced by EvtRhoPi::AmplitudeSquare().
00071 { 00072 double r1=_par[0]; 00073 double phase1=_par[1]; 00074 double r2=_par[2]; 00075 double phase2=_par[3]; 00076 double temp=0.0; 00077 int m; 00078 for(m=-1;m<=1;m+=2) temp=temp+pow(abs(r1*exp(EvtComplex(0,phase1))*amps1(m,0,1)+r2*exp(EvtComplex(0,phase2))*amps1(m,0,2)+amps1(m,1,2)),2.); 00079 return temp; 00080 }
EvtComplex VVS::amps1 | ( | int | m, | |
int | i, | |||
int | j | |||
) |
Definition at line 49 of file EvtVVS.cc.
References _pd, _res, EvtHelSys::Angles(), EvtVector4R::d3mag(), Djmn(), Fij(), EvtVector4R::mass2(), phi1, R00(), s, and subSeperate::temp.
Referenced by amps().
00049 { 00050 double mrho=_res[0]; 00051 double wrho=_res[1]; 00052 double s; 00053 EvtComplex img(0.0,1.0),dpro; 00054 EvtVector4R prho; 00055 prho=_pd[i]+_pd[j]; 00056 s=prho.mass2(); 00057 dpro=s-mrho*mrho+img*sqrt(s)*wrho; 00058 EvtHelSys angles(prho,_pd[i]),labAngles; 00059 double theta,phi,ct1,st1,phi1,st,ct; 00060 double rpp=angles.getHelAng(0); 00061 theta=angles.getHelAng(1); 00062 phi =angles.getHelAng(2); 00063 ct1 =labAngles.Angles(prho,1); 00064 phi1=labAngles.Angles(prho,2); 00065 int lamb; 00066 EvtComplex temp(0.0,0.0); 00067 for(lamb=-1;lamb<=1;lamb++) temp=temp+Fij(lamb,0,prho.d3mag())*Djmn(1,m,lamb,phi1,ct1,0.0)/dpro*R00(rpp)*Djmn(1,lamb,0,phi,theta,0.0); 00068 return temp; 00069 }
double VVS::Fij | ( | int | i, | |
int | j, | |||
double | r | |||
) |
Definition at line 40 of file EvtVVS.cc.
References _pd, mass, and subSeperate::temp.
Referenced by amps1().
00040 { 00041 double mpsi=(_pd[0]+_pd[1]+_pd[2]).mass(); 00042 double temp=mpsi*r; 00043 if(i==0&& j==0) return 0; 00044 if(i==1&& j==0) return temp; 00045 if(i==-1&& j==0) return -temp; 00046 00047 }
double VVS::R00 | ( | double | r | ) |
double VVS::_par[5] [private] |
EvtVector4R VVS::_pd[5] [private] |
double VVS::_res[5] [private] |