if highpass>0 : x.process_inplace("filter.highpass.gauss",{"cutoff_freq":1.0/highpass})
get_application().processEvents() // keeps the GUI happy
y=ptcl.process("misc.directional_sum",{"axis":"y","first":ptcl["nx"]/2-layers,"last":ptcl["nx"]/2+layers+1})
if lowpass>0 : y.process_inplace("filter.lowpass.gauss",{"cutoff_freq":1.0/lowpass})
if highpass>0 : y.process_inplace("filter.highpass.gauss",{"cutoff_freq":1.0/highpass})
get_application().processEvents()
After Change
if highpass>0 : y.process_inplace("filter.highpass.gauss",{"cutoff_freq":1.0/highpass})
get_application().processEvents()
z=ptcl.process("misc.directional_sum",{"axis":"z","first":first,"last":last})
if lowpass>0 : z.process_inplace("filter.lowpass.gauss",{"cutoff_freq":1.0/lowpass})
if highpass>0 : z.process_inplace("filter.highpass.gauss",{"cutoff_freq":1.0/highpass})
// different directions sometimes have vastly different standard deviations, independent normalization may help balance
