anchor.diff:  !  g4_physics = 0
   
     call g4_collimation_init(e0, rnd_seed, g4_recut, g4_aecut, g4_rcut, g4_rangecut_mm, g4_v0, trim(g4_phys_str), &
  -&                          g4_debug, g4_keep_stable, g4_edep)
  +    g4_debug, g4_keep_stable, g4_edep)
   #endif
   
     write (lout,"(a)") ""
  
positive.code:         ! first time particle hits collimator: entering jaw
          linside(j)=.true.
          if(dowrite_impact) then
             if ( tiltangle.gt.zero ) then
                x_Dump=(x+c_aperture/two+tiltangle*sp)*mirror+c_offset
             else
                x_Dump=(x+c_aperture/two+tiltangle*(sp-c_length))*mirror+c_offset
             end if
             xpDump=(xp+tiltangle)*mirror
             y_Dump=z
             ypDump=zp
             s_Dump=sp+real(j_slices-1,fPrec)*c_length
             write(profiling_jaw_unit,'(3(1x,i7),5(1x,e17.9),1x,i1)') &
                    icoll,iturn,name(j),x_Dump,xpDump,y_Dump,ypDump,s_Dump,1
          end if
        end if
  !JUNE2005
  
negative.code: 
    call f_close(outlun)
    call f_close(coll_gapsUnit)
    call f_close(coll_jawProfileUnit)
  
    if(dowritetracks) then
      call f_close(coll_tracksUnit)
  #ifdef HDF5
      if(h5_writeTracks2) call h5tr2_finalise
  #endif
    end if
  
    if(do_select) then
      call f_close(coll_ellipseUnit)
    end if
  
    if(dowrite_impact) then
      call f_close(coll_allImpactUnit)
  