anchor.diff:  
     if(myenom == zero) then
       write(lerr,"(a)") "COLL> ERROR Beam energy cannot be zero"
  +    call prror
  +  end if
  +
  +  if(cdb_fileName == " ") then
  +    write(lerr,"(a)") "COLL> ERROR No collimator database file specified"
       call prror
     end if
   
  
positive.code:               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)) + &
               ((xv2(j)*c1m3)**2   / (tbetay(ie)*myemity0_collgap)) .ge. sigsecut3)) ) &
               then
  
              xj  = (xv1(j)-torbx(ie))  /c1e3
              xpj = (yv1(j)-torbxp(ie)) /c1e3
  
negative.code:              other(j)     .eq. 4 .or. &
               scatterhit(j).eq. 8       ) .and. &
               (xv1(j).lt.99.0_fPrec .and. xv2(j).lt.99.0_fPrec) .and. &
               ((((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
  
            xj  = (xv1(j)-torbx(ie)) /c1e3
            xpj = (yv1(j)-torbxp(ie))/c1e3
  