HOME ESR EXPERIMENTS LABS etc
E121 E127 E132 E125 E143 laser_cooling_2021 E142 Ê128 E146 E0052 E0018
  proton-capture on 118Te, Page 1 of 29  Not logged in ELOG logo
ID Date Author Category Subjectdown Year
  38   Tue Mar 10 09:41:26 2020 JanDetectorsx-ray background from BaF 
The BaF2 detectors planned to be installed at the target induce an additional background in the x-ray detectors sitting below inside the target chamber. 
We made test measurements with an x-ray detector for estimating the background in the region of interest for our K-REC measurements (35 to 50 keV).

For the interval [35-50 keV]:
Normal background level: 0.065 (s*keV)^-1
BaF2 in 5cm distance:    0.078 (s*keV)^-1
This is an increase by about 22%

Assuming an increase by 100% induced by 6 BaF2 in the vincinity of the Ge-detectors this seems acceptable.

However additional BaF-induced peaks are also visible at ~30 keV and ~60 keV. While these are not in the region of interest for the K-REC, we need to be aware of them!

spectra are available on the Frankfurt exp-astro Server at /home/glorius/e127/
  241   Mon Mar 23 12:30:39 2020 laszloAnalysiswith and without scraper 
Attachment 1: 124Xe_with_scraper_3D.png
124Xe_with_scraper_3D.png
Attachment 2: 124Xe_with_scraper_2D.png
124Xe_with_scraper_2D.png
Attachment 3: 124Xe_without_scraper_2D.png
124Xe_without_scraper_2D.png
Attachment 4: 124Xe_without_scraper_3D.png
124Xe_without_scraper_3D.png
  62   Wed Mar 18 09:53:28 2020 Laszlo, JanGeneralvacuum after opening gate valve 
Gate valve to DSSSD setup was opened on Tue 17.03. ~9:00

The vacuum on Wed 18.03. ~9:00 was
1.6e-10 mbar at DSSSD setup
1.0e-11 mbar at ESR dipole
Attachment 1: vacuum_ESR1.jpg
vacuum_ESR1.jpg
Attachment 2: vacuum_setup1.jpg
vacuum_setup1.jpg
  277   Thu Apr 23 19:09:58 2020 LaszloCalibrationtheoretical K-REC cross sections  
Find attached Andrey Surzhykov's calculations for the theta angle in respect to the beam direction (in lab. frame) vs. cross section for 124Xe54+ and 118Te52+.
The calculations made for collision with two H atoms with the accuracy of 1%. There are no molecular corrections done, but these corrections are within 1%.

The photon-emission is symmetrical in the azimuthal phi angle, but asymetric in theta. The K-REC cross section is given for each integer theta angle. The problem is that our 90° and 145° xray detectors cover more integer theta-angles --> The disk shaped entrance window of the Xray-detectors is sliced for each covered theta angles, and the CS values are averaged together with the weights of the area of the corresponding slice. The 35° det had a non-disk shaped slit collimator! Was it aligned vertical or horizontal or random? I assumed that it had a vertical position --> only cs at theta=35° needs to taken into account

90°145°35°
weighted cs_Xe [barn/sr]128,44441,96745,550
weighted cs_Te [barn/sr]118,61638,91841,740


For the steradian values I calculated the area of the entrance window of the Xray detectors and I divided it with the area of the sphere with radius = distance between source and det. If r=radius of the entrance window, D=distance between source and det.: covered steradian = r^2*pi/(4*pi*D^2) [%], covered steradian = 4Pi*[r^2*pi/(4*pi*D^2)] [4pi]. This is valid for the 90° and 145° Xray detectors. The area of 35° det. was calculated individually.

90°145°35°
sr [%]0,002370,000760,00002
sr [4pi]0,02980,00960,00030


I converted barn/sr to barn in the following way: steradian[4pi]*weigthed_cs [barn/sr] = weigthed_cs [barn]. I am not 100% sure if I need here the steradian in [4pi] or in [%].

