Collaboration diagram for generators:

Modules
	generators data objects

	generators modules

Classes
class	AAFHInterface
	Class to inferface AAFH/DIAG36 Generator written in Fortran. More...

class	BabayagaNLO
	C++ Interface for the Fortran Bhabha and exclusive two photon generator BABAYAGA.NLO. More...

class	BBBrem
	Generator for low scattering angle radiative Bhabha events (Beam-Beam Bremsstrahlung). More...

class	BHWide
	C++ Interface for the Fortran Bhabha scattering generator BHWide. More...

class	SGCosmic
	Class to generate tracks in the cosmics generator and store them in a MCParticle graph. More...

class	CRY
	C++ Interface for the generator CRY. More...

class	EvtBSemiTauonicDecayRateCalculator
	Class to calculate the differential decay rate, R(D), R(D), polarizations of tau and D using BSTD model based on [M. More...

class	EvtGenInterface
	Class to interface EvtGen. More...

class	EvtGenModelRegister
	Class to keep a register of all Belle2 EvtDecayBases. More...

class	EvtB0toKsKK
	Register Decay model EvtB0toKsKK. More...

class	EvtBCL
	The class provides the form factors for orbitally excited semileptonic decays. More...

class	EvtBCLFF
	BCL Form Factors. More...

class	EvtBSemiTauonic
	The EvtGen model of semi-tauonic B decays including new physics effects based on [M. More...

class	EvtBSemiTauonic2HDMType2
	The EvtGen model of semi-tauonic B decays including the charged higgs effect of 2HDM Type2 based on [M. More...

class	EvtBSemiTauonicAmplitude
	The base class for the calculation of the spin dependent amplitudes for the BSemiTauonic model based on [M. More...

class	EvtBSemiTauonicHelicityAmplitudeCalculator
	The class calculates the helicity amplitude of semi-tauonic B decays including new physics effects based on [M. More...

class	EvtBSemiTauonicScalarMesonAmplitude
	The class calculates the spin dependent amplitudes of B->Dtaunu decays for the BSemiTauonic model based on [M. More...

class	EvtBSemiTauonicVectorMesonAmplitude
	The class calculates the spin dependent amplitudes of B->D*taunu decays for the BSemiTauonic model based on [M. More...

class	EvtBtoXsnunu_FERMI
	The evtgen model to produce non-resonant B-> Xs nu nubar decay sample. More...

class	EvtEtaFullDalitz
	Class for the simulation of the eta -> pi+pi-pi0 decay with an improved dalitz description. More...

class	EvtEtaPi0Dalitz
	Class for the simulation of the eta -> 3pi0 and eta'->3pi0 decays. More...

class	EvtEtaPrimeDalitz
	Class for the simulation of the eta' -> pi+ pi- eta and pi0 pi0 eta decays. More...

class	EvtKnunu
	The evtgen model to produce B-> K nu nubar decay sample. More...

class	EvtKstarnunu_REV
	The evtgen model to produce B-> Kstar nu nubar decay sample. More...

class	EvtPhspCP
	Register Decay model EvtPhspCP. More...

class	HepevtReader
	Class to read Hepevt files and store the content in a MCParticle graph. More...

class	HepMCReader
	Class to read HepMC files and store the content in a MCParticle graph. More...

class	KKGenInterface
	Interface class for using the KKMC generator. More...

class	KoralW
	C++ interface for the FORTRAN 4-fermion final state generator KoralW. More...

class	LHEReader
	Class to read LHE files and store the content in a MCParticle graph. More...

class	ParticleGun
	Class to generate tracks in the particle gun and store them in a MCParticle graph. More...

class	Phokhara
	C++ Interface for the Fortran generator PHOKHARA. More...

class	ReaderSAD
	Class to read files that have been created by SAD and store their content in a MCParticle graph. More...

class	Teegg
	C++ Interface for the Fortran generator TEEGG. More...

class	EvtBCLFFTest
	The fixture for testing the EvtBCLFF. More...

class	ParticleGunTest
	The fixture for testing the Particlegun. More...

class	TouschekReaderTURTLE
	Class to read Touschek files and store their content in a MCParticle graph. More...

struct	DualNumber
	Simple structure implementing dual numbers which are used for exact evaluation of the derivatives, see https://en.wikipedia.org/wiki/Automatic_differentiation#Automatic_differentiation_using_dual_numbers. More...

struct	GeneralVector< T >
	3-vector with members of arbitrary type, especially members can be dual numbers More...

class	InitialParticleGeneration
	Generate Collision. More...

Macros
#define	B2_EVTGEN_REGISTER_MODEL(classname)
	Class to register B2_EVTGEN_REGISTER_MODEL. More...

Functions
	B2_EVTGEN_REGISTER_MODEL (EvtB0toKsKK)
	register the model in EvtGen

double	s13_min (const double &s12, const double &M, const double &m1, const double &m2, const double &m3)
	minimum s13

double	s13_max (const double &s12, const double &M, const double &m1, const double &m2, const double &m3)
	maximum s13

	B2_EVTGEN_REGISTER_MODEL (EvtBCL)
	register the model in EvtGen

	B2_EVTGEN_REGISTER_MODEL (EvtBSemiTauonic)
	register the model in EvtGen

	B2_EVTGEN_REGISTER_MODEL (EvtBSemiTauonic2HDMType2)
	register the model in EvtGen

	B2_EVTGEN_REGISTER_MODEL (EvtBtoXsnunu_FERMI)
	register the model in EvtGen

	B2_EVTGEN_REGISTER_MODEL (EvtEtaFullDalitz)
	register the model in EvtGen

	B2_EVTGEN_REGISTER_MODEL (EvtEtaPi0Dalitz)
	register the model in EvtGen

	B2_EVTGEN_REGISTER_MODEL (EvtEtaPrimeDalitz)
	register the model in EvtGen

	B2_EVTGEN_REGISTER_MODEL (EvtKnunu)
	register the model in EvtGen

	B2_EVTGEN_REGISTER_MODEL (EvtKstarnunu_REV)
	register the model in EvtGen

	B2_EVTGEN_REGISTER_MODEL (EvtPhspCP)
	register the model in EvtGen

	TEST_F (EvtBCLFFTest, Scalar)
	Test calculation of scalar BCL FF.

	TEST_F (EvtBCLFFTest, Vector)
	Test calculation of vector BCL FF.

	TEST_F (ParticleGunTest, MomentumDistributions)
	Tests momentum generation parameters.

	TEST_F (ParticleGunTest, VertexDistributions)
	Tests vertex generation parameters.

	TEST_F (ParticleGunTest, AngularDistributions)
	Tests angular generation parameters.

double	sqrt (double a)
	sqrt for double

double	tan (double a)
	tan for double

double	atan (double a)
	atan for double

DualNumber	operator+ (DualNumber a, DualNumber b)
	addition of dual numbers

DualNumber	operator+ (DualNumber a, double b)
	addition of dual number and double

DualNumber	operator- (DualNumber a, double b)
	subtraction of dual number and double

DualNumber	operator- (double a, DualNumber b)
	subtraction of double and dual number

DualNumber	operator/ (double a, DualNumber b)
	division of double and dual number

DualNumber	operator/ (DualNumber a, DualNumber b)
	division of two dual numbers

DualNumber	operator- (DualNumber a, DualNumber b)
	subtraction of two dual numbers

DualNumber	operator* (DualNumber a, DualNumber b)
	multiplication of two dual numbers

DualNumber	operator* (double a, DualNumber b)
	multiplication of double and dual number

DualNumber	sqrt (DualNumber a)
	sqrt of dual number

DualNumber	atan (DualNumber a)
	atan of dual number

DualNumber	tan (DualNumber a)
	tan of dual number

template<typename T >
GeneralVector< T >	operator+ (GeneralVector< T > a, GeneralVector< T > b)
	addition of two general vectors

template<typename T >
T	dot (GeneralVector< T > a, GeneralVector< T > b)
	dot product of two general vectors

template<typename T >
GeneralVector< T >	operator* (T a, GeneralVector< T > b)
	product of scalar and general vector

template<typename T >
T	getEcms (GeneralVector< T > p1, GeneralVector< T > p2)
	get CMS energy from HER & LER momentum 3-vectors

template<typename T >
GeneralVector< T >	getBoost (GeneralVector< T > p1, GeneralVector< T > p2)
	get boost vector from HER & LER momentum 3-vectors

template<typename T >
std::pair< T, T >	getAnglesCMS (GeneralVector< T > p1, GeneralVector< T > p2)
	get collision axis angles in CM frame from HER & LER momentum 3-vectors

template<typename T >
GeneralVector< T >	getFourVector (T energy, T angleX, T angleY, bool isHER)
	get 4-momentum from energy and angles of beam

Eigen::MatrixXd	getGradientMatrix (double eH, double thXH, double thYH, double eL, double thXL, double thYL)
	get gradient matrix of (eH, thXH, thYH, eL, thXL, thYL) -> (eCMS, boostX, boostY, boostY, angleXZ, angleYZ) transformation

Eigen::VectorXd	getCentralValues (double eH, double thXH, double thYH, double eL, double thXL, double thYL)
	get vector including (Ecms, boostX, boostY, boostZ, angleXZ, angleYZ) from beam energies and angles

Eigen::MatrixXd	transformCov (const Eigen::MatrixXd &covHER, const Eigen::MatrixXd &covLER, const Eigen::MatrixXd &grad)
	transform covariance matrix to the variables (Ecms, boostX, boostY, boostZ, angleXZ, angleYZ)

double	getTotalEnergy (const std::vector< ROOT::Math::PxPyPzMVector > &ps)
	Main function scaleParticleEnergies in this header scales momenta in CMS of all particles in the final state by a constant factor such that the overall collision energy is slightly changed. More...

double	getScaleFactor (const std::vector< ROOT::Math::PxPyPzMVector > &ps, double EcmsTarget)
	Find approximative value of the momentum scaling factor which makes the total energy of the particles provided in the vector ps equal to EcmsTarget.

std::vector< ROOT::Math::PxPyPzMVector >	scaleMomenta (const std::vector< ROOT::Math::PxPyPzMVector > &ps, double C)
	Scale momenta of the particles in the vector ps with a factor C. More...

void	scaleParticleEnergies (MCParticleGraph &mpg, double EcmsTarget)
	Scale momenta of the particles by a constant factor such that total energy is changed to EcmsTarget. More...

static Eigen::MatrixXd	toEigen (TMatrixDSym m)
	transform matrix from ROOT to Eigen format

void	init ()
	Initializes the generator.

void	generateEvent (MCParticleGraph &mcGraph)
	Generates one single event. More...

void	term ()
	Terminates the generator. More...

std::string	getName ()
	Get function Name

EvtDecayBase *	clone ()
	Clone the decay of B0toKsKK.

void	init ()
	Initialize standard stream objects

void	initProbMax ()
	Initialize standard stream objects for probability function

void	decay (EvtParticle *p)
	Member of particle in EvtGen.

EvtVector4R	umu (const EvtVector4R &p4a, const EvtVector4R &p4b, const EvtVector4R &p4c)
	Function 4Vector umu.

EvtVector4R	Smu (const EvtVector4R &p4a, const EvtVector4R &p4b, const EvtVector4R &p4c)
	Function 4Vector Smu.

EvtVector4R	Lmu (const EvtVector4R &p4a, const EvtVector4R &p4b, const EvtVector4R &p4c)
	Function 4Vector Lmu.

EvtTensor4C	gmunu_tilde (const EvtVector4R &p4a, const EvtVector4R &p4b, const EvtVector4R &p4c)
	Function Tensor gmunu

EvtTensor4C	Tmunu (const EvtVector4R &p4a, const EvtVector4R &p4b, const EvtVector4R &p4c)
	Function Tensor Tmunu

EvtTensor4C	Multiply (const EvtTensor4C &t1, const EvtTensor4C &t2)
	Function Tensor Multiply

EvtTensor4C	RaiseIndices (const EvtTensor4C &t)
	Function RaiseIndices

void	RaiseIndex (EvtVector4R &vector)
	Member function RaiseIndices.

EvtTensor4C	Mmunu (const EvtVector4R &p4a, const EvtVector4R &p4b, const EvtVector4R &p4c)
	Function Tensor Mmunu.

double	BWBF (const double &q, const unsigned int &L)
	Meson radius

double	BWBF (const double &q, const double &q0, const unsigned int &L)
	Meson radius

EvtComplex	BreitWigner (const double &m, const double &m0, const double &Gamma0, const double &q, const double &q0, const unsigned int &L)
	BreitWigner Shape.

EvtVector4R	Boost (const EvtVector4R &p4, const EvtVector4R &boost)
	Parameter for boost frame

double	p (const double &mab, const double &M, const double &mc)
	Constants p

double	q (const double &mab, const double &ma, const double &mb)
	Constants q.

EvtComplex	Flatte_k (const double &s, const double &m_h)
	Constant Flatte_k.

EvtComplex	Flatte (const double &m, const double &m0)
	Constant Flatte.

EvtComplex	A_f0ks (const EvtVector4R &p4ks, const EvtVector4R &p4kp, const EvtVector4R &p4km)
	A_f0ks is amplitude of f0.

EvtComplex	A_phiks (const EvtVector4R &p4ks, const EvtVector4R &p4kp, const EvtVector4R &p4km)
	A_phiks is amplitude of phi.

EvtComplex	A_fxks (const EvtVector4R &p4ks, const EvtVector4R &p4kp, const EvtVector4R &p4km)
	A_fxks is amplitude of fxks.

EvtComplex	A_chic0ks (const EvtVector4R &p4ks, const EvtVector4R &p4kp, const EvtVector4R &p4km)
	A_chic0ks is amplitude of chic0ks.

EvtComplex	A_kknr (const EvtVector4R &p4k1, const EvtVector4R &p4k2, const double &alpha_kk)
	A_kknr is amplitude of kknr.

	EvtBCL ()
	Default constructor.

virtual	~EvtBCL ()
	virtual destructor

std::string	getName ()
	Returns name of module.

EvtDecayBase *	clone ()
	Clones module.

void	decay (EvtParticle *p)
	Creates a decay.

void	initProbMax ()
	Sets maximal probab.

void	init ()
	Initializes module.

	EvtBCLFF (int numarg, double *arglist)
	constructor

void	getscalarff (EvtId parent, EvtId daughter, double t, double, double fpf, double f0f)
	Scalar FF's. More...

void	getvectorff (EvtId parent, EvtId daughter, double t, double, double a1f, double a2f, double vf, double a0f)
	Vector FF's. More...

void	gettensorff (EvtId parent, EvtId daughter, double t, double, double hf, double kf, double bp, double bm)
	Not Implemented.

void	getbaryonff (EvtId, EvtId, double, double, double , double , double , double )
	Not Implemented.

void	getdiracff (EvtId, EvtId, double, double, double , double , double , double , double , double )
	Not Implemented.

void	getraritaff (EvtId, EvtId, double, double, double , double , double , double , double , double , double , double )
	Not Implemented.

	EvtBSemiTauonic ()
	The default constructor

virtual	~EvtBSemiTauonic ()
	The destructor

std::string	getName ()
	The function returns the model name. More...

EvtDecayBase *	clone ()
	The function makes a copy of an EvtBSTD object. More...

void	decay (EvtParticle *p)
	The function evaluates the decay amplitude of the parent particle. More...

void	initProbMax ()
	The function sets the maximum value of the probability.

void	init ()
	The function initializes the decay. More...

	EvtBSemiTauonic2HDMType2 ()
	The default constructor

virtual	~EvtBSemiTauonic2HDMType2 ()
	The destructor

std::string	getName ()
	The function returns the model name. More...

EvtDecayBase *	clone ()
	The function makes a copy of an EvtBSemiTauonic2HDMType2 object. More...

void	decay (EvtParticle *p)
	The function evaluates the decay amplitude of the parent particle. More...

void	initProbMax ()
	The function sets the maximum value of the probability. More...

void	init ()
	The function initializes the decay. More...

EvtSpinDensity	RotateToHelicityBasisInBoostedFrame (const EvtParticle *p, EvtVector4R p4boost)
	The function calculates the rotation matrix to convert the spin basis to the helicity basis in the boosted frame. More...

double	CalcMaxProb (EvtId parent, EvtId meson, EvtId lepton, EvtId nudaug, EvtBSemiTauonicHelicityAmplitudeCalculator *HelicityAmplitudeCalculator)
	The function calculates the maximum probability. More...

	EvtBSemiTauonicHelicityAmplitudeCalculator ()
	The default constructor. More...

	EvtBSemiTauonicHelicityAmplitudeCalculator (const double rho12, const double rhoA12, const double ffR11, const double ffR21, const double AS1, const double AR3, const double bottomMass, const double charmMass, const EvtComplex &CV1, const EvtComplex &CV2, const EvtComplex &CS1, const EvtComplex &CS2, const EvtComplex &CT, const double parentMass, const double DMass, const double DstarMass)
	The constructor with HQET form factor parameters, Wilson coefficients of new physics contributions and parent B, daughter D(*) meson masses. More...

EvtComplex	helAmp (double mtau, int tauhel, int Dhel, double w, double costau) const
	The function calculates the helicity amplitude. More...

EvtComplex	helAmp (const EvtComplex &CV1, const EvtComplex &CV2, const EvtComplex &CS1, const EvtComplex &CS2, const EvtComplex &CT, double mtau, int tauhel, int Dhel, double w, double costau) const
	The function calculates helicity amplitudes with given Wilson coefficients. More...

double	Lep (const double mtau, int tauhel, int whel, double q2, double costau) const
	The function to calculate the Leptonic Amplitudes for B->D*taunu decay of the vector type contribution. More...

double	Lep (const double mtau, int tauhel, double q2, double costau) const
	The function to calculate the Leptonic Amplitudes for B->Dtaunu decay of the scalar type contribution. More...

double	Lep (const double mtau, int tauhel, int whel1, int whel2, double q2, double costau) const
	The function to calculate the Leptonic Amplitudes for B->D*taunu decay of the tensor type contribution. More...

double	HadV1 (int Dhel, int whel, double w) const
	The function to calculate the Hadronic Amplitudes of left handed (V-A) type contribution. More...

double	HadV2 (int Dhel, int whel, double w) const
	The function to calculate the Hadronic Amplitudes of right handed (V+A) type contribution. More...

double	HadS1 (int Dhel, double w) const
	The function to calculate the Hadronic Amplitudes of scalar (S+P) type contribution. More...

double	HadS2 (int Dhel, double w) const
	The function to calculate the Hadronic Amplitudes of scalar (S-P) type contribution. More...

double	HadT (int Dhel, int whel1, int whel2, double w) const
	The function to calculate the Hadronic Amplitudes of tensor type contribution. More...

double	helampSM (double mtau, int tauhel, int Dhel, double w, double costau) const
	Helicity Amplitudes of SM (left handed) contribution. More...

double	helampV1 (double mtau, int tauhel, int Dhel, double w, double costau) const
	Helicity Amplitudes of left handed (V-A) contribution. More...

double	helampV2 (double mtau, int tauhel, int Dhel, double w, double costau) const
	Helicity Amplitudes of right handed (V+A) contribution. More...

double	helampS1 (double mtau, int tauhel, int Dhel, double w, double costau) const
	Helicity Amplitudes of scalar (S+P) type contribution. More...

double	helampS2 (double mtau, int tauhel, int Dhel, double w, double costau) const
	Helicity Amplitudes of scalar (S-P) type contribution. More...

double	helampT (double mtau, int tauhel, int Dhel, double w, double costau) const
	Helicity Amplitudes of tensor type contribution. More...

double	hp (double w) const
	HQET D vector form factor h_+(w). More...

double	hm (double w) const
	HQET D vector form factor h_-(w). More...

double	hA1 (double w) const
	HQET D* axial vector form factor h_{A1}(w). More...

double	hV (double w) const
	HQET D* axial vector form factor h_V(w). More...

double	hA2 (double w) const
	HQET D* axial vector form factor h_{A2}(w). More...

double	hA3 (double w) const
	HQET D* axial vector form factor h_{A3}(w). More...

double	hS (double w) const
	D scalar form factor h_S(w) in terms of vector form factors. More...

double	hP (double w) const
	D* pseudo scalar form factor h_P(w) in terms of axial vector form factors. More...

double	hT (double w) const
	D tensor form factor h_T(w) in terms of vector form factors. More...

double	hT1 (double w) const
	D* tensor form factor h_{T1}(w) in terms of axial vector form factors. More...

double	hT2 (double w) const
	D* tensor form factor h_{T2}(w). More...

double	hT3 (double w) const
	D* tensor form factor h_{T3}(w). More...

double	z (double w) const
	CLN form factor z. More...

double	ffV1 (double w) const
	CLN form factor V1. More...

double	ffS1 (double w) const
	CLN form factor S1. More...

double	ffA1 (double w) const
	CLN form factor A1. More...

double	ffR1 (double w) const
	CLN form factor R1. More...

double	ffR2 (double w) const
	CLN form factor R2. More...

double	ffR3 (double w) const
	CLN form factor R3. More...

double	dS1 (double w) const
	HQET correction factor for the scalar form factor for B->Dtaunu. More...

double	dR3 (double w) const
	HQET correction factor for the scalar form factor for B->D*taunu. More...

double	mD (int Dhel) const
	Daughter D(*) meson mass. More...

double	v (double mtau, double q2) const
	Function to calculate the tau velocity. More...

double	qh2 (int Dhel, double w) const
	Function to calculate the q^2 divided by the square of parent mass (m_B^2). More...

double	q2 (int Dhel, double w) const
	Function to calculate the q^2 of the decay (square of l+nu invariant mass). More...

bool	chkDhel (int Dhel) const
	sanity checkers More...

bool	chkwhel (int whel) const
	Function to check if whel is in the valid range. More...

bool	chktauhel (int tauhel) const
	Function to check if tauhel is in the valid range. More...

void	CalcAmp (EvtParticle parent, EvtAmp &amp, EvtBSemiTauonicHelicityAmplitudeCalculator HelicityAmplitudeCalculator) override
	The function calculates the spin dependent amplitude. More...

void	CalcAmp (EvtParticle parent, EvtAmp &amp, EvtBSemiTauonicHelicityAmplitudeCalculator HelicityAmplitudeCalculator) override
	The function calculates the spin dependent amplitude. More...

virtual	~EvtBtoXsnunu_FERMI ()
	Destructor.

std::string	getName ()
	The function which returns the name of the model.

EvtDecayBase *	clone ()
	The function which makes a copy of the model.

void	decay (EvtParticle *p)
	The function to determine kinematics of daughter particles based on dGdsb distribution.

void	initProbMax ()
	The function to sets a maximum probability. More...

void	init ()
	The function for an initialization.

double	dGdsbProb (double _sb)
	The function returns the probability density value depending on sb.

double	FermiMomentum (double pf)
	The function returns a momentum of b quark. More...

double	FermiMomentumProb (double pb, double pf)
	The function returns a probability based on the Fermi motion model. More...

virtual	~EvtEtaFullDalitz ()
	Default Destructor.

std::string	getName ()
	Returns the model name: ETA_FULLDALITZ.

EvtDecayBase *	clone ()
	Makes a copy of the pointer to the class.

void	init ()
	Checks that the number of input parameters are correct: More...

void	initProbMax ()
	Sets the Maximum probability for the PHSP reweight. More...

void	decay (EvtParticle *p)
	Function that implements the energy-dependent Dalitz.

virtual	~EvtEtaPi0Dalitz ()
	Default destructor.

std::string	getName ()
	Returns the name of the model: ETA_PI0DALITZ.

EvtDecayBase *	clone ()
	Makes a copy of the class object.

void	init ()
	Checks that the number of input parameters are correct: More...

void	initProbMax ()
	Sets the Maximum probability for the PHSP reweight. More...

void	decay (EvtParticle *p)
	Function that implements the energy-dependent Dalitz.

virtual	~EvtEtaPrimeDalitz ()
	Default destructor.

std::string	getName ()
	Returns the model name: ETAPRIME_DALITZ.

EvtDecayBase *	clone ()
	Returns a copy of the class object.

void	init ()
	Checks that the number of input parameters are correct: More...

void	initProbMax ()
	Sets the Maximum probability for the PHSP reweight. More...

void	decay (EvtParticle *p)
	Function that implements the energy-dependent Dalitz.

virtual	~EvtKnunu ()
	Destructor.

std::string	getName ()
	The function which returns the name of the model.

EvtDecayBase *	clone ()
	The function which makes a copy of the model.

void	decay (EvtParticle *p)
	The function to calculate a quark decay amplitude.

void	init ()
	The function for an initialization.

void	initProbMax ()
	The function to sets a maximum probability.

virtual	~EvtKstarnunu_REV ()
	Destructor.

std::string	getName ()
	The function which returns the name of the model.

EvtDecayBase *	clone ()
	The function which makes a copy of the model.

void	decay (EvtParticle *p)
	The function to calculate a quark decay amplitude.

void	init ()
	The function for an initialization.

void	initProbMax ()
	The function to sets a maximum probability.

std::string	getName ()
	Get function Name

EvtDecayBase *	clone ()
	Clone the decay.

void	init ()
	Initialize standard stream objects

void	initProbMax ()
	Initialize standard stream objects for probability function

void	decay (EvtParticle *p)
	Member of particle in EvtGen.

double	dGammadwdcostau (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, double mtau, int tauhel, int Dhel, double w, double costau)
	Function calculates the differential decay rate dGamma/dw/dcostau. More...

double	dGammadw (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, double mtau, int tauhel, int Dhel, double w)
	Function calculates the differential decay rate dGamma/dw, integrated for costau. More...

double	dGammadcostau (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, double mtau, int tauhel, int Dhel, double costau)
	Function calculates the differential decay rate dGamma/dcostau, integrated for w. More...

double	Gamma (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, double mtau, int tauhel, int Dhel)
	Function calculates the helicity dependent decay rate Gamma, integrated for w and costau. More...

double	GammaD (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, double mtau)
	Function calculates the decay rate Gamma for Dtaunu decay, integrated for w and costau and summed for helicities. More...

double	GammaDstar (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, double mtau)
	Function calculates the differential decay rate Gamma for D*taunu decay, integrated for w and costau and summed for helicities. More...

double	GammaSMD (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, const double mlep=0.0005110)
	Function calculates the SM decay rate Gamma for Dlnu decay, integrated for w and costau and summed for helicities. More...

double	GammaSMDstar (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, const double mlep=0.0005110)
	Function calculates the SM decay rate Gamma for D*lnu decay, integrated for w and costau and summed for helicities. More...

double	RGammaD (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, const double mtau, const double mlep=0.0005110)
	Function calculates the ratio of Br(B->Dtaunu)/Br(B->Dlnu), R(D). More...

double	RGammaDstar (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, const double mtau, const double mlep=0.0005110)
	Function calculates the ratio of Br(B->Dtaunu)/Br(B->Dlnu), R(D*). More...

double	PtauD (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, const double mtau)
	Function calculates the polarization of tau, (RH - LH)/(LH + RH), in B->Dtaunu decay. More...

double	PtauDstar (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, const double mtau)
	Function calculates the polarization of tau, (RH - LH)/(LH + RH), in B->D*taunu decay. More...

double	PDstar (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, const double mtau)
	Function calculates the polarization of D, longitudinal/(longitudinal + transverse), in B->Dtaunu decay. More...

double	pf (const EvtBSemiTauonicHelicityAmplitudeCalculator &BSTD, double mtau, int Dhel, double w)
	Phase space factor, which is multiplied to the helicity amplitude to calculate the decay rate. More...

double	EvaluateByW (double x, double param)
	Function used internally for numerical integration.

double	EvaluateByCostau (double x, double param)
	Function used internally for numerical integration.

double	EvaluateBy2D (double x, double param)
	Function used internally for numerical integration.

static EvtGenModelRegister &	getInstance ()
	Return reference to the instance. More...

static std::list< EvtDecayBase * >	getModels ()
	Return a list of models. More...

void	checkVariable (const std::string &name, const std::map< int, bool > &allowed, ParticleGun::EDistribution &dist, std::vector< double > &pars)
	check one of the variables given a list of allowed and excluded distributions

void	initialize ()
	function to be executed on initialize()

ROOT::Math::XYZVector	generateVertex (const ROOT::Math::XYZVector &initial, const TMatrixDSym &cov, MultivariateNormalGenerator &gen) const
	generate the vertex More...

TMatrixDSym	adjustCovMatrix (TMatrixDSym cov) const
	adjust smearing covariance matrix based on the generation flags

ConditionalGaussGenerator	initConditionalGenerator (const ROOT::Math::PxPyPzEVector &pHER, const ROOT::Math::PxPyPzEVector &pLER, const TMatrixDSym &covHER, const TMatrixDSym &covLER)
	Initialize the conditional generator using HER & LER 4-vectors and HER & LER covariance matrices describing spread.

ROOT::Math::PxPyPzEVector	generateBeam (const ROOT::Math::PxPyPzEVector &initial, const TMatrixDSym &cov, MultivariateNormalGenerator &gen) const
	generate 4 vector for one beam More...

MCInitialParticles &	generate ()
	Generate a new event.

MCInitialParticles &	generate (int allowedFlags)
	Generate a new event wit a particular set of allowed flags. More...

ROOT::Math::XYZVector	getVertexConditional ()
	Generate vertex position and possibly update the generator of Lorentz transformation.

ROOT::Math::XYZVector	updateVertex (bool force=false)
	Update the vertex position: More...

Variables
static const double	me = Const::electron.getMass()
	electron mass

static EvtGen *	s_evtgen = nullptr
	pointer to the evtgen instance

Detailed Description

Macro Definition Documentation

◆ B2_EVTGEN_REGISTER_MODEL

#define B2_EVTGEN_REGISTER_MODEL ( classname )

Value:

    namespace {\
    Belle2::EvtGenModelRegister::Factory<classname> EvtGenModelFactory_##classname; \
  }

Class to register B2_EVTGEN_REGISTER_MODEL.

Definition at line 75 of file EvtGenModelRegister.h.

Function Documentation

◆ CalcAmp() [1/2]

void CalcAmp	(	EvtParticle *	parent,
		EvtAmp &	amp,
		EvtBSemiTauonicHelicityAmplitudeCalculator *	HelicityAmplitudeCalculator
	)

overridevirtual

The function calculates the spin dependent amplitude.

Parameters

parent	a pointer to the parent particle.
amp	a pointer to fill the calculated spin dependent amplitude.
HelicityAmplitudeCalculator	a pointer to the calculator of the helicity dependent amplitude. The function calculate the spin dependent amplitude of the semi-tauonic decay to a scalar meson (D meson).

Implements EvtBSemiTauonicAmplitude.

Definition at line 25 of file EvtBSemiTauonicScalarMesonAmplitude.cc.

   {
  
     static EvtId EM = EvtPDL::getId("e-");
     static EvtId MUM = EvtPDL::getId("mu-");
     static EvtId TAUM = EvtPDL::getId("tau-");
 //  static EvtId EP = EvtPDL::getId("e+");
 //  static EvtId MUP = EvtPDL::getId("mu+");
 //  static EvtId TAUP = EvtPDL::getId("tau+");
  
     // calculate w and costau
  
     EvtVector4R p4d = p->getDaug(0)->getP4();
     EvtVector4R p4l = p->getDaug(1)->getP4();
     EvtVector4R p4n = p->getDaug(2)->getP4();
     EvtVector4R p4ln(p4l + p4n);
  
     EvtVector4R p4dln = boostTo(p4d, p4ln, true);
     EvtVector4R p4lln = boostTo(p4l, p4ln, true);
  
     const double gmB = p->getP4().mass();
     const double gmd = p4d.mass();
     const double gr = gmd / gmB;
  
     const double q2 = (p4l + p4n).mass2();
     const double w = (1. + gr * gr - q2 / gmB / gmB) / 2. / gr;
     // const double w=CalcHelAmp->wfunc(2,q2); avoid possible w<0 caused by the decay width
     const double costau = p4dln.dot(p4lln) / p4dln.d3mag() / p4lln.d3mag();
     const double ml = p4l.mass();
  
     // obtain helicity amplitudes
     EvtComplex helamp[2]; // tauhel={1,-1}
     helamp[0] = CalcHelAmp->helAmp(ml, 1, 2, w, costau); // note the parameter order is tauhel, Dhel
     helamp[1] = CalcHelAmp->helAmp(ml, -1, 2, w, costau); // Dhel=2 ==> D meson
  
     // lepton theta and phi in l+nu rest frame
     //const double l_theta=acos(p4lln.get(3)/p4lln.d3mag());
     //const double l_phi=atan2(p4lln.get(2),p4lln.get(1));
  
     // spin (in l rest frame) -> helicity (in l+nu rest frame) rotation matrix
     // ( sp0->hel0 , sp1->hel0 )
     // ( sp0->hel1 , sp1->hel1 )
     EvtSpinDensity l_HelFromSp = RotateToHelicityBasisInBoostedFrame(p->getDaug(1),
                                  p4ln);
 //                 l_phi,
 //                 l_theta,
 //                 -l_phi);
  
     // helicity (in l+nu rest frame) -> spin (in l rest frame) rotation matrix
     // ( hel0->sp0 , hel1->sp0 )
     // ( hel0->sp1 , hel1->sp1 )
     EvtComplex l_SpFromHel[2][2]; // {0,1} from {1,-1}
     EvtId l_num = p->getDaug(1)->getId();
 //  if (l_num == EM || l_num == MUM || l_num == TAUM) {
     l_SpFromHel[0][0] = conj(l_HelFromSp.get(0, 0));
     l_SpFromHel[0][1] = conj(l_HelFromSp.get(1, 0));
     l_SpFromHel[1][0] = conj(l_HelFromSp.get(0, 1));
     l_SpFromHel[1][1] = conj(l_HelFromSp.get(1, 1));
 //  } else {
 //    l_SpFromHel[0][1] = conj(l_HelFromSp.get(0, 0));
 //    l_SpFromHel[0][0] = conj(l_HelFromSp.get(1, 0));
 //    l_SpFromHel[1][1] = conj(l_HelFromSp.get(0, 1));
 //    l_SpFromHel[1][0] = conj(l_HelFromSp.get(1, 1));
 //  }
  
     // calculate spin amplitudes
     EvtComplex spinamp[2];
  
     for (int lsp = 0; lsp < 2; lsp++) {
       for (int lhel = 0; lhel < 2; lhel++) {
         // b -> l
         if (l_num == EM || l_num == MUM || l_num == TAUM) {
           spinamp[lsp] += l_SpFromHel[lsp][lhel] * helamp[lhel];
         }
         // b-bar -> anti-l
         else {
           spinamp[lsp] += l_SpFromHel[lsp][lhel] * (lhel == 0 ? +1 : -1) * conj(helamp[1 - lhel]);
         }
       }
     }
  
     amp.vertex(0, spinamp[0]);
     amp.vertex(1, spinamp[1]);
  
     // consistency check
     double helprob = abs2(helamp[0]) + abs2(helamp[1]);
     double spinprob = abs2(spinamp[0]) + abs2(spinamp[1]);
     if (fabs(helprob - spinprob) / helprob > 1E-6 || !finite(helprob) || !finite(spinprob)) {
       B2ERROR("EvtBSemiTauonicScalarMesonAmplitude total helicity prob does not match with total spin prob.");
       B2ERROR("helprob: " << helprob << " spinprob: " << spinprob);
       B2ERROR("w: " << w << " costau: " << costau << " hel probs: " << abs2(helamp[0])
               << "\t" << abs2(helamp[1])
               << "\t" << abs2(helamp[0]) + abs2(helamp[1]));
  
       B2ERROR("w: " << w << " costau: " << costau << " spin probs: " << abs2(spinamp[0])
               << "\t" << abs2(spinamp[1])
               << "\t" << abs2(spinamp[0]) + abs2(spinamp[1]));
  
 //    EvtGenReport(EVTGEN_ERROR, "EvtGen") <<
 //                                         "EvtBSemiTauonicScalarMesonAmplitude total helicity prob does not match with total spin prob."
 //                                         << std::endl;
 //    EvtGenReport(EVTGEN_ERROR, "EvtGen") << "helprob: "<<helprob<<" spinprob: "<<spinprob<<std::endl;
 //    EvtGenReport(EVTGEN_ERROR, "EvtGen") << "w: "<<w<<" costau: "<<costau
 //                                         <<" hel probs: "<<abs2(helamp[0])
 //                                         <<"\t"<<abs2(helamp[1])
 //                                   <<"\t"<<abs2(helamp[0])+abs2(helamp[1])<<std::endl;
 //
 //    EvtGenReport(EVTGEN_ERROR, "EvtGen") << "w: "<<w<<" costau: "<<costau
 //                                         <<" spin probs: "<<abs2(spinamp[0])
 //                                         <<"\t"<<abs2(spinamp[1])
 //                                   <<"\t"<<abs2(spinamp[0])+abs2(spinamp[1])<<std::endl;
       // abort();
     }
  
     return;
   }

◆ CalcAmp() [2/2]

void CalcAmp	(	EvtParticle *	parent,
		EvtAmp &	amp,
		EvtBSemiTauonicHelicityAmplitudeCalculator *	HelicityAmplitudeCalculator
	)

overridevirtual

The function calculates the spin dependent amplitude.

Parameters

parent	a pointer to the parent particle.
amp	a pointer to fill the calculated spin dependent amplitude.
HelicityAmplitudeCalculator	a pointer to the calculator of the helicity dependent amplitude. The function calculate the spin dependent amplitude of the semi-tauonic decay to a vector meson (D* meson).

Implements EvtBSemiTauonicAmplitude.

Definition at line 24 of file EvtBSemiTauonicVectorMesonAmplitude.cc.

   {
     static EvtId EM = EvtPDL::getId("e-");
     static EvtId MUM = EvtPDL::getId("mu-");
     static EvtId TAUM = EvtPDL::getId("tau-");
 //  static EvtId EP = EvtPDL::getId("e+");
 //  static EvtId MUP = EvtPDL::getId("mu+");
 //  static EvtId TAUP = EvtPDL::getId("tau+");
  
     // calculate w and costau
  
     EvtVector4R p4d = p->getDaug(0)->getP4();
     EvtVector4R p4l = p->getDaug(1)->getP4();
     EvtVector4R p4n = p->getDaug(2)->getP4();
     EvtVector4R p4ln(p4l + p4n);
  
     EvtVector4R p4dln = boostTo(p4d, p4ln, true);
     EvtVector4R p4lln = boostTo(p4l, p4ln, true);
  
     const double gmB = p->getP4().mass();
     const double gmd = p4d.mass();
     const double gr = gmd / gmB;
  
     const double q2 = (p4l + p4n).mass2();
     const double w = (1. + gr * gr - q2 / gmB / gmB) / 2. / gr;
     // const double w=CalcHelAmp->wfunc(1,q2); avoid w<1 caused by the D* width
     const double costau = p4dln.dot(p4lln) / p4dln.d3mag() / p4lln.d3mag();
     const double ml = p4l.mass();
  
     // Set D* mass from its momentum (to take into account fluctuation due to its width)
     const double orig_mD = CalcHelAmp->getMDst();
     CalcHelAmp->setMDst(gmd);
  
     // obtain helicity amplitudes
     EvtComplex helamp[3][2]; // Dhel={1,0,-1}, tauhel={1,-1}
     helamp[0][0] = CalcHelAmp->helAmp(ml, 1, 1, w, costau); // note the parameter order is tauhel, Dhel
     helamp[0][1] = CalcHelAmp->helAmp(ml, -1, 1, w, costau);
     helamp[1][0] = CalcHelAmp->helAmp(ml, 1, 0, w, costau);
     helamp[1][1] = CalcHelAmp->helAmp(ml, -1, 0, w, costau);
     helamp[2][0] = CalcHelAmp->helAmp(ml, 1, -1, w, costau);
     helamp[2][1] = CalcHelAmp->helAmp(ml, -1, -1, w, costau);
  
     // lepton theta and phi in l+nu rest frame
     // const double l_theta=acos(p4lln.get(3)/p4lln.d3mag());
     // const double l_phi=atan2(p4lln.get(2),p4lln.get(1));
  
     // spin (in l rest frame) -> helicity (in l+nu rest frame) rotation matrix
     // ( sp0->hel0 , sp1->hel0 )
     // ( sp0->hel1 , sp1->hel1 )
     EvtSpinDensity l_HelFromSp = RotateToHelicityBasisInBoostedFrame(p->getDaug(1),
                                  p4ln);
 //                 l_phi,
 //                 l_theta,
 //                 -l_phi);
  
     // helicity (in l+nu rest frame) -> spin (in l rest frame) rotation matrix
     // ( hel0->sp0 , hel1->sp0 )
     // ( hel0->sp1 , hel1->sp1 )
     EvtComplex l_SpFromHel[2][2]; // {0,1} from {1,-1}
     EvtId l_num = p->getDaug(1)->getId();
 //  if (l_num == EM || l_num == MUM || l_num == TAUM) {
     l_SpFromHel[0][0] = conj(l_HelFromSp.get(0, 0));
     l_SpFromHel[0][1] = conj(l_HelFromSp.get(1, 0));
     l_SpFromHel[1][0] = conj(l_HelFromSp.get(0, 1));
     l_SpFromHel[1][1] = conj(l_HelFromSp.get(1, 1));
 //  } else {
 //    l_SpFromHel[0][1] = conj(l_HelFromSp.get(0, 0));
 //    l_SpFromHel[0][0] = conj(l_HelFromSp.get(1, 0));
 //    l_SpFromHel[1][1] = conj(l_HelFromSp.get(0, 1));
 //    l_SpFromHel[1][0] = conj(l_HelFromSp.get(1, 1));
 //  }
  
     // meson spin state to helicity state
     const double D_theta = acos(p4d.get(3) / p4d.d3mag());
     const double D_phi = atan2(p4d.get(2), p4d.get(1));
     EvtSpinDensity D_HelFromSp = p->getDaug(0)->rotateToHelicityBasis(D_phi, D_theta, -D_phi);
  
     EvtComplex D_SpFromHel[3][3]; //  {0,1,2} from {1,0,-1}
     D_SpFromHel[0][0] = conj(D_HelFromSp.get(0, 0));
     D_SpFromHel[0][1] = conj(D_HelFromSp.get(1, 0));
     D_SpFromHel[0][2] = conj(D_HelFromSp.get(2, 0));
     D_SpFromHel[1][0] = conj(D_HelFromSp.get(0, 1));
     D_SpFromHel[1][1] = conj(D_HelFromSp.get(1, 1));
     D_SpFromHel[1][2] = conj(D_HelFromSp.get(2, 1));
     D_SpFromHel[2][0] = conj(D_HelFromSp.get(0, 2));
     D_SpFromHel[2][1] = conj(D_HelFromSp.get(1, 2));
     D_SpFromHel[2][2] = conj(D_HelFromSp.get(2, 2));
  
     // calculate spin amplitudes
     EvtComplex spinamp[3][2];
  
     for (int dsp = 0; dsp < 3; dsp++) {
       for (int lsp = 0; lsp < 2; lsp++) {
         for (int dhel = 0; dhel < 3; dhel++) {
           for (int lhel = 0; lhel < 2; lhel++) {
             // b -> l, D*+
             if (l_num == EM || l_num == MUM || l_num == TAUM) {
               spinamp[dsp][lsp] += l_SpFromHel[lsp][lhel] * D_SpFromHel[dsp][dhel] * helamp[dhel][lhel];
             }
             // b-bar -> anti-l, D*-
             else {
               spinamp[dsp][lsp] += l_SpFromHel[lsp][lhel] * D_SpFromHel[dsp][dhel] * (lhel == 0 ? +1 : -1) * (dhel == 1 ? +1 : -1)
                                    * conj(helamp[2 - dhel][1 - lhel]);
             }
           }
         }
       }
     }
  
     amp.vertex(0, 0, spinamp[0][0]);
     amp.vertex(0, 1, spinamp[0][1]);
     amp.vertex(1, 0, spinamp[1][0]);
     amp.vertex(1, 1, spinamp[1][1]);
     amp.vertex(2, 0, spinamp[2][0]);
     amp.vertex(2, 1, spinamp[2][1]);
  
     // Set D* mass to its original value
     CalcHelAmp->setMDst(orig_mD);
  
     // consistency check
     double helprob = abs2(helamp[0][0]) + abs2(helamp[0][1]) + abs2(helamp[1][0]) + abs2(helamp[1][1]) + abs2(helamp[2][0])
                      + abs2(helamp[2][1]);
     double spinprob = abs2(spinamp[0][0]) + abs2(spinamp[0][1]) + abs2(spinamp[1][0]) + abs2(spinamp[1][1])
                       + abs2(spinamp[2][0]) + abs2(spinamp[2][1]);
     if (fabs(helprob - spinprob) / helprob > 1E-6 || !finite(helprob) || !finite(spinprob)) {
       B2ERROR("EvtBSemiTauonicVectorMesonAmplitude total helicity prob does not match with total spin prob, or nan.");
       B2ERROR("helprob: " << helprob << " spinprob: " << spinprob);
       B2ERROR("w: " << w << " costau: " << costau << " hel probs: "
               << abs2(helamp[0][0]) << "\t" << abs2(helamp[0][1]) << "\t"
               << abs2(helamp[1][0]) << "\t" << abs2(helamp[1][1]) << "\t"
               << abs2(helamp[2][0]) << "\t" << abs2(helamp[2][1]) << "\t"
               << abs2(helamp[0][0]) + abs2(helamp[0][1]) + abs2(helamp[1][0]) + abs2(helamp[1][1]) + abs2(helamp[2][0]) + abs2(helamp[2][1])
              );
  
       B2ERROR("w: " << w << " costau: " << costau << " spin probs: "
               << abs2(spinamp[0][0]) << "\t" << abs2(spinamp[0][1]) << "\t"
               << abs2(spinamp[1][0]) << "\t" << abs2(spinamp[1][1]) << "\t"
               << abs2(spinamp[2][0]) << "\t" << abs2(spinamp[2][1]) << "\t"
               << abs2(spinamp[0][0]) + abs2(spinamp[0][1]) + abs2(spinamp[1][0]) + abs2(spinamp[1][1]) + abs2(spinamp[2][0]) + abs2(spinamp[2][1])
              );
  
 //    EvtGenReport(EVTGEN_ERROR, "EvtGen") <<
 //                                         "EvtBSemiTauonicVectorMesonAmplitude total helicity prob does not match with total spin prob, or nan."
 //                                         << std::endl;
 //    EvtGenReport(EVTGEN_ERROR, "EvtGen") "helprob: "<<helprob<<" spinprob: "<<spinprob<< std::endl;
 //    EvtGenReport(EVTGEN_ERROR, "EvtGen") "w: "<<w<<" costau: "<<costau<<" hel probs: "
 //            <<abs2(helamp[0][0])<<"\t"<<abs2(helamp[0][1])<<"\t"
 //            <<abs2(helamp[1][0])<<"\t"<<abs2(helamp[1][1])<<"\t"
 //            <<abs2(helamp[2][0])<<"\t"<<abs2(helamp[2][1])<<"\t"
 //            <<abs2(helamp[0][0]) + abs2(helamp[0][1]) + abs2(helamp[1][0]) + abs2(helamp[1][1]) + abs2(helamp[2][0]) + abs2(helamp[2][1])
 //      <<std::endl;
 //
 //    EvtGenReport(EVTGEN_ERROR, "EvtGen") "w: "<<w<<" costau: "<<costau<<" spin probs: "
 //            <<abs2(spinamp[0][0])<<"\t"<<abs2(spinamp[0][1])<<"\t"
 //            <<abs2(spinamp[1][0])<<"\t"<<abs2(spinamp[1][1])<<"\t"
 //            <<abs2(spinamp[2][0])<<"\t"<<abs2(spinamp[2][1])<<"\t"
 //            <<abs2(spinamp[0][0]) + abs2(spinamp[0][1]) + abs2(spinamp[1][0]) + abs2(spinamp[1][1]) + abs2(spinamp[2][0]) + abs2(spinamp[2][1])
 //      <<std::endl;
  
       // Debugging information
 //    fprintf(stderr, "q2 by Helamp: %g\n", CalcHelAmp->q2(1, w));
 //    fprintf(stderr, "helampSM by Helamp: %g\n", CalcHelAmp->helampSM(ml, 1, 1, w, costau));
 //    fprintf(stderr, "Lep by Helamp: %g\n", CalcHelAmp->Lep(ml, 1, 1, CalcHelAmp->q2(1, w), costau));
 //    fprintf(stderr, "HadV2 by Helamp: %g\n", CalcHelAmp->HadV2(1, 1, w));
 //    fprintf(stderr, "v by Helamp: %g\n", CalcHelAmp->v(ml, CalcHelAmp->v(ml, CalcHelAmp->q2(1, w))));
 //
 //    fprintf(stderr, "B mass : %g \t nominal %g\n", gmB, EvtPDL::getMeanMass(p->getId()));
 //    fprintf(stderr, "D mass : %g \t nominal %g\n", gmd, EvtPDL::getMeanMass(p->getDaug(0)->getId()));
 //    fprintf(stderr, "lepton mass : %g \t nominal %g\n", ml, EvtPDL::getMeanMass(p->getDaug(1)->getId()));
  
       // abort();
     }
  
     return;
   }

◆ CalcMaxProb()

double CalcMaxProb	(	EvtId	parent,
		EvtId	meson,
		EvtId	lepton,
		EvtId	nudaug,
		EvtBSemiTauonicHelicityAmplitudeCalculator *	HelicityAmplitudeCalculator
	)

The function calculates the maximum probability.

Parameters

parent	a ID of the parent meson.
meson	a ID of the daughter meson.
lepton	a ID of the daughter lepton.
nudaug	a ID of the daughter neutrino.
HelicityAmplitudeCalculator	a pointer to the calculator of the helicity dependent amplitude.

Returns: the maxmum probability multiplied by 1.1. The function scan the q^2 and angle between the daughter meson and lepton and search for the the maximum probability value.

Definition at line 77 of file EvtBSemiTauonicAmplitude.cc.

   {
     //This routine takes the arguements parent, meson, and lepton
     //number, and a form factor model, and returns a maximum
     //probability for this semileptonic form factor model.  A
     //brute force method is used.  The 2D cos theta lepton and
     //q2 phase space is probed.
  
     //Start by declaring a particle at rest.
  
     //It only makes sense to have a scalar parent.  For now.
     //This should be generalized later.
  
     EvtScalarParticle* scalar_part;
     EvtParticle* root_part;
  
     scalar_part = new EvtScalarParticle;
  
     //cludge to avoid generating random numbers!
     scalar_part->noLifeTime();
  
     EvtVector4R p_init;
  
     p_init.set(EvtPDL::getMass(parent), 0.0, 0.0, 0.0);
     scalar_part->init(parent, p_init);
     //root_part = (EvtParticle*)scalar_part;
     root_part = static_cast<EvtParticle*>(scalar_part);
  
     root_part->setDiagonalSpinDensity();
  
     EvtParticle* daughter, *lep, *trino;
  
     EvtAmp amp;
  
     EvtId listdaug[3];
     listdaug[0] = meson;
     listdaug[1] = lepton;
     listdaug[2] = nudaug;
  
     amp.init(parent, 3, listdaug);
  
     root_part->makeDaughters(3, listdaug);
     daughter = root_part->getDaug(0);
     lep = root_part->getDaug(1);
     trino = root_part->getDaug(2);
  
     //cludge to avoid generating random numbers!
     daughter->noLifeTime();
     lep->noLifeTime();
     trino->noLifeTime();
  
  
     //Initial particle is unpolarized, well it is a scalar so it is
     //trivial
     EvtSpinDensity rho;
     rho.setDiag(root_part->getSpinStates());
  
     double mass[3];
  
     double m = root_part->mass();
  
     EvtVector4R p4meson, p4lepton, p4nu, p4w;
  
     double q2, elepton, plepton;
     int i, j;
     double erho, prho, costl;
  
     double maxfoundprob = 0.0;
     double prob;
     int massiter;
  
     for (massiter = 0; massiter < 3; massiter++) {
  
       mass[0] = EvtPDL::getMeanMass(meson);
       mass[1] = EvtPDL::getMeanMass(lepton);
       mass[2] = EvtPDL::getMeanMass(nudaug);
       if (massiter == 1) {
         mass[0] = EvtPDL::getMinMass(meson);
       }
       if (massiter == 2) {
         mass[0] = EvtPDL::getMaxMass(meson);
         if ((mass[0] + mass[1] + mass[2]) > m) mass[0] = m - mass[1] - mass[2] - 0.00001;
       }
  
       double q2min = mass[1] * mass[1]; // limit to minimum=lepton mass square
       double q2max = (m - mass[0]) * (m - mass[0]);
       double dq2 = (q2max - q2min) / 25;
 //    std::cout<<"m: "<<m<<" mass[0]: "<<mass[0]<<" q2min: "<<q2min<<" q2max: "<<q2max<<std::endl;
  
       //loop over q2
  
       for (i = 0; i < 25; i++) {
         q2 = q2min + (i + 0.5) * dq2; // <-- !! not start from unphysical q2=0 !!
  
         erho = (m * m + mass[0] * mass[0] - q2) / (2.0 * m);
  
         prho = sqrt(erho * erho - mass[0] * mass[0]);
  
         p4meson.set(erho, 0.0, 0.0, -1.0 * prho);
         p4w.set(m - erho, 0.0, 0.0, prho);
  
 //      std::cout<<"q2: "<<q2<<std::endl;
 //      std::cout<<"p4meson: "<<p4meson<<std::endl;
  
         //This is in the W rest frame
         elepton = (q2 + mass[1] * mass[1]) / (2.0 * sqrt(q2));
         plepton = sqrt(elepton * elepton - mass[1] * mass[1]);
 //      std::cout<<"elepton: "<<elepton<<" plepton: "<<plepton<<std::endl;
  
         double probctl[3];
  
         for (j = 0; j < 3; j++) {
  
           costl = 0.99 * (j - 1.0);
  
           //These are in the W rest frame. Need to boost out into
           //the B frame.
           p4lepton.set(elepton, 0.0,
                        plepton * sqrt(1.0 - costl * costl), plepton * costl);
           p4nu.set(plepton, 0.0,
                    -1.0 * plepton * sqrt(1.0 - costl * costl), -1.0 * plepton * costl);
  
           EvtVector4R boost((m - erho), 0.0, 0.0, 1.0 * prho);
           p4lepton = boostTo(p4lepton, boost);
           p4nu = boostTo(p4nu, boost);
  
           //Now initialize the daughters...
  
           daughter->init(meson, p4meson);
           lep->init(lepton, p4lepton);
           trino->init(nudaug, p4nu);
  
           CalcAmp(root_part, amp, CalcHelAmp);
  
           //Now find the probability at this q2 and cos theta lepton point
           //and compare to maxfoundprob.
  
           //Do a little magic to get the probability!!
           prob = rho.normalizedProb(amp.getSpinDensity());
  
           probctl[j] = prob;
         }
  
         //probclt contains prob at ctl=-1,0,1.
         //prob=a+b*ctl+c*ctl^2
  
         double a = probctl[1];
         double b = 0.5 * (probctl[2] - probctl[0]);
         double c = 0.5 * (probctl[2] + probctl[0]) - probctl[1];
  
         prob = probctl[0];
         if (probctl[1] > prob) prob = probctl[1];
         if (probctl[2] > prob) prob = probctl[2];
  
         if (fabs(c) > 1e-20) {
           double ctlx = -0.5 * b / c;
           if (fabs(ctlx) < 1.0) {
             double probtmp = a + b * ctlx + c * ctlx * ctlx;
             if (probtmp > prob) prob = probtmp;
           }
  
         }
  
         if (prob > maxfoundprob) {
           maxfoundprob = prob;
         }
  
       }
       if (EvtPDL::getWidth(meson) <= 0.0) {
         //if the particle is narrow dont bother with changing the mass.
         massiter = 4;
       }
  
     }
     root_part->deleteTree();
  
     maxfoundprob *= 1.1;
     return maxfoundprob;
  
   }

◆ chkDhel()

bool chkDhel ( int Dhel ) const

private

sanity checkers

Function to check if Dhel is in the valid range.

Parameters

Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.

Definition at line 527 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

   {
     if (Dhel == -1 || Dhel == 0 || Dhel == 1 || Dhel == 2)return true;
     else return false;
   }

◆ chktauhel()

bool chktauhel ( int tauhel ) const

private

Function to check if tauhel is in the valid range.

Parameters

tauhel helicity of the lepton in the (l+nu) rest frame {+1,-1}.

Definition at line 539 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ chkwhel()

bool chkwhel ( int whel ) const

private

Function to check if whel is in the valid range.

Parameters

whel	helicity of the virtual vector boson {+1,0,1,2}.

Definition at line 533 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ clone() [1/2]

EvtDecayBase * clone ( )

The function makes a copy of an EvtBSTD object.

Returns: a pointer to a new uninitialized EvtBSemiTauonic object

Definition at line 47 of file EvtBSemiTauonic.cc.

   {
     return new EvtBSemiTauonic;
   }

◆ clone() [2/2]

EvtDecayBase * clone ( )

The function makes a copy of an EvtBSemiTauonic2HDMType2 object.

Returns: a pointer to a new uninitialized EvtBSemiTauonic2HDMType2 object

Definition at line 47 of file EvtBSemiTauonic2HDMType2.cc.

   {
     return new EvtBSemiTauonic2HDMType2;
   }

◆ decay() [1/2]

void decay ( EvtParticle * p )

The function evaluates the decay amplitude of the parent particle.

Parameters

p	a pointer to the parent particle

Definition at line 53 of file EvtBSemiTauonic.cc.

◆ decay() [2/2]

void decay ( EvtParticle * p )

The function evaluates the decay amplitude of the parent particle.

a pointer to the parent particle

Definition at line 53 of file EvtBSemiTauonic2HDMType2.cc.

◆ dGammadcostau()

double dGammadcostau	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		double	mtau,
		int	tauhel,
		int	Dhel,
		double	costau
	)

Function calculates the differential decay rate dGamma/dcostau, integrated for w.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame. Note that overall factor (GF/sqrt(2) Vcb)^2 is omitted.

Definition at line 59 of file EvtBSemiTauonicDecayRateCalculator.cc.

   {
     m_BSTD = &BSTD;
     TF1 f1("f1", this, &EvtBSemiTauonicDecayRateCalculator::EvaluateByW, BSTD.wmin(), BSTD.wmax(mtau, Dhel), 4,
            "EvtBSemiTauonicDecayRateCalculator", "EvaluateByW");
     f1.SetParameter(0, mtau);
     f1.SetParameter(1, tauhel);
     f1.SetParameter(2, Dhel);
     f1.SetParameter(3, costau);
     ROOT::Math::WrappedTF1 wf1(f1);
 //    ROOT::Math::GSLIntegrator ig;
     ROOT::Math::GaussIntegrator ig;
     ig.SetFunction(wf1);
     return ig.Integral(BSTD.wmin(), BSTD.wmax(mtau, Dhel));
   }

◆ dGammadw()

double dGammadw	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		double	mtau,
		int	tauhel,
		int	Dhel,
		double	w
	)

Function calculates the differential decay rate dGamma/dw, integrated for costau.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable. Note that overall factor (GF/sqrt(2) Vcb)^2 is omitted.

Definition at line 41 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ dGammadwdcostau()

double dGammadwdcostau	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		double	mtau,
		int	tauhel,
		int	Dhel,
		double	w,
		double	costau
	)

Function calculates the differential decay rate dGamma/dw/dcostau.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame. Note that overall factor (GF/sqrt(2) Vcb)^2 is omitted.

Definition at line 31 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ dR3()

double dR3 ( double w ) const

HQET correction factor for the scalar form factor for B->D*taunu.

Parameters

w	velocity transfer variable.

Returns: calculated factor.

Definition at line 493 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ dS1()

double dS1 ( double w ) const

HQET correction factor for the scalar form factor for B->Dtaunu.

Parameters

w	velocity transfer variable.

Returns: calculated factor.

Definition at line 488 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ EvtBSemiTauonicHelicityAmplitudeCalculator() [1/2]

EvtBSemiTauonicHelicityAmplitudeCalculator ( )

The default constructor.

Initializes with the default parameter values used by the aurhos of PRD87,034028.

Definition at line 21 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ EvtBSemiTauonicHelicityAmplitudeCalculator() [2/2]

EvtBSemiTauonicHelicityAmplitudeCalculator	(	const double	rho12,
		const double	rhoA12,
		const double	ffR11,
		const double	ffR21,
		const double	AS1,
		const double	AR3,
		const double	bottomMass,
		const double	charmMass,
		const EvtComplex &	CV1,
		const EvtComplex &	CV2,
		const EvtComplex &	CS1,
		const EvtComplex &	CS2,
		const EvtComplex &	CT,
		const double	parentMass,
		const double	DMass,
		const double	DstarMass
	)

The constructor with HQET form factor parameters, Wilson coefficients of new physics contributions and parent B, daughter D(*) meson masses.

Parameters

rho12	HQET form factor parameter rho_1^2 obtained by Dlnu decay data.
rhoA12	HQET form factor parameter rho_A1^2 obtained by D*lnu decay data.
ffR11	HQET form factor parameter R_1(1) obtained by D*lnu decay data.
ffR21	HQET form factor parameter R_2(1) obtained by D*lnu decay data.
AS1	a parameter to take into account the theoretical error of the scalar form factor for Dtaunu decay.
AR3	a parameter to take into account the theoretical error of the scalar form factor for D*taunu decay.
CV1	Wilson coeffcient of the left handed vector type NP contribution.
CV2	Wilson coeffcient of the right handed vector type NP contribution.
CS1	Wilson coeffcient of the scalar (S+P) type NP contribution.
CS2	Wilson coeffcient of the scalar (S-P) type NP contribution.
CT	Wilson coeffcient of the tensor type NP contribution.
parentMass	mass of the parent (B) meson.
DMass	mass of the scalar type daughter (D) meson.
DstarMass	mass of the vector type daughter (D*) meson.
bottomMass	mass of the bottom quark mass (running mass at the energy of bottom quark mass)
charmMass	mass of the charm quark mass (running mass at the energy of bottom quark mass) The constructor initializes the parameters of the decay. The recommended values of AS1 and AR3 by authors of PRD87,034028 are 1+/-1.

Definition at line 59 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ FermiMomentum()

double FermiMomentum ( double pf )

The function returns a momentum of b quark.

The distribution of the momentum is based on the Fermi motion model.

Definition at line 310 of file EvtBtoXsnunu_FERMI.cc.

   {
     // Pick a value for the b-quark Fermi motion momentum
     // according to Ali's Gaussian model
  
     // reference: Ali, Ahmed, et al. "Power corrections in the decay rate and distributions in B->Xs l+l- 2 in the standard model"
     // see (eq.57)
  
     double pb, pbmax;
     pb = 0.0;
     pbmax = 5.0 * pf;
  
     while (pb == 0.0) {
       double xbox = EvtRandom::Flat(pbmax);
       double ybox = EvtRandom::Flat();
       if (ybox < FermiMomentumProb(xbox, pf)) { pb = xbox; }
     }
  
     return pb;
   }

◆ FermiMomentumProb()

double FermiMomentumProb	(	double	pb,
		double	pf
	)

The function returns a probability based on the Fermi motion model.

reference: Ali, Ahmed, et al. "Power corrections in the decay rate and distributions in B->Xs l+l- 2 in the standard model" see (eq.57)

Parameters

pb	momentum of the b quark
pf	momentum width parameter in the Fermi motion model

Definition at line 331 of file EvtBtoXsnunu_FERMI.cc.

◆ ffA1()

double ffA1 ( double w ) const

CLN form factor A1.

Parameters

w	velocity transfer variable.

Returns: calculated form factor value.

Definition at line 466 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ ffR1()

double ffR1 ( double w ) const

CLN form factor R1.

Parameters

w	velocity transfer variable.

Returns: calculated form factor value.

Definition at line 473 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ ffR2()

double ffR2 ( double w ) const

CLN form factor R2.

Parameters

w	velocity transfer variable.

Returns: calculated form factor value.

Definition at line 478 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ ffR3()

double ffR3 ( double w ) const

CLN form factor R3.

Parameters

w	velocity transfer variable.

Returns: calculated form factor value.

Definition at line 483 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ ffS1()

double ffS1 ( double w ) const

CLN form factor S1.

Parameters

w	velocity transfer variable.

Returns: calculated form factor value.

Definition at line 461 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ ffV1()

double ffV1 ( double w ) const

CLN form factor V1.

Parameters

w	velocity transfer variable.

Returns: calculated form factor value.

Definition at line 454 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ Gamma()

double Gamma	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		double	mtau,
		int	tauhel,
		int	Dhel
	)

Function calculates the helicity dependent decay rate Gamma, integrated for w and costau.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D. Note that overall factor (GF/sqrt(2) Vcb)^2 is omitted.

Definition at line 77 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ GammaD()

double GammaD	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		double	mtau
	)

Function calculates the decay rate Gamma for Dtaunu decay, integrated for w and costau and summed for helicities.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter lepton mass. Note that overall factor (GF/sqrt(2) Vcb)^2 is omitted.

Definition at line 98 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ GammaDstar()

double GammaDstar	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		double	mtau
	)

Function calculates the differential decay rate Gamma for D*taunu decay, integrated for w and costau and summed for helicities.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter lepton mass. Note that overall factor (GF/sqrt(2) Vcb)^2 is omitted.

Definition at line 106 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ GammaSMD()

double GammaSMD	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		const double	mlep = `0.0005110`
	)

