Energy vs Time plots

I updated the monte carlo written by Alex to include time and also created new SRIM tables for 250 Torr rather than 300 Torr. I ran the Monte Carlo for a beam energy of 4.44MeV/u and plotted energy vs time for each detector. It should be noted that there are 2 codes, one for scattering 21Ne(p,p)and one for the reaction 21Ne(p,a) and the number of scattered particles as opposed to the number of reaction products from (p,a) are not what is expected in the data but the energies and timing are.

Attachment 1 shows energy vs time in S2-1 (downstream-heavy ion detector) from the data runs at 4.44Mev/u. Beam pulses are 85ns apart. The cuts on this include: equal energy cuts between the front and back of each detector removing interstrip hits, coincidence hits between S2-1 and S2-2, must be a hit in the front and back of each detector, phi cut so that particles must hit S2-1 and S2-2 at 180 degrees to one another. A feature not shown in the monte carlo is the second weaker line of particles, hitting the detector just after the first stronger line. I am not sure what this is, could be beam scattering off the entrance window or from N gas left over in the chamber?

Attachment 2 shows the same plot (energy vs time) simulated with the monte carlo for scattering reactions and for the (p,a) reaction. This shows that the reaction is dominated by beam scattering from the target H2 gas. The 18F extends to lower energies as shown in the monte carlo and there are some lower energy points in the data  (Att 1) which could correspond to this, however I am not sure whether this is F or beam scattering from the entrance window which also produces an extended range of energies. 

Attachment 3 shows energy vs time data plot in S2-2(upstream - light ion detection), at same energy and with same cuts as for s2-1 in Att 1. Again comparing this to attachment 4 which shows the Monte Carlo calculated energy vs time plot for S2-2 the data is dominated by scattering events. The small number of events which happen after 8.5 Mev (protons will punch through the detector after this) could be alphas. The Monte Carlo shows the alpha particles have a wider energy range than the protons. 

Other plots

I ran a 'gas out' data set to see how much scattering was coming from the entrance window to the chamber. There is no gas out run for 4.44 MeV/u so I used the one at 4.11MeV/u. This shows that there was some scattering from the window getting through the cuts described above, in S2-2 they go no higher than 8 MeV which makes me think it is scattering with protons because of the punch through energy. To check it was from the window I checked a run where there was no target and no window in the chamber and saw nothing after the cuts above. 

Other plots I am looking at include: time(S2-1) vs time (S2-2), energy(S2-1) vs energy(S2-2), energy vs strip for each detector, time of flight difference between S2-1 and S2-2 vs energy difference between two detectors, difference in theta (front strip) vs sum energy of two detectors, time of flight vs energy difference. 

eol

 I have attached some single strip energy vs time plots to show that the structure seen in the collective E vs t plots is in each strip. 

Attachments 1 and 2 show two different strips in S2-2 with no cuts. You can see that the 3 loci are different to one another in the single strips as well as in the combined S2-2 E vs t plot shown in attachment 3. Attachement 3 is plotted after applying the following cuts: equal energy in the front and back of each detector, low energy cut, must be a hit in each detector, multiplicity must be greater than 3. The 3 different loci are a beam pulse apart. The energy cal is 40keV/chan. 

Attachments 4 and 5 show energy vs time for single strips in S2-1 with no cuts. Attachment 6 shows the combined front strips in S2-1 after the same cuts as above for the combined S2-2 plot have been applied. Now there are 2 loci in the same beam pulse which are still there after the cuts. Not sure what the second line is. The energy cal in attachment 6 is 80keV/chan.

See attached.

Ran all night OK, currently 14epA, H2 stable, DAQ stable

ISSUES:

1. Analysis of S4200 in offline sort lead to conclusion that specttrum is dominated by protons, 8MeV punchthrough visible??.Not many events at high energies, where Alex MC predicts them to be.

2.Noticed that 3 channels in particular in S2-2 have very high rates. All in same bank - can we raise thresholds?  Strips 34,36,46 in S2-2, we think (from mapping in sort code)

3.  The loci in the E vs T spectrum are (a) still aparrently upside down from what expected and (b) very wide in time. There is also a suspiciously sharp edge to the timing peak. Look at eny E vs T loci in the sort

4.Been sorting over night - nearly caught up.

We looked a run (54) from the first H2 run at 4.4 MeV/u and compared with the current data. The The S4200 spectrum and the E vs T loci look exactly the same.

 OK - summary of events during this shift.

1.TDCs. We have been sorting offline, and we have now established that between run 66 and run 79, the TDC mode was incorrectly set (not in common stop mode) after a DAQ crash causing broad loci in the E vs T spectra and times of the wrong sign. This was corrected from Run 81 onwards. 

2. The offline sort is behaving very badly and many of the sorts run since yesterday to now are wrong. Make sure you (a) sort on a different machine than the DAQ and (b) make sure that the sort is killed each time you start a new sort and (c) make sure that it is sorting the righrt run number.

3.  We have changed the trigger to a coincidence (S2-1 plus S2-2) OR S2-2 on its own OR PD.

4. We have increased the beam rate to about 120 epA and we stay on this. Trigger rate should be about 7K with 4K live time.

5. We will not change energy until early morning tomorrow, at the earliest. we will make a decision during the next two shifts based on what we see in the sort.

Mike

1.TDCs. We have been sorting offline, and we have now established that between run 66 and run 79, the TDC mode
was incorrectly set (not in common stop mode) after a DAQ crash causing broad loci in the E vs T spectra and
times of the wrong sign. This was corrected from Run 81 onwards.

TDC channel 0 (reference channel) should have been recorded in the EbyE data - common stop TDC values
can therefore be calculated offline by switching the egid2logical logical variable csm from TRUE to FALSE.
  
2. The offline sort is behaving very badly and many of the sorts run since yesterday to now are wrong. Make sure
you (a) sort on a different machine than the DAQ and (b) make sure that the sort is killed each time you start a
new sort and (c) make sure that it is sorting the righrt run number.

Amen

 We have been running smoothly all night. Only problem is that the DAQ occasionally does not start properly after a stop-go. The dead-time goes to nearly 100%. This happened twice in the night,. and a second stop-go did the trick. Needs looking at.

We are going to stop now and (a) check the tune properly and (b) retune if necessary and then carry on the same energy. There is not really evidence in the data of the structure we are looking for, so we need to carry on at this energy for a while. This process will take 6 hours or so

MAB

 Update: 

W decided at 8.00 to stop, check the tune, retune and get back going. This has nearly, finished but it is a very long process indeed - note for future - this takes AT LEAST 12 hours, if one is systematic about it.

We tested the tune as it was - it was OK, (90% transmission or so acdording to a straight ratio of the cups) but there are a fair bit on the 4mm collimator. 

Marco then tuned, and improved it with a high intensity tune, reduced the rate on the collimator by at least a factor of 5 - maybe more. Good position on scintillator - see photo below.

We have since pumped out and run several tests with the window in, but no gas: (1) 10mm collimator on first target ladder, nothing on second - transmission was about the same, measured the same way. Rate similar to same test we did a couple of days ago. (2) 5mm collimator and blank second target - result about the same. (3) 5mm collimator and scintillator in - photo below. Beam spot blew up a bit compared with no gas. See below. (4) 5mm collimator and PD2 (5) 5mm collimator and blank. Ran this to disk. and (6) CH2 target. This made no sense and we got a lot of horrible ate of this so we stopped it.

We saw noise in S2-2 when we had no gas. It was bad - 1200Hz with no beam - all from S2-2. When we put gas in, it all went away! We suspect that it was a gauge that we disconnected, perhaps.

Waiting for beam, as buncher has now tripped.

Note that full cycle (pump, flush, vent, tune, flush, pump H2) - with all the tests as you go along, took 12 hours!

The pictures below are tunes without and with the window. From Mikes iphone. !

Buncher problem was resolved around 22:00 (driver module needed to be replaced).
Checked current on FC0 and trigger rate, around 4000 presented. We then started a run, only got 20 acquired,
stoped and restarted a run, then all relevant numbers seemed normal.
The pressure inside TUDA is slowly creeping down (~0.3 torr/h), maybe a small leak in the window. It's not the
hall cooling down, the gas temperature actually rose in the meantime.

Summary at start of run

Read through the runbook and end of run summar and the situation appears to be:

1. The check of the tune showed it was OK, although a bit of beam (or halo) clipping 4mm collimator

2. There are some changes to the gas filling scheme - Jamie should be present for the next one.

3. Evidence that the gas pressure is dropping slowly - need to watch this

4. Trigger is Coinc.OR.SS_2.OR.PD.  Rate at this is 200/s at 4epA (50/50 Coinc and S2_2)

5. Suggestion we could increase the beam a factor of two to get back to the same coinc rate.

Analysis.  Looking at the rates I see:

Running at 66 pnA and see about 4,000/s triggers and 3,200/s accepts. Scaler 9 over 40,000/s.

Looking back in run book the max trigger accept rate we see is about 3,800/s.  Could double beam and gain only

modest increase on accepted rate, but at cost of faster damage of detector.  Suggest we run

like this for a bit and make sure all is well.

Took a look through the spectra.  All ADCs OK, but no TDC signals.  Checked past saved runs.  There up

to run 98.  The next saved on disc is run 104 and they aren't there.  Nor in any saved runs after this.

Looked at electronics and reason clear as the module providing the TDC trigger (lower section of a

429 unit) has no input.  Eventually determined that it would make sense if this unit was operating

in 2 x 8 mode rather than what it is currently set on (4x4).  Switched this over and now working OK.

Assume this was changed inadvertently during all the work on the chamber yesterday.

0300.  Running again and all ADC and TDC spectra make sense.  Off for a cup of tea.

Running status

Everything seems to be running OK.  ADC and TDC spectra look OK.  Chamber pressure has stabalised.  Detector leakages stabalised

Beam now up to 60pnA with trigger rate about 4,000/s and accepts about 3,000/s

Have spent some time looking at sorted spectra and put notes in the runbook, along with a copy of an email from Tom with some suggestions.  Need someone who can change the sort programme. The key is probably in getting the plot of Tdiff vs Esum calibrated [Tdiff is time difference between S2-1 and S2-2 signals and Esum is the sum of their energies].  The problem is that the TDC signals have not yet been lined up.  We also need to understand the features we see in the E vs t plots, including the puzzeling feature that one of the loci has its "blob" of events at a different energy than the others.

My feeling is that we should continue to run ast this energy, while continuing to refine the sort to try and confirm we can select out the (a,p) events of interest.

