diff --git a/fermi_stacking/analyze_results/AnalyzeResults.py b/fermi_stacking/analyze_results/AnalyzeResults.py
index 969594b..85d85f6 100755
--- a/fermi_stacking/analyze_results/AnalyzeResults.py
+++ b/fermi_stacking/analyze_results/AnalyzeResults.py
@@ -539,7 +539,7 @@ def calc_upper_limit(self,srcname,ul_emin,ul_emax,comp_list=[0,1,2,3],mult_lt=Fa
         ul_emin : float
             Lower energy bound for UL calculation in MeV.
         ul_emax : float
-            Upper energy bound for UL calculation.
+            Upper energy bound for UL calculation in MeV.
         comp_list : list, optional
             List of components to add to the SummedLikelihood object 
             for the JLA (default is for typical JLA with 4 components). 
@@ -549,7 +549,11 @@ def calc_upper_limit(self,srcname,ul_emin,ul_emax,comp_list=[0,1,2,3],mult_lt=Fa
         mult_lt : bool, optional
             If using lt cubes for each component, set to True (default is False, 
             for single lt cube. 
-	
+
+        Note 
+        ----
+        If getting 'IndexError', try moving the lower energy bound slightly below the bin. 
+
         Returns
         -------
         float, float
@@ -567,7 +571,7 @@ def calc_upper_limit(self,srcname,ul_emin,ul_emax,comp_list=[0,1,2,3],mult_lt=Fa
                 if each > 9:
                     print("***ERROR***")
                     print("Looks like you have a complex analysis!")
-                    print("More than 9 components is not currently supported.")
+                    print("More than 10 components is not currently supported.")
                     print("You'll need to add more definitions to the source code.")
                     print()
                     sys.exit()
@@ -591,8 +595,8 @@ def calc_upper_limit(self,srcname,ul_emin,ul_emax,comp_list=[0,1,2,3],mult_lt=Fa
         my_expCube = self.ltcube
         
         # Analysis selections (fixed for now):
-        irfs = "P8R3_SOURCE_V2"
-        optimizer = "Minuit" # Minuit or NewMinuit
+        irfs = self.irfs
+        optimizer = "NewMinuit" # Minuit or NewMinuit
         conf = BinnedConfig(edisp_bins=-1) #need to account for energy dispersion!
         
         if self.JLA == True:
diff --git a/setup.py b/setup.py
index 829f771..b9c2207 100644
--- a/setup.py
+++ b/setup.py
@@ -4,7 +4,7 @@
 # Setup:
 setup(
     name='fermi_stacking',
-    version="v0.1.2-alpha",
+    version="v0.1.3-alpha",
     url='https://github.com/ckarwin/Fermi_Stacking_Analysis.git',
     author='Chris Karwin',
     author_email='christopher.m.karwin@nasa.gov',