anchor.diff:  
   ! Log output energy + nucleons as per the FLUKA coupling
       nnuc1       = nnuc1 + naa(j)                          ! outcoming nucleons
  -    ien1        = ien1  + rcp(j) * c1e3                   ! outcoming energy
  +    ien1        = ien1  + (rcp(j) * c1e3)                 ! outcoming energy
   
   ! Fix hits
   ! if(part_hit_pos(j) .eq.ie .and. part_hit_turn(j).eq.iturn)
  
positive.code:       call root_EnergyDeposition(icoll, nnuc0-nnuc1,c1m3*(ien0-ien1))
      end if
  #endif
      write(unit208,"(2(i5,1x),e24.16)") icoll, (nnuc0-nnuc1), c1m3*(ien0-ien1)
      flush(unit208)
    end if
  
  
negative.code:   call k2coll_merlinInit
  #endif
  
    ! Open the edep file
    if(unit208 == -1) then
      call f_requestUnit(fort208,unit208)
      call f_open(unit=unit208,file=fort208,formatted=.true.,mode="w",status="replace")
    end if
  
  #ifdef G4COLLIMATION
  
    if(n_slices /= 0) then
      write(lerr,"(a)") "COLL> ERROR Cannot use jaw fit in G4COLLIMATION version of SixTrack"
      call prror
    end if
  
    ! This function lives in the G4Interface.cpp file in the g4collimat folder
  