code: def _pangledAnglet(self,t,angleperp,dangle,smallest):
          
          if smallest: eigIndx= 0
          else: eigIndx= 1      
          if isinstance(angleperp,(int,float,numpy.float32,numpy.float64)):
              angleperp= numpy.array([angleperp])
              t= numpy.array([t])
          out= numpy.zeros_like(angleperp)
          tindx= t < self._tdisrupt
          out[tindx]=\
              numpy.exp(-0.5*angleperp[tindx]**2.\
                             /(t[tindx]**2.*self._sortedSigOEig[eigIndx]+self._sigangle2))/\
                             numpy.sqrt(t[tindx]**2.*self._sortedSigOEig[eigIndx]+self._sigangle2)\
                             *self.ptdAngle(t[t < self._tdisrupt],dangle)
          return out
docs: p(angle_perp|angle_par,time)