[u/mrichter/AliRoot.git] / TEvtGen / EvtGenModels / EvtSLBKPoleFF.cxx

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//      Copyright (C) 1998      Caltech, UCSB
//
// Module: EvtSLBKPoleFF.cc
//
// Description: Routine to implement semileptonic form factors
//              according to the model SLBKPoles
//
// Modification history:
//
//    liheng       October 20,2005       Module created
//
//------------------------------------------------------------------------
// 
#include "EvtGenBase/EvtPatches.hh"
#include "EvtGenBase/EvtReport.hh"
#include "EvtGenModels/EvtSLBKPoleFF.hh"//modified
#include <string>
#include "EvtGenBase/EvtPDL.hh"
#include <math.h>
#include <stdlib.h>

EvtSLBKPoleFF::EvtSLBKPoleFF(int numarg, double *arglist) {//modified
   numSLBKPoleargs = numarg;//modified
   for (int i=0; i<numarg; i++) {
      SLBKPoleargs[i] = arglist[i]; }//modified

   return;
}


void EvtSLBKPoleFF::getscalarff(EvtId parent,EvtId daught,
                       double t, double mass, double *fpf,
			    double *f0f ) {

// Form factors have a general form, with parameters passed in
// from the arguments.

   if ( numSLBKPoleargs != 4 ) {//modified
     report(ERROR,"EvtGen") << "Problem in EvtSLBKPoleFF::getscalarff\n";
     report(ERROR,"EvtGen") << "wrong number of arguments!\n";
     report(ERROR,"EvtGen") << "number args:"<<numSLBKPoleargs<<" (expected 4)\n";
     report(ERROR,"EvtGen") << "Parent:"<<EvtPDL::name(parent)<<"\n";
     report(ERROR,"EvtGen") << "Daughter:"<<EvtPDL::name(daught)<<"\n";

   }


  double f0,af,powf;

  //double a_0, a_1, a_2, a_3, a_4, a_5, a_6, a_7;

  f0 = SLBKPoleargs[0];//f0
  af = SLBKPoleargs[1];//alpha
  //bf = SLBKPoleargs[2];
  double mass_star2 = SLBKPoleargs[3]*SLBKPoleargs[3];
  powf = 1.0;
  *fpf = f0/(pow( 1.0 - (1.0+af)*(t/mass_star2) + (af*((t/mass_star2)*(t/mass_star2))),powf));//modified

  f0 = SLBKPoleargs[0];//f0
  af = SLBKPoleargs[2];//beta
  //bf = SLBKPoleargs[6];
  powf = 1.0;

  *f0f = f0/(pow( 1.0 - (t/mass_star2/af),powf));//modified

  return;
}

void EvtSLBKPoleFF::getvectorff(EvtId parent,EvtId daught,
                       double t, double mass, double *a1f,
			     double *a2f, double *vf, double *a0f ){

   if ( numSLBKPoleargs !=8 ) {//modified
     report(ERROR,"EvtGen") << "Problem in EvtSLBKPoleFF::getvectorff\n";//modified
     report(ERROR,"EvtGen") << "wrong number of arguements!!!\n";
     report(ERROR,"EvtGen") << numSLBKPoleargs<<"\n";//modified
//     printf("\n*********************%d*********************",numSLBKPoleargs);
  }


   report(INFO,"EvtGen")<<"Check the implementation of EvtSLBKPoleFF::getvectorff()!\n";


  double mb=EvtPDL::getMeanMass(parent);
  double mb2 = mb*mb;

//modified-begin
  static EvtId B0=EvtPDL::getId("B0");
  static EvtId B0B=EvtPDL::getId("anti-B0");  
  static EvtId BP=EvtPDL::getId("B+");
  static EvtId BM=EvtPDL::getId("B-");
  static EvtId BS0=EvtPDL::getId("B_s0");

  static EvtId B0S=EvtPDL::getId("B*0");
  static EvtId BPMS=EvtPDL::getId("B*+");  
  static EvtId BS0S=EvtPDL::getId("B_s*0");
  
  static EvtId D0=EvtPDL::getId("D0");
  static EvtId D0B=EvtPDL::getId("anti-D0");
  static EvtId DP=EvtPDL::getId("D+");
  static EvtId DM=EvtPDL::getId("D-");
  static EvtId DSP=EvtPDL::getId("D_s+");
  static EvtId DSM=EvtPDL::getId("D_s-");

  static EvtId D0S=EvtPDL::getId("D*0");
  static EvtId DPMS=EvtPDL::getId("D*+");
  static EvtId DSPMS=EvtPDL::getId("D_s*+");

  double mass_star;
  double mass_star2;
  if(parent==B0||parent==B0B){
     mass_star=EvtPDL::getMeanMass(B0S);
     mass_star2=mass_star*mass_star;
  }
  if(parent==BP||parent==BM){
     mass_star=EvtPDL::getMeanMass(BPMS);
     mass_star2=mass_star*mass_star;
  }
  if(parent==BS0){
     mass_star=EvtPDL::getMeanMass(BS0S);
     mass_star2=mass_star*mass_star;
  }

  if(parent==D0||parent==D0B){
     mass_star=EvtPDL::getMeanMass(D0S);
     mass_star2=mass_star*mass_star;
  }
  if(parent==DP||parent==DM){
     mass_star=EvtPDL::getMeanMass(DPMS);
     mass_star2=mass_star*mass_star;
  }
  if(parent==DSP||parent==DSM){
     mass_star=EvtPDL::getMeanMass(DSPMS);
     mass_star2=mass_star*mass_star;
  }
//modified-end

  double f0,af,bf,powf;

  f0 = SLBKPoleargs[2];//A1
  af = SLBKPoleargs[6];//b'
  bf = 0;//0
  powf = 1.0;//1.0
  *a1f = f0/(pow( 1.0 - af*t/mass_star2,powf));//modified

  f0 = SLBKPoleargs[3];//A2
  af = SLBKPoleargs[6];//b'
  bf = SLBKPoleargs[7];//b''==0
  powf = 1.0;//1.0

  *a2f = f0/(pow(1.0 - (af+bf)*(t/mass_star2) + (af*bf)*((t/mass_star2)*(t/mass_star2)),powf));//modified

  f0 = SLBKPoleargs[0];//V0
  af = SLBKPoleargs[4];//a
  bf = 0;//0
  powf = 1.0;//1.0

  *vf = f0/(pow( 1.0 - (1.0+af)*(t/mass_star2) + af*(t/mass_star2)*(t/mass_star2),powf));//modified

  f0 = SLBKPoleargs[1];//A0
  af = SLBKPoleargs[5];//a'
  bf = 0;//0
  powf = 1.0;//1.0

  *a0f = f0/(pow( 1.0 - (1.0+af)*(t/mb2) + af*((t/mb2)*(t/mb2)),powf));//modified
  return;
 }


void EvtSLBKPoleFF::gettensorff(EvtId parent,EvtId daught,
                       double t, double mass, double *hf,
			     double *kf, double *bpf, double *bmf ){

  if ( numSLBKPoleargs !=16 ) {
     report(ERROR,"EvtGen") << "Problem in EvtSLBKPoleFF::gettensorff\n";
     report(ERROR,"EvtGen") << "wrong number of arguements!!!\n";
  }

  report(INFO,"EvtGen")<<"Check the implementation of EvtSLBKPoleFF::gettensorff()!\n";

  double mb=EvtPDL::getMeanMass(parent);
  double mb2 = mb*mb;

  double f0,af,bf,powf;

  f0 = SLBKPoleargs[0];
  af = SLBKPoleargs[1];
  bf = SLBKPoleargs[2];
  powf = SLBKPoleargs[3];
  *hf = f0/(pow( 1.0 + (af*t/mb2) + (bf*((t/mb2)*(t/mb2))),powf));

  f0 = SLBKPoleargs[4];
  af = SLBKPoleargs[5];
  bf = SLBKPoleargs[6];
  powf = SLBKPoleargs[7];

  *kf = f0/(pow( 1.0 + (af*t/mb2) + (bf*((t/mb2)*(t/mb2))),powf));

  f0 = SLBKPoleargs[8];
  af = SLBKPoleargs[9];
  bf = SLBKPoleargs[10];
  powf = SLBKPoleargs[11];

  *bpf = f0/(pow( 1.0 + (af*t/mb2) + (bf*((t/mb2)*(t/mb2))),powf));

  f0 = SLBKPoleargs[12];
  af = SLBKPoleargs[13];
  bf = SLBKPoleargs[14];
  powf = SLBKPoleargs[15];

  *bmf = f0/(pow( 1.0 + (af*t/mb2) + (bf*((t/mb2)*(t/mb2))),powf));
  return;
 }

