anchor.diff:                scatterhit(j).eq.8         ) .and. &
                (xv1(j).lt.99.0_fPrec .and. xv2(j).lt.99.0_fPrec).and.&
   !GRD HERE WE APPLY THE SAME KIND OF CUT THAN THE SIGSECUT PARAMETER
  -             ((((xv1(j)*c1m3)**2 / (tbetax(ie)*myemitx0_collgap)) .ge. real(sigsecut2,fPrec)) .or. &
  -             (((xv2(j)*c1m3)**2  / (tbetay(ie)*myemity0_collgap)) .ge. real(sigsecut2,fPrec)) .or. &
  +             ((((xv1(j)*c1m3)**2 / (tbetax(ie)*myemitx0_collgap)) .ge. sigsecut2) .or. &
  +             (((xv2(j)*c1m3)**2  / (tbetay(ie)*myemity0_collgap)) .ge. sigsecut2) .or. &
                (((xv1(j)*c1m3)**2  / (tbetax(ie)*myemitx0_collgap)) + &
  -             ((xv2(j)*c1m3)**2   / (tbetay(ie)*myemity0_collgap)) .ge. sigsecut3)) ) &
  -             then
  +             ((xv2(j)*c1m3)**2   / (tbetay(ie)*myemity0_collgap)) .ge. sigsecut3)) ) then
   
               xj  = (xv1(j)-torbx(ie))  /c1e3
               xpj = (yv1(j)-torbxp(ie)) /c1e3
  
positive.code: subroutine collimate_postInput(gammar)
  
    use crcoall
  
    real(kind=fPrec), intent(in) :: gammar
  
  
negative.code:   case("DO_RADIAL")
      if(nSplit /= 4) then
        write(lerr,"(a,i0)") "COLL> ERROR DO_RADIAL expects 3 values, got ",nSplit-1
        write(lerr,"(a)")    "COLL>       DO_RADIAL true|false amp smear"
        iErr = .true.
        return
      end if
  