anchor.diff:    character(len=:), allocatable, public, save :: cdb_famName(:)     ! Family name
     real(kind=fPrec), allocatable, public, save :: cdb_famNSig(:)     ! Family sigma
     real(kind=fPrec), allocatable, public, save :: cdb_famNSigOrig(:) ! Family sigma (original value from DB)
  -  integer,          allocatable, public, save :: cdb_famStage(:)    ! Family type
  +  integer,          allocatable, public, save :: cdb_famStage(:)    ! Family stage
   
     ! Element Map
     integer,          allocatable, public, save :: cdb_elemMap(:)     ! Map from single elements to DB
  
positive.code:   character(len=cdb_fNameLen) famName
    character(len=mNameLen) collDummy
    logical cErr, fExists
    integer j, dbUnit, ioStat, iLine, famID, matID, collType
  
    cErr = .false.
  
  
negative.code: !  Updated: 2019-10-16
  !  Checks the stage of the collimator based on LHC naming convention (for support with old format)
  ! ================================================================================================ !
  subroutine cdb_getCollStage(inElem, collStage)
  
    use string_tools
  
    character(len=mNameLen), intent(in)  :: inElem
    integer,                 intent(out) :: collStage
  
    character(len=mNameLen) elemName
  
    collStage = cdb_stgOther
    elemName  = chr_toLower(inElem)
  
    if(elemName(1:3) == "tcp") then
      collStage = cdb_stgPrimary
    else if(elemName(1:3) == "tcs") then
      collStage = cdb_stgSecondary
    else if(elemName(1:3) == "tcl" .or. elemName(1:3) == "tct" .or. &
            elemName(1:3) == "tcd" .or. elemName(1:3) == "tdi") then
      collStage = cdb_stgTertiary
    end if
  
  end subroutine cdb_getCollStage
  