Fabio Iocco / ctools

Commit abb14b993d7203747229c72cc341b384d718893d

Authored by Jürgen Knödlseder
1 parent d150779c

Pre-select events in cslightcrv

Inline Side-by-side
Showing 1 changed file with 26 additions and 1 deletions
cscripts/cslightcrv.py
  View file @abb14b9
... ... @@ -511,7 +511,32 @@ class cslightcrv(ctools.cscript):
511 511 # Write observation into logger
512 512 if self._logTerse():
513 513 self._log('\n')
514   - self._log.header1(gammalib.number('Observation',len(self._obs)))
  514 + self._log.header1(gammalib.number('Input observation',len(self._obs)))
  515 + self._log(str(self._obs))
  516 + self._log('\n')
  517 +
  518 + # Get time boundaries
  519 + tmin = self._tbins.tstart(0)
  520 + tmax = self._tbins.tstop(self._tbins.size()-1)
  521 +
  522 + # Select events
  523 + select = ctools.ctselect(self._obs)
  524 + select['emin'] = self['emin'].real()
  525 + select['emax'] = self['emax'].real()
  526 + select['tmin'] = tmin.convert(self._time_reference())
  527 + select['tmax'] = tmax.convert(self._time_reference())
  528 + select['rad'] = 'UNDEFINED'
  529 + select['ra'] = 'UNDEFINED'
  530 + select['dec'] = 'UNDEFINED'
  531 + select.run()
  532 +
  533 + # Extract observations
  534 + self._obs = select.obs().copy()
  535 +
  536 + # Write observation into logger
  537 + if self._logTerse():
  538 + self._log('\n')
  539 + self._log.header1(gammalib.number('Selected observation',len(self._obs)))
515 540 self._log(str(self._obs))
516 541 self._log('\n')
517 542  
... ...