void EvtSLBKPoleFF::getbaryonff(EvtId, EvtId, double, double, double*, 
				double*, double*, double*){
  
  report(ERROR,"EvtGen") << "Not implemented :getbaryonff in EvtSLBKPoleFF.\n";  
  ::abort();

}
Commit	Line	Data
da0e9ce3	1	//--------------------------------------------------------------------------
	2	//
	3	// Environment:
	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.
	7	//
	8	// Copyright Information: See EvtGen/COPYRIGHT
	9	// Copyright (C) 1998 Caltech, UCSB
	10	//
	11	// Module: EvtSLBKPoleFF.cc
	12	//
	13	// Description: Routine to implement semileptonic form factors
	14	// according to the model SLBKPoles
	15	//
	16	// Modification history:
	17	//
	18	// liheng October 20,2005 Module created
	19	//
	20	//------------------------------------------------------------------------
	21	//
	22	#include "EvtGenBase/EvtPatches.hh"
	23	#include "EvtGenBase/EvtReport.hh"
	24	#include "EvtGenModels/EvtSLBKPoleFF.hh"//modified
	25	#include <string>
	26	#include "EvtGenBase/EvtPDL.hh"
	27	#include <math.h>
	28	#include <stdlib.h>
	29
	30	EvtSLBKPoleFF::EvtSLBKPoleFF(int numarg, double *arglist) {//modified
	31	numSLBKPoleargs = numarg;//modified
	32	for (int i=0; i<numarg; i++) {
	33	SLBKPoleargs[i] = arglist[i]; }//modified
	34
	35	return;
	36	}
	37
	38
	39	void EvtSLBKPoleFF::getscalarff(EvtId parent,EvtId daught,
	40	double t, double mass, double *fpf,
	41	double *f0f ) {
	42
	43	// Form factors have a general form, with parameters passed in
	44	// from the arguments.
	45
	46	if ( numSLBKPoleargs != 4 ) {//modified
	47	report(ERROR,"EvtGen") << "Problem in EvtSLBKPoleFF::getscalarff\n";
	48	report(ERROR,"EvtGen") << "wrong number of arguments!\n";
	49	report(ERROR,"EvtGen") << "number args:"<<numSLBKPoleargs<<" (expected 4)\n";
	50	report(ERROR,"EvtGen") << "Parent:"<<EvtPDL::name(parent)<<"\n";
	51	report(ERROR,"EvtGen") << "Daughter:"<<EvtPDL::name(daught)<<"\n";
	52
	53	}
	54
	55
	56	double f0,af,powf;
	57
	58	//double a_0, a_1, a_2, a_3, a_4, a_5, a_6, a_7;
	59
	60	f0 = SLBKPoleargs[0];//f0
	61	af = SLBKPoleargs[1];//alpha
	62	//bf = SLBKPoleargs[2];
	63	double mass_star2 = SLBKPoleargs[3]*SLBKPoleargs[3];
	64	powf = 1.0;
65	fpf = f0/(pow( 1.0 - (1.0+af)(t/mass_star2) + (af((t/mass_star2)(t/mass_star2))),powf));//modified
66
67	f0 = SLBKPoleargs[0];//f0
68	af = SLBKPoleargs[2];//beta
69	//bf = SLBKPoleargs[6];
70	powf = 1.0;
71
72	*f0f = f0/(pow( 1.0 - (t/mass_star2/af),powf));//modified
73
74	return;
75	}
76
77	void EvtSLBKPoleFF::getvectorff(EvtId parent,EvtId daught,
78	double t, double mass, double *a1f,
79	double a2f, double vf, double *a0f ){
80
81	if ( numSLBKPoleargs !=8 ) {//modified
82	report(ERROR,"EvtGen") << "Problem in EvtSLBKPoleFF::getvectorff\n";//modified
83	report(ERROR,"EvtGen") << "wrong number of arguements!!!\n";
84	report(ERROR,"EvtGen") << numSLBKPoleargs<<"\n";//modified
85	// printf("\n*******************%d*******************",numSLBKPoleargs);
86	}
87
88
89	report(INFO,"EvtGen")<<"Check the implementation of EvtSLBKPoleFF::getvectorff()!\n";
90
91
92	double mb=EvtPDL::getMeanMass(parent);
93	double mb2 = mb*mb;
94
95	//modified-begin
96	static EvtId B0=EvtPDL::getId("B0");
97	static EvtId B0B=EvtPDL::getId("anti-B0");
98	static EvtId BP=EvtPDL::getId("B+");
99	static EvtId BM=EvtPDL::getId("B-");
100	static EvtId BS0=EvtPDL::getId("B_s0");
101
102	static EvtId B0S=EvtPDL::getId("B*0");
103	static EvtId BPMS=EvtPDL::getId("B*+");
104	static EvtId BS0S=EvtPDL::getId("B_s*0");
105
106	static EvtId D0=EvtPDL::getId("D0");
107	static EvtId D0B=EvtPDL::getId("anti-D0");
108	static EvtId DP=EvtPDL::getId("D+");
109	static EvtId DM=EvtPDL::getId("D-");
110	static EvtId DSP=EvtPDL::getId("D_s+");
111	static EvtId DSM=EvtPDL::getId("D_s-");
112
113	static EvtId D0S=EvtPDL::getId("D*0");
114	static EvtId DPMS=EvtPDL::getId("D*+");
115	static EvtId DSPMS=EvtPDL::getId("D_s*+");
116
117	double mass_star;
118	double mass_star2;
119	if(parent==B0\|\|parent==B0B){
120	mass_star=EvtPDL::getMeanMass(B0S);
121	mass_star2=mass_star*mass_star;
122	}
123	if(parent==BP\|\|parent==BM){
124	mass_star=EvtPDL::getMeanMass(BPMS);
125	mass_star2=mass_star*mass_star;
126	}
127	if(parent==BS0){
128	mass_star=EvtPDL::getMeanMass(BS0S);
129	mass_star2=mass_star*mass_star;
130	}
131
132	if(parent==D0\|\|parent==D0B){
133	mass_star=EvtPDL::getMeanMass(D0S);
134	mass_star2=mass_star*mass_star;
135	}
136	if(parent==DP\|\|parent==DM){
137	mass_star=EvtPDL::getMeanMass(DPMS);
138	mass_star2=mass_star*mass_star;
139	}
140	if(parent==DSP\|\|parent==DSM){
141	mass_star=EvtPDL::getMeanMass(DSPMS);
142	mass_star2=mass_star*mass_star;
143	}
144	//modified-end
145
146	double f0,af,bf,powf;
147
148	f0 = SLBKPoleargs[2];//A1
149	af = SLBKPoleargs[6];//b'
150	bf = 0;//0
151	powf = 1.0;//1.0
152	a1f = f0/(pow( 1.0 - aft/mass_star2,powf));//modified
153
154	f0 = SLBKPoleargs[3];//A2
155	af = SLBKPoleargs[6];//b'
156	bf = SLBKPoleargs[7];//b''==0
157	powf = 1.0;//1.0
158
159	a2f = f0/(pow(1.0 - (af+bf)(t/mass_star2) + (afbf)((t/mass_star2)*(t/mass_star2)),powf));//modified
160
161	f0 = SLBKPoleargs[0];//V0
162	af = SLBKPoleargs[4];//a
163	bf = 0;//0
164	powf = 1.0;//1.0
165
166	vf = f0/(pow( 1.0 - (1.0+af)(t/mass_star2) + af(t/mass_star2)(t/mass_star2),powf));//modified
167
168	f0 = SLBKPoleargs[1];//A0
169	af = SLBKPoleargs[5];//a'
170	bf = 0;//0
171	powf = 1.0;//1.0
172
173	a0f = f0/(pow( 1.0 - (1.0+af)(t/mb2) + af((t/mb2)(t/mb2)),powf));//modified
174	return;
175	}
176
177
178
179	void EvtSLBKPoleFF::gettensorff(EvtId parent,EvtId daught,
180	double t, double mass, double *hf,
181	double kf, double bpf, double *bmf ){
182
183	if ( numSLBKPoleargs !=16 ) {
184	report(ERROR,"EvtGen") << "Problem in EvtSLBKPoleFF::gettensorff\n";
185	report(ERROR,"EvtGen") << "wrong number of arguements!!!\n";
186	}
187
188	report(INFO,"EvtGen")<<"Check the implementation of EvtSLBKPoleFF::gettensorff()!\n";
189
190	double mb=EvtPDL::getMeanMass(parent);
191	double mb2 = mb*mb;
192
193	double f0,af,bf,powf;
194
195	f0 = SLBKPoleargs[0];
196	af = SLBKPoleargs[1];
197	bf = SLBKPoleargs[2];
198	powf = SLBKPoleargs[3];
199	hf = f0/(pow( 1.0 + (aft/mb2) + (bf((t/mb2)(t/mb2))),powf));
200
201	f0 = SLBKPoleargs[4];
202	af = SLBKPoleargs[5];
203	bf = SLBKPoleargs[6];
204	powf = SLBKPoleargs[7];
205
206	kf = f0/(pow( 1.0 + (aft/mb2) + (bf((t/mb2)(t/mb2))),powf));
207
208	f0 = SLBKPoleargs[8];
209	af = SLBKPoleargs[9];
210	bf = SLBKPoleargs[10];
211	powf = SLBKPoleargs[11];
212
213	bpf = f0/(pow( 1.0 + (aft/mb2) + (bf((t/mb2)(t/mb2))),powf));
214
215	f0 = SLBKPoleargs[12];
216	af = SLBKPoleargs[13];
217	bf = SLBKPoleargs[14];
218	powf = SLBKPoleargs[15];
219
220	bmf = f0/(pow( 1.0 + (aft/mb2) + (bf((t/mb2)(t/mb2))),powf));
221	return;
222	}
223
224	void EvtSLBKPoleFF::getbaryonff(EvtId, EvtId, double, double, double*,
225	double, double, double*){
226
227	report(ERROR,"EvtGen") << "Not implemented :getbaryonff in EvtSLBKPoleFF.\n";
228	::abort();
229
230	}
231
232