Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F10881345
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
2 KB
Subscribers
None
View Options
diff --git a/current_generator/jW_juno.frm b/current_generator/jW_juno.frm
index 3288077..84410f1 100644
--- a/current_generator/jW_juno.frm
+++ b/current_generator/jW_juno.frm
@@ -1,62 +1,65 @@
*/**
* \brief Contraction of W current with unordered current
*
* \authors The HEJ collaboration (see AUTHORS for details)
* \date 2020
* \copyright GPLv2 or later
*/
#include- include/helspin.frm
#include- include/write.frm
v pb,p2,pa,p1,pg,pl,plbar,pW,q2,q3,p,pr;
i mu,nu,sigma,alpha;
cf jW,epsg,m2,sqrt;
s h;
#do h2={+,-}
#do hg={+,-}
* eq:juno in developer manual with pa -> pb, p1 -> pg
* contracted with jW
l [U1_jW `h2'`hg'] = (
+ Current(`h2'1, p2, nu, pg)*Current(`h2'1, pg, mu, pb)
+ 2*p2(nu)*Current(`h2'1, p2, mu, pb)
)/m2(p2 + pg)*jW(-1, pa, p1, plbar, pl)*epsg(`hg'1, nu, pr);
l [U2_jW `h2'`hg'] = (
+ 2*Current(`h2'1, p2, mu, pb)*pb(nu)
- Current(`h2'1, p2, mu, pg)*Current(`h2'1, pg, nu, pb)
)/m2(pb - pg)*jW(-1, pa, p1, plbar, pl)*epsg(`hg'1, nu, pr);
l [L_jW `h2'`hg'] = (
- (q2(nu) + q3(nu))*Current(`h2'1, p2, mu, pb)
- 2*pg(mu)*Current(`h2'1, p2, nu, pb)
+ 2*Current(`h2'1, p2, pg, pb)*d_(mu, nu)
+ m2(q2)*Current(`h2'1, p2, mu, pb)*(p1(nu)/m2(p1+pg) + pa(nu)/m2(pa+pg))
)/m2(q3)*jW(-1, pa, p1, plbar, pl)*epsg(`hg'1, nu, pr);
#enddo
#enddo
+* choice of auxiliary vector
+id epsg(h?, mu?, pr)*Current(h?, ?a) = epsg(h, mu, p2)*Current(h, ?a);
+also epsg(h?, mu?, pr) = epsg(h, mu, pb);
+
* eq:Weffcur1 in developer manual
id jW(h?, pa?, p1?, plbar?, pl?) = Current(-1, pl, alpha, plbar)*(
+ Current(h,p1,alpha,p1+pW,mu,pa)/m2(p1+pW)
+ Current(h,p1,mu,pa-pW,alpha,pa)/m2(pa-pW)
);
id epsg(-1, mu?, pr?) = sqrt(2)/2*SpinorChain(pr, mu, pg)/AngleChain(pg,pr);
id epsg(+1, mu?, pr?) = sqrt(2)/2*SpinorChain(pg, mu, pr)/SquareChain(pg,pr);
multiply replace_(
pW,pl+plbar,
q2,p2+pg-pb,
- q3,p2-pb,
- pr,pb
+ q3,p2-pb
);
.sort
#call ContractCurrents
.sort
format O4;
format c;
#call WriteHeader(`OUTPUT')
#call WriteOptimised(`OUTPUT',U1_jW,2,pa,p1,pb,p2,pg,pl,plbar)
#call WriteOptimised(`OUTPUT',U2_jW,2,pa,p1,pb,p2,pg,pl,plbar)
#call WriteOptimised(`OUTPUT',L_jW,2,pa,p1,pb,p2,pg,pl,plbar)
#call WriteFooter(`OUTPUT')
.end
File Metadata
Details
Attached
Mime Type
text/x-diff
Expires
Sat, May 3, 6:13 AM (1 d, 7 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4980228
Default Alt Text
(2 KB)
Attached To
rHEJ HEJ
Event Timeline
Log In to Comment