oyvsyo / moed_task1

Phisical spectrum modulation program on python and ROOT

Task1

Task1 is the open source project for spectrum modulation designed by Dvoiak Stepan,
Kiev National University student for MOED course.
This program can modulate peaks, compton and background by parameters.

The code and this instruction also available at project site

Getting Started

Clone this repo to your local machine:

$ git clone https://github.com/oyvsyo/moed_task1.git

To start nodulation jus run this command in terminal (bash shell or mac shell) in project directory:

$ python2 task1.py

You will see two windows (left - for config window, rigth - for drawing spectrums).
In config window you can set parameters of spectrums, load configurations, save configuration, save spectrums in txt file and eps pictures.

Default parameters

To draw simple default specters just press red button Draw.
You got spectrums on canvas (drawing window)
If you want to clear all entries - press pink button Clear All.
If you want to set defaul parameters again - press Default button.

Your own configuration

The best way to set you own spectrum parameters - it's change defaul!
All rules and sintaxys diclarated in next paragraf.

Also, you can write config file manualy (syntaxis - like json
or standart python dict)
To load your file - simply open menu File->Open config

The rules of config file

Concept of config file is python dict - structure of "key": "value"
pairs. There are 6 keys of parameters:

• interval - the python list of two values: minimum and maximum
  of Energy range of investigated spectrum [Emin, Emax] (unsigned float)
• num_bin - number channels (unsigned interer)
• peaks - it's a python list of three dimentional lists with peaks
  parameters values [Intensity(1/s), Energy(keV), FWMH(keV)] (unsigned float)
• opt1 - optional chose for FWHM simulation (integer 1 or 0), if 0 - FWHM gets
  from last peak parameters, if 1 - FWHM is calculated by formula: a0+a1*sqrt(E)
• FWHM_pars - the parameters of FWHM formula [a0, a1] (float numbers)
• background - parameters of background simulation formula:
  a*exp(b* E) + c+d* E (float numbers)

Notice! Its not recomendet to set num_bin parameter over 1000.
Calculating can take too much time (4 seconds on my computer).
Operations numbers increase like (num_bin)^2 times.

Also you can comment some text with "<" and ">" symbols like:
<comment here>.

Notice! You can't use "<" and ">" symbols inside comment.

An example of default confuguration file:

<spectrum simulation>
<config file example>
<after square brackets "[..]" always put ",">
<all parameters must be set>
{
    <Energy(keV) range - must be positive:>
    'interval': [50, 3000],
    <number of bins (put "," after bin number value):>
    <must be positive int number>
    'num_bin': 300,

    <peaks parameters sets in square brackets>
    <in defined order: [Intensity(1/s), Energy(keV), FWMH(keV)]>
    <peaks parameters must be positive>
    'peaks':[
        [15000, 2000.89, 150.8],
        [14000, 1400.5, 135.0],
        [13000, 1000, 131.0],
        [17000, 1200, 75.0],
        [12000, 2999, 50.0],
    ],

    <optional chose for FWMH simulation>
    < 0 - FWMH is set from peak parameters>
    < 1 - calculating by formula: FWHM = a0+a1*sqrt(E)>
    <after opt1 value put ",">
    "opt1": 1,

    <set the FWMH params for formula FWHM = a0+a1*sqrt(E)>
    <if opt1=1 [a0, a1]:>
    "FWHM_pars": [0.5, 1],

    <background parameters [a, b, c, d] : a*exp(b*E) + c+d*E>
    <parameter 'a' must be positive>
    'background': [40, -0.01, 50, -0.005],
}

Notice! The config file must have .txt exestention

Using GUI

The main logic based on config file syntax.
The simplest way - just edit default parameters!

To set parameters with GUI you need to know some rules:

1. All entries are required.
2. Enter the vaule of "key":"value" pairs from config file.
3. If key value is list - write parameters in list format (like [a, b, c]).
4. For FWMH simulation by formula: FWHM = a0+a1*sqrt(E) just check
   FWHM optin calculate checkbuttun
5. Other parameters sets like vaule vor keys in config file dict.

Notice! GUI does not support comments.

Example for peak entry:

[[14000, 2700.89, 110.8]
 [16000, 700.5, 135.0]
 [17000, 1600, 75.0]]

Ntice: for peaks its not nessecery to write "," after "]"

To save your comfiguration press File->Save config as..
and enter name.

Saving results of simulation

To save a pictures and .txt files of simulated spectrum press
File->Save output to.. and write a path to directory that
will be created for yuor output files (dir name - name_out).
You can find 4 eps pictures and 4 txt spectrum files. (txt files - its two columns of float numbers - first - energy, second - counts).

Requirments

This software is require python 2.7* versions, some standart for ubuntu 14.04 python librarys:

• math
• re
• os
• ast
• sys
• Tkinter
• tkFileDialog
• tkMessageBox

There are one non standart library:

• numpy 1.8.2

And tha main core of simulation - ROOT.

Notice! For using this program you must have ROOT classes avaliable in python scope. Read more about PyRoot

Notice! ROOT 5.* versions are required (probably must works on 6.* versions, but not tested)

Tested on ubuntu 14.04 whith python 2.7.6 and ROOT 5.34/36, gcc 4.8.4

Authors

Dvoiak Stepan - oyvsyo

License

This project is licensed under the GNU GENERAL PUBLIC LICENSE - see the LICENSE file for details