Function calculates the SM decay rate Gamma for Dlnu decay, integrated for w and costau and summed for helicities.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mlep	daughter lepton mass (default is the electron mass). Note that overall factor (GF/sqrt(2) Vcb)^2 is omitted.

Definition at line 117 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ GammaSMDstar()

double GammaSMDstar	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		const double	mlep = `0.0005110`
	)

Function calculates the SM decay rate Gamma for D*lnu decay, integrated for w and costau and summed for helicities.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mlep	daughter lepton mass (default is the electron mass). Note that overall factor (GF/sqrt(2) Vcb)^2 is omitted.

Definition at line 132 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ generate()

MCInitialParticles & generate ( int allowedFlags )

private

Generate a new event wit a particular set of allowed flags.

Parameters

[in] allowedFlags Allowed flags.

Definition at line 145 of file InitialParticleGeneration.cc.

   {
     if (!m_event) {
       m_event.create();
     }
     if (!m_beamParams.isValid()) {
       B2FATAL("Cannot generate beam without valid BeamParameters");
     }
     if (m_beamParams.hasChanged()) {
       m_generateHER.reset();
       m_generateLER.reset();
       m_generateVertex.reset();
     }
     m_event->setGenerationFlags(m_beamParams->getGenerationFlags() & allowedFlags);
     ROOT::Math::PxPyPzEVector her = generateBeam(m_beamParams->getHER(), m_beamParams->getCovHER(), m_generateHER);
     ROOT::Math::PxPyPzEVector ler = generateBeam(m_beamParams->getLER(), m_beamParams->getCovLER(), m_generateLER);
     ROOT::Math::XYZVector vtx = generateVertex(m_beamParams->getVertex(), m_beamParams->getCovVertex(), m_generateVertex);
     m_event->set(her, ler, vtx);
     //Check if we want to go to CMS, if so boost both
     if (m_beamParams->hasGenerationFlags(BeamParameters::c_generateCMS)) {
       m_event->setGenerationFlags(0);
       her = m_event->getLabToCMS() * her;
       ler = m_event->getLabToCMS() * ler;
       m_event->set(her, ler, vtx);
       m_event->setGenerationFlags(m_beamParams->getGenerationFlags() & allowedFlags);
     }
     return *m_event;
   }

