00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022 #include "EvtGenBase/EvtPatches.hh"
00023 #include "EvtGenModels/EvtPFermi.hh"
00024
00025
00026
00027
00028 #include <math.h>
00029
00030
00031
00032
00033
00034 EvtPFermi::EvtPFermi(const double &a, const double &mB, const double &mb)
00035 {
00036 _a = a;
00037 _mb = mb;
00038 _mB = mB;
00039 }
00040
00041
00042
00043
00044
00045 EvtPFermi::~EvtPFermi( )
00046 {
00047 }
00048
00049
00050
00051
00052
00053 double EvtPFermi::getFPFermi(const double &kplus)
00054 {
00055 double FKplus;
00056 double x = kplus/(_mB-_mb);
00057
00058 if ( x >= 1) return 0;
00059 if ( kplus <= -_mb) return 0;
00060
00061 FKplus = pow(1-x,_a)*exp((1+_a)*x);
00062
00063 return FKplus;
00064 }