diff options
author | Aaron LI <aaronly.me@outlook.com> | 2017-03-04 17:30:04 +0800 |
---|---|---|
committer | Aaron LI <aaronly.me@outlook.com> | 2017-03-04 17:30:04 +0800 |
commit | 6ce0d5aa7c3d9280bde5e5b35606b3b9a56cafc3 (patch) | |
tree | 29043d40d84317249eb76757317050c63b0c95dc /acispy | |
parent | 0dd39d4d94d75341bfeb5184a4032f037dc7175f (diff) | |
download | chandra-acis-analysis-6ce0d5aa7c3d9280bde5e5b35606b3b9a56cafc3.tar.bz2 |
acispy/spectrum.py: Add "verbose" parameter
Diffstat (limited to 'acispy')
-rw-r--r-- | acispy/spectrum.py | 9 |
1 files changed, 6 insertions, 3 deletions
diff --git a/acispy/spectrum.py b/acispy/spectrum.py index 6b4481c..3fa0216 100644 --- a/acispy/spectrum.py +++ b/acispy/spectrum.py @@ -27,7 +27,7 @@ class Spectrum: self.EXPOSURE = self.header.get("EXPOSURE") self.BACKSCAL = self.header.get("BACKSCAL") - def calc_flux(self, elow, ehigh): + def calc_flux(self, elow, ehigh, verbose=False): """ Calculate the flux: flux = counts / exposure / area @@ -40,11 +40,14 @@ class Spectrum: chlow = ACIS.energy2channel(elow) chhigh = ACIS.energy2channel(ehigh) counts = self.counts[(chlow-1):chhigh].sum() + if verbose: + print("counts / exposure / backscale :: %d / %.1f / %.5g" % + (counts, self.EXPOSURE, self.BACKSCAL)) flux = counts / self.EXPOSURE / self.BACKSCAL return flux - def calc_pb_flux(self, elow=9500, ehigh=12000): + def calc_pb_flux(self, elow=9500, ehigh=12000, verbose=False): """ Calculate the particle background (default: 9.5-12 keV) flux. """ - return self.calc_flux(elow=elow, ehigh=ehigh) + return self.calc_flux(elow=elow, ehigh=ehigh, verbose=verbose) |