% % Standard Model - unitary and t'Hooft-Feynman gauges. % option chepPDWidth=200. keys gauge_fixing=Feynman. keys CKMdim=1. do_if gauge_fixing==Feynman. model 'SM+DM-CONT-FDM'/5. do_else_if gauge_fixing==unitary. model 'SM+DM-CONT(UG)'/15. do_else. write('Error: the key "gauge" should be either "Feynman" or "unitary".'). quit. end_if. read sm. prtcproperty pdg:(n1=12, n2=14, n3=16, e1=11, e2=13, e3=15). %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%% EFFECTIVE Fermion DM Interactions %%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% spinor '~fdm'/'~Fdm': (FDM, mass Mfdm=100, pdg 9000002, aux 0). let sig^mu^nu=i*(gamma^mu*gamma^nu-gamma^nu*gamma^mu)/2. parameter Lam=1000. parameter gD1,gD1Q,gD2,gD3,gD4,gD1T,gD2T,gD3T,gD4T,gD5,gD6,gD7,gD8,gD9,gD10. %%%%%%%%%%%%%%%%%%%%% lterm gD1/Lam**2*anti('~fdm')*'~fdm'*anti(q)*q where q=u; q=d ; q=s; q=c. lterm gD1Q*Mq/Lam**3*anti('~fdm')*'~fdm'*anti(psi)*psi where psi=u,Mq=Mu; psi=d,Mq=Md; psi=s,Mq=Ms; psi=c,Mq=Mcp; psi=b,Mq=Mbp; psi=t,Mq=Mtop. lterm i*gD2/Lam**2*(anti('~fdm')*gamma5*'~fdm')*(anti(q)*q) where q=u; q=d ; q=s; q=c. lterm i*gD3/Lam**2*(anti('~fdm')*'~fdm')*(anti(q)*gamma5*q) where q=u; q=d ; q=s; q=c. lterm gD4/Lam**2*(anti('~fdm')*gamma5*'~fdm')*(anti(q)*gamma5*q) where q=u; q=d ; q=s; q=c. %%%%%%%%%%%%%%%%%%%%% let aprt_u=AuxPrt(0,3,1). let aprt_d=AuxPrt(0,3,1). let aprt_s=AuxPrt(0,3,1). let aprt_c=AuxPrt(0,3,1). parameter cf_D1T=sqrt(gD1T)/Lam. lterm cf_D1T*(anti(q)*'~fdm'*aprt + anti('~fdm')*q*anti(aprt)) where aprt=aprt_u, q=u; aprt=aprt_d, q=d; aprt=aprt_c, q=c; aprt=aprt_s, q=s. %%%%%%% let aprt_u=AuxPrt(0,3,1). let aprt_d=AuxPrt(0,3,1). let aprt_s=AuxPrt(0,3,1). let aprt_c=AuxPrt(0,3,1). parameter cf_D2T=sqrt(gD2T)/Lam/sqrt(2). lterm cf_D2T*(anti(q)*'~fdm'*aprt + i*anti('~fdm')*gamma5*q*anti(aprt)) where aprt=aprt_u, q=u; aprt=aprt_d, q=d; aprt=aprt_c, q=c; aprt=aprt_s, q=s. let aprt_u=AuxPrt(0,3,1). let aprt_d=AuxPrt(0,3,1). let aprt_s=AuxPrt(0,3,1). let aprt_c=AuxPrt(0,3,1). lterm cf_D2T*(i*anti(q)*gamma5*'~fdm'*aprt + anti('~fdm')*q*anti(aprt)) where aprt=aprt_u, q=u; aprt=aprt_d, q=d; aprt=aprt_c, q=c; aprt=aprt_s, q=s. %%%%%%%% let aprt_u=AuxPrt(0,3,1). let aprt_d=AuxPrt(0,3,1). let aprt_s=AuxPrt(0,3,1). let aprt_c=AuxPrt(0,3,1). parameter cf_D3T=sqrt(gD3T)/Lam/sqrt(2). lterm cf_D3T*(anti(q)*'~fdm'*aprt + anti('~fdm')*gamma5*q*anti(aprt)) where aprt=aprt_u, q=u; aprt=aprt_d, q=d; aprt=aprt_c, q=c; aprt=aprt_s, q=s. let aprt_u=AuxPrt(0,3,1). let aprt_d=AuxPrt(0,3,1). let aprt_s=AuxPrt(0,3,1). let aprt_c=AuxPrt(0,3,1). lterm cf_D3T*(anti(q)*gamma5*'~fdm'*aprt - anti('~fdm')*q*anti(aprt)) where aprt=aprt_u, q=u; aprt=aprt_d, q=d; aprt=aprt_c, q=c; aprt=aprt_s, q=s. %%%%%%%%% let aprt_u=AuxPrt(0,3,1). let aprt_d=AuxPrt(0,3,1). let aprt_s=AuxPrt(0,3,1). let aprt_c=AuxPrt(0,3,1). parameter cf_D4T=sqrt(gD4T)/Lam. lterm cf_D4T*(anti(q)*gamma5*'~fdm'*aprt + anti('~fdm')*gamma5*q*anti(aprt)) where aprt=aprt_u, q=u; aprt=aprt_d, q=d; aprt=aprt_c, q=c; aprt=aprt_s, q=s. %%%%%%%%%%%%%%%%%%%%% lterm gD5/Lam**2*(anti('~fdm')*gamma^mu*'~fdm')*(anti(q)*gamma^mu * q) where q=u; q=d ; q=s; q=c. lterm gD6/Lam**2*(anti('~fdm')*gamma^mu*gamma5*'~fdm')*(anti(q)*gamma^mu *q) where q=u; q=d ; q=s; q=c. lterm gD7/Lam**2*(anti('~fdm')*gamma^mu*'~fdm')*(anti(q)*gamma^mu *gamma5*q) where q=u; q=d ; q=s; q=c. lterm gD8/Lam**2*(anti('~fdm')*gamma^mu*gamma5*'~fdm')*(anti(q)*gamma^mu *gamma5*q) where q=u; q=d ; q=s; q=c. lterm gD9/Lam**2*(anti('~fdm')*i*(gamma^mu*gamma^nu-gamma^nu*gamma^mu)/2*'~fdm')*(anti(q)*i*(gamma^mu*gamma^nu-gamma^nu*gamma^mu)/2*q) where q=u; q=d ; q=s; q=c. lterm i*gD10/Lam**2*(anti('~fdm')*i*(gamma^mu*gamma^nu-gamma^nu*gamma^mu)/2*gamma5*'~fdm')*(anti(q)*i*(gamma^mu*gamma^nu-gamma^nu*gamma^mu)/2*q) where q=u; q=d ; q=s; q=c. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% SetAngle(1-SW**2=CW**2). SetEM(A,EE). CheckHerm.