[u/mrichter/AliRoot.git] / TEvtGen / EvtGenModels / EvtVubBLNP.hh


//////////////////////////////////////////////////////////////////////
//
// Module: EvtVubBLNP.hh
//
// Description: Modeled on Riccardo Faccini's EvtVubNLO module
// author: sheila  date: October 2005
//
// tripleDiff from BLNP's notebook, based on hep-ph/0504071
//
//////////////////////////////////////////////////////////////////


#ifndef EVTVUBBLNP_HH
#define EVTVUBBLNP_HH

#include <vector>
#include "EvtGenBase/EvtDecayIncoherent.hh"

class EvtParticle;

class EvtVubBLNP:public  EvtDecayIncoherent  {

public:
  
  EvtVubBLNP() {}
  virtual ~EvtVubBLNP();

  std::string getName();

  EvtDecayBase* clone();

  void initProbMax();

  void init();

  void decay(EvtParticle *Bmeson); 

private:

  // Input parameters
  double mBB;
  double lambda2;

  // Shape function parameters
  double b;
  double Lambda;
  double Ecut;
  double wzero;

  // SF and SSF modes
  int itype;
  double dtype;
  int isubl;

  // flags
  int flag1;
  int flag2;
  int flag3;

  // Quark mass
  double mb;

  // Matching scales
  double muh;
  double mui;
  double mubar;

  // Perturbative quantities
  double CF;
  double CA;

  double beta0;
  double beta1;
  double beta2;

  double zeta3;

  double Gamma0;
  double Gamma1;
  double Gamma2;

  double gp0;
  double gp1;

  double Lbar;
  double mupisq;
  double moment2;

  int flagpower;
  int flag2loop;

  int maxLoop;
  double precision;

  std::vector<double> gvars;
  
  double rate3(double Pp, double Pl, double Pm);
  double F1(double Pp, double Pm, double muh, double mui, double mubar, double doneJS, double done1);
  double F2(double Pp, double Pm, double muh, double mui, double mubar, double done3);
  double F3(double Pp, double Pm, double muh, double mui, double mubar, double done2);
  double DoneJS(double Pp, double Pm, double mui);
  double Done1(double Pp, double Pm, double mui);
  double Done2(double Pp, double Pm, double mui);
  double Done3(double Pp, double Pm, double mui);
  static double IntJS(double what, const std::vector<double> &vars);
  static double Int1(double what, const std::vector<double> &vars);
  static double Int2(double what, const std::vector<double> &vars);
  static double Int3(double what, const std::vector<double> &vars);
  static double g1(double w, const std::vector<double> &vars);
  static double g2(double w, const std::vector<double> &vars);
  static double g3(double w, const std::vector<double> &vars);
  static double Shat(double w, const std::vector<double> &vars);
  static double Mzero(double muf, double mu, double mupisq, const std::vector<double> &vars);
  double wS(double w);
  double t(double w);
  double u(double w);
  double v(double w);
  double myfunction(double w, double Lbar, double mom2);
  double myfunctionBIK(double w, double Lbar, double mom2);
  double dU1nlo(double muh, double mui);
  double U1lo(double muh, double mui);
  double Sfun(double mu1, double mu2, double epsilon);
  double S0(double a1, double r);
  double S1(double a1, double r);
  double S2(double a1, double r);
  double aGamma(double mu1, double mu2, double epsilon);
  double agp(double mu1, double mu2, double epsilon);
  double alo(double muh, double mui);
  double anlo(double muh, double mui);   // d/depsilon of aGamma
  static double alphas(double mu, const std::vector<double> &vars);
  double PolyLog(double v, double z);
  static double Gamma(double z);
  static double Gamma(double a, double x);
  static double gamser(double a, double x, double LogGamma);
  static double gammcf(double a, double x, double LogGamma);
  double findBLNPWhat();
  std::vector<double> _pf;
};

