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

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

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

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

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

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

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)

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.

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. 

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 

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.

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. 

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.

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     

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.

# 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

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