Hide keyboard shortcuts

Hot-keys on this page

r m x p   toggle line displays

j k   next/prev highlighted chunk

0   (zero) top of page

1   (one) first highlighted chunk

1import astropy.units as u 

2from lsst.pipe.base import Struct, Task 

3from lsst.pex.config import Config 

4from lsst.verify import Measurement 

5 

6__all__ = ("StarFracTask", ) 

7 

8 

9class StarFracTask(Task): 

10 ConfigClass = Config 

11 _DefaultName = "starFracTask" 

12 

13 def run(self, catalog, metric_name, vIds): 

14 self.log.info(f"Measuring {metric_name}") 

15 if not catalog.isContiguous(): 

16 catalog = catalog.copy(deep=True) 

17 extended = catalog.get('base_ClassificationExtendedness_value') 

18 good_extended = extended[~catalog.get('base_ClassificationExtendedness_flag')] 

19 n_gals = sum(good_extended) 

20 frac = 100*(len(good_extended) - n_gals)/len(good_extended) 

21 meas = Measurement("starFrac", frac * u.percent) 

22 return Struct(measurement=meas)