#endif
Commit	Line	Data
da0e9ce3	1
	2	//////////////////////////////////////////////////////////////////////
	3	//
	4	// Module: EvtVubBLNP.hh
	5	//
	6	// Description: Modeled on Riccardo Faccini's EvtVubNLO module
	7	// author: sheila date: October 2005
	8	//
	9	// tripleDiff from BLNP's notebook, based on hep-ph/0504071
	10	//
	11	//////////////////////////////////////////////////////////////////
	12
	13
	14	#ifndef EVTVUBBLNP_HH
	15	#define EVTVUBBLNP_HH
	16
	17	#include <vector>
	18	#include "EvtGenBase/EvtDecayIncoherent.hh"
	19
	20	class EvtParticle;
	21
	22	class EvtVubBLNP:public EvtDecayIncoherent {
	23
	24	public:
	25
	26	EvtVubBLNP() {}
	27	virtual ~EvtVubBLNP();
	28
	29	std::string getName();
	30
	31	EvtDecayBase* clone();
	32
	33	void initProbMax();
	34
	35	void init();
	36
	37	void decay(EvtParticle *Bmeson);
	38
	39	private:
	40
	41	// Input parameters
	42	double mBB;
	43	double lambda2;
	44
	45	// Shape function parameters
	46	double b;
	47	double Lambda;
	48	double Ecut;
	49	double wzero;
	50
	51	// SF and SSF modes
	52	int itype;
	53	double dtype;
	54	int isubl;
	55
	56	// flags
	57	int flag1;
	58	int flag2;
	59	int flag3;
	60
	61	// Quark mass
	62	double mb;
	63
	64	// Matching scales
65	double muh;
66	double mui;
67	double mubar;
68
69	// Perturbative quantities
70	double CF;
71	double CA;
72
73	double beta0;
74	double beta1;
75	double beta2;
76
77	double zeta3;
78
79	double Gamma0;
80	double Gamma1;
81	double Gamma2;
82
83	double gp0;
84	double gp1;
85
86	double Lbar;
87	double mupisq;
88	double moment2;
89
90	int flagpower;
91	int flag2loop;
92
93	int maxLoop;
94	double precision;
95
96	std::vector<double> gvars;
97
98	double rate3(double Pp, double Pl, double Pm);
99	double F1(double Pp, double Pm, double muh, double mui, double mubar, double doneJS, double done1);
100	double F2(double Pp, double Pm, double muh, double mui, double mubar, double done3);
101	double F3(double Pp, double Pm, double muh, double mui, double mubar, double done2);
102	double DoneJS(double Pp, double Pm, double mui);
103	double Done1(double Pp, double Pm, double mui);
104	double Done2(double Pp, double Pm, double mui);
105	double Done3(double Pp, double Pm, double mui);
106	static double IntJS(double what, const std::vector<double> &vars);
107	static double Int1(double what, const std::vector<double> &vars);
108	static double Int2(double what, const std::vector<double> &vars);
109	static double Int3(double what, const std::vector<double> &vars);
110	static double g1(double w, const std::vector<double> &vars);
111	static double g2(double w, const std::vector<double> &vars);
112	static double g3(double w, const std::vector<double> &vars);
113	static double Shat(double w, const std::vector<double> &vars);
114	static double Mzero(double muf, double mu, double mupisq, const std::vector<double> &vars);
115	double wS(double w);
116	double t(double w);
117	double u(double w);
118	double v(double w);
119	double myfunction(double w, double Lbar, double mom2);
120	double myfunctionBIK(double w, double Lbar, double mom2);
121	double dU1nlo(double muh, double mui);
122	double U1lo(double muh, double mui);
123	double Sfun(double mu1, double mu2, double epsilon);
124	double S0(double a1, double r);
125	double S1(double a1, double r);
126	double S2(double a1, double r);
127	double aGamma(double mu1, double mu2, double epsilon);
128	double agp(double mu1, double mu2, double epsilon);
129	double alo(double muh, double mui);
130	double anlo(double muh, double mui); // d/depsilon of aGamma
131	static double alphas(double mu, const std::vector<double> &vars);
132	double PolyLog(double v, double z);
133	static double Gamma(double z);
134	static double Gamma(double a, double x);
135	static double gamser(double a, double x, double LogGamma);
136	static double gammcf(double a, double x, double LogGamma);
137	double findBLNPWhat();
138	std::vector<double> _pf;
139	};
140
141	#endif
142