90°145°35°
weighted cs_Xe [barn]3,82470,40120,0137
weighted cs_Te [barn]3,53210,37200,0125
Attachment 1: CS_theory_all.png
CS_theory_all.png
Attachment 2: CS_theory_all.ps
CS_theory_all.ps
Attachment 3: CS_theory_90Xe.png
CS_theory_90Xe.png
Attachment 4: CS_theory_90Xe.ps
CS_theory_90Xe.ps
Attachment 5: CS_theory_90Te.png
CS_theory_90Te.png
Attachment 6: CS_theory_90Te.ps
CS_theory_90Te.ps
Attachment 7: DCS_Z_54_Tp_10MeV.dat
0. 0.8764393262751978
1. 0.9184409606576515
2. 1.0443879628122903
3. 1.2541067241116186
4. 1.5473082111252325
5. 1.9235884362763038
6. 2.3824291152290664
7. 2.9231985094824493
8. 3.5451524522204467
9. 4.2474355550499725
10. 5.029082592847121
11. 5.889020063533556
12. 6.826067919219877
13. 7.83894146477025
14. 8.926253419493086
15. 10.086516137307964
16. 11.318143980418107
17. 12.619455841203488
18. 13.988677806759991
19. 15.423945960238417
20. 16.92330931288666
21. 18.484732860465083
22. 20.106100757502197
23. 21.785219602666896
24. 23.51982182837228
25. 25.307569187584907
26. 27.146056330695448
27. 29.03281446521163
28. 30.965315090961862
29.000000000000004 32.94097380344995
30. 34.95715415797087
31. 37.011171587095994
32. 39.100297364149334
33. 41.2217626053315
34. 43.37276230320916
35. 45.55045938435645
36. 47.75198878403283
37. 49.97446153088989
38. 52.21496883482496
39.00000000000001 54.470586171241614
40. 56.7383773551312
41. 59.01539859855422
42. 61.29870254528142
43.00000000000001 63.58534227653866
44. 65.87237528199884
45. 68.15686739036538
46. 70.43589665410244
47.00000000000001 72.70655718308338
48. 74.96596292214365
49. 77.2112513677509
50. 79.43958721922418
51. 81.64816596016036
52. 83.83421736595126
53. 85.99500893349557
54. 88.1278492294348
55.00000000000001 90.23009115346127
56. 92.29913511346422
57. 94.33243210949429
58.00000000000001 96.32748672373847
59. 98.28186001390355
60. 100.19317230761206
61. 102.05910589560905
62. 103.87740762177793
63.00000000000001 105.64589136814853
64. 107.36244043326852
65.00000000000001 109.0250098024874
66. 110.63162830887903
67.00000000000001 112.18040068369685
68. 113.6695094954226
69.00000000000001 115.09721697662884
70. 116.46186673802805
71. 117.76188536923578
72. 118.99578392591334
73.00000000000001 120.16215930310027
74. 121.25969549467807
75.00000000000001 122.28716473903538
76. 123.24342855113079
77. 124.12743864126509
78.00000000000001 124.93823772099097
79. 125.67496019669385
80. 126.33683275148397
81. 126.9231748161348
82. 127.43339892989992
83.00000000000001 127.86701099212635
84. 128.223610405669
85. 128.50289011318924
86.00000000000001 128.70463652749714
87. 128.82872935716867
88. 128.8751413287365
89.00000000000001 128.84393780681756
90. 128.7352763136031
91. 128.54940594919196
92. 128.2866667143067
93.00000000000001 127.94748873698224
94.00000000000001 127.53239140486892
95. 127.04198240483984
96. 126.47695667163981
97. 125.83809524735832
98. 125.12626405355418
99.00000000000001 124.3424125779007
100. 123.4875724772668
101. 122.56285609918638
102. 121.56945492371455
103. 120.50863792770586
104. 119.38174987359533
105.00000000000001 118.19020952479832
106. 116.93550778989012
107. 115.61920579776323
108. 114.24293290600282
109. 112.80838464475747
110.00000000000001 111.31732059842409
111. 109.77156222750445
112. 108.17299063302771
113. 106.52354426597255
114. 104.82521658415874
114.99999999999999 103.08005365911461
116.00000000000001 101.29015173546108
117. 99.45765474538705
118. 97.58475178082271
119. 95.67367452595023
120. 93.72669465271724
121.00000000000001 91.746121182048
122. 89.73429781347248
123.00000000000001 87.69360022591539
124. 85.62643335241032
125. 83.53522863152084
126.00000000000001 81.42244123826835
127.00000000000003 79.29054729738027
128. 77.14204108168238
129. 74.97943219846945
130.00000000000003 72.80524276669398
131. 70.62200458781662
132. 68.43225631316218
133. 66.23854061062511
134.00000000000003 64.04340133356263
135. 61.84938069470726
136. 59.65901644792196
137. 57.474839080607154
138.00000000000003 55.29936901955791
139. 53.13511385304818
140. 50.984565571903104
141. 48.85019783229762
142. 46.73446324299586
143. 44.63979067971726
144. 42.56858262929029
145.00000000000003 40.523212566222384
146.00000000000003 38.50602236427895
147. 36.51931974563394
148. 34.56537577011455
149. 32.64642236702124
150.00000000000003 30.76464991196859
151. 28.92220485114267
152. 27.121187375328226
153. 25.363649146012833
154. 23.65159107582154
155. 21.986961165486573
156.00000000000003 20.371652399503976
157. 18.80750070256989
158. 17.296282958834453
159. 15.839715095951185
160. 14.439450235837072
161.00000000000003 13.097076913995858
162. 11.81411736919071
163. 10.592025905185897
164. 9.432187326207346
165. 8.335915447699179
166.00000000000003 7.30445168388348
166.99999999999997 6.338963713551622
168. 5.440544225441671
169. 4.610209744476686
170. 3.848899540058887
171. 3.157474617532845
172.00000000000003 2.536716793847944
173. 1.9873278583654939
174. 1.5099288196705962
175.00000000000003 1.1050592391610732
176. 0.7731766520971258
176.99999999999997 0.514656076708857
178.00000000000003 0.32978961186464684
179. 0.21878612371563544
180. 0.1817710216397566
Attachment 8: DCS_Z_52_Tp_10MeV.dat
0. 0.6838809884595187
1. 0.7224229853331035
2. 0.8379964473235315
3. 1.030443872137824
4. 1.2995030344282912
5. 1.64480740285091
6. 2.065886722645698
7. 2.562167762440693
8. 3.1329752236222754
9. 3.7775328102535015
10. 4.494964457174398
11. 5.284295713575055
12. 6.144455279002064
13. 7.074276688431106
14. 8.072500142735743
15. 9.13777448057661
16. 10.268659287456684
17. 11.463627137414845
18. 12.721065962577109
19. 14.039281545545448
20. 15.416500129385255
21. 16.850871139764088
22. 18.340470013614155
23. 19.88330112852013
24. 21.477300826887046
25. 23.12034052881444
26. 24.81022992749396
27. 26.54472026085621
28. 28.321507653123813
29.000000000000004 30.138236519876482
30. 31.99250303019941
31. 33.88185861947699
32. 35.803813546395155
33. 37.75584048773818
34. 39.73537816460896
35. 41.73983499375132
36. 43.76659275772978
37. 45.81301028780408
38. 47.87642715343563
39.00000000000001 49.954167352475665
40. 52.043542996208096
41. 54.14185798355203
42. 56.24641165887502
43.00000000000001 58.35450244801732
44. 60.46343146728837
45. 62.5705061003606
46. 64.67304353815759
47.00000000000001 66.76837427700832
48. 68.85384557051732
49. 70.92682483078326
50. 72.98470297478
51. 75.02489771189985
52. 77.0448567688431
53. 79.0420610482222
54. 81.01402771743378
55.00000000000001 82.95831322453394
56. 84.87251623803525
57. 86.75428050772088
58.00000000000001 88.6012976437524
59. 90.4113098115192
60. 92.18211233985382
61. 93.91155624040339
62. 95.59755063611817
63.00000000000001 97.23806509697859
64. 98.8311318812456
65.00000000000001 100.37484808067475
66. 101.86737766828935
67.00000000000001 103.30695344745868
68. 104.69187890117422
69.00000000000001 106.02052994056153
70. 107.29135655180244
71. 108.5028843407837
72. 109.65371597491449
73.00000000000001 110.74253252168937
74. 111.76809468369437
75.00000000000001 112.72924392987584
76. 113.62490352300958
77. 114.454079443417
78.00000000000001 115.21586120908894
79. 115.90942259247545
80. 116.53402223430722
81. 117.08900415490466
82. 117.57379816352693
83.00000000000001 117.98792016639949
84. 118.33097237414394
85. 118.60264340941445
86.00000000000001 118.80270831562267
87. 118.93102846770623
88. 118.98755138596732
89.00000000000001 118.97231045407382
90. 118.88542454238262
91. 118.72709753780394
92. 118.49761778148746
93.00000000000001 118.1973574156676
94.00000000000001 117.8267716410606
95. 117.38639788626082
96. 116.87685489063509
97. 116.29884170226605
98. 115.65313659254451
99.00000000000001 114.9405958890601
100. 114.16215272848856
101. 113.3188157312197
102. 112.41166759952044
103. 111.44186364106933
104. 110.41063021974982
105.00000000000001 109.31926313563206
106. 108.16912593611966
107. 106.96164816028295
108. 105.69832351844508
109. 104.3807080091302
110.00000000000001 103.01041797552814
111. 101.58912810367194
112. 100.11856936456695
113. 98.60052690255144
114. 97.03683787220862
114.99999999999999 95.42938922619064
116.00000000000001 93.78011545634986
117. 92.0909962906127
118. 90.36405434806198
119. 88.60135275473215
120. 86.80499272264713
121.00000000000001 84.97711109466407
122. 83.1198778577115
123.00000000000001 81.23549362703562
124. 79.32618710409213
125. 77.39421251073887
126.00000000000001 75.44184700240568
127.00000000000003 73.47138806293053
128. 71.48515088376489
129. 69.4854657302612
130.00000000000003 67.47467529776368
131. 65.45513206022716
132. 63.429195614092336
133. 61.3992300201446
134.00000000000003 59.36760114608176
135. 57.33667401250889
136. 55.308810145072414
137. 53.28636493543301
138.00000000000003 51.27168501376619
139. 49.26710563546159
140. 47.27494808467541
141. 45.29751709737059
142. 43.337098306456504
143. 41.39595571161337
144. 39.476329176362796
145.00000000000003 37.580431954915426
146.00000000000003 35.71044825129239
147. 33.86853081318914
148. 32.056798563011185
149. 30.277334268472416
150.00000000000003 28.532182255112247
151. 26.82334616304045
152. 25.152786750176944
153. 23.522419744211014
154. 21.934113745450986
155. 20.389688182690108
156.00000000000003 18.890911324162325
157. 17.43949834560509
158. 16.037109457394024
159. 14.685348092655238
160. 13.385759158201713
161.00000000000003 12.139827350079225
162. 10.948975535442912
163. 9.81456320242211
164. 8.737884979563095
165. 7.720169226369559
166.00000000000003 6.762576696392841
166.99999999999997 5.866199274248659
168. 5.032058787863961
169. 4.261105897183412
170. 3.554219060485478
171. 2.9122035793804977
172.00000000000003 2.335790723483083
173. 1.8256369356687627
174. 1.382323118743823
175.00000000000003 1.0063540042709131
176. 0.6981576042097696
176.99999999999997 0.45808474594681375
178.00000000000003 0.28640869119862783
179. 0.18332483918885642
180. 0.14895051440953486
Attachment 9: source_vs_XrayDet_sketch.pdf
source_vs_XrayDet_sketch.pdf
  48   Sun Mar 15 17:13:21 2020 Laszlo, JanCalibrationthe x-ray gate issue 
We found a strange behaviour of our ADC gates for the x-ray detectors.
The setting as shown in the foto below corresponds to a delay of 5 us and a gate length of 15 us.
This is the only setting that produces nice spectra, even though the signals and gate are not aligned.
The setting is extremely sensitive to the delay 4 us still works but 3 or 7 us give already reduced energy values and beyond these values the spectra disapears.

For now it was decided to operate with this not-understood settings because the result is pretty good...
However, further investigation to understand whats happening would be desirable.

Additional note:
This effect appeared after using the new ortec SpecAmps. Before this we used old Tennelec Amps which delivered horrible resolution, but the gates and signals could be matched as expected at this time.
Attachment 1: IMG_20200315_172129.jpg
IMG_20200315_172129.jpg
  291   Thu Apr 22 10:56:20 2021 JanDetectorstarget-chamber det. distances2021
This is the lookup table for the distance of each x-ray detector to the target.
See attached foto for explanation.

It is the same values as in 2020: https://elog.gsi.de/esr/E127/36

Detector A (dA) [mm] B (dB) [mm] C (dC) [mm] calculate result [mm]
35° 136.0 (1.0) 20.0 (0.2) 450.0 (0.5) C + B - A 334.0 (?)
90° 76.0 (0.5) 20.0 (0.2) 447.0/2 = 223.5 (1.0) C + B - A 167.5 (?)
145° 165.0 (1.0) 20.0 (0.2) 450.0 (0.5) C + B - A 305.0 (?)

