Commit 9dfda7afd4372d679ca5534ab6f04020d7c64fbc

Authored by Michael Mayer
Committed by Jürgen Knödlseder
1 parent 49bccaea

Added seed parameter to cspull

cscripts/cspull.par
@@ -44,6 +44,8 @@ rad, r, h, 5.0,0,180, "Radius of ROI (deg)" @@ -44,6 +44,8 @@ rad, r, h, 5.0,0,180, "Radius of ROI (deg)"
44 pattern, s, h, single,single|four,, "Observation pattern" 44 pattern, s, h, single,single|four,, "Observation pattern"
45 offset, r, h, 1.5,,, "Observation pattern offset (deg)" 45 offset, r, h, 1.5,,, "Observation pattern offset (deg)"
46 profile, b, h, no,,, "Use likelihood profile method for errors?" 46 profile, b, h, no,,, "Use likelihood profile method for errors?"
  47 +seed, i, h, 1,,, "Initial random number generator seed"
  48 +
47 49
48 # 50 #
49 # Standard parameters 51 # Standard parameters
cscripts/cspull.py
@@ -64,6 +64,7 @@ class cspull(ctools.cscript): @@ -64,6 +64,7 @@ class cspull(ctools.cscript):
64 self._binsz = 0.0 64 self._binsz = 0.0
65 self._pattern = "single" 65 self._pattern = "single"
66 self._enumbins = 0 66 self._enumbins = 0
  67 + self._seed = 1
67 self._chatter = 2 68 self._chatter = 2
68 69
69 # Initialise observation container from constructor arguments. 70 # Initialise observation container from constructor arguments.
@@ -116,6 +117,7 @@ class cspull(ctools.cscript): @@ -116,6 +117,7 @@ class cspull(ctools.cscript):
116 self._ntrials = self["ntrials"].integer() 117 self._ntrials = self["ntrials"].integer()
117 self._edisp = self["edisp"].boolean() 118 self._edisp = self["edisp"].boolean()
118 self._offset = self["offset"].real() 119 self._offset = self["offset"].real()
  120 + self._seed = self["seed"].integer()
119 self._chatter = self["chatter"].integer() 121 self._chatter = self["chatter"].integer()
120 122
121 # Query some parameters 123 # Query some parameters
@@ -410,8 +412,8 @@ class cspull(ctools.cscript): @@ -410,8 +412,8 @@ class cspull(ctools.cscript):
410 # Loop over trials 412 # Loop over trials
411 for seed in range(self._ntrials): 413 for seed in range(self._ntrials):
412 414
413 - # Make a trial  
414 - result = self._trial(seed) 415 + # Make a trial and add initial seed
  416 + result = self._trial(seed + self._seed)
415 417
416 # Write out result immediately 418 # Write out result immediately
417 if seed == 0: 419 if seed == 0: