<p>I want to produce samples from a double Schechter luminosity function.</p>
<p>Ideally, I am looking for a python implementation of this problem but am open to other tools if necessary.</p>
<p>Do you know of any such package that will facilitate this?</p>

CC BY-SA 4.0
2020-09-23T00:14:15.600Z
Do you mean you want to draw a random luminosity from a PDF given by a double Schechter? The rather old package [CosmoloPy](https://github.com/roban/CosmoloPy) can return the cumulative distribution of a [regular Schechter and other LFs](https://roban.github.io/CosmoloPy/docAPI/cosmolopy.luminosityfunction-pysrc.html), so perhaps you can modify and put in the extra term?
pela
CC BY-SA 4.0
2020-09-23T17:18:33.267Z
That is what I meant, yes. Thank you for the link - looks promising!
jm22b
<galaxy><luminosity><resource-request>
Samples from the double Schechter luminosity funcion
jm22b
# Samples from the double Schechter luminosity funcion

(Asked by: username_0 on 2020-09-22)

I want to produce samples from a double Schechter luminosity function.

Ideally, I am looking for a python implementation of this problem but am open to other tools if necessary.

Do you know of any such package that will facilitate this?

**Tags:** galaxy, luminosity, resource-request

> **Comments:**
> **username_1**: Do you mean you want to draw a random luminosity from a PDF given by a double Schechter? The rather old package [CosmoloPy](https://github.com/roban/CosmoloPy) can return the cumulative distribution of a [regular Schechter and other LFs](https://roban.github.io/CosmoloPy/docAPI/cosmolopy.luminosityfunction-pysrc.html), so perhaps you can modify and put in the extra term?  
> **username_0**: That is what I meant, yes. Thank you for the link - looks promising!  

---

## 0 answers