◆ generateBeam()

ROOT::Math::PxPyPzEVector generateBeam	(	const ROOT::Math::PxPyPzEVector &	initial,
		const TMatrixDSym &	cov,
		MultivariateNormalGenerator &	gen
	)		const

private

generate 4 vector for one beam

Parameters

initial	beam
cov	covariance of the beam momentum (E, theta_x, theta_y)
gen	multivariate normal generator to be used

Definition at line 101 of file InitialParticleGeneration.cc.

◆ generateEvent()

void generateEvent ( MCParticleGraph & mcGraph )

Generates one single event.

Parameters

mcGraph Reference to the MonteCarlo graph into which the generated particles will be stored.

Definition at line 105 of file CRY.cc.

   {
     bool eventInAcceptance = 0;
     static std::vector<CRYParticle*> ev;
  
     //force at least particle in acceptance box
     for (int iTrial = 0; iTrial < m_maxTrials; ++iTrial) {
       m_totalTrials++;
       // Generate one event
       ev.clear();
       m_cryGenerator->genEvent(&ev);
       // check all particles
       for (auto* p : ev) {
         const int pdg = p->PDGid();
         const double kineticEnergy = p->ke() * Unit::MeV;
         if (kineticEnergy < m_kineticEnergyThreshold) continue;
  
         const double mass = TDatabasePDG::Instance()->GetParticle(pdg)->Mass();
         const double etot = kineticEnergy + mass;
         // since etot is at least mass this cannot be negative
         const double ptot = sqrt(etot * etot - mass * mass);
  
         // Momentum
         // We have x horizontal,  y up and z along the beam. So uvw -> zxy, xc yc zc -> zxy
         const double px = ptot * p->v();
         const double py = ptot * p->w();
         const double pz = ptot * p->u();
         // Vertex
         const double vx = p->y() * Unit::m;
         const double vy = p->z() * Unit::m;
         const double vz = p->x() * Unit::m;
  
         // Time
         /* In basf2, it is assumed that t = 0 when an event was produced,
            For the cosmic case, we set t = 0 when particle cross y=0 plane;
            The output time from CRY (p->t()) is too large (order of second) and  also
         increase as simulated time, so it is impossible to handle in basf2.
         For event which has more then one particle, the difference between their production
         times is also too large (> micro-second) to keep in basf2, so the time relation
         between particles in each event is also reset. Production of each particle in event is set t=0 at y=0.
         if one need production from CRY for a special study, you have to find a way to handle it...*/
         double time = 0;
  
         vecgeom::Vector3D<double> pos(vx, vy, vz);
         vecgeom::LorentzVector<double> mom(px, py, pz, etot);
         const double speed = (mass == 0 ? 1 : mom.Beta()) * Const::speedOfLight;
  
         // Project on the boundary of the world box
         auto inside = m_world->Inside(pos);
         if (inside == vecgeom::kInside) {
           // inside the world volume, go backwards in time to the world box
           const auto dir = -mom.vect().Unit();
           double dist = m_world->DistanceToOut(pos, dir);
           pos += dist * dir;
           time -= dist / speed;
         } else if (inside == vecgeom::kOutside) {
           // outside the world volume, go forwards in time to the world box
           // this should not happen but better safe then sorry ...
           const auto dir = mom.vect().Unit();
           double dist = m_world->DistanceToIn(pos, dir);
           if (dist == vecgeom::InfinityLength<double>()) continue;
           pos += dist * dir;
           time += dist / speed;
         }
         // Intersect with the acceptance box
         double dist = m_acceptance->DistanceToIn(pos, mom.vect().Unit());
         if (dist == vecgeom::InfinityLength<double>()) continue;
  
         // We want to keep this one
         auto& particle = mcGraph.addParticle();
         // all particle of a generator are primary
         particle.addStatus(MCParticle::c_PrimaryParticle);
         // all particle of CRY are stable
         particle.addStatus(MCParticle::c_StableInGenerator);
         particle.setPDG(pdg);
         particle.setFirstDaughter(0);
         particle.setLastDaughter(0);
         particle.setMomentum(ROOT::Math::XYZVector(mom.x(), mom.y(), mom.z()));
         particle.setMass(mass);
         particle.setEnergy(mom.e());
         particle.setProductionVertex(ROOT::Math::XYZVector(pos.x(), pos.y(), pos.z()));
         particle.setProductionTime(time);
         eventInAcceptance = true;
       }
       // clean up CRY event
       for (auto* p : ev) delete p;
       // and if we have something in the acceptance then we're done
       if (eventInAcceptance) {
         return;
       }
  
     }
     B2ERROR("Number of trials exceeds maxTrials increase number of maxTrials" << LogVar("maxTrials", m_maxTrials));
   }

