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