old_contents: """
  ==================
  PMT Positions.
  ==================
  
  Plot the PMT directions for a given DOM.
  
  """
  # Author: Tamas Gal <tgal@km3net.de>
  # License: BSD-3
  from mpl_toolkits.mplot3d import axes3d
  import matplotlib.pyplot as plt
  from matplotlib import rcParams
  
  import km3pipe as kp
  import km3pipe.style
  km3pipe.style.use('km3pipe')
  km3pipe.style.use('talk')
  
  DETX = 'data/orca_115strings_av23min20mhorizontal_18OMs_alt9mvertical_v1.detx'
  
  det = kp.hardware.Detector(DETX)
  
  #####################################################
  # Alternatively, you can use the `det_id` to retrieve the geometry from the DB.
  
  # det = kp.hardware.Detector(det_id=29)
  
  #####################################################
  # Let's take the first DOM ID
  dom_id = det.dom_ids[0]
  
  #####################################################
  # ...and get the table of the PMTs in the chosen DOM:
  
  pmts = det.pmts[det.pmts.dom_id == dom_id]
  
  #####################################################
  # Some cosmetics for the plot...
  rcParams['axes.titlepad'] = 25
  rcParams['axes.labelpad'] = 20
  
  #####################################################
  # The `quiver` function can directly plot the PMT data, since those are
  # stored as numpy arrays.
  fig = plt.figure()
  ax = fig.gca(projection='3d')
  
  ax.quiver(pmts.pos_x, pmts.pos_y, pmts.pos_z,
            pmts.dir_x, pmts.dir_y, pmts.dir_z,
            length=0.1)
  
  ax.set_xlabel("x [m]")
  ax.set_ylabel("y [m]")
  ax.set_zlabel("z [m]")
  
  ax.set_title("PMT directions on DOM {}".format(dom_id))
  
new_contents: """
  ==================
  PMT Directions
  ==================
  
  Plot the PMT directions for a given DOM.
  
  """
  # Author: Tamas Gal <tgal@km3net.de>
  # License: BSD-3
  from mpl_toolkits.mplot3d import axes3d
  import matplotlib.pyplot as plt
  from matplotlib import rcParams
  
  import km3pipe as kp
  import km3pipe.style
  km3pipe.style.use('km3pipe')
  km3pipe.style.use('talk')
  
  DETX = 'data/orca_115strings_av23min20mhorizontal_18OMs_alt9mvertical_v1.detx'
  
  det = kp.hardware.Detector(DETX)
  
  #####################################################
  # Alternatively, you can use the `det_id` to retrieve the geometry from the DB.
  
  # det = kp.hardware.Detector(det_id=29)
  
  #####################################################
  # Let's take the first DOM ID
  dom_id = det.dom_ids[0]
  
  #####################################################
  # ...and get the table of the PMTs in the chosen DOM:
  
  pmts = det.pmts[det.pmts.dom_id == dom_id]
  
  #####################################################
  # Some cosmetics for the plot...
  rcParams['axes.titlepad'] = 25
  rcParams['axes.labelpad'] = 20
  
  #####################################################
  # The `quiver` function can directly plot the PMT data, since those are
  # stored as numpy arrays.
  fig = plt.figure()
  ax = fig.gca(projection='3d')
  
  ax.quiver(pmts.pos_x, pmts.pos_y, pmts.pos_z,
            pmts.dir_x, pmts.dir_y, pmts.dir_z,
            length=0.1)
  
  ax.set_xlabel("x [m]")
  ax.set_ylabel("y [m]")
  ax.set_zlabel("z [m]")
  
  ax.set_title("PMT directions on DOM {}".format(dom_id))
  
subject: Fix title