In [1]:
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

from electronfactors import (
    new_reading, calc_and_display, energy_to_reference_depth, 
    initialise)

In [2]:
data = dict()

12 MeV


In [3]:
energy_to_reference_depth(12)


Out[3]:
array([25])

In [4]:
# Standard insert 10app 12MeV 100MU
np.mean([1.054, 1.055])


Out[4]:
1.0545

In [5]:
key = '6.7x12cm_12MeV_v1_col0'

data = initialise(
    key=key, data=data,
    reference=1.0545,
    energy=12
)

data = new_reading(
    key=key, data=data,
    ionisation=np.array([1.065, 1.065]),
    depth=25
)

data[key]['factor'] = calc_and_display(**data[key])


Cutout factor = 1.010 | 1.0%
Inverse factor = 0.990 | -1.0%

In [6]:
key = '6.7x12cm_12MeV_v1_col90'

data = initialise(
    key=key, data=data,
    reference=1.0545,
    energy=12
)

data = new_reading(
    key=key, data=data,
    ionisation=np.array([1.065, 1.065]),
    depth=25
)

data[key]['factor'] = calc_and_display(**data[key])


Cutout factor = 1.010 | 1.0%
Inverse factor = 0.990 | -1.0%

In [7]:
key = '6.7x12cm_12MeV_v2_col0'

data = initialise(
    key=key, data=data,
    reference=1.055,
    energy=12
)

data = new_reading(
    key=key, data=data,
    ionisation=np.array([1.054, 1.055]),
    depth=25
)

data[key]['factor'] = calc_and_display(**data[key])


Cutout factor = 1.000 | -0.0%
Inverse factor = 1.000 | 0.0%

In [8]:
key = '6.7x12cm_12MeV_v2_col0'

data = initialise(
    key=key, data=data,
    reference=1.055,
    energy=12
)

data = new_reading(
    key=key, data=data,
    ionisation=np.array([1.054, 1.055]),
    depth=25
)

data[key]['factor'] = calc_and_display(**data[key])


Cutout factor = 1.000 | -0.0%
Inverse factor = 1.000 | 0.0%

In [9]:
# Standard insert 10app 6MeV 100MU
np.mean([1.055, 1.055])


Out[9]:
1.0549999999999999

In [ ]:


In [ ]:


In [ ]: