# Dose simulation

**URL:** <https://openmc.discourse.group/t/dose-simulation/3297>\
**Category:** User Support\
**Created:** [August 10, 2023, 6:15am UTC](https://openmc.discourse.group/t/dose-simulation/3297 "2023-08-10T06:15:12Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![ZENGJIE](https://avatars.discourse-cdn.com/v4/letter/z/898d66/32.png) [@ZENGJIE](https://openmc.discourse.group/u/ZENGJIE)\
**Post date:** [August 10, 2023, 6:15am UTC](https://openmc.discourse.group/t/dose-simulation/3297/1 "2023-08-10T06:15:12Z")

</div>

I’m new to openmc。  
I am currently doing dose simulation of BNCT, which needs to simulate neutron dose, gamma dose, boron dose, hydrogen dose and so on.

On the forum, you see that you use the.dose\_Coefficients function to simulate the dose. Such as:

```auto
particle_filte = openmc.ParticleFilter('neutron')
filter_cell = openmc.CellFilter((cell1))

tally1 = openmc.Tally()
energy, dose = openmc.data.dose_coefficients('neutron','AP' )
dose_filter = openmc.EnergyFunctionFilter(energy, dose)
tally1.filters = [filter_cell, particle_filte, dose_filter]
tally1.scores = ['flux']

```

AP refers to frontal incidence and has a corresponding kerma factor。  
MOW，I want to replace the kerma factor for AP with the kerma factor for boron dose in order to calculate the boron dose。  
I want to know, what do I do?  
thanks！

---

<div class="post-metadata">

**Author:** ![Shimwell](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/shimwell/32/3876_2.png) [@Shimwell](https://openmc.discourse.group/u/Shimwell)\
**Post date:** [August 11, 2023, 11:56am UTC](https://openmc.discourse.group/t/dose-simulation/3297/2 "2023-08-11T11:56:08Z")

</div>

Hi Zengjie

The dose\_coefficients returns [this](https://github.com/openmc-dev/openmc/blob/develop/openmc/data/effective_dose/neutrons.txt) data which is list of energy and dose values. Perhaps you can replace these lists with your own list of values.