

In [1]:

    
from chxanalys.chx_packages import *
%matplotlib notebook
plt.rcParams.update({'figure.max_open_warning': 0})
#%reset -f  #for clean up things in the memory









    



/opt/conda_envs/analysis/lib/python3.5/site-packages/filestore/retrieve.py:15: UserWarning: Do not import filestore.retrieve, import filestore.api instead
  warnings.warn("Do not import filestore.retrieve, "



In [2]:

    
from chxanalys.Two_Time_Correlation_Function import *

Attenuator =1, do a discriminator sweep to find the voltage peak at different applied HV voltage



In [39]:

    
np.linspace( 700, 2000, 14)









    Out[39]:





array([  700.,   800.,   900.,  1000.,  1100.,  1200.,  1300.,  1400.,
        1500.,  1600.,  1700.,  1800.,  1900.,  2000.])



In [15]:

    
T=1
HV1 = np.linspace( 700, 2000, 14)
#         array([  700.,   800., 900.,1000.,1100.,1200.,1300.,1400.,1500.,  1600.,  1700.,  1800.,  1900.,  2000.])
vol1 = np.array(  [0.09, 0.14, 0.18, 0.24, 0.29, 0.32, 0.34,  0.38, 0.42,  0.47,   0.51,    0.56,   0.61,   0.66]  )
T= 0.5
HVp5 = np.linspace( 700, 2000, 14)
volp5 = np.array(  [0.05, 0.14, 0.2, 0.25, 0.29, 0.35, 0.4, 0.46, 0.51, 0.54, 0.58, 0.62, 0.65, 0.66]  )



In [37]:

    
fig, ax = plt.subplots()

plot1D( y=vol1, x=HV1, ax=ax, m='o', c='k',xlabel='HV, (V)', ylabel='peak_voltage, (V)', legend='peak@T=1')
plot1D( y=volp5, x=HVp5, ax=ax, m='s', c='r', xlabel='HV, (V)', ylabel='peak_voltage, (V)', 
       legend='peak@T=0.5',  title='NaI Scitillation Counter')

At different applied HV voltage, check count near the finding peak voltage against attenuator

the voltage range centered at peak position with width as 0.02
integration time 1000 us
dead time 20 ns



In [32]:

    
vol = 700
att_n = np.array(       [1,     .9,     0.7,   0.5,    0.4,    0.3,       0.2,    0.1,   ]                                           )
att_700 =   np.array(  [ 1,   0.812,   0.66,  0.54,   0.44,   0.288,      0.19,   0.1      ])
count_700 = np.array([ 430,  440,  305,   60,    0,      0,   0,     0       ]   )

count_800 = np.array([ 490,  470,  500, 540,  540, 80, 0, 0    ]   )
count_800_2 = np.array([ 400,  400,  400, 460,  440, 80, 0, 0    ]   )

att_1500 =   np.array(  [ 1,   0.812,   0.66,  0.54,   0.44,   0.288,      0.19,   
                         0.1, 0.0828,  0.055, 0.036, 0.019, 0.01, 0.0084, 0.0056, 0.0037,  0.002, 0.001    ])
count_1500 = np.array([ 285,  2, 0, 0,  0, 0, 10, 120, 340, 560, 700, 720, 190, 75,  20, 4, 1, 0     ]   )



In [36]:

    
fig, ax = plt.subplots()
plot1D( x=att_700, y=count_700, ax=ax, m='o', c='k',xlabel='attenuation', ylabel='count', legend='@HV=-700')
plot1D( x=att_700, y=count_800, ax=ax, m='s', c='r',xlabel='att', ylabel='count', legend='@HV=-800')
plot1D( x=att_700, y=count_800_2, ax=ax, m='s', c='b',xlabel='att', ylabel='count', legend='@HV=-800_2')
#plot1D( volp5, HVp5, ax=ax, m='s', c='r', xlabel='HV, (V)', ylabel='peak_voltage, (V)', legend='peak@T=0.5')
plot1D( x=att_1500, y=count_1500, ax=ax, m='D', c='g',xlabel='attenuation', ylabel='count', 
       legend='@HV=-1500', logx=True, title='NaI Scitillation Counter')



In [ ]:



In [ ]:



In [ ]: