HOME ESR EXPERIMENTS LABS etc
E121 E127 E132 E125 E143 laser_cooling_2021 E142 Ê128 E146 E0052 E0018 E0028 E0038
  proton-capture on 118Te, Page 16 of 29  Not logged in ELOG logo
Entry  Fri Apr 5 06:05:55 2019, Jan, DAQ, making changes to the DAQ ,  
If you need to add or change a module of the VME crate you need to follow this proceedure:

1. stop DRASI
2. add/change module in esrdaq_2018/r4l-58/main.cfg 
   a. add module name and VME address: CAEN_V830 (0x00A00000) {}
   b. use BARRIER keyword between the modules
3. restart DRASI (e.g. ./pulser.sh)
4. check empty unpacker data output
   > empty --stream=r4l-58 --print --data
   should run without error and now have your new module event data
5. adjust unpacker in ~/e127/unpackexps/esr
   a. make sure you have a .spec file defining readout of you module, e.g. "vme_caen_v775.spec" 
   b. edit main_user_esr.spec
      - add module as e.g. "tdc[0] = VME_CAEN_V775(geom=0x2, crate=0);"
      - select "geom=0x1" value according to ascending order of all modules (0x1, 0x2, 0x3, ...)
   c. edit det_mapping.hh 
      - add SIGNAL() keyword for your module
      - e.g.: SIGNAL( TDC_1, esr.tdc[0].data[0],TDC_32, esr.tdc[0].data[31], DATA12); 
      - choose data format, e.g. "DATA12", as specified in module .spec-file
   d. compile your unpacker: 
      > make -j
   e. try it out:
      > ./esr --stream=r4l-58 --print --data
      - should run without errors 
      > ./esr --stream=r4l-58 --ntuple=RAW,/dir/file.root
      - root-file should have your new branches with your module data
Entry  Wed Jun 5 11:29:13 2019, Jan, Detectors, position/channel calib.,  
This is the first position calibration done with the detector installed in the holder. The goal is to assign rough positions (top/bottom/left/right) to PreAmp- and Shaper-channels.


PreAmp Channel polarity position
1-1 (Y-1) neg. bottom
1-16 (Y-16) neg. top
2-1 (X-1) pos. left (close to beam)
2-16 (X-16) pos. right (far from beam)
Entry  Fri Jun 7 10:18:23 2019, Jan, Runs, prerun001 - MixedAlpha,  alpha_centered.jpg
save point: lxg1275:/datalocal1/e127/predata/prerun001-root

~2h run

DAQ setup:
32 ADC channels
 ch. 1-16 (section 0) >> y-strips
 ch. 17-32 (section1) >> x-strips
32 TDC channels
 ch. 1-16 (section 0) >> y-strips
 ch. 17-32 (section1) >> x-strips
32 scaler channels (not used)
16 MDPP channels (not used)

Detector setup:
bias -90V
MixedAlpha source ~in central position
Entry  Fri Jun 7 13:52:32 2019, Jan, Runs, prerun002 - MixedAlpha,  alpha_top_left.jpg
save point: lxg1275:/datalocal1/e127/predata/prerun002

~20min run

DAQ setup:
32 ADC channels
 ch. 1-16 (section 0) >> y-strips
 ch. 17-32 (section1) >> x-strips
32 TDC channels
 ch. 1-16 (section 0) >> y-strips
 ch. 17-32 (section1) >> x-strips
32 scaler channles 
 (not used)

Detector setup:
bias -90V
MixedAlpha source in upper left corner
Entry  Fri Jun 7 14:28:09 2019, Jan, Runs, prerun003 - MixedAlpha,  alpha_bottom_right.jpg
save point: lxg1275:/datalocal1/e127/predata/prerun003

~20min run

DAQ setup:
32 ADC channels
 ch. 1-16 (section 0) >> y-strips
 ch. 17-32 (section1) >> x-strips
32 TDC channels
 ch. 1-16 (section 0) >> y-strips
 ch. 17-32 (section1) >> x-strips
32 scaler channles 
 (not used)

Detector setup:
bias -90V
MixedAlpha source in lower right corner
Entry  Fri Jun 7 14:39:36 2019, Jan, Detectors, Pt100 installation & readout,  det_connector.pngPt100_readout.JPG
The new detector holder incorporates a Pt100 temperature sensor, which is UHV compatible. The Pt-resistor is read out via 4 wires. These wires are integrated into the D-SUB connector of the detector. The new wires use the pins as shown in the attached picture
"det_connector.png".