◆ generateVertex()

ROOT::Math::XYZVector generateVertex	(	const ROOT::Math::XYZVector &	initial,
		const TMatrixDSym &	cov,
		MultivariateNormalGenerator &	gen
	)		const

private

generate the vertex

Parameters

initial	nominal vertex position
cov	covariance of the vertex position
gen	multivariate normal generator to be used

Definition at line 27 of file InitialParticleGeneration.cc.

◆ getInstance()

EvtGenModelRegister & getInstance ( )

staticprivate

Return reference to the instance.

This class behaves like a purely static class but we need a singleton pattern to avoid initialisation hell.

Definition at line 21 of file EvtGenModelRegister.cc.

   {
     static unique_ptr<EvtGenModelRegister> instance(new EvtGenModelRegister());
     return *instance;
   }

◆ getModels()

list< EvtDecayBase * > getModels ( )

static

Return a list of models.

Caller takes responsibility to free the instances when no longer needed

Returns: list with pointers to instances of all registered models.

Definition at line 27 of file EvtGenModelRegister.cc.

◆ getName() [1/2]

std::string getName ( )

The function returns the model name.

Returns: name of the model "BSTD"

Definition at line 40 of file EvtBSemiTauonic.cc.

◆ getName() [2/2]

std::string getName ( )