Errors are estimated after measurement, final errorbars have to be double checked!
Attachment 1: E127_xray_distances.jpg
E127_xray_distances.jpg
  130   Sat Mar 21 14:52:05 2020 LaszloGeneraltarget optimization 
previous position of the beam at target was -3mm (for the night position)
max. intensity position of the beam at target 1s -1mm. --> now it is set to this position

now Ulli also increasing the density of the jet. Right now Ulli is warming the noozle. --> takes 1hour roughly.
Earlier: 1.4E13
Now: 
Attachment 1: IMG_20200321_145105.jpg
IMG_20200321_145105.jpg
Attachment 2: IMG_20200321_145059.jpg
IMG_20200321_145059.jpg
  68   Fri Mar 20 01:21:00 2020 LaszloGeneraltarget is switched on event based 
  147   Sat Mar 21 21:36:06 2020 JanRunstarget incident 
in run 66 or 67 we had a single ESR cycle for which the target didn't switch off for new injection and deceleration etc.
Very likely the DAQ also didn't get an OFF signal and collected OFF data to the TRIGGER=1 (ON data) branch.
As a consequence we collected about ~100 counts on the Si in one shot, which was about ~20 counts before this shot.

attached:
1.picture of the target density and beam current at this time.
Attachment 1: run67_target_incident.png
run67_target_incident.png
  201   Sun Mar 22 18:35:43 2020 Laszlo target dropped aagain 
