search

NuSTAR / nustar-idl

Repository for NuSTAR IDL-based analysis and support scripts
3 stars 4 forks source link

Flux resolved spectra. To create spectra for different flux levels (low, high, very high) #21

Open samuzalbarua opened 6 years ago

samuzalbarua commented 6 years ago

Generally the code nustar_filter_lightcurve.pro produces GTI file for a particular rate cutoff (say 0 to 0.4 or 0 to 0.5), which is clearly justified from line no. 17 and 23 in this script. If I want to create GTI for any range of flux (say 0.4 - 0.5), then how do I change the script. Please kindly help me.

bwgref commented 6 years ago

I don't think I'd checked in here since the Spring, my apologies.

That's not really what the tool is meant for. This is really in place to assist people in filtering out periods with high count rate due to solar flares. But in principle we could add a minimum flux filter to the script.

However, this script has several safety checks to make sure that you're not throwing away a lot of data. I would recommend starting a new script, removing the safety checks, and then add a new "minrate" keyword (which defaults to 0 counts / sec). Line 18 then becomes:

goodones = where(lc.rate LT rate and lc.rate GT minrate, ngood)

...with the rest of the script left in place.

This is a relatively minor change to the script, but if it's useful I can add a new version as described above. But I would guess that you've found any number of other work arounds by now.