Amorphous fix

hexrd/wppf/WPPF.py

@@ -1573,18 +1573,28 @@ def _set_params_vals_to_class(self, params, init=False, skip_phases=False):
                 self.amorphous_model.center = center
                 self.amorphous_model.fwhm = fwhm
 
+    @property
+    def total_area(self):
+        tth, intensity   = self.spectrum_sim.data
+        _, background    = self.background.data
+        total_intensity = intensity-background
+        '''put some guard rails around the total intensity
+        to protect against nans in the values
+        '''
+        mask = np.isnan(total_intensity)
+        sum_area = np.trapz(total_intensity[~mask], tth[~mask])
+        return sum_area
+
+
     @property
     def DOC(self):
         if self.amorphous_model is None:
             return 1.
         else:
-            tth, intensity   = self.spectrum_sim.data
-            _, background    = self.background.data
-            total_intensity = intensity-background
             amorphous_area = \
                 self.amorphous_model.integrated_area
-            return 1. - amorphous_area/np.trapz(
-                        total_intensity, tth)
+            total_area = self.total_area
+            return 1. - amorphous_area/total_area
 
 
 def _nm(x):

hexrd/wppf/amorphous.py

@@ -312,8 +312,8 @@ def fwhm(self, val):
     def amorphous_lineout(self):
         if self.model_type == "experimental":
 
+            lo = np.zeros_like(self.tth_list)
             for key in self.shift:
-                lo = np.zeros_like(self.tth_list)
                 smooth_model_data = gaussian_filter(
                                    self.model_data[key],
                                    self.smoothing