

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]:

    
# Standard insert 14app 12MeV 100MU
np.mean([1.039, 1.039])









    Out[5]:





1.0389999999999999



In [6]:

    
key = 'arbit'

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

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

data = new_reading(
    key=key, data=data,
    ionisation=np.array([1.050, 1.050]),
    depth=24
)

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









    












    



Cutout factor = 1.009 | 0.9%
Inverse factor = 0.991 | -0.9%



In [7]:

    
# Standard insert 14app 12MeV 100MU
np.mean([1.039, 1.039])









    Out[7]:





1.0389999999999999



In [8]:

    
# Standard insert 10app 12MeV 100MU
np.mean([1.053, 1.053])









    Out[8]:





1.0529999999999999



In [ ]: