We are now using new SpecAmps for the x-ray detectors from ORTEC, they deliver much improved resolution but have other strange effects as discribed in the next entry.

for some reason the event number of the target ON and OFF was misaligned with the pattern.
This means that the last few runs (run53, run52...) measured not in the target on phase but one before (at 10 instead of 11). Therefore, with these 
runs we measured the TargetOFF phase instead TargetON (that is why the pg peak seem to disappear...)

The measurements were taken by using a caliper ruler, with precision of +/-0.05mm

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

1. before using stoch. cooling in ESR, boot and update the windows PC
2. better diagnostics for low energy/intensity beam-target overlapp in ESR is needed (either a clever analysis or a new detector)
more to come

I have made also the inverse square law fits. We have data only for 90degre with the 241Am source, but both for the 1. and 2. calibrations. The 1. and 2. calibration data sets treated separate.
2 peaks were investigated, 59.5keV and 26.3keV, at two distances 184.8mm and 217.3mm. These distances are the sum of 4 distances:
a = width of plastic head. uncert.: +/- 0.05mm, measured with caliper.
b = width of the brass collimator. uncert.: +/- 0.05mm, measured with caliper.
c = width of protector plastic ring. uncert.: +/- 0.05mm, measured with caliper.
d= distance of the paper head. uncert: +/- 0.5mm, judged by the eye :)

The Be window and the dead layer of the Ge detector is not taken into account.

The distances calculated as D=a+b+c+d.

The conclusion is that the uncertainty coming from the distance measurements are negligible compared to the other uncertainties. The data obey the inverse square law.