there was a drop again from 8.6e13 to 7e13 in the target density
Attachment 1: target_density_after_118Te_afternoon_runs.png
target_density_after_118Te_afternoon_runs.png
  361   Fri May 21 10:00:51 2021 JanAcceleratortarget bump2021
The target overlap was optimized on xray and si rates.
The setting now is -10mm in SC9
  198   Sun Mar 22 17:45:07 2020 LaszloGeneralswitching to primary beam 
We have finished the measurement with the fragments.

Now Yury will setup up a primary beam through the FRS, but without target in FRS.
We are hoping for increased intensity and stable gas target.
The plan is to make a measurement with and without the scraper in order to see the effect of scraping the Rutherford scattering.
  64   Thu Mar 19 15:09:29 2020 JanGeneralstarting off with DSSSD 
here is my short guide how to start with the measurement as soon as we have a beam available.

0. ESR timing in DAQ (maybe Uwe can also assist here)
    - when the cycle is running already, even for tuning, you can try to make the cabling for ESR timing (jet ON/OFF)
    - replace the gate generator signals by the ESR timing signals in the Messhütte and check UDP reader (E127_rates)
    - the settings for ESR timing (event & machine) depend on the cycle, talk to Sergey or so, its the same settings as for the target in HKR
    - with the current cycle this would be: ON-{machine 11, event 55}, OFF-{machine 13, event 55}
    - the event for end of this process is alway 55
    - if it works all tpat rates should be switching from upper to lower panel with the ESR cycle in the rate monitor (as they do now with the gate generator)

Once the tuning is finished:

1. optimize target ON/measurement time
    - switch target to event and use the same machines/events as above
    - check number of ions in ESR at ON and OFF and lifetime of beam
    - we should measure for a drop of a factor of 3-4 roughly i would say
    - e.g. start with 2e5 and stop with 5e4
   - if we start with more, we can measure longer

1. without target,  put in Si detector (GE01DD4AS in device controle) to max-in position
    - see if it disturbs the ESR cycle (discuss optimization with ESR guys if necessary)
    - we should also check vacuum in ring (i think this can be done remotely, ask Markus, Sergey or so)

2. make Si settings with beam (call me i can do things remotely)
    - voltage to 60 volts, use E127_epics (i would expect the current at ~1uA, if no light source is close by)
    - check picoscope for signals during entire cycle
    - picoscope: (ch1: E_x ,ch2: Ex, ch4: gate)
    - select central strip in y and 1st or 2nd x-strip as monitor output with epics (E127_epics > MSCF Si 1 & 2)
    - if you don't get any at all, switch on target in measurement cycle (or permantely)
    - still no nice signals? -> polarity is probably wrong, we need to switch MSCF of X and Y
    - during our measurement period, the signals should be well below 10V (goal for reaction ions ~6V)
    - adjust shaping time, gain, pole zero, threshold (epics) and Si-gate delay and length (e127.trlo) until you get nice histos
    - if you see large rates at injection or decelleration, we should try to shield detector with the scraper
    - put it as close as possible without disturbing the cycle, i.e. loosing more beam than without scraper
    - we can also move it by event (e.g. move in or out on start or end of machines 4,6,11,12,13) depending on when it disturbs

3. find beam position at detector
    - stop the cycle in SC 11
    - try to scrape it away by moving the detector our with small steps
    - you will need several injections, because if the short liftime, can get tedious
    - follow instructions and numbers here: https://elog.gsi.de/esr/E127/37

4. last things
    - ignore BaF for now
    - check that x-rays are working, spectra okay? resolution?
    - check rates during ON (x-rays: some 100 Hz each, Si: low, probably below 100Hz) and OFF (below 100Hz for all but BaF)
    - check deadtime, do we need downscaling for any detector? (i would not expect this, BaF is downscaled already, don't worry)

5. first runs
     - measure with and without scraper to confirm that Rutherford is under controle
     - how long? > depends on the statistics
     - try to see a clean effect of scraping, e.g. by the Si-plot E_ion vs x_position (Laszlo)
  22   Wed Sep 25 13:03:46 2019 LaszloAnalysissimulations on 118Te+p 
I have made some simulations at 8AMeV, 7AMeV and 4AMeV energies for the main 118Te + p reaction channels. For the
Rutherford simulations credits to Yuanming! For each energy there are 3 simulations: without the scraper, with the
scraper (online), scraper (online) + E-truncation (offline). The scraper is treated as an "ideal scraper" meaning
that there is no scattering at the edges (maybe worth a GEANT4?). The slit position is at -3cm from the beam, the vertical
length is 7cm (centralized). The following channels are combined in the
simulation:

-8AMeV:
  Rutherford
  pg channel, 5cascade model. Photon emission is treated isotropically
  pn channel, including: -->gs, -->1-->gs, -->2-->gs, -->3-->gs, -->3-->1-->gs, -->4-->1-->gs, -->5-->3-->gs,
-->5-->3-->1-->gs, -->5-->4-->1-->gs decay chains with their weights. Neutron and photon emission is treated isotropically.

-7AMeV:
  Rutherford
  pg channel, 4cascade model. Photon emission is treated isotropically

-4AMeV:
  Rutherford
  pg channel, 3cascade model. Photon emission is treated isotropically


The detector position is -2.5cm from the beam in the radial direction and centered vertically with a 45° tilt. 
The cross section values for pg and pn, the pn channel mixing, and the pg cascade number is based on TALYS simulations.
The other two input parameters are the luminosity and the measurement time. For the simulation, a little
pessimistic (or realistic?) scenario is taken: L = 10^24 cm-2s-1, and t = 10^4 s. From this values, from the N=CS*L*t
equation, the pg counts (based on the TALYS
cross sections) are the following: 8AMeV ~267, 7AMeV ~163, 4AMeV ~5! This means that to reach the Gamow-window
energies is challenging, but with a well-working scraper it is not impossible to reach.


As it is visible, the 8AMeV 118Te(p,g) case is very similar to the 7AMeV 124Xe case. The p,n threshold is
somewhere ~7.5MeV.
Attachment 1: 118Te_cs_energy_line_.png
118Te_cs_energy_line_.png
Attachment 2: 8AMeV_combined.png
8AMeV_combined.png
Attachment 3: 8AMeV_combined_spectrum_withoutTruncation.root
Attachment 4: 8AMeV_combined_spectrum_Truncation.root
Attachment 5: 8AMeV_combined_spectrum_ETruncation.root
Attachment 6: 7AMeV_combined_spectrum_withoutTruncation.root
Attachment 7: 7AMeV_combined_spectrum_Truncation.root
Attachment 8: 7AMeV_combined_spectrum_ETruncation.root
Attachment 9: 4AMeV_combined_spectrum_withoutTruncation.root
Attachment 10: 4AMeV_combined_spectrum_Truncation.root
Attachment 11: 4AMeV_combined_spectrum_ETruncation.root
  360   Fri May 21 01:17:37 2021 LaszloAnalysissimple analysis code2021
//a simple code (template) for offline analysis
//made by Laszlo, serves as a simple demonstration for enthusiastic shifters
//it creates a "no double counting Si 2D pos" histo
//usage:
//
//save the file as eg. "simple_code.c"
//root -l
//.L simple_code.c++
//run()
//when counter finished: ".q"

#define INPUT1 "input.root"//first unpack the lmd, then give the path of the unpacked .root file.

#define OUTPUT "./"//folder of the output. minimum input: "./"

#define ROOT_NAME "dummy.root"//name of the output

#include <cmath>
#include <string>
#include <sstream>
#include <cstdlib>
#include <cstdio>
#include <ctime>
#include <fstream>
#include <iostream>
#include <stdint.h>
#include "TROOT.h"
#include "TAttText.h"
#include "TAxis.h"
#include "TCanvas.h"
#include "TChain.h"
#include "TCut.h"
#include "TF1.h"
#include "TFile.h"
#include "TGraph.h"
#include "TGraphAsymmErrors.h"
#include "TGraphErrors.h"
#include "TH1.h"
#include "TH2.h"
#include "THistPainter.h"
#include "TKey.h"
#include "TLatex.h"
#include "TLegend.h"
#include "TMath.h"
#include "TMatrixD.h"
#include "TMinuit.h"
#include "TMultiGraph.h"
#include "TNtuple.h"
#include "TPave.h"
#include "TPaveText.h"
#include "TPoint.h"
#include "TRandom.h"
#include "TRint.h"
#include "TStyle.h"
#include "TString.h"
#include "TTree.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TSystem.h"
#include "TProfile.h"
#include "TVirtualFitter.h"
#include "TCanvas.h"
#include "TLegend.h"
#include "TColor.h"
#include <time.h>

using namespace std;

inline bool exists_test0 (const std::string& name) {
    ifstream f(name.c_str());
    return f.good();
}

///////////////////////////////////////////////////////////////////////////////////////////////

void loop(TChain *fChain){
	
	//setting pedestal values
	int PEDESTAL_LOW=400;
	int PEDESTAL_HIGH=8000;
	

	//this normally should be in a separate header, branches are defined.
	UInt_t          trigger;
	fChain->SetBranchAddress("TRIGGER",&trigger);
	UInt_t          E_SiY[17];
	UInt_t          E_SiX[17];
	UInt_t          t_SiY[17];
	UInt_t          t_SiX[17];
	for(int a=1;a<17;a++){
		fChain->SetBranchAddress(Form("E_SiY%d",a),&E_SiY[a]);
		fChain->SetBranchAddress(Form("E_SiX%d",a),&E_SiX[a]);
		fChain->SetBranchAddress(Form("t_SiY%d",a),&t_SiY[a]);
		fChain->SetBranchAddress(Form("t_SiX%d",a),&t_SiX[a]);
	}

	//creating histos
	TH2D *h_pos_si_xy=new TH2D("h_pos_si_xy", "h_pos_si_xy",16,0.5,16.5,16,0.5,16.5);

	//creating and initializing some variables used in the event loop (for "no double counting")
	int r_pos_x=0,r_pos_y=0;
	int dc_Ex_max=-999;
	int dc_Ey_max=-999;
	
	
	Long64_t nentries = fChain->GetEntries();
	Long64_t nbytes = 0;

	//starting the entry loop
	for (Long64_t i=0; i<nentries;i++){

		nbytes += fChain->GetEntry(i);
		//event countdown
		if ((float(i)/100000.)==int(i/100000)){cout << "event: " << i << " \tof " << nentries << endl;}


		if(trigger==1){//trigger 1 = TargetON

			dc_Ex_max=-999;
			dc_Ey_max=-999;

			for(int i_x=1;i_x<17;i_x++){
				for(int i_y=1;i_y<17;i_y++){

					if( ((int)t_SiX[i_x])>0 && ((int)t_SiY[i_y])>0){
						if( PEDESTAL_LOW<((int)E_SiX[i_x]) && PEDESTAL_HIGH>((int)E_SiX[i_x]) &&
						    PEDESTAL_LOW<((int)E_SiY[i_y]) && PEDESTAL_HIGH>((int)E_SiY[i_y])){

							//assign the hit to StripX and StripY where the most energy is deposited (rel. Ecalibration is needed, but roughly ok)
							if(dc_Ex_max<((int)E_SiX[i_x]) && dc_Ey_max<((int)E_SiY[i_y])){
								r_pos_x=i_x;
								r_pos_y=i_y;
								dc_Ex_max=E_SiX[i_x];
								dc_Ey_max=E_SiY[i_y];
							}
						
						}
					}
					
				}
			}

			//Filling pos. histo
			if(dc_Ex_max!=-999 && dc_Ey_max!=-999){h_pos_si_xy	-> Fill(r_pos_x,r_pos_y);}
			
		}//trigger==1

	}//entry loop
	

	//writing out the root output file
	TFile *graphfile = TFile::Open((OUTPUT + (string)("") + ROOT_NAME).c_str(), "RECREATE");
	graphfile -> mkdir("map");
	graphfile -> cd("map");
	h_pos_si_xy -> Write();
	graphfile -> Close();

}//loop

//////////////////////////////////////////////////////////////////////////////////////////////////////////////

void run(){

    const char *command = new char[1000];

    char filename[100];
    TChain *fChain = new TChain("h101");

    sprintf(filename,INPUT1);

    if(exists_test0(INPUT1) && exists_test0(OUTPUT)){

    cout<<"\033[0;37m//loading run: "<<filename << "\033[0m" <<endl;
    fChain->Add(filename);

    loop(fChain);
    }
    else{cout<<"\033[0;31mError 404: non-existing INPUT1 or OUTPUT file!\033[0m" <<endl;}

	command = "rm  *.so";
	gSystem->Exec(command);
	command = "rm  *.d";
	gSystem->Exec(command);

	return;
}

int main(){
    run();
    return(0);
}
  30   Thu Dec 5 10:57:57 2019 Jan, LaszloDetectorssetting up the DSSD with beam 
This is a brief how-to about getting the DSSSD set up with stored beam and target in the ESR.

The main goals
1. adjust timing signals to TDC-range
2. adjust energy signals to ADC-range
3. optimize energy resolution

The signal chain is the following:

det [32ch.] > preamp [2x16ch] > 2x MSCF shapers [16ch]
  244   Mon Mar 23 13:26:59 2020 JanAnalysisruns_90to99 SCRAPER RESULTS 
With scraper, we have roughly 150-200 (p,g) counts now.

X-rays look fine!
Attachment 1: 124Xe_Si_map.png
124Xe_Si_map.png
Attachment 2: 124Xe_x-ray_35.png
124Xe_x-ray_35.png
Attachment 3: 124Xe_x-ray_90.png
124Xe_x-ray_90.png
Attachment 4: 124Xe_x-ray_145.png
124Xe_x-ray_145.png
Attachment 5: 124Xe_Si_projection.png
124Xe_Si_projection.png
  251   Mon Mar 23 14:32:03 2020 JanAnalysisruns_100to104 NO-SCRAPER RESULTS 
Without scraper, we have very few (p,g) counts, hard to get numbers quickly, maybe ~50.

X-rays look fine!
Attachment 1: 124Xe_Si_map.png
124Xe_Si_map.png
Attachment 2: 124Xe_Si_projection.png
124Xe_Si_projection.png
Attachment 3: 124Xe_x-ray_35.png
124Xe_x-ray_35.png
Attachment 4: 124Xe_x-ray_90.png
124Xe_x-ray_90.png
Attachment 5: 124Xe_x-ray_145.png
124Xe_x-ray_145.png
  5   Fri Apr 5 05:46:20 2019 JanDAQrunning DRASI for (p,g) 
We now have a drasi DAQ running for our (p,g) experiments.
The RIO we have it installed right now is: r4l-58
It shall be used with the user <litv-exp>

The DAQ-controle shall be done from lxg1275 inside a SCREEN SESSION named "daq"
tab 0 >> drasi 
tab 1 >> logger
tab 2 >> rate
tab 3 >> stream server
tab 4 >> file writing (optional)

To check and control running daq:

login on lxg1275
> ssh litv-exp@lxg1275

enter screen session
> screen -x daq

DRASI
go to tab 0
> ctrl-a 0 
should be quite, no errors showing, no terminal outputs
full restart proceedure
> ssh r4l-58
> cd esrdaq_2018/r4l-58
> ./start.sh (for det. triggered daq)
> ./pulser.sh (for internal pulser trigger)
> ctrl-c (for killing daq)

LOGGER
go to tab 1
> ctrl-a 1
should be quite, no errors showing, no terminal outputs
full restart proceedure
> ssh r4l-58
> cd esrdaq_2018/r4l-58
> ./logger.sh

RATE MONITOR
go to tab 2
> ctrl-a 2
rates should update each second
"incr" should be event rate per sec 
"send" column should be reasonable value if data is written
full restart proceedure
> ssh r4l-58
> cd esrdaq_2018/r4l-58
> ./rate.sh

EMPTY STREAM SERVER
go to tab 3
> ctrl-a 3 
empty stream should be running  
number of processed events should keep increasing
you should see the number of clients connected
full restart proceedure (on lxg1275)
> empty stream://r4l-58 --server=stream:6002

WRITE A FILE
can be done by litv-exp from any lgx-machine
ROOT FILE:
> e127_unpack --stream=lxg1275 --ntuple=RAW,/data.local1/e127/pre_data/file.root
LMD FILE:
> e127_unpack --stream=lxg1275 --output=/data.local1/e127/pre_data/test.lmd
 
  123   Sat Mar 21 12:06:50 2020 LaszloGeneralrun61 - SCRAPER moved out 
we want to see any change on the Silicon count rate, because we will have now also Rutherford scattering. 
ELOG V3.1.5-fc6679b