MatrixProduct / CKM

[AntonK]

Dear Florian,

I got a question regardless of MatrixProduct function. As an example, let us consider a provided THDM model. For W+- you get:

Code: Select all

Vertex[{bar[Fd], Fu, VWm}]

{{-((I g2 Delta[ct1, ct2] sum[j1, 1, 3, conj[Vu[gt2, j1]] Vd[gt1, j1]])/Sqrt[2]),  LorentzProduct[gamma[lt3], PL]},
 {0, PR}}

You can specify CKM matrix by:

Code: Select all

MatrixProduct -> {Vd,Vu}

and thus end up getting:

Code: Select all

{{-((I g2 Delta[ct1, ct2] Vckm[gt2, gt1])/Sqrt[2]),  LorentzProduct[gamma[lt3], PL]}, {0, PR}}

On the other hand, for the charged Higgs:

Code: Select all

 Vertex[{ bar[Fd], Fu, Hm }]

you get:

Code: Select all

{{-I Delta[ct1, ct2] sum[j2, 1, 3,  conj[Vu[gt2, j2]] sum[j1, 1, 3, conj[Ud[gt1, j1]] Yd[j1, j2]]] ZP[ gt3, 2],  PL},
{I Delta[ct1, ct2] sum[j2, 1, 3, sum[j1, 1, 3, conj[Yu[j1, j2]] Uu[gt2, j1]] Vd[gt1, j2]] ZP[gt3, 2], PR}}

In the documentation it is written that function

Code: Select all

{X, {MatrixProduct-> {A,B} }} 

repalces A^\dagger B by X but how to deal with the aforementioned case, i.e.
conj[Vu] conj[Ud] -> V_CKM <- Uu Vd

P.S. if it is relevant, I am still trying to figure out the previous issue with negative BR in CalcHEP and was curious if it is possible to specify V_CKM in this case as I was doing that by hand.

With best regards,
Anton

[FStaub]

Hi,

I'm not sure if I got the question. Did you check the SM and MSSM. There are CKM versions included in SARAH.

Cheers
Florian