> For all 3 detector the calibration data sets were combined to include the systematics in the fit results directly. Combining means not a weighted average, just simple all data points were included into the fit --> doubled efficiency value for most of the energies.
> The g(x) = a * (1-exp(-(x-c)/b)) * exp(-x/d) function was used to describe the behavior of the Germanium detectors for the whole range of energies (global behavior). From these fits the 80keV outlier point was excluded. This is very strange that it doesn't follow the trend, it would be nice to find out why not.
> Between 40-75keV a linear fit was carried out as well. This can also approximate quite well this local energy range, what we need for the K-REC peaks.
> l(x) = m*x+e
> 
> All the fits were done by gnuplot, but it was also confirmed that ROOT gives us the same parameters + errors + chisquare. One just need to choose well the starting values :)
> 
> 
> 90° fits:
> degrees of freedom    (FIT_NDF)                        : 8
> rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.419915
> variance of residuals (reduced chisquare) = WSSR/ndf   : 0.176329
> p-value of the Chisq distribution (FIT_P)              : 0.994094
> 
> Final set of parameters            Asymptotic Standard Error
> =======================            ==========================
> a               = 0.00308376       +/- 0.0005144    (16.68%)
> c               = 15.6259          +/- 2.035        (13.03%)
> b               = 9.36888          +/- 3.16         (33.73%)
> d               = 177.141          +/- 84.95        (47.95%)
> 
> degrees of freedom    (FIT_NDF)                        : 4
> rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.283372
> variance of residuals (reduced chisquare) = WSSR/ndf   : 0.0802998
> p-value of the Chisq distribution (FIT_P)              : 0.988405
> 
> Final set of parameters            Asymptotic Standard Error
> =======================            ==========================
> m               = -6.22672e-06     +/- 1.836e-06    (29.48%)
> e               = 0.00256426       +/- 9.873e-05    (3.85%)
> 
> 
> 
> 145° fits:
> degrees of freedom    (FIT_NDF)                        : 6
> rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.379472
> variance of residuals (reduced chisquare) = WSSR/ndf   : 0.143999
> p-value of the Chisq distribution (FIT_P)              : 0.990247
> 
> Final set of parameters            Asymptotic Standard Error
> =======================            ==========================
> a               = 0.0012346        +/- 0.0006794    (55.03%)
> c               = 11.3754          +/- 4.075        (35.82%)
> b               = 15.6961          +/- 13.2         (84.1%)
> d               = 116.084          +/- 99.83        (86%)
> 
> degrees of freedom    (FIT_NDF)                        : 2
> rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.327728
> variance of residuals (reduced chisquare) = WSSR/ndf   : 0.107405
> p-value of the Chisq distribution (FIT_P)              : 0.898161
> 
> Final set of parameters            Asymptotic Standard Error
> =======================            ==========================
> m               = -2.39246e-06     +/- 1.278e-06    (53.42%)
> e               = 0.000848859      +/- 7.329e-05    (8.634%)
> 
> 
> 
> 35° fit:
> degrees of freedom    (FIT_NDF)                        : 7
> rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.40676
> variance of residuals (reduced chisquare) = WSSR/ndf   : 1.97896
> p-value of the Chisq distribution (FIT_P)              : 0.0538635
> 
> Final set of parameters            Asymptotic Standard Error
> =======================            ==========================
> a               = 0.000114396      +/- 0.0002367    (206.9%)
> c               = 14.3178          +/- 8.322        (58.13%)
> b               = 16.7454          +/- 43.21        (258%)
> d               = 79.7991          +/- 173.2        (217.1%)
> 
> degrees of freedom    (FIT_NDF)                        : 3
> rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.21233
> variance of residuals (reduced chisquare) = WSSR/ndf   : 1.46975
> p-value of the Chisq distribution (FIT_P)              : 0.220529
> 
> Final set of parameters            Asymptotic Standard Error
> =======================            ==========================
> m               = -3.7879e-07      +/- 3.143e-07    (82.99%)
> e               = 7.31447e-05      +/- 1.788e-05    (24.45%)
> 
> 
> 
> > For the efficiency vs E fit of the 90degree Xray detector I have used the following phenomenological funciton:
> > 
> > f(x) = a * (1-exp(-(x-c)/b)) * exp(-x/d)
> > 
> > Here the first exponent member is a saturation curve. This part describes the passing through of the two Be windows (chamber + before detector) and through the dead layer of Ge crystal. One needs a minimum energy to enter to the detecting Ge crystal = C parameter. b parameter = characteristic absorbtion E of these nondetecting layers.
> > The second exponent is an exponential decrease of the detector efficiency. Photons with higher energy are less detectable by the germaniums. The d parameter is the characteristic E for hard Xray and gamma (>40keV) detectability. 
> > 
> > https://www.amptek.com/internal-products/si-pin-vs-cdte-comparison
> > 
> > //Jan's comment: the tail of this function should more or less follow a linear trend a bit above than 40 keV.
> > 
> > 
> > In the attachment there is an example fit for 90 degree with combined 1. and 2. (before and after beamtime) calibration datasets. 
> > I made the fit with gnuplot:
> > 
> > degrees of freedom    (FIT_NDF)                        : 8
> > rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.419915
> > variance of residuals (reduced chisquare) = WSSR/ndf   : 0.176329
> > p-value of the Chisq distribution (FIT_P)              : 0.994094
> > 
> > Final set of parameters            Asymptotic Standard Error
> > =======================            ==========================
> > a               = 0.00308376       +/- 0.0005144    (16.68%)
> > c               = 15.6259          +/- 2.035        (13.03%)
> > b               = 9.36888          +/- 3.16         (33.73%)
> > d               = 177.141          +/- 84.95        (47.95%)
> > 
> > Laszlo's out.

We are going into the ESR to check the electronics.
Use calibration sources (133Ba etc).

Yesterday evening we have found out that the energy signals from the silicon are negative --> wrong polarity of the shapers was chosen while using the jumpers.
Cables coming from the preamp are switched now. --> in the recorded data the X and Y coordinates of the Si strip detector is the other way around!
Hopefully, the left and right side we can judge well...

