Privacy and Security Notice
replay_errors
Replay
errors 03/14/06
After inserting the satcorr parameterization into the Analyzer, there
is a very common error message in the output from replay. Not
sure what this is about. Some runs are filled with them (49941)
and
others have none. I tracked the error message to
TRACKING/total_eloss.f (line 685), but it is not clear what went
wrong. When
I take the new satcorr parameterization out of the Analyzer, the errors
go away.
loss: eloss<=0!
electron?
ztgt atgt
thick dens
velocity e_loss
1.00000 7.32000 14.68000
0.01800 0.00121
inf 0.00000
velocity
beta pmass denscorr
inf nan 0.00051
10.84877
betagamma log10bg
tau gamma icon_ev hnup (eV)
inf
inf
inf inf
94.840 0.708
Ok, I figured it out. The
parameterization gives VERY large numbers outside of the
acceptance. This can make the sos momentum go negative, which
caused the errors. By limiting the ssxp_tar used in the
parameterization, the errors went away. Here are the new lines in
s_satcorr.f.
if (spcentral.gt.0.96296) then
deltacorr = 46.729-137.21*spcentral+181.15*spcentral**2
>
-76.089*spcentral**3
else
deltacorr = 14.6369
endif
c Extra correction added for
pionct 3/14/06
benxp = ssxp_tar
if (benxp.gt.0.041) then
benxp=0.041
endif
if (benxp.lt.-0.041) then
benxp=-0.041
endif
ssdelta = ssdelta +
> (0.0199-0.34515*benxp+12.8196*benxp**2+212.49192*benxp**3)
> +(-0.0160165-0.277808*benxp-14.74545*benxp**2)*spcentral
>
+(-174.44034*exp(138.8557*(benxp-0.07550238))+0.1112637+6.1343792*benxp)
>
*(1.077034-1.56109*spcentral+0.387612*spcentral**3)*exp(0.475572*spcentral)
endif