

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()

9 MeV



In [3]:

    
energy_to_reference_depth(9)









    Out[3]:





array([20])



In [4]:

    
# Standard insert 6app 9MeV 150MU
np.mean([1.0171, 1.0186, 1.0187])









    Out[4]:





1.0181333333333333



In [5]:

    
key = 'arbit'

data = initialise(
    key=key, data=data,
    reference=1.018,
    energy=9
)

data = new_reading(
    key=key, data=data,
    ionisation=np.array([0.89902, 0.89983, 0.89983]),
    depth=13
)

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









    



Cutout factor = 0.868 | -13.2%
Inverse factor = 1.152 | 15.2%



In [6]:

    
# Standard insert 6app 9MeV 150MU
np.mean([1.0173, 1.0181, 1.0177])









    Out[6]:





1.0177



In [ ]: