Question
Rpi can handle at most 4 DS18B20s, how come so few? How can greedy me do as many
as 60?
Short Answer
Update 2019nov20hkt1014

Update 2019jul26hkt2030
I have upgraded my Rpi3B+ stretch to Rpi4B buster.  I test the same Rpi3B+ program for Rpi4B /buster and found everything work smoothly.  The program v71 with updated documentation is listed in Appendix G.  I have deleted the old program in Appendix F, because this answer is again getting too long, hitting the forum answer's 32k words limit.
Update 2019jul17hkt1301
I connected 12 DS18B20 to GPIO2, and another 3 to GPIO3, and found
measurements stable.  I am using the I2C pins for now.  I hope to test
other pins after I finished testing Rpi4 and buster.

*** Begin testGetPrintSensorIdTempTimeList() **********************************************************
Time Testing = 2019-07-17 12:49 tlfong01 
Getting sensor temperatures, ...............
Begin printSensorIdTempTimeList(), ...
  ----------------------------------------------------------
  Sensor         Sensor    Temp                         Time
     Num             ID   deg C                             
  ----------------------------------------------------------
       1   01161cd97fee   28.12   2019-07-17 12:49:22.365111
       2   01161cda20ee   28.12   2019-07-17 12:49:23.325028
       3   0316a2795c11   27.81   2019-07-17 12:49:24.284948
       4   021313a052aa   28.56   2019-07-17 12:49:25.244955
       5   011319f99467   28.38   2019-07-17 12:49:26.204966
       6   000001d8e47e   28.44   2019-07-17 12:49:27.164934
       7   0316a27926c3   28.25   2019-07-17 12:49:28.124953
       8   01131a39c3e9   28.56   2019-07-17 12:49:29.085026
       9   0316a279d285   28.25   2019-07-17 12:49:30.044939
      10   0316a27907c7   27.69   2019-07-17 12:49:31.004974
      11   02131261bcaa   28.38   2019-07-17 12:49:31.965575
      12   01131a2edfe9   28.50   2019-07-17 12:49:32.924948
      13   0316a27938eb   27.69   2019-07-17 12:49:34.845093
      14   01131bb23b09   28.19   2019-07-17 12:49:35.804940
      15   000001d8e502   28.44   2019-07-17 12:49:36.764982
  ----------------------------------------------------------
End   printSensorIdTempTimeList().
Number of sensors tested = 15
Seconds lapsed           = 16.69
Seconds per sensor       = 1.11
*** End   testGetPrintSensorIdTempTimeList() **********************************************************

Long Answer to be shortened or removed
Now I am testing 12 DS18B20 devices, all connected to one GPIO pin 2, configuring as star.  Now I am finding that the syatem is not reliable, in that it takes longer to made one measurements, from about 1 second to up ta 1.5 seconds.  I need to find if using daisy chain configuration can improve the performance.
Testing 11 sensors OK.  All 11 sensors connected in parallel, power
3V3, all output pins connected to GPIO2 (Similarly to shown in the following post .  Test program list and sample output shown in Appendix B.

Now I am going to test 9 sensor modules, suing 2 or 3 GPIO pins.

Now I am testing 11 DS18B20 (7 with waterproof probes, two modules, two through hole chips), all on the same bus connected to GPIO2.  I found measurements stable.  See sample output and program listing in the appendix below.

Long winded Research Notes to shorten or delete later
I guess, DS18B20, using the 1-Wire method, is like the trouble making I2C, cannot have wires too long (preferably less than 30cm), pull up not too strong, and not too weak (total around 4k7)  capacitance less than 400pF etc.
For I2C, I had a very hard time and could achieve only placing 8 devices on the bus, but system still very unstable, often getting I/O error (due to impedance overloading, noise etc).
So for this 1-Wire DS18B20, I was humble and my goal was to get at most 8
devices on the 1-Wire bus. I started with onle one DS18B20, and added one
after another, until I used up all my 7 devices.  Surprisingly I found
measurements very stable.
Now let me show off my setup, results, and program.
Setup

I think I confused everybody by using I2C data pin for DS18B20.
Actually I can use any other pin.  The reason that I choose I2C Data
pin (GPIO 2) is that I know I2C CLK and DATA pins have a very strong
pull up of 1k8.  All other pins (only 90% sure) have pullups of 50k.
In other words, I don't need any extra 4k7.
References
DS18B20 Programmable Resolution 1-Wire Digital Thermometer Datasheet - Maxim
 
Guidelines for Reliable Long Line 1-Wire Networks, Tutorial 148 - Maxim
 
w1thermsensor DS18B20 Python Module 1.1.2 - pypi
 
Multiple DS18B20 1-Wire Temperature Sensors - rpi.org forum 2016dec09
 
DS18B20 Wiring Diagram - raspberry pi forum discussion
 
DS18B20 MAX DISTANCE -rpi forum geekinthesticks 2013maro6
 
Joining Multiple DS18B20 Sensors - rpi forum Darranking 2019jan14
 
1-Wire - Wikipedia
 
SunFounder DS18B20 Temperature Sensor Learning Kit
 
Rpi3B+ Python DS2482 I2C to 1-Wire Bridge Connecting 10+ DS18B20 Temperature Sensors Using 30 Metres Long CAT5e Cables
 
Appendices
Appendix A - Test Program 2019jul05
Appendix B - DS18B20 Test Program 2019jul11
Appendix C - DS18B20 12 devices testing notes
Appendix D/E - DS18B20 Test Program V50 2019jul14hkt1128
Listings too long hitting the StackExchange answer's 32k word limit, therefore removed.
Appendix F - Program listing of version 62  2019jul17hkt1259
Appendix G - Version 71 Rpi4B Throny python 3.7.3
Update 2019nov20hkt1031
(1) DS18B20 test71 was written for Rpi3B+ stretch 9 and IDLE python 3.5.3.  It is still compatible to Rpi4B buster python 3.7.3, without any modification.
(2) The program was developed in an almost pure functional programming and declarative programming paradigm, using lisp/scheme/LUA programming constructs, and OOP style without using OP constructs, only FP.  It automatically scales up for any number of DS18B20 sensors. For example, there is no C++ like switch/case, nested if-then-else, numerical for loops.  The same program can entertain more or less sensors without modifying any for loop range parameter, ie, "looping for list", not for "range(length ...) etc.
# DS18B20_test71 tlfong01 2019jul26hkt204
# Rpi3B+ Raspbian  9 stretch 2019jul15, IDLE   Python 3.5.3
# Rpi4B  Raspbian 10 buster  2019jul25, Throny Python 3.7.3
# python3-w1ThermSensor-1.0.5
*** Setup Notes ***
1. "$ sudo nano /boot/config.txt" to config w1thermsensor GPIO dtoverlay
2. GPIO 2/3 (I2C DATA/CLK) with built in 1k8 pull up is used.  For other
pins with weaker pullups, a 4k7 resistor is usually recommended.
3. dtoverlay=w1-gpio,gpiopin=2, dtoverlay=w1-gpio,gpiopin=3
4. For two or more GPIO pins, w1thermsensor auto detects connected sensors,
but might hang if GPIO is specified but no senors are connected.
And remember to reboot after making changes to dtoverlay settings.
6. Terminal mode sensor detection command = "$ ls /sys/bus/w1/devices"
*** Testing Results ***
1. Program V62 has been tested OK for 15 sensors (12 connected to GPIO2
and 3 to GPIO3.  All connecting cables are less than 1 metre long.
2. To test later - Rpi4, buster, long CAT5 cables, 5V power and logical
level shifter.
Program Summary and Testing Notes
1. The three datetime functions in this sections are mainly for demo and testing.  They
can be removed after development is completed.
2. The data structure "sensorIdList" contains one member for each sensor,
Each member consists of the id string of the sensor.
3. The get/print/getPrint SensorIdList is mainly used for detecting sensors which is
much quicker than the very slow sensor measurement process of around 1 second long.
4. The data structure "sensorIdTempTimeList", contains one member for each sensor.
Each member consists of the id, temperature read, time read of the sensor.
5. The getSensorIdTempTime fucntion is very slow, about 1 second to execute.
6. The main tests can be executed all at once, or comment/uncomment to select lists.
from time          import sleep
from datetime      import datetime
from datetime      import timedelta
from w1thermsensor import W1ThermSensor
*** Date Time Functions ***
def printDateTimeNow():
    timeNowLong  = datetime.now()
    timeNowShort = str(timeNowLong)[0:16]
    print('\n    timeNowLong         =',  timeNowLong)
    print('    timeNowShort        =', timeNowShort)
    return
def printElapsedTime():
    startTime = datetime.now()
    sleep(2)
    finishTime = datetime.now()
timeElapsed    = finishTime - startTime
secondsElapsed = timeElapsed.seconds
microSecondsElapsed = (secondsElapsed * 1000000) + timeElapsed.microseconds
milliSecondsElapsed = microSecondsElapsed / 1000
print('\n    Start Time          =', startTime)
print('    Finish Time         =', finishTime)
print('\n    timeElapsed         =', timeElapsed)
print('    secondsElapsed      =', secondsElapsed)
print('    milliSecondsElapsed =', '%d' %milliSecondsElapsed)    
print('    microSecondsElapsed =', microSecondsElapsed, '\n')
def testDateTime():
    print('\n*** Begin Test python datetime module, ... *****************************')
    printDateTimeNow()
    printElapsedTime()
    print('* End Testing python datetime module. **********************************')
    return
*** DS18B20 Functions ***
*** SensorIdList Functions ***
def getSensorIdList():
    sensorIdList = []
    for sensor in W1ThermSensor.get_available_sensors():
        sensorIdList.append(sensor.id)
    return sensorIdList
def printSensorIdList(sensorIdList):
    print('    ------------------------------------------------------------------------')
    print('    Sensor        ', '   Sensor')
    print('       Num        ', '       ID')
    print('    ------------------------------------------------------------------------')
    sensorNum = 1
    for sensorId in sensorIdList:
        if sensorNum < 10:
            print('        ', sensorNum, '    ', sensorId)
        else:
            print('       ', sensorNum, '    ', sensorId)
        sensorNum = sensorNum + 1
    print('    ------------------------------------------------------------------------')
def getPrintSensorIdList():
    sensorIdList = getSensorIdList()
    printSensorIdList(sensorIdList)
    return
def testGetPrintSensorIdList():
    print('\n*** Begin testGetPrintSensorIdList() ****************************************************************\n')
    getPrintSensorIdList()
    print('\n* End   testGetPrintSensorIdList() ******************************************************************\n')
    return
*** SensorIdTempTimeList Functions ***
def getSensorIdTempTimeList():
    sensorIdTempTimeList = []
    print('    Getting sensor temperatures, ', end = '')
    for sensor in W1ThermSensor.get_available_sensors():
        idTempTimeList = (sensor.id, sensor.get_temperature(), datetime.now())
        sensorIdTempTimeList.append(idTempTimeList)
        print('.', end = '')
    return sensorIdTempTimeList
def printSensorIdTempTimeList(sensorIdTempTimeList):
    print('\n\n    Begin printSensorIdTempTimeList(), ...\n')
    print('      ----------------------------------------------------------')
    print('      Sensor', '        Sensor', '   Temp', '                        Time')
    print('         Num', '            ID', '  deg C', '                            ')
    print('      ----------------------------------------------------------')
sensorNum = 1
for idTempTimeList in sensorIdTempTimeList:
    if sensorNum < 10:
        print('          ', sensorNum, end = '')
    else:
        print('         ',  sensorNum, end = '')  
    print('  ', idTempTimeList[0], end = '  ')
    print('', '%.2f' %(idTempTimeList[1]), end = '  ')
    print('', idTempTimeList[2])
    sensorNum = sensorNum + 1
print('      ----------------------------------------------------------')
print('\n    End   printSensorIdTempTimeList().')
return    
def testGetPrintSensorIdTempTimeList():
    print('\n*** Begin testGetPrintSensorIdTempTimeList() **********************************************************\n')
    startTime = datetime.now()
    print('    Time Testing =', str(startTime)[0:16], 'tlfong01 \n')
sensorIdTempTimeList = getSensorIdTempTimeList()
printSensorIdTempTimeList(sensorIdTempTimeList)
finishTime = datetime.now()
timeElapsed    = finishTime - startTime
secondsElapsed = timeElapsed.seconds
microSecondsElapsed = (secondsElapsed * 1000000) + timeElapsed.microseconds
sensorNum = len(sensorIdTempTimeList)
print('\n    Number of sensors tested =', sensorNum)
print('    Seconds lapsed           =', '%.2f' %(microSecondsElapsed / 1000000))
print('    Seconds per sensor       =', '%.2f' %((microSecondsElapsed / 1000000) / sensorNum))
print('\n*** End   testGetPrintSensorIdTempTimeList() **********************************************************\n')
return  
*** Main Test ***
testDateTime()
testGetPrintSensorIdList()
testGetPrintSensorIdTempTimeList()
'''
Python 3.7.3 (/usr/bin/python3)
>>> %cd /home/pi/Python_Programs/test1321/ds18b20
>>> %Run ds18b20_test70_2019jul2501.py
*** Begin Test python datetime module, ... *******************************
timeNowLong         = 2019-07-26 20:01:22.402122
timeNowShort        = 2019-07-26 20:01
Start Time          = 2019-07-26 20:01:22.402724
Finish Time         = 2019-07-26 20:01:24.404806
timeElapsed         = 0:00:02.002082
secondsElapsed      = 2
milliSecondsElapsed = 2002
microSecondsElapsed = 2002082 
*** End Testing python datetime module. **********************************
*** Begin testGetPrintSensorIdList() ******************************************************************
------------------------------------------------------------------------
Sensor            Sensor
   Num                ID
------------------------------------------------------------------------
     1      01161cd97fee
     2      01161cda20ee
     3      0316a2795c11
     4      021313a052aa
     5      011319f99467
     6      000001d8e47e
     7      0316a27926c3
     8      01131a39c3e9
     9      0316a279d285
    10      0316a27907c7
    11      02131261bcaa
    12      01131a2edfe9
    13      0316a27938eb
    14      01131bb23b09
    15      000001d8e502
------------------------------------------------------------------------
*** End   testGetPrintSensorIdList() ******************************************************************
*** Begin testGetPrintSensorIdTempTimeList() **********************************************************
Time Testing = 2019-07-26 20:01 tlfong01 
Getting sensor temperatures, ...............
Begin printSensorIdTempTimeList(), ...
  ----------------------------------------------------------
  Sensor         Sensor    Temp                         Time
     Num             ID   deg C                             
  ----------------------------------------------------------
       1   01161cd97fee   29.81   2019-07-26 20:01:25.373759
       2   01161cda20ee   29.75   2019-07-26 20:01:26.253930
       3   0316a2795c11   29.38   2019-07-26 20:01:27.133707
       4   021313a052aa   30.00   2019-07-26 20:01:28.014356
       5   011319f99467   29.88   2019-07-26 20:01:28.903601
       6   000001d8e47e   29.81   2019-07-26 20:01:29.853169
       7   0316a27926c3   29.88   2019-07-26 20:01:30.733717
       8   01131a39c3e9   29.94   2019-07-26 20:01:31.613679
       9   0316a279d285   29.62   2019-07-26 20:01:32.493157
      10   0316a27907c7   29.12   2019-07-26 20:01:33.373696
      11   02131261bcaa   29.88   2019-07-26 20:01:34.313676
      12   01131a2edfe9   29.62   2019-07-26 20:01:35.213715
      13   0316a27938eb   29.50   2019-07-26 20:01:36.093166
      14   01131bb23b09   29.44   2019-07-26 20:01:36.973689
      15   000001d8e502   29.94   2019-07-26 20:01:37.853752
  ----------------------------------------------------------
End   printSensorIdTempTimeList().
Number of sensors tested = 15
Seconds lapsed           = 13.48
Seconds per sensor       = 0.90
*** End   testGetPrintSensorIdTempTimeList() **********************************************************
>>> 
'''
*** End of Sample Output ***
*** End of Program Documentaton ***