anchor.diff:  !                  endif
   !CB
   
  -      if(n_slices.gt.one .and. totals.gt.smin_slices .and. totals.lt.smax_slices .and. &
  +      if(n_slices.gt.1 .and. totals.gt.smin_slices .and. totals.lt.smax_slices .and. &
    &      (cdb_cNameUC(icoll)(1:4).eq.'TCSG' .or. cdb_cNameUC(icoll)(1:3).eq.'TCP' .or. cdb_cNameUC(icoll)(1:4).eq.'TCLA'.or. &
    &       cdb_cNameUC(icoll)(1:3).eq.'TCT' .or. cdb_cNameUC(icoll)(1:4).eq.'TCLI'.or. cdb_cNameUC(icoll)(1:4).eq.'TCL.'.or.  &
   !     RB: added slicing of TCRYO as well
  
positive.code: #endif
  
    case(16)
      if(nSplit > 0)  coll_db = lnSPlit(1)
      if(nSplit > 1)  call chr_cast(lnSPlit(2), ibeam, iErr)
  
    case(17)
      if(nSplit > 0)  call chr_cast(lnSPlit(1), dowritetracks,iErr)
      if(nSplit > 1)  call chr_cast(lnSPlit(2), cern,         iErr)
      if(nSplit > 2)  castordir = lnSPlit(3)
      if(nSplit > 3)  call chr_cast(lnSPlit(4), jobnumber,    iErr)
      if(nSplit > 4)  call chr_cast(lnSPlit(5), sigsecut2,    iErr)
      if(nSplit > 5)  call chr_cast(lnSPlit(6), sigsecut3,    iErr)
  
    case default
      write(lout,"(a,i0,a)") "COLL> ERROR Unexpected line ",iLine," encountered."
  
negative.code:         do jjj=1,n_slices+1
            x_sl(jjj) = (jjj-1) * c_length / real(n_slices,fPrec)
  
            y1_sl(jjj) = jaw_fit(1,1) + jaw_fit(1,2)*x_sl(jjj) + jaw_fit(1,3)/c_length*(x_sl(jjj)**2) +           &
   &                           jaw_fit(1,4)*(x_sl(jjj)**3) + jaw_fit(1,5)*(x_sl(jjj)**4) + jaw_fit(1,6)*(x_sl(jjj)**5)
  
            y2_sl(jjj) = -one * (jaw_fit(2,1) + jaw_fit(2,2)*x_sl(jjj) + jaw_fit(2,3)/c_length*(x_sl(jjj)**2) +   &
   &                           jaw_fit(2,4)*(x_sl(jjj)**3) + jaw_fit(2,5)*(x_sl(jjj)**4) + jaw_fit(2,6)*(x_sl(jjj)**5))
          end do
  
  !       Apply the slicing scaling factors (ssf's):
  !       CB:10-2007 coordinates rotated of the tilt
          do jjj=1,n_slices+1
            y1_sl(jjj) = jaw_ssf(1) * y1_sl(jjj)
            y2_sl(jjj) = jaw_ssf(2) * y2_sl(jjj)
  ! CB code
            x1_sl(jjj) = x_sl(jjj) *cos_mb(db_tilt(icoll,1))-y1_sl(jjj)*sin_mb(db_tilt(icoll,1))
            x2_sl(jjj) = x_sl(jjj) *cos_mb(db_tilt(icoll,2))-y2_sl(jjj)*sin_mb(db_tilt(icoll,2))
  