GBARpy 0.0.3

Libraries and software to analyse data (the MCP's picture and CsI signals) in the framework of the GBAR experiment

How to use the library

This library is made to work with python 3.

Table of Contents

• Installation
• MCPpicture library
  • Basic code for an analysis
  • More examples
  • BeamSpot class
    • Attributes
    • Methods

Installation

MCPpicture library

Basic code for an analysis

The examples corresponds to the python scripts MCP_example_basic.py and MCP_example_small_functions.py.

# Import the library
import GBARpy.MCPPicture as mcp
import matplotlib.pyplot as plt

Let's see how to import a beam spot picture (you can try with the file IMG0008.bmp):

# reshape analyse the beam spot
pic = mcp.BeamSpot("IMG0008.bmp",reshape=[1250,1000,600])
# the reshape array is a made of 3 elements: positionX, positionY, length (in pixel)
# if you don't want to reshape
pic = mcp.BeamSpot("IMG0008.bmp")

With this minimal code, the picture is analysed and the parameters of the fit can be obtained usign print(img). If it then possible to see and save the pictures as in the following example:

#plot the image
fig1 = plt.figure(figsize=(5,5))
plt.imshow(pic.img)
fig1.savefig("fig_example_1.pdf")

Even if it can be written manually, there are line codes to plot the intgrals along the x-axis and the y-axis:

#plot the fit
fig2 = plt.figure(figsize=(5,5))
pic.plot_X_int()
pic.plot_Y_int()
fig2.savefig("fig_example_2.pdf")

or to plot a summary of the fit:

#plot all
pic.plot("fig_example_3.pdf")

More examples

To import the required librairies for the following examples:

# Import the library
import GBARpy.MCPPicture as mcp
import matplotlib.pyplot as plt

For some reasons, you might desire to import the picture without analysing it:

### Import the Picture
img = mcp.import_image("IMG0008.bmp")
fig4 = plt.figure(figsize=(5,5))
plt.imshow(img)
fig4.savefig("fig_example_4.png")

### Import the Picture and reshape
img = mcp.import_image("IMG0008.bmp",reshape=[1250,1000,600])
fig5 = plt.figure(figsize=(5,5))
plt.imshow(img)
fig5.savefig("fig_example_5.png")

Once the pictures imported as a 2D array, it it possible to get the integrals along the x or y axis

### Integrals along the X and Y axis
Px,Ix = mcp.integrate_picture_along_X(img)
Py,Iy = mcp.integrate_picture_along_Y(img)
fig6 = plt.figure(figsize=(10,5))
plt.subplot(121)
plt.imshow(img)
plt.subplot(122)
plt.plot(Px,Ix)
plt.plot(Py,Iy)
fig6.savefig("fig_example_6.png")

and then, using the fit function defined in the library

### Fit of the integrals
poptX,perrX=mcp.fit_gaussian_offset_filtered(Px,Ix)
poptY,perrY=mcp.fit_gaussian_offset_filtered(Py,Iy)
fig7 = plt.figure(figsize=(5,5))
plt.plot(Px,Ix,'.',color='tab:red',ms=1)
plt.plot(Px,mcp.gaussian_offset(Px,*poptX),color='tab:red')
plt.plot(Py,Iy,'.',color='tab:blue',ms=1)
plt.plot(Py,mcp.gaussian_offset(Py,*poptY),color='tab:blue')
fig7.savefig("fig_example_7.png")

BeamSpot class

Attributes

• BeamSpot.fname: string, file name of the picture
• BeamSpot.img: 2D array, picture as an array
• BeamSpot.pix: the pixels along the x axis
• BeamSpot.piy: the pixels along the y axis
• BeamSpot.Ix: array of floats, integral along the x axis
• BeamSpot.Iy: array of floatt, integral along the y axis
• BeamSpot.Ax: float, Amplitude, fit along the x axis
• BeamSpot.Ay: float, Amplitude, fit along the y axis
• BeamSpot.sigx: float, Sigma, fit along the x axis
• BeamSpot.sigy: float, Sigma, fit along the x axis
• BeamSpot.r0x: float, Center, fit along the x axis
• BeamSpot.r0y: float, Center, fit along the x axis
• BeamSpot.offsetx: float, offset, fit along the x axis
• BeamSpot.offsety: float, offset, fit along the x axis
• BeamSpot.poptx: array of floats, the parameters of the fit along the x-axis
• BeamSpot.perrx: array of floats, errors on the parameters of the fit along the x-axis
• BeamSpot.popty: array of floats, the parameters of the fit along the y-axis
• BeamSpot.perry: array of floats, errors on the parameters of the fit along the y-axis
• BeamSpot.reshape: array of int, the parameters to reshape, see help(import_image)

Methods

__init__(self,fname,reshape=[]): Constructor of the class

Parameters

• fname: string, file name of the picture, the accepted file format ["tif","jpg","jpeg","png","asc","bmp"]
• reshape: array of 3 integers (optional), to reshape the pictures (square): x,y,length

Example

import GBARpy.MCPPicture as mcp
bs = mcp.BeamSpot("name.tif")

__repr__(self): To represent the object as a string

Returns

• a string variable

Example

import GBARpy.MCPPicture as mcp
bs = mcp.BeamSpot("name.tif")
repr = bs.__repr__()
#or to print it in the python console
print(bs)

plot_Y_int(self,label=""): To plot the integral of the picture along the "y" axis

Parameters

• label: (optional) a string

Example

import GBARpy.MCPPicture as mcp
bs = mcp.BeamSpot("name.tif")
bs.plot_Y_int("Integral along the y-axis")

plot_X_int(self,label=""): To plot the integral of the picture along the "x" axis

Parameters

• label: (optional) a string

Example

import GBARpy.MCPPicture as mcp
bs = mcp.BeamSpot("name.tif")
bs.plot_X_int("Integral along the x-axis")

plot_X_int_revert(self): To plot the integral of the picture along the "x" axis and reverse the picture

Example

import GBARpy.MCPPicture as mcp
bs = mcp.BeamSpot("name.tif")
bs.plot_X_int("Integral along the x-axis")

plot(self,fname="",figsize=(12,10),fontsize=12,ftsizeticks=12): To plot the picture and the analysis

Parameters

• fname: string (optional), the name of the file to save the plot
• figsize: tuple (size in inch X, Y) (optional), size of the figure
• fontsize: int (optional), size of the font
• ftsizeticks: int (optional), size of the ticks' font

Returns

• fig: a matplotlib.pyplot.figure

Example

import GBARpy.MCPPicture as mcp
bs = mcp.BeamSpot("name.tif")
fig = bs.plot("analysis.pdf")
# or
fig = bs.plot()
fig.savefig("analysis.pdf")