The function returns the model name.

Returns: name of the model "BSTD_2HDMTYPE2"

Definition at line 40 of file EvtBSemiTauonic2HDMType2.cc.

◆ getscalarff()

void getscalarff	(	EvtId	parent,
		EvtId	daughter,
		double	t,
		double	,
		double *	fpf,
		double *	f0f
	)

Scalar FF's.

Implementation follows arXiv:1509.06938v3.

For the scalar FF, the arglist in the constructor should contain 8 expansion parameters: b+_0, b+_1, b+_2, b+_3, b0_0, b0_1, b0_2, b0_3

Parameters

parent
daughter
t	Momentum transfer, also called q2. q2 = (p_B - p_M)^2
fpf	f_+(q2).
f0f	f_0(q2).

Definition at line 30 of file EvtBCLFF.cc.

   {
  
     if (m_numBCLFFCoefficients != 8) {
       EvtGenReport(EVTGEN_ERROR, "EvtGen") << "Wrong number of arguments for EvtBCLFF::getscalarff!\n";
     }
  
     const auto mB = EvtPDL::getMeanMass(parent);
     const auto mM = EvtPDL::getMeanMass(daughter);
  
     const auto tplus = (mB + mM) * (mB + mM);
     const auto tzero = (mB + mM) * (std::sqrt(mB) - std::sqrt(mM)) * (std::sqrt(mB) - std::sqrt(mM));
  
     const auto mR2 = m_resonance1Minus * m_resonance1Minus;
     const auto pole = 1 / (1 - t / mR2);
  
     const std::vector<double> bplus = {m_BCLFFCoefficients[0], m_BCLFFCoefficients[1], m_BCLFFCoefficients[2], m_BCLFFCoefficients[3]};
     const std::vector<double> bzero = {m_BCLFFCoefficients[4], m_BCLFFCoefficients[5], m_BCLFFCoefficients[6], m_BCLFFCoefficients[7]};
  
     const auto N_fpf = bplus.size();
     const auto N_f0f = bzero.size();
  
     auto z = [tplus, tzero](decltype(t) q2) {
       const auto term1 = std::sqrt(tplus - q2);
       const auto term2 = std::sqrt(tplus - tzero);
       return (term1 - term2) / (term1 + term2);
     };
  
     double sum_fpf = 0;
     for (unsigned int n = 0; n < N_fpf; ++n) {
       sum_fpf += bplus[n] * (std::pow(z(t), n) - std::pow(-1, n - N_fpf) * n / N_fpf * std::pow(z(t), N_fpf));
     }
     *fpf = pole * sum_fpf;
  
     double sum_f0f = 0;
     for (unsigned int n = 0; n < N_f0f; ++n) {
       sum_f0f += bzero[n] * std::pow(z(t), n);
     }
     *f0f = sum_f0f;
   }

