diff --git a/doc/source/reference_manual/cspull.rst b/doc/source/reference_manual/cspull.rst
index 7d0d69e..8e96bf2 100644
--- a/doc/source/reference_manual/cspull.rst
+++ b/doc/source/reference_manual/cspull.rst
@@ -26,9 +26,15 @@ General parameters
 ``outfile [file]``
     ASCII file containing the individual pull values.
  	 	 
-``ntrials [integer]``
-    Number of samples for generating the pull distribution.
- 	 	 
+``expcube [file]``
+    Exposure cube file (only needed for stacked analysis).
+
+``psfcube [file]``
+    PSF cube file (only needed for stacked analysis).
+
+``bkgcube [file]``
+    Background cube file (only needed for stacked analysis).
+
 ``caldb [string]``
     Calibration database.
  	 	 
@@ -38,6 +44,9 @@ General parameters
 ``(edisp = no) [boolean]``
     Apply energy dispersion to response computation.
 
+``ntrials [integer]``
+    Number of samples for generating the pull distribution.
+ 	 	 
 ``ra [real]``
     Right Ascension of CTA pointing (J2000, in degrees).
  	 	 
diff --git a/scripts/cspull.py b/scripts/cspull.py
index f6cb8fb..9abf3ae 100755
--- a/scripts/cspull.py
+++ b/scripts/cspull.py
@@ -98,8 +98,11 @@ class cspull(ctools.cscript):
             pars.append(gammalib.GApplicationPar("outfile","f","a","pull.dat","","","Output file name"))
             pars.append(gammalib.GApplicationPar("expcube","s","a","NONE","","","Exposure cube file (only needed for stacked analysis)"))
             pars.append(gammalib.GApplicationPar("psfcube","s","a","NONE","","","PSF cube file (only needed for stacked analysis)"))
+            pars.append(gammalib.GApplicationPar("bkgcube","s","a","NONE","","","Background cube file (only needed for stacked analysis)"))
             pars.append(gammalib.GApplicationPar("caldb","s","a","dummy","","","Calibration database"))
             pars.append(gammalib.GApplicationPar("irf","s","a","cta_dummy_irf","","","Instrument response function"))
+            pars.append(gammalib.GApplicationPar("edisp","b","h","no","","","Apply energy dispersion?"))
+            pars.append(gammalib.GApplicationPar("ntrials","i","a","10","","","Number of trials"))
             pars.append(gammalib.GApplicationPar("ra","r","a","83.6331","0","360","RA of pointing (deg)"))
             pars.append(gammalib.GApplicationPar("dec","r","a","22.0145","-90","90","Dec of pointing (deg)"))
             pars.append(gammalib.GApplicationPar("emin","r","a","0.1","","","Lower energy limit (TeV)"))
@@ -113,8 +116,6 @@ class cspull(ctools.cscript):
             pars.append(gammalib.GApplicationPar("rad","r","h","5.0","0","180","Radius of ROI (deg)"))         
             pars.append(gammalib.GApplicationPar("pattern","s","h","single","","","Observation pattern (single/four)"))
             pars.append(gammalib.GApplicationPar("offset","r","h","1.5","","","Observation pattern offset (deg)"))
-            pars.append(gammalib.GApplicationPar("ntrials","i","a","10","","","Number of trials"))
-            pars.append(gammalib.GApplicationPar("edisp","b","h","no","","","Apply energy dispersion?"))
             pars.append_standard()
             pars.append(gammalib.GApplicationPar("logfile","f","h","cspull.log","","","Log filename"))
             pars.save(parfile)