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:     call f_requestUnit(coll_impactFile,    coll_impactUnit)
      call f_requestUnit(coll_flukImpFile,   coll_flukImpUnit)
      call f_requestUnit(coll_flukImpAllFile,coll_flukImpAllUnit)
  
      call f_open(unit=coll_allImpactUnit, file=coll_allImpactFile, formatted=.true.,mode="w")
      call f_open(unit=coll_allAbsorbUnit, file=coll_allAbsorbFile, formatted=.true.,mode="w")
  
negative.code:   character(len=13), parameter :: coll_twissLikeFile  = "twisslike.out"
    character(len=17), parameter :: coll_sigmaSetFile   = "sigmasettings.out"
    character(len=16), parameter :: coll_settingsFile   = "collsettings.dat"
    character(len=16), parameter :: coll_jawProfileFile = "jaw_profiles.dat"
  
    integer, private, save :: coll_survivalUnit   = -1
    integer, private, save :: coll_gapsUnit       = -1
  