def contract_triple_charmed_omega(prop_charm,
spin_matrix,
pol_matrix,
diquark=None):
if diquark is None:
diquark = quark_contract_13(gpt.eval(prop_charm * spin_matrix),
gpt.eval(spin_matrix * prop_charm))
return baryon_decuplet_base_contraction(prop_charm, prop_charm,
[diquark for _ in range(3)],
pol_matrix)
def contract_charmed_omega_star(prop_strange,
prop_charm,
spin_matrix,
pol_matrix,
diquarks=None):
if diquarks is None:
diquarks = [
quark_contract_13(gpt.eval(prop_charm * spin_matrix),
gpt.eval(spin_matrix * prop_strange)),
quark_contract_13(gpt.eval(prop_strange * spin_matrix),
gpt.eval(spin_matrix * prop_charm)),
quark_contract_13(gpt.eval(prop_strange * spin_matrix),
gpt.eval(spin_matrix * prop_strange))
]
return baryon_decuplet_base_contraction(prop_charm, prop_strange, diquarks,
pol_matrix)
def contract_double_charmed_xi_plus_star(prop_down,
prop_charm,
spin_matrix,
pol_matrix,
diquarks=None):
if diquarks is None:
diquarks = [
quark_contract_13(gpt.eval(prop_down * spin_matrix),
gpt.eval(spin_matrix * prop_charm)),
quark_contract_13(gpt.eval(prop_charm * spin_matrix),
gpt.eval(spin_matrix * prop_down)),
quark_contract_13(gpt.eval(prop_charm * spin_matrix),
gpt.eval(spin_matrix * prop_charm))
]
return baryon_decuplet_base_contraction(prop_down, prop_charm, diquarks,
pol_matrix)