Track doesn't radiate in quadrupole

[ldeniau]

I suspect that pt is not updated properly in Track...

x0 = 0.002;
px0=-0.001;

!el: sbend, l=1, angle=6.3e-4, e1=0.1, e2=0.1;
el: quadrupole, l=1, k1=0.2;

ss: sequence, l=1;
el, at=0.5;
endsequence;

beam,particle="positron",energy=100,radiate=true; ! FCC-ee
use,sequence=ss;

twiss,x=x0,px=px0,betx=1,bety=1;
twiss_x =table(twiss,ss$end, x);
twiss_y =table(twiss,ss$end, y);
twiss_px=table(twiss,ss$end,px);
twiss_py=table(twiss,ss$end,py);
twiss_pt=table(twiss,ss$end,pt);
write, table=twiss, file="twx.tfs";

track,onepass=true;
start,x=x0,px=px0;
run,turns=1;
endtrack;
track_x =table(tracksumm, x,2);
track_y =table(tracksumm, y,2);
track_px=table(tracksumm,px,2);
track_py=table(tracksumm,py,2);
track_pt=table(tracksumm,pt,2);
write, table=tracksumm, file="tkx.tfs";

option, -echo;
value,twiss_x,twiss_px,twiss_y,twiss_py,twiss_pt;
value,track_x,track_px,track_y,track_py,track_pt;

Gives the output (I added PTC for comparison)

twiss_x            =    0.0008363125922 ;
twiss_px           =    -0.001288453569 ;
twiss_y            =                  0 ;
twiss_py           =                  0 ;
twiss_pt           =   -1.323288621e-06 ;
track_x            =    0.0008363127732 ;
track_px           =    -0.001288454962 ;
track_y            =                  0 ;
track_py           =                  0 ;
track_pt           =                  0 ;
ptc_x              =    0.0008363117858 ;
ptc_px             =    -0.001288453481 ;
ptc_y              =                  0 ;
ptc_py             =                  0 ;
ptc_pt             =   -1.323288746e-06 ;