00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020 #include "EvtGenBase/EvtPatches.hh"
00021
00022 #include "EvtGenBase/EvtPatches.hh"
00023 #include "EvtGenBase/EvtPDL.hh"
00024 #include "EvtGenModels/EvtbTosllAliFF.hh"
00025 #include <math.h>
00026
00027 EvtbTosllAliFF::EvtbTosllAliFF(){}
00028
00029
00030 void EvtbTosllAliFF::getScalarFF(EvtId parent, EvtId daught,
00031 double t, double mass,
00032 double& fp,double& f0,double& ft){
00033
00034 double m=EvtPDL::getMeanMass(parent);
00035
00036
00037 double shat=t/(m*m);
00038 double shat2=shat*shat;
00039 double shat3=shat2*shat;
00040
00041 fp = 0.278*exp(1.568*shat+0.470*shat2+0.885*shat3);
00042 f0 = 0.278*exp(0.740*shat+0.080*shat2+0.425*shat3);
00043 ft = 0.300*exp(1.600*shat+0.501*shat2+0.796*shat3);
00044
00045 }
00046
00047
00048 void EvtbTosllAliFF::getVectorFF(EvtId parent, EvtId daught,
00049 double t, double mass,
00050 double& a1,double& a2,double& a0, double& v,
00051 double& t1, double& t2, double& t3 ){
00052
00053
00054 double m=EvtPDL::getMeanMass(parent);
00055
00056
00057 double shat=t/(m*m);
00058 double shat2=shat*shat;
00059
00060
00061 a1=0.294*exp(0.656*shat+0.456*shat2);
00062 a2=0.246*exp(1.237*shat+0.822*shat2);
00063 a0=0.412*exp(1.543*shat+0.954*shat2);
00064 v=0.399*exp(1.537*shat+1.123*shat2);
00065
00066 t1=0.334*exp(1.575*shat+1.140*shat2);
00067 t2=0.334*exp(0.562*shat+0.481*shat2);
00068 t3=0.234*exp(1.230*shat+1.089*shat2);
00069
00070 }
00071
00072