Additionally, a new cable for air-sided connection to the feedthrough is made. This one has an additional outgoing branch for the 4-wire Pt-100 readout, which can be connected to the readout device.

The readout device is an Omega DP32Pt, which provides pre-calibrated conversion to temperature values for many sensor types.
https://www.omega.de/temperature/pdf/DPPT_SERIES.pdf
https://assets.omega.com/manuals/M5460_DE.pdf
https://data2.manualslib.com/pdf6/133/13210/1320926-omega/cn32pt.pdf?0f771b894ebb783ac172a013ea140d1d
Entry  Fri Jun 7 15:57:41 2019, Jan, DAQ, DAQ with & without MDPP-16,  
The DAQ located in 

../esrdaq_2018/r4l-58/ 
currently includes the readout of
MADC
Caen 775 TDC
Caen 830 Scaler
MDPP-16

../esrdaq_2018/r4l-58_rewind/
is without the MDPP-16:
MADC
Caen 775 TDC
Caen 830 Scaler
Entry  Tue Jun 11 10:13:56 2019, Jan, Runs, prerun003 - MixedAlpha,  
save point: lxg1275:/datalocal1/e127/predata/prerun004
~10min run

DAQ setup:
32 ADC channels
 ch. 1-16 (section 0) >> y-strips
 ch. 17-32 (section1) >> x-strips
32 TDC channels
 ch. 1-16 (section 0) >> y-strips
 ch. 17-32 (section1) >> x-strips
32 scaler channles 
 (not used)

Detector setup:
bias -90V
MixedAlpha source in center
new air-side cable used (with temp. sensor)
Entry  Tue Jun 25 21:19:50 2019, Laszlo, Simulations, MOCADI input parameters,  template.pdfsteck2004.pdf
initial BEAM input parameters for MOCADI are based on:
- https://web-docs.gsi.de/~weick/mocadi/download/esr-exl-test.in --> twist parameters (ratio of X/A and B/Y)
- M. Steck et al, Electron cooling experiments at the ESR, NIM A 532 (2004) 357-365 --> epsilon_x = 5·10^-7
m*rad, dp/p = 10^-4

Example input for 6AMeV 111Sn beam:

BEAM
1000000
6 , 0 , 110.8803121305 , 50
4
0.27386, 0.182573, 0, 0, 0
4
0.131477, 0.383186, 0, 0, 0
1
0.02, 0, 0, 0, 0


For more explanation please see the attached pdf.
Entry  Fri Jun 28 14:06:28 2019, Jan, General, CAD of (p,g) setup,  
 
Entry  Mon Jul 8 14:57:07 2019, Yuanming Xing, Detectors, The measurement of the tilting angle,  TiltingAngle.png
As shown in the figure, the detector is tilted by a certain angle which is about 45 degree. 
The tilted angle is defined as angle B.
We measured angle A using the vertical laser line (red) and papers with a straight line edge. So we can marked
the line denoted by the laser line when the paper is put on the holder.
In the beginning, the angle A is measured as 43.3±0.5 degree. So angle B= 46.7±0.5 degree
After the adjustment of the detector (to make it in the center of the chamber), the angle A is measured as
43.1±0.5 degree with another paper.  So angle B= 46.9±0.5 degree. 
Entry  Tue Jul 9 15:26:59 2019, Yuanming Xing, Simulations, Simulated Rutherford scattering distribution on detector,  6MeV.png7MeV.png8MeV.png
Simulated Rutherford scattering distribution on detector

   Beam: 111Sn
   Luminosity: 3.08 /mb/s
   
   Distance between the detector and central orbit: 2.80cm

   Particle Number in the simulation: 5 Million

   Beam quality used in MOCADI input file:
4,
0.27386, 0.182573, 0, -4.0, 0
4,
0.131477, 0.383186, 0, 0.0,  0
1,
0.02, 0, 0, 0.0, 0


1. Beam Energy: 6MeV/u
   
2. Beam Energy: 7MeV/u 

