From a7d7e0b4ba4b69c7a87690d1f197cc005fa957aa Mon Sep 17 00:00:00 2001 From: Karl Berry Date: Fri, 28 Sep 2018 22:31:34 +0000 Subject: pst-feyn (27sep18) git-svn-id: svn://tug.org/texlive/trunk@48781 c570f23f-e606-0410-a88d-b1316a301751 --- Master/texmf-dist/dvips/pst-feyn/pst-feyn.pro | 671 ++++++++++++++++++++++++++ 1 file changed, 671 insertions(+) create mode 100644 Master/texmf-dist/dvips/pst-feyn/pst-feyn.pro (limited to 'Master/texmf-dist/dvips/pst-feyn') diff --git a/Master/texmf-dist/dvips/pst-feyn/pst-feyn.pro b/Master/texmf-dist/dvips/pst-feyn/pst-feyn.pro new file mode 100644 index 00000000000..e49c6a47fee --- /dev/null +++ b/Master/texmf-dist/dvips/pst-feyn/pst-feyn.pro @@ -0,0 +1,671 @@ +%% $Id: pst-feyn.pro 826 2018-09-27 09:21:43Z herbert $ +%% +%% This is file `pst-feyn.pro', +%% It is based on the old axodraw package from Jos Vermaseren +%% +%% IMPORTANT NOTICE: +%% +%% Package `pst-feyn.tex' +%% +%% Herbert Voss +%% +%% This program can be redistributed and/or modified under the terms +%% of the LaTeX Project Public License Distributed from CTAN archives +%% in directory macros/latex/base/lppl.txt. +%% +%% DESCRIPTION: +%% `pst-feyn' is a PSTricks package to draw 3d curves and graphical objects +%% +%% +%% version 0.01 / 2018-09-26 Herbert Voss +%% +% +/tx@FeynDict 100 dict def +tx@FeynDict begin +% +/ED { exch def } def +/SLW { setlinewidth } def +/p1 {/y1 ED /x1 ED } def +/p2 {/y2 ED /x2 ED } def +/abox{ newpath x1 y1 moveto x1 y2 lineto x2 y2 lineto x2 y1 lineto closepath } def +% +% +/arrowdown{ + /nwidth { CLW 1 add } def + newpath + 0 nwidth 1.5 mul neg moveto % Arrow is a triangle + nwidth 1.2 mul nwidth 3 mul rlineto + nwidth 2.4 mul neg 0 rlineto + nwidth 1.2 mul nwidth 3 mul neg rlineto + closepath fill % and it is filled +} def +% +/arrowup{ + /nwidth { CLW 1 add } def + newpath + 0 nwidth 1.5 mul moveto % Arrow is a triangle + nwidth 1.2 mul nwidth 3 mul neg rlineto + nwidth 2.4 mul neg 0 rlineto + nwidth 1.2 mul nwidth 3 mul rlineto + closepath fill % and it is filled +} def +% +/arrowright{ + /nwidth { CLW 1 add } def + newpath + nwidth 1.5 mul 0 moveto % Arrow is a triangle + nwidth 3 mul neg nwidth 1.2 mul rlineto + 0 nwidth 2.4 mul neg rlineto + nwidth 3 mul nwidth 1.2 mul rlineto + closepath fill % and it is filled +} def +% +/gluon{ + %dup scale /width ED width SLW + /num ED /ampi ED /y2 ED /x2 ED /y1 ED /x1 ED + + /dy y2 y1 sub def /dx x2 x1 sub def +% +% We have a 'head' and a 'tail' and inbetween the 'body' +% The head + tail is 2 windings. The body is num-1 windings. +% + /dr dx dx mul dy dy mul add sqrt def +% + x1 y1 translate dy dx atan rotate +% + /num num 0.5 sub round def + /inc dr num 2 mul 2 add div def % increment per half winding + /amp8 ampi 0.9 mul def + amp8 0 lt {/amp8 amp8 neg def} if +% + /x1 inc 2 mul def +% + newpath + 0 0 moveto + inc 0.1 mul ampi 0.5 mul inc 0.5 mul ampi inc 1.4 mul ampi curveto + x1 amp8 add dup ampi exch ampi neg dup x1 exch curveto +% + 2 1 num { + pop + x1 amp8 sub dup ampi neg exch ampi dup x1 inc add exch curveto + /x1 x1 inc dup add add def + x1 amp8 add dup ampi exch ampi neg dup x1 exch curveto + } for +% + x1 amp8 sub dup ampi neg exch ampi dup x1 inc 0.6 mul add exch curveto + x1 inc 1.5 mul add ampi dr inc 0.1 mul sub ampi 0.5 mul dr 0 curveto +% stroke +% +} def +% +/photon{ +% +% Draws a photon from x1,y1 to x2,y2 with amplitude A and n wiggles +% +% gsave dup scale /width ED width SLW + /num ED /ampi ED p2 p1 + /dy y2 y1 sub def /dx x2 x1 sub def + /dr dx dx mul dy dy mul add sqrt def +% + x1 y1 translate + dy dx atan rotate + /num num 2 mul 0.5 sub round def + /x2 dr num div def + /pi 3.141592 def + /sign 1 def + 1 1 num { + pop + newpath + 0 0 moveto + 4 3 div x2 mul pi div dup neg x2 add + 4 3 div ampi sign mul mul dup 3 1 roll + x2 0 curveto + stroke + /sign sign neg def + x2 0 translate + } for +% +% grestore +} def +% +/zigzag{ +% +% Draws a zigzag line from x1,y1 to x2,y2 with amplitude A and n zigzags +% +% gsave dup scale /width ED width SLW + /num ED /ampi ED p2 p1 + /dy y2 y1 sub def /dx x2 x1 sub def + /dr dx dx mul dy dy mul add sqrt def +% + x1 y1 translate + dy dx atan rotate + /num num 2 mul 0.5 sub round def + /x2 dr num div def + /pi 3.141592 def + /sign 1 def + 1 1 num { + pop + newpath + 0 0 moveto + x2 2 div ampi sign mul lineto + x2 0 lineto + stroke + /sign sign neg def + x2 0 translate + } for +% +% grestore +} def +% +/photonarc{ +% +% Draws a photonarc center at x1,y1, radius arcstart,arcend, amplitude +% number of wiggles, width, scale +% +% gsave dup scale /width ED width SLW + /num ED /ampli ED /arcend ED /arcstart ED /radius ED +% + translate % The center of the circle is now the origin +% + /num num 2 mul round def % number of half wiggles + arcend arcstart lt { /arcend arcend 360 add def } if + /arcend arcend arcstart sub num div def % phi + arcstart rotate + /arcstart arcend 2 div def % phi/2 + /cp arcend cos def + /sp arcend sin def + /cp2 arcstart cos def + /sp2 arcstart sin def +% + newpath + 1 1 num { + pop + radius 0 moveto + /beta radius arcend mul 180 ampli mul div def + /tt sp cp beta mul sub cp sp beta mul add div def + /amp1 radius ampli add 8 mul beta cp2 mul sp2 sub mul beta 4 cp add mul + tt cp mul 3 mul sp 4 mul sub add radius mul sub + beta tt sub 3 mul div def % this is x2 + radius ampli add 8 mul cp2 mul 1 cp add radius mul sub 3 div amp1 sub + dup radius sub beta mul % x1,y1 + amp1 dup radius cp mul sub tt mul radius sp mul add % x2,y2 + radius cp mul radius sp mul % x3 y3 + curveto + /ampli ampli neg def + arcend rotate + } for +% stroke +% +% grestore +} def +% +/gluearc{ +% +% Draws a gluon on an arcsegment +% radius,stat_angle,end_angle,gluon_radius,num +% in which num is the number of windings of the gluon. +% Method: +% 1: compute length of arc. +% 2: generate gluon in x and y as if the arc is a straight line +% 3: x' = (radius+y)*cos(x*const) +% y' = (radius+y)*sin(x*const) +% + gsave %dup scale /width ED width SLW + /num ED + /ampi ED + /arcend ED + /arcstart ED + /radius ED +% + arcstart rotate % segment starts at zero + /darc arcend arcstart sub def % argsegment + /dr darc 180 div 3.141592 mul radius mul def % length of segment. + /const darc dr div def % conversion constant +% + /num num 0.5 sub round def + /inc dr num 2 mul 2 add div def % increment per half winding +% + /amp8 ampi 0.9 mul def + /amp1 radius ampi add def + /amp2 radius ampi sub def + /amp3 radius ampi 2 div add def + /amp4 amp1 inc amp8 add const mul cos div def + /amp5 amp2 amp8 const mul cos div def + /amp6 amp1 inc 0.6 mul amp8 add const mul cos div def + /amp7 amp1 inc 0.9 mul const mul cos div def + amp8 0 lt {/amp8 amp8 neg def} if +% + /x1 inc 2 mul def +% + newpath + radius 0 moveto +% + inc 0.1 mul const mul dup cos amp3 mul exch sin amp3 mul + inc 0.5 mul const mul dup cos amp7 mul exch sin amp7 mul + inc 1.4 mul const mul dup cos amp1 mul exch sin amp1 mul + curveto + x1 amp8 add const mul dup cos amp6 mul exch sin amp6 mul + x1 amp8 add const mul dup cos amp5 mul exch sin amp5 mul + x1 const mul dup cos amp2 mul exch sin amp2 mul + curveto +% + 2 1 num { + pop + x1 amp8 sub const mul dup cos amp5 mul exch sin amp5 mul + x1 amp8 sub const mul dup cos amp4 mul exch sin amp4 mul + x1 inc add const mul dup cos amp1 mul exch sin amp1 mul + curveto + /x1 x1 inc dup add add def + x1 amp8 add const mul dup cos amp4 mul exch sin amp4 mul + x1 amp8 add const mul dup cos amp5 mul exch sin amp5 mul + x1 const mul dup cos amp2 mul exch sin amp2 mul + curveto + } for +% + x1 amp8 sub const mul dup cos amp5 mul exch sin amp5 mul + x1 amp8 sub const mul dup cos amp6 mul exch sin amp6 mul + x1 inc 0.6 mul add const mul dup cos amp1 mul exch sin amp1 mul + curveto + x1 inc 1.5 mul add const mul dup cos amp7 mul exch sin amp7 mul + dr inc 0.1 mul sub const mul dup cos amp3 mul exch sin amp3 mul + dr const mul dup cos radius mul exch sin radius mul + curveto + stroke +% + grestore +} def +% +/arrowarc{ +% +% Draws an anticlockwise arc with an arrow in the middle +% The arc is x_center, y_center, radius, start_angle, end_angle +% +% gsave dup scale /width ED width SLW + /arcend ED /arcstart ED /radius ED +% + translate % x and y are still on stack + newpath 0 0 radius arcstart arcend arc stroke + arcstart arcend gt { + /arcend arcend 360 add def } if + arcstart arcend add 2 div rotate % middle of arc + radius 0 translate % move to it + arrowup +% grestore +} def +% +/longarrowarc{ +% +% Draws an anticlockwise arc with an arrow at the end +% The arc is x_center, y_center, radius, start_angle, end_angle +% +% gsave dup scale /width ED width SLW + /arcend ED /arcstart ED /radius ED +% + translate % x and y are still on stack + arcstart arcend gt { + /arcend arcend 360 add def } if + /arcmid 540 CLW 1 add mul 3.14159 div radius div def + % discount for arrow +% newpath 0 0 radius arcstart arcend arcmid sub arc stroke + arcend arcmid 2 div sub rotate % middle of arrow + radius 0 translate % move to it + arrowup +% grestore +} def +% +/dasharrowarc{ +% +% Draws a dashed anticlockwise arc with an arrow in the middle +% The arc is x_center, y_center, radius, start_angle, end_angle dsize +% + gsave dup scale /width ED width SLW /dsize ED /arcend1 ED /arcstart1 ED /radius ED +% + translate % x and y are still on stack +% + arcend1 arcstart1 lt { /arcend1 arcend1 360 add def } if + /arcmid1 arcend1 arcstart1 add 2 div def +% + 0 0 radius arcstart1 arcmid1 dsize width 1 dashcarc + 0 0 radius arcmid1 arcend1 dsize width 1 dashcarc + arcmid1 rotate + radius 0 translate + arrowup + grestore +} def +% +/arrowarcn{ +% +% Draws a clockwise arc with an arrow in the middle +% The arc is x_center, y_center, radius, start_angle, end_angle +% +% gsave dup scale /width ED width SLW + /arcend ED /arcstart ED /radius ED +% + translate % x and y are still on stack + newpath 0 0 radius arcstart arcend arcn stroke + arcstart arcend lt { + /arcstart arcstart 360 add def } if + arcstart arcend add 2 div rotate % middle of arc + radius 0 translate % move to it + arrowdown +% grestore +} def +% +/longarrowarcn{ +% +% Draws a clockwise arc with an arrow in the end +% The arc is x_center, y_center, radius, start_angle, end_angle +% + gsave dup scale /width ED width SLW /arcend ED /arcstart ED /radius ED +% + translate % x and y are still on stack + arcstart arcend lt { + /arcstart arcstart 360 add def } if + /arcmid 540 width 1 add mul 3.14159 div radius div def + % correction for arrow + newpath 0 0 radius arcstart arcend arcmid add arcn stroke + arcend arcmid 2 div add rotate % middle of arrow + radius 0 translate % move to it + arrowdown + grestore +} def +% +/dasharrowarcn{ +% +% Draws a dashed clockwise arc with an arrow in the middle +% The arc is x_center, y_center, radius, start_angle, end_angle +% + gsave dup scale /width ED width SLW /dsize ED /arcend1 ED /arcstart1 ED /radius ED +% + translate % x and y are still on stack + arcstart1 arcend1 lt { + /arcstart1 arcstart1 360 add def } if + /arcmid1 arcstart1 arcend1 add 2 div def + 0 0 radius arcmid1 arcstart1 dsize width 1 dashcarc + 0 0 radius arcend1 arcmid1 dsize width 1 dashcarc + arcmid1 rotate + radius 0 translate + arrowdown + grestore +} def +% +/arrowline{ +% +% Draws a straight line with an arrow in the middle +% x1,y1,x2,y2 +% + %gsave dup scale /width ED width SLW + p2 p1 + /dx x2 x1 sub def /dy y2 y1 sub def + /dr dx dx mul dy dy mul add sqrt def +% + x1 y1 translate +% newpath +% 0 0 moveto +% dx dy rlineto +% stroke + dy dx atan rotate + dr 2.0 div 0 translate + arrowright +% grestore +} def +% +/longarrow{ +% +% Draws a straight line with an arrow at the end +% x1,y1,x2,y2 +% + gsave dup scale /width ED width SLW p2 p1 + /dx x2 x1 sub def /dy y2 y1 sub def + /dr dx dx mul dy dy mul add sqrt def +% + x1 y1 translate + dy dx atan rotate + newpath + 0 0 moveto + dr width 3 mul sub 0 rlineto + stroke + dr width 1.5 mul sub 0 translate + arrowright + grestore +} def +% +/dasharrowline{ +% +% Draws a straight dashed line with an arrow in the middle +% x1,y1,x2,y2 +% +% The pattern is ideally [dsize dsize] 0 setdash +% but we want to have (2*n+1)/2 patterns, so dsize must be rounded +% Actually we want the center to be black too so that the arrow +% fits in nice. This means that n must be odd. So +% r = dsize*(4*m+3) +% + gsave dup scale /width ED width SLW /dsize ED p2 p1 + /dx x2 x1 sub def /dy y2 y1 sub def + /dr dx dx mul dy dy mul add sqrt 2 div def +% + x1 y1 translate + dy dx atan rotate +% + 0 0 dr 0 dsize width 1 dashline + dr 0 translate + 0 0 dr 0 dsize width 1 dashline + arrowright + grestore +} def +% +/carc{ +% +% Draws an anti-clockwise arc segment: +% x_center, y_center, radius, start_angle, end_angle +% + %dup scale /width ED width SLW + /arcend ED /arcstart ED /radius ED +% + translate % x and y are still on stack + newpath 0 0 radius arcstart arcend arc +% stroke +% grestore +} def +% +/dashcarc{ +% +% Draws an anti-clockwise arc segment: +% x_center, y_center, radius, start_angle, end_angle, dsize +% + gsave dup scale /width ED width SLW /dsize ED /arcend ED /arcstart ED /radius ED +% + translate % x and y are still on stack +% +% Compute the length of the line +% + /dr arcend arcstart sub dup 0 lt { 360 add } if + 3.14159 mul 180 div radius mul def + /dsize dr dsize 2 mul div 0.5 sub round dup 0 le { pop 0 } if 2 mul 1 add + dr exch div def + [dsize dsize] 0 setdash +% + newpath 0 0 radius arcstart arcend arc stroke + grestore +} def +% +/bcirc{ +% +% Draws an anti-clockwise blanked circle: +% x_center, y_center, radius +% + gsave dup scale /width ED width SLW /radius ED +% + translate % x and y are still on stack +% + gsave + 1 setgray + newpath 0 0 radius 0 360 arc fill + grestore + newpath 0 0 radius 0 360 arc stroke + grestore +} def +% +/gcirc{ +% +% Draws an anti-clockwise blanked gray circle: +% x_center, y_center, radius, grayscale +% + gsave dup scale /width ED width SLW /gcolor ED /radius ED +% + translate % x and y are still on stack +% + 1 setgray + newpath 0 0 radius 0 360 arc fill + gcolor setgray + newpath 0 0 radius 0 360 arc fill + 0 setgray + newpath 0 0 radius 0 360 arc stroke + grestore +} def +% +/ccirc1{ +% +% Draws an anti-clockwise circle : +% x_center, y_center, radius +% Part 1: the contents in background color +% + gsave dup scale /width ED width SLW /radius ED +% + translate % x and y are still on stack +% + newpath 0 0 radius 0 360 arc fill + grestore +} def +% +/ccirc2{ +% +% Draws an anti-clockwise circle : +% x_center, y_center, radius +% Part 1: the contents in background color +% + gsave dup scale /width ED width SLW /radius ED +% + translate % x and y are still on stack +% + newpath 0 0 radius 0 360 arc stroke + grestore +} def +% +% +/pathlength{ + flattenpath + /dist 0 def + { /yfirst ED /xfirst ED /ymoveto yfirst def /xmoveto xfirst def } + { /ynext ED /xnext ED /dist dist ynext yfirst sub dup mul + xnext xfirst sub dup mul add sqrt add def + /yfirst ynext def /xfirst xnext def } + {} + {/ynext ymoveto def /xnext xmoveto def + /dist ynext yfirst sub dup mul + xnext xfirst sub dup mul add sqrt add def + /yfirst ynext def /xfirst xnext def } + pathforall + dist +} def +% +/centerdash{ + /pathlen pathlength def + /jj pathlen dsize div 2.0 div cvi def + /ddsize pathlen jj 2.0 mul div def + [ddsize] ddsize 2 div setdash +} def +% +/logaxis{ +% +% Draws an axis from x1,y1 to x2,y2 with nl log divisions +% size of the hashes hs, offset F +% and width W. The stack looks like +% x1,y1,x2,y2,nl,hs,F,W,scale +% After the rotation the hash marks are on top if nl is positive and +% on the bottom if nl is negative +% +% gsave dup scale /width ED width SLW + /offset ED /hashsize ED /nlogs ED p2 p1 + x1 y1 translate + /y2 y2 y1 sub def /x2 x2 x1 sub def + y2 x2 atan rotate + /rr x2 dup mul y2 dup mul add sqrt def + offset 0 ne { /offset offset ln 10 ln div def } if + /offset offset dup cvi sub def + newpath + 0 0 moveto + rr 0 lineto + /lsize rr nlogs div def + 0 1 nlogs { /x2 ED + x2 offset ge { + /y2 x2 offset sub lsize mul def + y2 rr le { + y2 0 moveto + y2 hashsize 1.2 mul lineto + } if + } if + } for + stroke + CLW 0.6 mul setlinewidth + newpath + 0 1 nlogs { /x2 ED + 2 1 9 { + ln 10 ln div x2 add + /xx2 ED + xx2 offset ge { + /y2 xx2 offset sub lsize mul def + y2 rr le { + y2 0 moveto + y2 hashsize 0.8 mul lineto + } if + } if + } for + } for +% stroke +% grestore +} def +% +/linaxis{ +% +% x1,y1,x2,y2,num_decs,per_dec,hashsize,offset,width,scale +% + %gsave dup scale /width ED width SLW + /offset ED /hashsize ED /perdec ED /numdec ED p2 p1 + x1 y1 translate + /y2 y2 y1 sub def /x2 x2 x1 sub def + y2 x2 atan rotate + /rr x2 dup mul y2 dup mul add sqrt def + newpath + 0 0 moveto + rr 0 lineto + /x1 rr numdec perdec mul div def + /y1 rr numdec div def + /offset offset x1 mul def + 0 1 numdec { y1 mul offset sub + dup 0 ge { + dup rr le { + dup 0 moveto + hashsize 1.2 mul lineto + } if + } if + } for + stroke + CLW 0.6 mul setlinewidth + newpath + offset cvi 1 numdec perdec mul offset add { + x1 mul offset sub + dup 0 ge { + dup rr le { + dup 0 moveto + hashsize 0.8 mul lineto + } if + } if + } for + stroke + %grestore +} def +% +end % tx@FeynDict +% \ No newline at end of file -- cgit v1.2.3