code: def get_flux(self, energies):
  
          
  
          results = [component.shape(energies) for component in self.components.values()]
  
          return numpy.sum(results, 0)
docs: Get the total flux of this particle source at the given energies (summed over the components)