1 //--------------------------------------------------------------------------
4 // This software is part of the EvtGen package developed jointly
5 // for the BaBar and CLEO collaborations. If you use all or part
6 // of it, please give an appropriate acknowledgement.
8 // Copyright Information: See EvtGen/COPYRIGHT
9 // Copyright (C) 1998 Caltech, UCSB
11 // Module: EvtMelikhovFF.cc
13 // Description: form factors for B->D*lnu according to HQET
15 // Modification history:
17 // DJL April 17, 1998 Module created
19 //------------------------------------------------------------------------
21 #include "EvtGenBase/EvtPatches.hh"
22 #include "EvtGenBase/EvtPatches.hh"
23 #include "EvtGenBase/EvtReport.hh"
24 #include "EvtGenModels/EvtMelikhovFF.hh"
25 #include "EvtGenBase/EvtId.hh"
28 #include "EvtGenBase/EvtPDL.hh"
31 EvtMelikhovFF::EvtMelikhovFF(double arg1) {
33 whichfit = int(arg1+0.1);
37 void EvtMelikhovFF::getvectorff(EvtId parent,EvtId,
38 double t, double mass, double *a1f,
39 double *a2f, double *vf, double *a0f ){
41 double ma1(0.0),ra1(0.0),na1(0.0);
42 double ma2(0.0),ra2(0.0),na2(0.0);
43 double mv(0.0),rv(0.0),nv(0.0);
90 double mb=EvtPDL::getMeanMass(parent);
91 //double w = ((mb*mb)+(mass*mass)-t)/(2.0*mb*mass);
93 double melr = mass/mb;
94 double mely = t/(mb*mb);
97 *a1f = ((1.0+melr*melr-mely)/(1+melr))*ra1/pow(1.0-(t/(ma1*ma1)),na1);
98 *a2f = (1+melr)*((1.0-melr*melr-mely)/((1+melr)*(1+melr)-mely))
99 *ra2/pow(1.0-(t/(ma2*ma2)),na2);
100 *vf = (1+melr)*rv/pow(1.0-(t/(mv*mv)),nv);
107 void EvtMelikhovFF::getscalarff(EvtId, EvtId, double, double, double*,
110 report(ERROR,"EvtGen") << "Not implemented :getvectorff in EvtMelikhovFF.\n";
117 void EvtMelikhovFF::gettensorff(EvtId, EvtId, double, double, double*,
118 double*, double*, double*){
120 report(ERROR,"EvtGen") << "Not implemented :gettensorff in EvtMelikhovFF.\n";
127 void EvtMelikhovFF::getbaryonff(EvtId, EvtId, double, double, double*,
128 double*, double*, double*){
130 report(ERROR,"EvtGen") << "Not implemented :getbaryonff in EvtMelikhovFF.\n";