Should the sources have a method get_neutron_rate ?

[RemDelaporteMathurin]

The sources should maybe have a method to compute the neutron rate in neutron/s .

usage could be

from openmc_plasma_source import FusionPointSource

my_source = FusionPointSource(.....)

neutron_rate = my_source.get_neutron_rate(fusion_power=1e9)  # for 1 GW

[shimwell]

Sounds useful for steady state machines, not sure how this would work for pulsed machines

[RemDelaporteMathurin]

Ah yes I guess the neutron rate would be in n/s or n/shot.

Anyway the calculation is the same isn't it? for SS machines, the input is the fusion power in J/s (W) for pulsed machines it's in J/pulse

[shimwell]

I'm just not sure if we need to separate keywords or not, happy to go along with your vision here

[RemDelaporteMathurin]

I think it will just have to be described in the docstrings. It's still a fusion power in both cases, the unit of time changes that's all.