Worth thinking as to whether we should take more beam (it is available.

 All ran smoothly.

As a test, ran for a short period at 1200 midday to run without S2-2 in the trigger. Finished that at 1339.

After discussion with Alison, Alex, Brian, Chris - we decided to change energy to 3.98 MeV/u. This started at 1340.

Pumped out H2, Marco steered with low intensity beam. He optimized on our PD2 with the 3mm collimator, then he observed the spot on the scintillator. He believed (as do I) that it looked as good as it did after the high-intensity tune yesterday. Beam spot blown up a bit, as before, but very well centred. picture from iphone below.

Did a run without gas, then filled, then got going on H2. H2 run started 1730. Thus, 4 hours for beam change and low-intensity tune.

Did analysis of rates/triggers with and without gas - all with 100 epA

                     WITHOUT GAS        WITH GAS

S2-2 all           3000                        3800

PD                   1500                        1900

Coincidence  10                            1800 

S2-1 all           13000                     57000

The trigger (7000 ish) is a sum of the first three of course. Interesting that all the coincidences are gas-"related" - which is good I think. It confirms to me that we should consider, at some point, running on only coincidences in the trigger to cut down the dead time (currently 60% live, 40% dead)

that's it. Jamie checking the spectra, if nothing else in this e-log then they are ok.

Sort started on new data, will run it overnight and compare in the morning with the old sort.

Mike

For run 123:

Beam outage for ~15 min but not stopping run.

FC0 = 120pA when we got beam back.

We'll let the run go 15 min longer.

Run 125 stopped early because of strange rate values on scalar 1.

Sorted Runs:

Finished 121

122

123

124

Run139

Started 19:05:03
Stopped 20:36:36

Rates were reaching well above 8500, ops were measuring 160 epa on FC0 so they decreased it to 120 epa before starting next run.

Run 140

Started 20:42:12    - Scalers were acting up, stopped the run
Stopped 20:55:50

Run 141

Started 20:56:14
Stopped at 22:50:02

Rates were too low (~4289), ops were measuring 70 epa on FC0, got ops to increase current to 120 epa before starting next run.

Run 142

Started 22:56:34
Stopped 23:32:13

Rates were reaching above 9000 again, ops were reading 150 epa on FC0.  Brought it back down to 115 epa before starting next run.

Run 143

Started 23:36:38 - scalers were acting up, stopped the run
Stopped 23:37:46

Run 144

Started 23:40:17


****submitting summary early in case browser crashes again.  will update with any changes.

Sarah informed me that Jamie was performing the sorting and saving of spectra differently and didn't know what exactly to do so I will not be performing any sorting.

01:40 Run 144 Stopped

01:45 Run 145 Started

03:43 Run 145 Stopped

Stripper foil changed, still seeing fluctuations so the problem could be in the SC-LINAC

03:57 Run 146 Started

Rate durring run 146 has been a little low, will ask for more at the run change.

05:50 Run 146 Stopped

FC0 averaged 65pA for run146.  Got it back up to 100pA for the next run.

05:54 Run 147 Started

06:15 No beam, lost a SC cavity.  Run 147 stopped

06:42 Beam back.  Run 148 started and stopped due to low scaler 1 reading

06:43 Run 149 Started

07:10 Rate dropped back down to ~4200, FC0 = 60pA.  Stopped run 149, saved and zeroed the histograms

FC0 back up to 100pA

07:13 Run 150 Started

No major problems presented, A pluser walk was performed at the begining of the shift run161 S2-1 front, run162 S2-2 front, run163 S2-1 back and run165 S2-2 back.

Runs 166, 168, 169 are void due to bad scaler reading

Two times the control room reported the dipole magnet tripped but beam was recovered within few minutes. FC0 readings range 150-120epA

S2-1 L. current 17.60uA @ 16:00 and 18.85uA @ 23:10

Beam current dropped to 85 epA at the start of the shift, operator checked for reasons and decided to change
stripper foil. That brought the current back up, but it started decreasing again, so from run to run slits were
gradually opened further to keep FC0 current at about 100 epA.
I forgot the part about saving and clearing the spectra after every run, so runs 171-177 were summed together.

 Energy calibration for S2-2 was done using the pulser data (set 6 of the pulser runs) to calculate the offsets and alpha run R186 to calculate the gains and constant. For S2-1 pulser walkthrough 6 was also used to calculate the offsets and the gains calculated using run 31 which is 21Ne beam scattering from Au/C target. The constants were then calculated by fitting a least squares fit to the alpha peaks from run 186 and the gold scattering peak from run 31 energy vs peak position plot.

TDCs were calibrated with different sources for each detector. For S2-1 the offsets were calculated using the TDC spectra for the 21Ne on Au/C (Run 31) data. For S2-2 the run 23 was used which is a pulser run.

Currently applying equal energy and coincidence cuts to the data and plotting using Monte Carlo to compare to what I am seeing in the data. Plotting things such as Energy in S2-1 vs Energy in S2-2, back strip in S2-1 vs back strip in s2-2, energy sum.

I have sorted the 4.11MeV/u data runs (the ones which had the common stop mode on the TDCs) and applied the following cuts:

1. There must be a hit registered in each side of the detectors.
2. equal energy must be deposited in both sides of the detectors.
3. Low energy cut is applied to each detector.
4. Must have a hit in the TDC for the front strips.
5. Must have multiplicity greater than 3.
6. Must have hits in both S2-1 and S2-2.
7. Particles must hit at 180 degrees to in the detectors. (back spoke cut)

Looking at the back spoke in S2-1 vs back spoke in S2-2 plot cut number 7 should be changed to +/- 3 back strips to ensure no particles are being cut out due to position straggling in the gas.

Plotting S2-2 energy vs time with these cuts gives attachment 1. I need to redo the timing calibration for S2-2 as strip 96 has the wrong offset and the resolution needs to be improved. The protons can be seen clearly as they punch through the detector and around 8MeV. I have also attached as attachment 2 the monte carlo plot for S2-2 at this energy. There are two monte carlo codes one for 21Ne beam scattering and one for the reaction which means the number of each type of event are not correctly proportional, in the data there is far more scattering. Here you can see the protons from ~4-8MeV lying underneath the alpha line. The time spread in the actual data of the proton/alpha line is about 25ns whereas in the monte carlo it is about 10 ns for combined alpha and protons. 

Because the protons will punch through the detector at approx 8MeV I believe that the points after this could be alphas. Attachment 3 shows the projection for the polygon which selects particles after the protons in the energy vs time plot. As can be seen from the energy projection the energy dips where you would expect it to for alphas, however there are some higher energy points and I am not sure what these might be.

I then plotted Energy vs time in S2-1 gating on energy in S2-2 greater than 8.2MeV. This plot is shown in attachment 6. Attachment 4 shows the Monte Carlo simulated energy vs time for 18F in S2-1, with attachment 5 shows simulated energy vs time for scattered 21Ne in S2-1. I plotted these separately as they have the similar timing and so cant be seen separately on the same plot. Here you can see the 21Ne clearly in the data, so some of the protons are depositing more than 8.2 MeV in S2-2. If I change the cut so that the energy in S2-2 must be greater than 9MeV and plot energy vs time for S2-1 again the 21Ne is no longer visible (will not let me attach this file due to capacity so I will attach it in a new elog entry straight after this one). The energy spread of the points in the S2-1 energy vs time plot cut on what could be alphas in S2-2 gives the correct energy spread for 18F in the monte carlo simulated data.

Next:

• Redo timing calibration for S2-2
• Look at gas out run
• Work out how to separate lower energy alphas and protons
• Work out what higher energy points in S2-2 could be.

 Attachment 1 shows S2-1 energy vs time with a gate on energy in S2-2 being greater than 9MeV i.e. See previous elog also posted today for explanation.