| 1 | import numpy as np |
|---|
| 2 | #import matplotlib.pyplot as plt |
|---|
| 3 | from tube import calibrate |
|---|
| 4 | from tube_spec import TubeSpec |
|---|
| 5 | |
|---|
| 6 | #high intensity Version 30x30mm beam |
|---|
| 7 | w258=Load("258") |
|---|
| 8 | for i in range(259,266): |
|---|
| 9 | print i |
|---|
| 10 | wtemp=Load(str(i)) |
|---|
| 11 | w258=w258+wtemp |
|---|
| 12 | |
|---|
| 13 | w258int=Integration(w258,20000.0,100000.0) |
|---|
| 14 | ws2=Integration(w258,45000.0,100000.0) |
|---|
| 15 | |
|---|
| 16 | known_pos1 = np.array([-0.297,-0.2,-0.1,0.0,0.1,0.2]) |
|---|
| 17 | known_pos1 = known_pos1 +(0.6165-0.571) |
|---|
| 18 | print known_pos1 |
|---|
| 19 | |
|---|
| 20 | peaks_form1 = 6*[1] # all the peaks are gaussian peaks |
|---|
| 21 | tubes1=['LARMOR/LARMORSANSDetector/tube%d'%(i) for i in range(1,81)] |
|---|
| 22 | print tubes1 |
|---|
| 23 | |
|---|
| 24 | calibrationTable = calibrate(ws2,tubes1,known_pos1, peaks_form1,margin=20) |
|---|
| 25 | |
|---|
| 26 | |
|---|
