anchor.diff:              a_tmp2 = y2_sl(jjj+1)
             end if
   
  -!!     Write down the information on slice centre and offset
  -!                     if (firstrun) then
  -!                        write(*,*) 'Processing slice number ',jjj,
  -!     &                       ' of ',n_slices,' for the collimator ',
  -!     &                       cdb_cNameUC(icoll)
  -!                        write(*,*) 'Aperture [m]= ',
  -!     &                       a_tmp1 - a_tmp2
  -!                        write(*,*) 'Offset [m]  = ',
  -!     &                       0.5 * ( a_tmp1 + a_tmp2 )
  -!                     endif
  -!!
   !     Be careful! the initial tilt must be added!
   !     We leave it like this for the moment (no initial tilt)
   !         c_tilt(1) = c_tilt(1) + angle1(jjj)
  
positive.code:      &                 onesided, flukaname, secondary, 1, nabs_type)
  #endif
        end if !if (n_slices.gt.one .and.
      end if !if(cdb_cNameUC(icoll)(1:4).eq.'COLM') then
    end if !if (found) then
  end subroutine collimate_do_collimator
  
  !>
  
negative.code:   use hdf5_output
  #endif
  
    implicit none
  
    logical onesided,hit
  ! integer nprim,filel,mat,nev,j,nabs,nhit,np,icoll,nabs_tmp
  