◆ getTotalEnergy()

double Belle2::getTotalEnergy ( const std::vector< ROOT::Math::PxPyPzMVector > & ps )

inline

Main function scaleParticleEnergies in this header scales momenta in CMS of all particles in the final state by a constant factor such that the overall collision energy is slightly changed.

It is used for MC generators which do not provide option to generate events with smeared energy of initial particles. There is an assumption that matrix element of the process does not change much when Ecms is varied by 5MeV/10580MeV = 0.05% (typical Ecms smearing at Belle II). This is typically the case if the cross section does not have the resonance peak around collision energy. Get total energy of all particles in the provided vector

Definition at line 34 of file scaleParticleEnergies.h.

   {
     double eTot = 0;
     for (auto p : ps)
       eTot += p.E();
     return eTot;
   }

◆ getvectorff()

void getvectorff	(	EvtId	parent,
		EvtId	daughter,
		double	t,
		double	,
		double *	a1f,
		double *	a2f,
		double *	vf,
		double *	a0f
	)

Vector FF's.

Implementation follows arXiv:1503.05534v3. It is assumed that each expansion has three terms (hardcoded). However, this can be easily expanded or generalized. It is not done, because this way we can check if the number of arguments in the decay file is the correct one.

For the vector FF, the arglist in the constructor should contain 11 expansion parameters: A0_1, A0_2, A1_0, A1_1, A1_2, A12_0, A12_1, A12_2, V_0, V_1, V_2 Nota bene: A0_0 is correlated to A12_0.

