本文整理汇总了C++中Flavour::IsoWeak方法的典型用法代码示例。如果您正苦于以下问题:C++ Flavour::IsoWeak方法的具体用法?C++ Flavour::IsoWeak怎么用?C++ Flavour::IsoWeak使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Flavour
的用法示例。
在下文中一共展示了Flavour::IsoWeak方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: c_FFV
void Interaction_Model_SM_Zprime::c_FFV(vector<Single_Vertex>& vertex,int& vanz)
{
// create the vertices for the standard model
Interaction_Model_SM::c_FFV(vertex,vanz);
// create FFV vertices with Z' if it's on
kf_code kfZp=32;
Flavour flZprime(kfZp);
if (flZprime.IsOn()) {
// parse through all fermions than couple to Z' and create vertices
int PossibleFermions[12] = {1,2,3,4,5,6,11,12,13,14,15,16};
for (int i=0; i<12; i++) {
// initialize the currently parsed fermion
int FermionNumber = PossibleFermions[i];
Flavour flFermion = Flavour((kf_code)(FermionNumber));
Kabbala B = Kabbala(string("B_{")+flFermion.TexName()+string("}"),
flFermion.BaryonNumber());
Kabbala L = Kabbala(string("L_{")+ flFermion.TexName()+string("}"),
flFermion.LeptonNumber());
Kabbala Y3R = Kabbala(string("YR_{")+flFermion.TexName()+string("}"),
flFermion.IsoWeak());
if (flFermion.IsOn()) {
// create the vertex for that particular fermion and a Z'.
// Right-handed neutrinos will not take part in any interaction.
Kabbala kcpl0;
if ((FermionNumber==12)||(FermionNumber==14)||(FermionNumber==16))
{kcpl0 = Kabbala("0.0", 0.);}
else {kcpl0 = -M_I * g2 * (Y3R * alphaLR + (L-B)/(alphaLR*2));};
Kabbala kcpl1 = -M_I * g2 * (L-B) / (alphaLR*2);
// set couplings and particle info for current vertex
vertex[vanz].in[0] = flFermion;
vertex[vanz].in[1] = flZprime;
vertex[vanz].in[2] = Flavour((kf_code)(FermionNumber));
vertex[vanz].cpl[0] = kcpl0;
vertex[vanz].cpl[1] = kcpl1;
vertex[vanz].Str = (kcpl0*PR+kcpl1*PL).String();
// Color Function for vertex
if (flFermion.Strong()) {
vertex[vanz].Color.push_back(Color_Function(cf::D));;
vertex[vanz].Color.back().SetParticleArg(0,2);
vertex[vanz].Color.back().SetStringArg('0','2');
}
else
vertex[vanz].Color.push_back(Color_Function(cf::None));;
// Lorenz function for vertex
vertex[vanz].Lorentz.push_back(LF_Getter::GetObject("Gamma",LF_Key()));
vertex[vanz].Lorentz.back()->SetParticleArg(1);
vertex[vanz].on = 1;
vertex.push_back(Single_Vertex());vanz++;
};
};
};
}