General remarks to the MOCADI simulations:

1. scraper size vertically > (p,g) spot at the scraping position --> (p,g) can be separated by truncating the backscattered Rutherford events in the energy.


2. replacing the atomic masses to nuclear masses does not change significantly the absolute positions of the (p,g) and the Rutherford. However, the (p,g) spot in the simulation seems bigger using the nuclear masses, than in the 124Xe experiment with a factor ~1.6.


3. going for lighter elements, the separation between (p,g) and Rutherford getting better, than slit can be placed even more far from the beam in x (radial) direction. For 124Xe the minimal distance from the beam axis is <4cm, for 73As it is <7.5cm.


4. in the simulations the scraping is mostly until x=-inf. However, if the scraping is incomplete in x, we can end up with underlying background events below the p,g peak (see the last slides)! The Rutherford cone gets bigger by going down with A,Z. For 91Nb the minimal width of the scraper: x>6cm, for 73As x>7cm. Therefore, I suggest to have a scraper width x=9cm. For the y, y=6cm should be safe.


5. for lighter elements, the (p,g) spot size increases: in the simulation for 73As it reaches the detector size. However, the (p,g) spot size might be overestimated, please read the 2. point.

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.

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 

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.

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

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)

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

save point: lxg1275:/datalocal1/e127/predata/prerun007

start: Fr 25.10.19 - 16:35 
stopp: 

DAQ setup:
64 ADC channels
 ch. 01-64 (ADC0) >> empty
 ch. 33-48 (ADC1 - section0) >> x-strips
 ch. 49-64 (ADC1 - section1) >> y-strips
32 TDC channels
 ch. 01-64 (TDC0) >> empty
 ch. 33-48 (TDC1 - section0) >> x-strips
 ch. 49-64 (TDC1 - section1) >> y-strips
64 scaler channels 
 (not used)

Detector setup:
bias -30V
MixedAlpha source in center

beam: 118Te52+
energy: 7.4 MeV/u

purpose: data with TARGET ON

Detector position (Si): xxx mm

run start at 14:21 , first file: run0xx_0001.lmd
run stop  at 17:11 , last file:  run088_0178.lmd

_______________________________________________________________________________

time: 14:45

ON Rates
Si_X:       29 Hz
Si_Y:       31 Hz
Xray_35:   592 Hz
Xray_90:   168 Hz
Xray_145:   xx Hz
BaF_OR:    xx Hz

Si voltage:          x.x V  
Si leakage current:  x.x uA


SIS particles before:           5e9
ESR particles at injection:     1e6
ESR particles after decel.:     1e5
Target ON density:              5.5e13


copy and repeat the above (below the line) every 30 minutes

Setting for the machines.

beam: 124Xe54+
energy: 10 MeV/u

purpose: data with TARGET ON

Detector position (Si): -25 mm

run start at 04:31 , first file: run046_0001.lmd
run stop  at xx:xx , last file:  run0xx_0xxx.lmd

_______________________________________________________________________________

time: 04:31:44

ON Rates
Si_X:       1 Hz
Si_Y:       1 Hz
Xray_35:   3 Hz
Xray_90:   33 Hz
Xray_145:   23 Hz
BaF_OR:    2463 Hz

Si voltage:          60.02 V  
Si leakage current:  3.93 uA


SIS particles before:           2e9
ESR particles at injection:     5e7
ESR particles after decel.:     2e6
Target ON density:              1.8e13


copy and repeat the above (below the line) every 30 minutes

beam: 124Xe54+
energy: 10 MeV/u

purpose: data with TARGET ON

Detector position (Si): -25 mm

run start at 04:31 , first file: run046_0001.lmd
run stop  at xx:xx , last file:  run0xx_0xxx.lmd

_______________________________________________________________________________

time: 05:00

ON Rates
Si_X:       1 Hz
Si_Y:       1 Hz
Xray_35:   3 Hz
Xray_90:   28 Hz
Xray_145:   24 Hz
BaF_OR:    2452 Hz

Si voltage:          60.02 V  
Si leakage current:  3.93 uA