new event server is set up at lxg1299 on port 6003
It runs as:
/data.local1/gstore/mrevserv64 lxg1275

the attached efficiency curve for the 90 deg. xray detector is made with the runs
e127b_run0008.lmd
e127b_run0009.lmd
e127b_run0010.lmd

All of them were taken in the lab (not in ESR).

This is a very preliminary analysis!
E.g. deadtime has been estimated and was not cleanly determined.

the attached efficiency curve for the 35 deg. xray detector is made with the runs
e127b_run0015.lmd
e127b_run0016.lmd
e127b_run0017.lmd + e127b_run0019.lmd

All of them were taken in the lab (not in ESR).

This is a very preliminary analysis!
E.g. deadtime has been estimated and was not cleanly determined.

For all 3 detector the calibration data sets were combined to include the systematics in the fit results directly. Combining means not a weighted average, just simple all data points were included into the fit --> doubled efficiency value for most of the energies.
The g(x) = a * (1-exp(-(x-c)/b)) * exp(-x/d) function was used to describe the behavior of the Germanium detectors for the whole range of energies (global behavior). From these fits the 80keV outlier point was excluded. This is very strange that it doesn't follow the trend, it would be nice to find out why not.
Between 40-75keV a linear fit was carried out as well. This can also approximate quite well this local energy range, what we need for the K-REC peaks.
l(x) = m*x+e

All the fits were done by gnuplot, but it was also confirmed that ROOT gives us the same parameters + errors + chisquare. One just need to choose well the starting values :)


90° fits:
degrees of freedom    (FIT_NDF)                        : 8
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.419915
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.176329
p-value of the Chisq distribution (FIT_P)              : 0.994094

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 0.00308376       +/- 0.0005144    (16.68%)
c               = 15.6259          +/- 2.035        (13.03%)
b               = 9.36888          +/- 3.16         (33.73%)
d               = 177.141          +/- 84.95        (47.95%)

degrees of freedom    (FIT_NDF)                        : 4
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.283372
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.0802998
p-value of the Chisq distribution (FIT_P)              : 0.988405

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
m               = -6.22672e-06     +/- 1.836e-06    (29.48%)
e               = 0.00256426       +/- 9.873e-05    (3.85%)



145° fits:
degrees of freedom    (FIT_NDF)                        : 6
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.379472
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.143999
p-value of the Chisq distribution (FIT_P)              : 0.990247

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 0.0012346        +/- 0.0006794    (55.03%)
c               = 11.3754          +/- 4.075        (35.82%)
b               = 15.6961          +/- 13.2         (84.1%)
d               = 116.084          +/- 99.83        (86%)

degrees of freedom    (FIT_NDF)                        : 2
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.327728
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.107405
p-value of the Chisq distribution (FIT_P)              : 0.898161

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
m               = -2.39246e-06     +/- 1.278e-06    (53.42%)
e               = 0.000848859      +/- 7.329e-05    (8.634%)



35° fit:
degrees of freedom    (FIT_NDF)                        : 7
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.40676
variance of residuals (reduced chisquare) = WSSR/ndf   : 1.97896
p-value of the Chisq distribution (FIT_P)              : 0.0538635

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 0.000114396      +/- 0.0002367    (206.9%)
c               = 14.3178          +/- 8.322        (58.13%)
b               = 16.7454          +/- 43.21        (258%)
d               = 79.7991          +/- 173.2        (217.1%)

degrees of freedom    (FIT_NDF)                        : 3
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.21233
variance of residuals (reduced chisquare) = WSSR/ndf   : 1.46975
p-value of the Chisq distribution (FIT_P)              : 0.220529

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
m               = -3.7879e-07      +/- 3.143e-07    (82.99%)
e               = 7.31447e-05      +/- 1.788e-05    (24.45%)



