The obsutils.sim() function does not correctly handle stacked observations. It appends stacked response cubes to a counts cube in case that multiple observations have been simulated and stacked analysis was requested.
The stacked response cube computation is done using the obsutils.get_stacked_response() function.
The function obsutils.set_ts_model() has been added to setup a model for TS determination. This is used by some of the cscripts, and avoids code duplication.
Finally, the cstsdist and cspull scripts have been refactored.