3. Beam Energy: 8MeV/u 
Entry  Tue Jul 23 10:46:29 2019, Laszlo, Analysis, 124Xe pg peak fit using mixed cascade simulation,  7x
This entry is the continuation of the cascade effects on the pg peak shape entry:

https://exp-astro.physik.uni-frankfurt.de:8080/E108/475

The main point is that having more cascade makes the impulse carried by one photon smaller --> recoil cone is
smaller --> pg peak gets more centered

When one wants to make an "all-inklusive" simulation the problem comes, that there are too many populated states
after CN* decay (many are just theoretical), also many cascades afterwards as well.
Instead, one can make a "statistical approach":
 1, take an average 1st populated state after the gamma emission of the CN*.
 2, take an average of the gamma energies after the p-capture. (it includes also the photons from the CM* to the
1st state after g emission with energy of [(E_CM - Q) - E_1st_state] --> should these be excluded from the avg.
gamma energy calculation?)
 3, the excited compound nucleus first decays to 1, then cascading down with the avg energy of 2,

    number_of_cascades = E_1st_state/E_avg_gamma. 

    For the 124Xe at 7AMeV using the TALYS code the above values are:
    E_1st_state = 6.74 MeV
    E_avg_gamma = 2.75 MeV
    number_of_cascades = 2.45

To reproduce the resulting avg. pg peak, I took a 2 and a 3 cascade simulations produced by MOCADI and mixed
them as 55% of the 2-cascade and 45% of the 3-cascade for this naiv, statistical model.
Then, to test the peakshape, I fit this mixed distribution to a part our measurement data at 7AMeV after a
sloppy Rutherford background removal (the background removal can be more improved, I used this as a quick method
for now). The calculated Chi^2/NDF ~ 3, which is not a super good value, but taking into account that this is a
strongly naiv model, and the background removal can be also improved, the chi^2/NDF value can be also
interpreted as promising.
Entry  Fri Aug 9 13:14:09 2019, Laszlo, , pn channel simulations for 111Sn(p,g)112Sb experiment,  
The p,n channel opening for the 111Sn+p reaction is at 5.8895 MeV/u beam energy. Therefore, when choosing as 6MeV/u, 7MeV/u, 8MeV/u beam energies, there is not much kinetic energy left for the outgoing neutron --> there are "less" nuclear levels populated above groundstate in 111Sb, they can be handled explicitly. My simulations are based on TALYS results. The following steps did I made:

1, choose an energy for the 111Sn beam. for this example let it be 8AMeV. 
Entry  Fri Aug 9 13:43:15 2019, Laszlo, Simulations, pn channel simulations for 111Sn(p,g)112Sb experiment,  12x
The p,n channel opening for 111Sn beam is at 5.8895 MeV/u beam energy. Therefore, when choosing 6MeV/u, 7MeV/u,
8MeV/u beam energies, there is not much kinetic energy left for the excitation of the 111Sb nucleus --> there
can be "less" nuclear levels populated above groundstate of 111Sb. These "less" states can be handled explicitly
(no continuum levels) with the TALYS code. I make the following steps to simulate the (p,n) products for a given
energy:

1, choose an energy for the 111Sn beam. for this example let it be 8AMeV. Then I use my energy calculator:

/u/lvarga/public/lab2cm_updated 111Sn 8

this gives the last output: " equvivalent (E_lab)^rel for p in talys = 8.06689 MeV ". This energy is the value
what i should put into TALYS as the p projectile energy. In TALYS we have only normal kinematics.

2, make the calculation with the TALYS code for outpopulation, outgamdis and outangle. For the first two cases,
an extra command is given: "maxlevelsbin n 30" to get discrete levels.

3, based on the outpopulation and outgamdis I create an excellsheet which tells me the probability of decay
happening after the n-emission.
   from outpopul, I can tell which level is populated in which % after the n-emission. (using the cross sections
and normalizing to the summed cross section)
   from outgamdis, I can tell that for a specific level after the n-emission what is the probability of the
decay afterwards. (using the cross sections given). 

4, regarding outangle, the cross sections are given for each 2. degree angle only. To get a continuous
dependency, I fit these with a 6. order polinomial function for each level. Later on, I use this polinomial
function for angular corrections.

5, for each decay scheme, I make a MOCADI code. For the target, instead of Hydrogen mass, I am using the proton
mass (otherwise, TALYS and MOCADI will not be compatible, ~511keV gap will remain in the CM energy, which of
course highly influencing the kinematics). For the nuclear masses, my calculator can be used:

