def q_prior(q, m=<NUM_LIT:1>, gamma=<NUM_LIT>, qmin=<NUM_LIT:0.1>):
if q < qmin or q > <NUM_LIT:1>:<EOL><INDENT>return <NUM_LIT:0><EOL><DEDENT>C = <NUM_LIT:1>/(<NUM_LIT:1>/(gamma+<NUM_LIT:1>)*(<NUM_LIT:1> - qmin**(gamma+<NUM_LIT:1>)))<EOL>return C*q**gamma<EOL>
Default prior on mass ratio q ~ q^gamma