> For the efficiency vs E fit of the 90degree Xray detector I have used the following phenomenological funciton:
> 
> f(x) = a * (1-exp(-(x-c)/b)) * exp(-x/d)
> 
> Here the first exponent member is a saturation curve. This part describes the passing through of the two Be windows (chamber + before detector) and through the dead layer of Ge crystal. One needs a minimum energy to enter to the detecting Ge crystal = C parameter. b parameter = characteristic absorbtion E of these nondetecting layers.
> The second exponent is an exponential decrease of the detector efficiency. Photons with higher energy are less detectable by the germaniums. The d parameter is the characteristic E for hard Xray and gamma (>40keV) detectability. 
> 
> https://www.amptek.com/internal-products/si-pin-vs-cdte-comparison
> 
> //Jan's comment: the tail of this function should more or less follow a linear trend a bit above than 40 keV.
> 
> 
> In the attachment there is an example fit for 90 degree with combined 1. and 2. (before and after beamtime) calibration datasets. 
> I made the fit with gnuplot:
> 
> degrees of freedom    (FIT_NDF)                        : 8
> rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.419915
> variance of residuals (reduced chisquare) = WSSR/ndf   : 0.176329
> p-value of the Chisq distribution (FIT_P)              : 0.994094
> 
> Final set of parameters            Asymptotic Standard Error
> =======================            ==========================
> a               = 0.00308376       +/- 0.0005144    (16.68%)
> c               = 15.6259          +/- 2.035        (13.03%)
> b               = 9.36888          +/- 3.16         (33.73%)
> d               = 177.141          +/- 84.95        (47.95%)
> 
> Laszlo's out.

For the efficiency vs E fit of the 90degree Xray detector I have used the following phenomenological funciton:

f(x) = a * (1-exp(-(x-c)/b)) * exp(-x/d)

Here the first exponent member is a saturation curve. This part describes the passing through of the two Be windows (chamber + before detector) and through the dead layer of Ge crystal. One needs a minimum energy to enter to the detecting Ge crystal = C parameter. b parameter = characteristic absorbtion E of these nondetecting layers.
The second exponent is an exponential decrease of the detector efficiency. Photons with higher energy are less detectable by the germaniums. The d parameter is the characteristic E for hard Xray and gamma (>40keV) detectability. 

https://www.amptek.com/internal-products/si-pin-vs-cdte-comparison

//Jan's comment: the tail of this function should more or less follow a linear trend a bit above than 40 keV.


In the attachment there is an example fit for 90 degree with combined 1. and 2. (before and after beamtime) calibration datasets. 
I made the fit with gnuplot:

degrees of freedom    (FIT_NDF)                        : 8
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.419915
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.176329
p-value of the Chisq distribution (FIT_P)              : 0.994094

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 0.00308376       +/- 0.0005144    (16.68%)
c               = 15.6259          +/- 2.035        (13.03%)
b               = 9.36888          +/- 3.16         (33.73%)
d               = 177.141          +/- 84.95        (47.95%)

Laszlo's out.

after like 20min run, the density of the gasjet already dropped from 10^14 to 8.7*10^13. Now it seems to remain in that value

Below in arbitrary units one can see the trend of the gas jet density starting from 10^14 to 8.7*10^13

For the 124Xe primary beam measurement:
-with the DSSSD we scraped the beam at position -40
-position of the detector is set to -25. (1.5cm away from the beam)
-with the scraper (at Eggelhof 1 position) we scraped the beam at ~(-15)-(-13)
-position of the scraper is set to +20 (3.5cm away from beam axis (sollbahn))

For the 118Te fragment measurement:
-with the DSSSD we scraped the beam at position -39
-position of the detector is set to -24 (1.5cm away from the beam)
-with the scraper (at Eggelhof 1 position) we scraped the beam at ~ -14
-position of the scraper is set to +19 (3.5cm away from beam axis (sollbahn))

measured beam diagnosis with 124Xe48+ at 326 MeV/u