SIS particles before:           2e9
ESR particles at injection:     5e7
ESR particles after decel.:     2e6
Target ON density:              1.8e13


copy and repeat the above (below the line) every 30 minutes

beam: 124Xe54+
energy: 10 MeV/u

purpose: data with TARGET ON

Detector position (Si): -25 mm

run start at 05:30 , first file: run047_0001.lmd
run stop  at xx:xx , last file:  run0xx_0xxx.lmd

_______________________________________________________________________________

time: 05:30

ON Rates
Si_X:       0 Hz
Si_Y:       0 Hz
Xray_35:   2 Hz
Xray_90:   27 Hz
Xray_145:   20 Hz
BaF_OR:    2379 Hz

Si voltage:          60.02 V  
Si leakage current:  3.93 uA


SIS particles before:           2e9
ESR particles at injection:     5e7
ESR particles after decel.:     5e6
Target ON density:              1.97e13


copy and repeat the above (below the line) every 30 minutes

beam: 124Xe54+
energy: 10 MeV/u

purpose: data with TARGET ON

Detector position (Si): -25 mm

run start at 05:30 , first file: run047_0001.lmd
run stop  at xx:xx , last file:  run0xx_0xxx.lmd

_______________________________________________________________________________

time: 06:00

ON Rates
Si_X:      1 Hz
Si_Y:       1 Hz
Xray_35:   3 Hz
Xray_90:   33 Hz
Xray_145:   23 Hz
BaF_OR:    2512 Hz

Si voltage:          60.02 V  
Si leakage current:  3.94 uA


SIS particles before:           2e9
ESR particles at injection:     5.3e7
ESR particles after decel.:     5.8e6
Target ON density:              1.97e13


copy and repeat the above (below the line) every 30 minutes

beam: 124Xe54+
energy: 10 MeV/u

purpose: data with TARGET ON

Detector position (Si): -25 mm

run start at 06:33 , first file: run048_0001.lmd
run stop  at xx:xx , last file:  run0xx_0xxx.lmd

_______________________________________________________________________________

time: 06:36

ON Rates
Si_X:      0 Hz
Si_Y:       0 Hz
Xray_35:   0 Hz
Xray_90:   39 Hz
Xray_145:   16 Hz
BaF_OR:    2590 Hz

Si voltage:          60.02 V  
Si leakage current:  4.0 uA


SIS particles before:           2e9
ESR particles at injection:     5.8e7
ESR particles after decel.:     5.2e6
Target ON density:              1.55e13


copy and repeat the above (below the line) every 30 minutes

beam: 124Xe54+
energy: 10 MeV/u

purpose: data with TARGET ON

Detector position (Si): -25 mm

run start at 06:33 , first file: run048_0001.lmd
run stop  at xx:xx , last file:  run0xx_0xxx.lmd

_______________________________________________________________________________

time: 07:00

ON Rates
Si_X:      0 Hz
Si_Y:       0 Hz
Xray_35:  1 Hz
Xray_90:   39 Hz
Xray_145:   38 Hz
BaF_OR:    2498 Hz

Si voltage:          60.02 V  
Si leakage current:  4.03 uA


SIS particles before:           2e9
ESR particles at injection:     5.8e7
ESR particles after decel.:     5.4e6
Target ON density:              1.57e13


copy and repeat the above (below the line) every 30 minutes

beam: 124Xe54+
energy: 10 MeV/u

purpose: data with TARGET ON

Detector position (Si): -25 mm

run start at 07:31 , first file: run049_0001.lmd
run stop  at xx:xx , last file:  run0xx_0xxx.lmd

_______________________________________________________________________________

time: 07:30

ON Rates
Si_X:      0 Hz
Si_Y:       0 Hz
Xray_35:  0 Hz
Xray_90:   35 Hz
Xray_145:   23 Hz
BaF_OR:    2483 Hz

Si voltage:          60.02 V  
Si leakage current:  4.02 uA


SIS particles before:           2.3e9
ESR particles at injection:     5.8e7
ESR particles after decel.:     5.8e6
Target ON density:              1.35e13


copy and repeat the above (below the line) every 30 minutes