/u/lvarga/public/atomic2nuclearMass 111Sn

It is important, that after the first n-emission I put a save point (which is the SAVE #1)for the later angular
corrections.

6, I mix the root file outputs from MOCADI based on the probabilities from my excell table. Also angular
corrections can be given using the (180-tof[1]) angle. There is not much difference however, but the computation
time increases dramatically.
//in the simulations at the PIN-diode position there is a scraping edge 3cm away from the beam.
Entry  Thu Aug 29 10:02:22 2019, Jan, Detectors, bakeout at ESR setup,  T_set_bakeout.png
The bakeout for our detector setup has started on the evening of 27.08.2019 
We have a temperature sensor in the vacuum next to the detector.
The goal is to maximize bakeout temperature, while not going above 125°C at the detector/sensor.

Here is a daily protocol of temperatures and vacuum pressure
T_12 is the set value of the 1st and 2nd controller device
T_34 is the set value of the 3rd and 4th controller device 
T_56 is the set value of the 5th and 6th controller device
T_sens is the temperature measured by the sensor close to the detector

day   date        P_setup    T_sens  T_12   T_34   T_56   comment
1     tue 28.08.  1e-6 mbar  ~100°   ~120°  ~120°  ~120°  carefull start
2     wed 29.08.  3e-7 mbar  ~121°   ~180°  ~150°  ~150°  T_sens reached ~126° during NEG-pump conditioning
3     thu 30.08.  6e-7 mbar  ~124°   ~185   ~150°  ~175°  readjusted T_56, outgassing a bit increased
7     mon 02.09.  5e-9 mbar   ~80°                        during cooldown
9     wed 04.09.  2e-10mbar   ~30°                        after cooldown, NexTorr Pump active     
Entry  Wed Sep 25 13:03:46 2019, Laszlo, Analysis, simulations on 118Te+p,  11x
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.
Entry  Wed Sep 25 13:58:24 2019, Laszlo, General, Gamow window energies for 118Te(p,g),  PhysRevC.81.045807.pdf
# Z   A   Iso Re Rc   T9 upper width mxpos  shift
 52 118 te118 pg 20  0.5  1.39  0.58  1.06  -0.01
 52 118 te118 pg 20  1.0  2.27  1.04  1.69  -0.01
 52 118 te118 pg 20  1.5  3.04  1.46  2.21  -0.01
 52 118 te118 pg 20  2.0  3.71  1.82  2.68  -0.01
 52 118 te118 pg 20  2.5  4.32  2.16  3.10  -0.02
 52 118 te118 pg 20  3.0  4.85  2.44  3.47  -0.06
 52 118 te118 pg 20  3.5  5.35  2.73  3.79  -0.12
 52 118 te118 pg 20  4.0  5.79  2.97  4.08  -0.19
 52 118 te118 pg 20  5.0  6.55  3.38  4.61  -0.35

source: https://journals.aps.org/prc/supplemental/10.1103/PhysRevC.81.045807
Entry  Wed Oct 23 11:40:29 2019, Jan, Runs, prerun005 - MixedAlpha,  
save point: lxg1275:/datalocal1/e127/predata/prerun005

start: 23.10.19 - 11:40 
stopp: 23.10.19 - ~13:30

DAQ setup:
32 ADC channels
 ch. 1-16 (section 0) >> x-strips
 ch. 17-32 (section1) >> y-strips
32 TDC channels
 ch. 1-16 (section 0) >> x-strips
 ch. 17-32 (section1) >> y-strips
32 scaler channels 
 (not used)

Detector setup:
bias -30V
MixedAlpha source in center
Entry  Wed Oct 23 14:51:43 2019, Jan, Laszlo, DAQ, VME modules & main.cfg,  main.cfg
Here is the current list of VME modules in our DRASI daq including HW addresses:

Slot <module> [HW addr.]

1 <RIO4>
2 <Vulom4b> [0500]
3 <ENV3>
5 <scaler V830> [00E1]
7 <scaler V830> [00E2]
15 <TDC V775> [00C1]
17 <TDC V775> [00C2]
19 <MADC> [00A1]
21 <MADC> [00A2]

See also main.cfg attached
ELOG V3.1.5-fc6679b