Parameters

parent
daughter
t	Momentum transfer, also called q2. q2 = (p_B - p_M)^2
a1f	A1(q2)
a2f	A2(q2)
vf	V(q2)
a0f	A0(q2)

Definition at line 71 of file EvtBCLFF.cc.

◆ hA1()

double hA1 ( double w ) const

HQET D* axial vector form factor h_{A1}(w).

Definition at line 367 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hA2()

double hA2 ( double w ) const

HQET D* axial vector form factor h_{A2}(w).

Definition at line 379 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hA3()

double hA3 ( double w ) const

HQET D* axial vector form factor h_{A3}(w).

Definition at line 385 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ HadS1()

double HadS1	(	int	Dhel,
		double	w
	)		const

The function to calculate the Hadronic Amplitudes of scalar (S+P) type contribution.

Parameters

Dhel	helicity of the daughter D() meson {+1,0,1} for D and 2 for D.
w	velocity transfer variable.

Returns: calculated amplitude value.

Definition at line 242 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ HadS2()

double HadS2	(	int	Dhel,
		double	w
	)		const

The function to calculate the Hadronic Amplitudes of scalar (S-P) type contribution.

Parameters

Dhel	helicity of the daughter D() meson {+1,0,1} for D and 2 for D.
w	velocity transfer variable.

Returns: calculated amplitude value.

Definition at line 254 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ HadT()

double HadT	(	int	Dhel,
		int	whel1,
		int	whel2,
		double	w
	)		const

The function to calculate the Hadronic Amplitudes of tensor type contribution.

Parameters

Dhel	helicity of the daughter D() meson {+1,0,1} for D and 2 for D.
whel1	helicity of the one virtual vector boson {+1,0,1,2}.
whel2	helicity of the another virtual vector boson {+1,0,1,2}.
w	velocity transfer variable.

Returns: calculated amplitude value. The tensor contribution is described by two vector contributions.

Definition at line 264 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ HadV1()

double HadV1	(	int	Dhel,
		int	whel,
		double	w
	)		const

The function to calculate the Hadronic Amplitudes of left handed (V-A) type contribution.

Parameters

Dhel	helicity of the daughter D() meson {+1,0,1} for D and 2 for D.
whel	helicity of the virtual vector boson {+1,0,1,2}.
w	velocity transfer variable.

Returns: calculated amplitude value.

Definition at line 197 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ HadV2()

double HadV2	(	int	Dhel,
		int	whel,
		double	w
	)		const

The function to calculate the Hadronic Amplitudes of right handed (V+A) type contribution.

Parameters

Dhel	helicity of the daughter D() meson {+1,0,1} for D and 2 for D.
whel	helicity of the virtual vector boson {+1,0,1,2}.
w	velocity transfer variable.

Returns: calculated amplitude value.

Definition at line 229 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ helAmp() [1/2]

EvtComplex helAmp	(	const EvtComplex &	CV1,
		const EvtComplex &	CV2,
		const EvtComplex &	CS1,
		const EvtComplex &	CS2,
		const EvtComplex &	CT,
		double	mtau,
		int	tauhel,
		int	Dhel,
		double	w,
		double	costau
	)		const

The function calculates helicity amplitudes with given Wilson coefficients.

Parameters

CV1	Wilson coeffcient of the left handed vector type NP contribution.
CV2	Wilson coeffcient of the right handed vector type NP contribution.
CS1	Wilson coeffcient of the scalar (S+P) type NP contribution.
CS2	Wilson coeffcient of the scalar (S-P) type NP contribution.
CT	Wilson coeffcient of the tensor type NP contribution.
mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame. The overall factor GF/sqrt(2) Vcb omitted because it does not change the distribution.

Definition at line 108 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ helAmp() [2/2]

EvtComplex helAmp	(	double	mtau,
		int	tauhel,
		int	Dhel,
		double	w,
		double	costau
	)		const

The function calculates the helicity amplitude.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame. The overall factor GF/sqrt(2) Vcb omitted because it does not change the distribution.

Definition at line 101 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ helampS1()

double helampS1	(	double	mtau,
		int	tauhel,
		int	Dhel,
		double	w,
		double	costau
	)		const

Helicity Amplitudes of scalar (S+P) type contribution.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame.

Returns: calculated amplitude value. Overall factor GF/sqrt(2) Vcb omitted. Wilson coefficients CXX ommited.

Definition at line 325 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ helampS2()

double helampS2	(	double	mtau,
		int	tauhel,
		int	Dhel,
		double	w,
		double	costau
	)		const

Helicity Amplitudes of scalar (S-P) type contribution.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame.

Returns: calculated amplitude value. Overall factor GF/sqrt(2) Vcb omitted. Wilson coefficients CXX ommited.

Definition at line 331 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ helampSM()

double helampSM	(	double	mtau,
		int	tauhel,
		int	Dhel,
		double	w,
		double	costau
	)		const

Helicity Amplitudes of SM (left handed) contribution.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame. Overall factor GF/sqrt(2) Vcb omitted. Wilson coefficients CXX ommited.

Definition at line 297 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ helampT()

double helampT	(	double	mtau,
		int	tauhel,
		int	Dhel,
		double	w,
		double	costau
	)		const

Helicity Amplitudes of tensor type contribution.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame.

Returns: calculated amplitude value. Overall factor GF/sqrt(2) Vcb omitted. Wilson coefficients CXX ommited.

Definition at line 337 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ helampV1()

double helampV1	(	double	mtau,
		int	tauhel,
		int	Dhel,
		double	w,
		double	costau
	)		const

Helicity Amplitudes of left handed (V-A) contribution.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame.

Returns: calculated amplitude value. Overall factor GF/sqrt(2) Vcb omitted. Wilson coefficients CXX ommited.

Definition at line 308 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ helampV2()

double helampV2	(	double	mtau,
		int	tauhel,
		int	Dhel,
		double	w,
		double	costau
	)		const

Helicity Amplitudes of right handed (V+A) contribution.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame.

Returns: calculated amplitude value. Overall factor GF/sqrt(2) Vcb omitted. Wilson coefficients CXX ommited.

Definition at line 314 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hm()

double hm ( double w ) const

HQET D vector form factor h_-(w).

Definition at line 361 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hp()

double hp ( double w ) const

HQET D vector form factor h_+(w).

Definition at line 354 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hP()

double hP ( double w ) const

D* pseudo scalar form factor h_P(w) in terms of axial vector form factors.

Definition at line 404 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hS()

double hS ( double w ) const

D scalar form factor h_S(w) in terms of vector form factors.

Definition at line 392 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hT()

double hT ( double w ) const

D tensor form factor h_T(w) in terms of vector form factors.

Definition at line 416 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hT1()

double hT1 ( double w ) const

D* tensor form factor h_{T1}(w) in terms of axial vector form factors.

Definition at line 423 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hT2()

double hT2 ( double w ) const

D* tensor form factor h_{T2}(w).

Definition at line 432 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hT3()

double hT3 ( double w ) const

D* tensor form factor h_{T3}(w).

Definition at line 440 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ hV()

double hV ( double w ) const

HQET D* axial vector form factor h_V(w).

Definition at line 373 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ init() [1/5]

void init ( )

The function initializes the decay.

The function reads paramters from the decay file and initializes the parameters for the decay amplitude calculator.

Definition at line 78 of file EvtBSemiTauonic.cc.

◆ init() [2/5]

void init ( )

The function initializes the decay.

