Loading examples/macro_fit_and_mask 0 → 100644 +35 −0 Original line number Diff line number Diff line macro ! ===================================== ! test for masking tau and phase offset ! ===================================== clear data set ns 1e-9 set A 1e-10 bins pix nbins 16 bins tof custom 4 14 25 38 histo tau nbins 30 min 0.05*ns max 50*ns log histo q nbins 2 min 0.03/A max 0.07/A datapath /SNS/NSE/IPTS-25952 read s10190.echo as res ! 8A Q=0.050 match all fit res read s10556.echo as sam ! 8A Q=0.050! match all ! here masking produces warning mask tau set 2 run 10556 fit sample flag offset phase_offset 180.0 !collect ! here masking works ! mask tau set 2 run 10556 !collect sources/fit_utils.f90 +1 −1 Original line number Diff line number Diff line Loading @@ -373,7 +373,7 @@ contains enddo dtxpd if(ssq_norpix == 0d0 .or. ssq_norpix >= HUGE(1d0)) then call msg_error('get_phase_offset','phase offset determination failed ('& call msg_warning('get_phase_offset','phase offset determination failed ('& ! //trim(msg_fmt("('ssq normalization: ', g9.5)", ssq_norpix))//')' & //'Check: r-parameters r.ave_min, r.fqt_min or r.fqt_maxsigma !)') djoffset = 0 Loading Loading
examples/macro_fit_and_mask 0 → 100644 +35 −0 Original line number Diff line number Diff line macro ! ===================================== ! test for masking tau and phase offset ! ===================================== clear data set ns 1e-9 set A 1e-10 bins pix nbins 16 bins tof custom 4 14 25 38 histo tau nbins 30 min 0.05*ns max 50*ns log histo q nbins 2 min 0.03/A max 0.07/A datapath /SNS/NSE/IPTS-25952 read s10190.echo as res ! 8A Q=0.050 match all fit res read s10556.echo as sam ! 8A Q=0.050! match all ! here masking produces warning mask tau set 2 run 10556 fit sample flag offset phase_offset 180.0 !collect ! here masking works ! mask tau set 2 run 10556 !collect
sources/fit_utils.f90 +1 −1 Original line number Diff line number Diff line Loading @@ -373,7 +373,7 @@ contains enddo dtxpd if(ssq_norpix == 0d0 .or. ssq_norpix >= HUGE(1d0)) then call msg_error('get_phase_offset','phase offset determination failed ('& call msg_warning('get_phase_offset','phase offset determination failed ('& ! //trim(msg_fmt("('ssq normalization: ', g9.5)", ssq_norpix))//')' & //'Check: r-parameters r.ave_min, r.fqt_min or r.fqt_maxsigma !)') djoffset = 0 Loading