The function reads paramters from the decay file and initializes the parameters for the decay amplitude calculator.

Definition at line 78 of file EvtBSemiTauonic2HDMType2.cc.

◆ init() [3/5]

void init ( )

Checks that the number of input parameters are correct:

3 scalar particles (check on spin but not on charge)
6 real parameters (no check on the actual values)

Definition at line 41 of file EvtEtaFullDalitz.cc.

   {
  
     // check that there is are arguments
     checkNArg(6);
     checkNDaug(3);
  
  
     checkSpinParent(EvtSpinType::SCALAR);
  
     checkSpinDaughter(0, EvtSpinType::SCALAR);
     checkSpinDaughter(1, EvtSpinType::SCALAR);
     checkSpinDaughter(2, EvtSpinType::SCALAR);
   }

◆ init() [4/5]

void init ( )

Checks that the number of input parameters are correct:

3 scalar particles (check on spin but not on charge)
1 real parameters (no check on the actual values)

Definition at line 56 of file EvtEtaPi0Dalitz.cc.

   {
  
     // check that there is 1 argument
     checkNArg(1);
     checkNDaug(3);
  
  
     checkSpinParent(EvtSpinType::SCALAR);
  
     checkSpinDaughter(0, EvtSpinType::SCALAR);
     checkSpinDaughter(1, EvtSpinType::SCALAR);
     checkSpinDaughter(2, EvtSpinType::SCALAR);
   }

◆ init() [5/5]

void init ( )

Checks that the number of input parameters are correct:

3 scalar particles (check on spin but not on charge)
4 real parameters (no check on the actual values)

Definition at line 58 of file EvtEtaPrimeDalitz.cc.

   {
  
     // check that there is are arguments
     checkNArg(4);
     checkNDaug(3);
  
  
     checkSpinParent(EvtSpinType::SCALAR);
  
     checkSpinDaughter(0, EvtSpinType::SCALAR);
     checkSpinDaughter(1, EvtSpinType::SCALAR);
     checkSpinDaughter(2, EvtSpinType::SCALAR);
   }

◆ initProbMax() [1/5]

void initProbMax ( )

The function sets the maximum value of the probability.

Definition at line 59 of file EvtBSemiTauonic2HDMType2.cc.

◆ initProbMax() [2/5]

void initProbMax ( )

The function to sets a maximum probability.

In this model, noProbMax() is called because maximum probability is determined by m_dGdsbProbMax

Definition at line 161 of file EvtBtoXsnunu_FERMI.cc.

◆ initProbMax() [3/5]

void initProbMax ( )

Sets the Maximum probability for the PHSP reweight.

Maximum value is hardcoded and inherited from the EvtEtaDalitz class.

Definition at line 57 of file EvtEtaFullDalitz.cc.

◆ initProbMax() [4/5]

void initProbMax ( )

Sets the Maximum probability for the PHSP reweight.

Maximum value is hardcoded and inherited from the EvtEtaDalitz class.

Definition at line 72 of file EvtEtaPi0Dalitz.cc.

◆ initProbMax() [5/5]

void initProbMax ( )

Sets the Maximum probability for the PHSP reweight.

Maximum value is hardcoded and inherited from the EvtEtaDalitz class.

Definition at line 74 of file EvtEtaPrimeDalitz.cc.

◆ Lep() [1/3]

double Lep	(	const double	mtau,
		int	tauhel,
		double	q2,
		double	costau
	)		const

The function to calculate the Leptonic Amplitudes for B->Dtaunu decay of the scalar type contribution.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
q2	q^2 of the decay (square of l+nu invariant mass).
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame.

Returns: calculated amplitude value.

Definition at line 150 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ Lep() [2/3]

double Lep	(	const double	mtau,
		int	tauhel,
		int	whel,
		double	q2,
		double	costau
	)		const

The function to calculate the Leptonic Amplitudes for B->D*taunu decay of the vector type contribution.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
whel	helicity of the virtual vector boson {+1,0,1,2}.
q2	q^2 of the decay (square of l+nu invariant mass).
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame.

Returns: calculated amplitude value.

Definition at line 127 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ Lep() [3/3]

double Lep	(	const double	mtau,
		int	tauhel,
		int	whel1,
		int	whel2,
		double	q2,
		double	costau
	)		const

The function to calculate the Leptonic Amplitudes for B->D*taunu decay of the tensor type contribution.

Parameters

mtau	daughter lepton mass.
tauhel	helicity of the lepton in the (l+nu) rest frame {+1,-1}.
whel1	helicity of the one virtual vector boson {+1,0,1,2}.
whel2	helicity of the another virtual vector boson {+1,0,1,2}.
q2	q^2 of the decay (square of l+nu invariant mass).
costau	cosine of the angle between D(*) meson and the lepton in the (l+nu) rest frame.

Returns: calculated amplitude value. The tensor contribution is described by two vector contributions.

Definition at line 165 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ mD()

double mD ( int Dhel ) const

Daughter D(*) meson mass.

kinematics

Parameters

Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.

Returns: daughter D(*) meson mass corresponding to the Dhel.

Definition at line 499 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ PDstar()

double PDstar	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		const double	mtau
	)

Function calculates the polarization of D*, longitudinal/(longitudinal + transverse), in B->D*taunu decay.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter tau lepton mass for B->Dtaunu decay.

Definition at line 187 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ pf()

double pf	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		double	mtau,
		int	Dhel,
		double	w
	)

Phase space factor, which is multiplied to the helicity amplitude to calculate the decay rate.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter tau lepton mass for B->Dtaunu decay.
Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.

Definition at line 199 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ PtauD()

double PtauD	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		const double	mtau
	)

Function calculates the polarization of tau, (RH - LH)/(LH + RH), in B->Dtaunu decay.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter tau lepton mass for B->Dtaunu decay.

Definition at line 170 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ PtauDstar()

double PtauDstar	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		const double	mtau
	)

Function calculates the polarization of tau, (RH - LH)/(LH + RH), in B->D*taunu decay.

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter tau lepton mass for B->Dtaunu decay.

Definition at line 177 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ q2()

double q2	(	int	Dhel,
		double	w
	)		const

Function to calculate the q^2 of the decay (square of l+nu invariant mass).

Parameters

Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.

Returns: calculated q^2.

Definition at line 521 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ qh2()

double qh2	(	int	Dhel,
		double	w
	)		const

Function to calculate the q^2 divided by the square of parent mass (m_B^2).

Parameters

Dhel	helicity of the D() meson in the rest frame of the parent meson {+1,0,-1} for D and 2 for D.
w	velocity transfer variable.

Returns: calculated q^2/m_B^2.

Definition at line 515 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ RGammaD()

double RGammaD	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		const double	mtau,
		const double	mlep = `0.0005110`
	)

Function calculates the ratio of Br(B->Dtaunu)/Br(B->Dlnu), R(D).

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter tau lepton mass for B->Dtaunu decay.
mlep	daughter lepton mass for B->Dlnu decay (default is the electron mass). Note that overall factor (GF/sqrt(2) Vcb)^2 is omitted.

Definition at line 149 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ RGammaDstar()

double RGammaDstar	(	const EvtBSemiTauonicHelicityAmplitudeCalculator &	BSTD,
		const double	mtau,
		const double	mlep = `0.0005110`
	)

Function calculates the ratio of Br(B->Dtaunu)/Br(B->Dlnu), R(D*).

Parameters

BSTD	helicity ampliutude calculator initialized properly by user.
mtau	daughter tau lepton mass for B->Dtaunu decay.
mlep	daughter lepton mass for B->Dlnu decay (default is the electron mass).

Definition at line 158 of file EvtBSemiTauonicDecayRateCalculator.cc.

◆ RotateToHelicityBasisInBoostedFrame()

EvtSpinDensity RotateToHelicityBasisInBoostedFrame	(	const EvtParticle *	p,
		EvtVector4R	p4boost
	)

The function calculates the rotation matrix to convert the spin basis to the helicity basis in the boosted frame.

Parameters

p	a pointer to the particle
p4boost	a 4 vector to specify the rest frame to boost to. The function calculate the rotation matrix from the spin basis defined in the p rest frame to the helicity basis in the rest frame of p4boost.

Definition at line 27 of file EvtBSemiTauonicAmplitude.cc.

◆ scaleMomenta()

std::vector<ROOT::Math::PxPyPzMVector> Belle2::scaleMomenta	(	const std::vector< ROOT::Math::PxPyPzMVector > &	ps,
		double	C
	)

inline

Scale momenta of the particles in the vector ps with a factor C.

Result is returned as returned as a vector with the scaled momenta

Definition at line 58 of file scaleParticleEnergies.h.

◆ scaleParticleEnergies()

void Belle2::scaleParticleEnergies	(	MCParticleGraph &	mpg,
		double	EcmsTarget
	)

inline

Scale momenta of the particles by a constant factor such that total energy is changed to EcmsTarget.

It also changes momenta of the incoming particles such that energy is conserved. Should be called in CM reference frame.

Definition at line 71 of file scaleParticleEnergies.h.

◆ term()

void term ( )

Terminates the generator.

Closes the internal generator.

Definition at line 200 of file CRY.cc.

◆ updateVertex()

ROOT::Math::XYZVector updateVertex ( bool force = false )

Update the vertex position:

If there is no initial particles object generate a new one with nominal values without smearing
If the BeamParameters disallow smearing of the vertex it does nothing
If initial particles already exist check if the vertex smearing has already been applied. If not, apply vertex smearing if allowed.
Return the shift in vertex to the previous value (or the origin if there was no previous value).

This function does not update the energies as this would possibly introduce inconsistency between values used by the generator and the values contained in the initial particles. But it is useful to smear the vertex after generation.

Parameters

force if true the vertex will be regenerated even if vertex smearing was already applied.

Definition at line 202 of file InitialParticleGeneration.cc.

◆ v()

double v	(	double	mtau,
		double	q2
	)		const

Function to calculate the tau velocity.

Parameters

mtau	daughter lepton mass.
q2	q^2 of the decay (square of l+nu invariant mass).

Returns: calculated tau velocity.

Definition at line 509 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

◆ z()

double z ( double w ) const

CLN form factor z.

Parameters

w	velocity transfer variable.

Returns: calculated form factor value.

Definition at line 449 of file EvtBSemiTauonicHelicityAmplitudeCalculator.cc.

Modules

Classes

Macros

Functions

Variables

Detailed Description

Macro Definition Documentation

◆ B2_EVTGEN_REGISTER_MODEL

Function Documentation

◆ CalcAmp() [1/2]

◆ CalcAmp() [2/2]

◆ CalcMaxProb()

◆ chkDhel()

◆ chktauhel()

◆ chkwhel()

◆ clone() [1/2]

◆ clone() [2/2]

◆ decay() [1/2]

◆ decay() [2/2]

◆ dGammadcostau()

◆ dGammadw()

◆ dGammadwdcostau()

◆ dR3()

◆ dS1()

◆ EvtBSemiTauonicHelicityAmplitudeCalculator() [1/2]

◆ EvtBSemiTauonicHelicityAmplitudeCalculator() [2/2]

◆ FermiMomentum()

◆ FermiMomentumProb()

◆ ffA1()

◆ ffR1()

◆ ffR2()

◆ ffR3()

◆ ffS1()

◆ ffV1()

◆ Gamma()

◆ GammaD()

◆ GammaDstar()

◆ GammaSMD()

◆ GammaSMDstar()

◆ generate()

◆ generateBeam()

◆ generateEvent()

◆ generateVertex()

◆ getInstance()

◆ getModels()

◆ getName() [1/2]

◆ getName() [2/2]

◆ getscalarff()

◆ getTotalEnergy()

◆ getvectorff()

◆ hA1()

◆ hA2()

◆ hA3()

◆ HadS1()

◆ HadS2()

◆ HadT()

◆ HadV1()

◆ HadV2()

◆ helAmp() [1/2]

◆ helAmp() [2/2]

◆ helampS1()

◆ helampS2()

◆ helampSM()

◆ helampT()

◆ helampV1()

◆ helampV2()

◆ hm()

◆ hp()

◆ hP()

◆ hS()

◆ hT()

◆ hT1()

◆ hT2()

◆ hT3()

◆ hV()

◆ init() [1/5]

◆ init() [2/5]

◆ init() [3/5]

◆ init() [4/5]

◆ init() [5/5]