Difference between revisions of "Faraday Cup Analysis"

From New IAC Wiki
Jump to navigation Jump to search
 
(59 intermediate revisions by one other user not shown)
Line 1: Line 1:
[http://wiki.iac.isu.edu/index.php/PhotoFission_with_Polarized_Photons_from_HRRL Go Back]
+
[https://wiki.iac.isu.edu/index.php/Roman_calculation Go Back]
  
 
=Kicker magnet calibration=
 
=Kicker magnet calibration=
Line 21: Line 21:
 
1. x axis: ADC mean value for one pulse<br>
 
1. x axis: ADC mean value for one pulse<br>
 
2. y axis: number of pulse w/ that ADC mean value<br>
 
2. y axis: number of pulse w/ that ADC mean value<br>
[[File:1477_1.png]][[File:1473_1.png]]<br>
+
[[File:1477_1.png|200 px]][[File:1473_1.png|200 px]]<br>
[[File:1461_1.png]][[File:1465_1.png]]<br><br>
+
[[File:1461_1.png|200 px]][[File:1465_1.png|200 px]]<br><br>
  
 
==by using ADC channel charge distribution==
 
==by using ADC channel charge distribution==
Line 39: Line 39:
 
[[File:Fc_poisson_error.png]]
 
[[File:Fc_poisson_error.png]]
 
<br><br>
 
<br><br>
 
  
 
=Faraday Cup 3D plot=
 
=Faraday Cup 3D plot=
Line 46: Line 45:
 
[[File:Far.jpg]]
 
[[File:Far.jpg]]
  
=What about moving beam?=
+
=What about moving beam? splitting beam?=
=="rain" plots==
+
==in time charge distribution plots (Run 1477 -5A)==
 +
[[File:ADC6.png|200 px]][[File:ADC7.png|200 px]][[File:ADC8.png|200 px]]<br>
 +
[[File:ADC9.png|200 px]][[File:ADC10.png|200 px]][[File:ADC11.png|200 px]][[File:ADC12.png|200 px]]<br><br>
 +
As we can see from plots above:<br>
 +
1. There are no beam on blocks 6 and 12.<br>
 +
2. There are beam on blocks (7, 8, 9, 10)<br>
 +
3. The total charge for blocks (7, 8, 9, 10, 11) is strongly fluctuates with time.<br>
  
As we can see from plots above, the total charge for one block is strongly fluctuates with time.<br>
+
Also it is interesting to note the switching in mean value of charge distribution from "up" (blocks 7,8) to "down" (blocks (9,10,11)<br>
Question: Is it moving beam or is it current fluctuation in beam pulse? Take in, for example 1477.dat(-5A)<br>
+
This switching in mean value position is easy to see if we plot the histogram of charge distribution. (don't mind about that, wrong observation)<br><br>
 +
[[File:ADC6h.png|200 px]][[File:ADC7h.png|200 px]][[File:ADC8h.png|200 px]]<br>
 +
[[File:ADC9h.png|200 px]][[File:ADC10h.png|200 px]][[File:ADC11h.png|200 px]][[File:ADC12h.png|200 px]]<br><br>
  
==correlation plots==
+
'''Question:'''<br>
 +
'''1. Is it moving beam?'''<br>
 +
'''2. Or is it maybe splitting of beam into two beams with very unstable current fluctuation?'''<br><br>
 +
 
 +
==correlation plots (Run 1477 -5A)==
 
[[File:ADC6-ADC8.png|200 px]][[File:ADC6-ADC9.png|200 px]][[File:ADC6-ADC10.png|200 px]]<br>
 
[[File:ADC6-ADC8.png|200 px]][[File:ADC6-ADC9.png|200 px]][[File:ADC6-ADC10.png|200 px]]<br>
 
[[File:ADC7-ADC8.png|200 px]][[File:ADC7-ADC9.png|200 px]][[File:ADC7-ADC10.png|200 px]][[File:ADC7-ADC11.png|200 px]][[File:ADC7-ADC12.png|200 px]]<br>
 
[[File:ADC7-ADC8.png|200 px]][[File:ADC7-ADC9.png|200 px]][[File:ADC7-ADC10.png|200 px]][[File:ADC7-ADC11.png|200 px]][[File:ADC7-ADC12.png|200 px]]<br>
 
[[File:empty.png|200 px]][[File:ADC8-ADC9.png|200 px]][[File:ADC8-ADC10.png|200 px]][[File:ADC8-ADC11.png|200 px]][[File:ADC8-ADC12.png|200 px]]<br>
 
[[File:empty.png|200 px]][[File:ADC8-ADC9.png|200 px]][[File:ADC8-ADC10.png|200 px]][[File:ADC8-ADC11.png|200 px]][[File:ADC8-ADC12.png|200 px]]<br>
 
[[File:empty.png|200 px]][[File:empty.png|200 px]][[File:ADC9-ADC10.png|200 px]][[File:ADC9-ADC11.png|200 px]][[File:ADC9-ADC12.png|200 px]]<br>
 
[[File:empty.png|200 px]][[File:empty.png|200 px]][[File:ADC9-ADC10.png|200 px]][[File:ADC9-ADC11.png|200 px]][[File:ADC9-ADC12.png|200 px]]<br>
[[File:empty.png|200 px]][[File:empty.png|200 px]][[File:empty.png|200 px]][[File:ADC10-ADC11.png|200 px]][[File:ADC10-ADC12.png|200 px]]<br>
+
[[File:empty.png|200 px]][[File:empty.png|200 px]][[File:empty.png|200 px]][[File:ADC10-ADC11.png|200 px]][[File:ADC10-ADC12.png|200 px]]<br><br>
 +
 
 +
As we can see from plots above, there are strong correlations between blocks (7, 8, 10, 11)<br>
 +
There are weak correlation between blocks (9) and (7, 8, 10, 11)<br>
 +
And there are no correlation between blocks (6, 12) and all other blocks.<br>
  
As we can see for case above above, there are strong correlations between blocks 7, 8, 9 and 11<br>
+
''For moving beam I would expect the negative correlation between blocks (7, 8, 9, 10), if somewhere decrease, somewhere increase''<br>
 +
'''So because the positive correlation between blocks (7, 8, 10, 11)'''<br>
 +
'''we have beam splitting into two beams with very unstable current?'''<br>
  
==Conclusion:==We have possibly moving up and down beam
+
This structure can be seen if I plot histogram for charge distribution for all ADC blocks<br>
 +
[[File:ADC channel distribution.png|200 px]]<br>
 +
Here is interesting to note that the total charge for first peak (blocks 7, 8) is about 3 times more then for second peak (blocks 10, 11)<br>
 +
It would be interesting to see this ratio for all other kicker magnet current values. Is it constant?<br>
 +
Also it would be interesting to see the structure of the splitting for different current values.<br>
 +
Let see how it looks<br><br>
  
 
= Error Analysis=
 
= Error Analysis=
Line 114: Line 136:
 
Let see what the total charge looks like for 3 cases.
 
Let see what the total charge looks like for 3 cases.
  
[http://wiki.iac.isu.edu/index.php/FC_Analysis Go Up]
+
[http://wiki.iac.isu.edu/index.php/Faraday_Cup_Analysis Go Up]
 
[http://wiki.iac.isu.edu/index.php/PhotoFission_with_Polarized_Photons_from_HRRL Go Back]
 
[http://wiki.iac.isu.edu/index.php/PhotoFission_with_Polarized_Photons_from_HRRL Go Back]

Latest revision as of 19:03, 24 May 2012

Go Back

Kicker magnet calibration

by using pulse by pulse ADC channel mean value distribution

For each beam pulse:

[math] X^{pulse}_{avg}=\frac{\sum_{i=1}^{16}{Q_{i}*i}}{\sum_{i=1}^{16}{Q_{i}}}[/math] (average beam position for one pulse)

For distribution over all beam pulses:

[math] X_{avg}=\frac{\sum_{i=1}^{pulses}{X_{avg}^{pulse}}}{pulses}[/math] (average beam position over all pulses)
[math] X_{err}={\sqrt{\frac{1}{pulses}\sum_{i=1}^{pulses}{\left(X^{pulse}_{avg}-X_{avg}\right)^{2}}}}[/math] (absolute error)

Here is:
1. ADC# = bridge#
2. Pulse# = ReadOut# = Entry# = Event#

FC data 23.png

FC plot 2 4.png



Some examples of ADC mean value distribution. Here are:
1. x axis: ADC mean value for one pulse
2. y axis: number of pulse w/ that ADC mean value
1477 1.png1473 1.png
1461 1.png1465 1.png

by using ADC channel charge distribution

For each ADC channel:

[math] Q_{channel}=\sum_{i=1}^{pulses}{Q^{pulse}_{channel}}[/math]     (total charge per one ADC channel)
[math] q_{channel}=\sum_{i=1}^{pulses}{\sqrt{Q^{pulse}_{channel}}}[/math] (absolute error in charge distribution per one ADC channel assuming Poisson)
For each ADC channel we have (Q +- q)

For distribution over all ADC channel:

[math] X_{avg}=\frac{\sum_{i=1}^{16}{iQ_{i}}}{\sum_{i=1}^{16}{Q_{i}}}[/math]     (average beam position over all pulses)
[math] X^2_{err}=\frac{\sum_{i=1}^{16}{i^2q^2_{i}}}{\sum_{i=1}^{16}{Q^2_{i}}}+\frac{(\sum_{i=1}^{16}{iQ_{i}})^2\sum_{i=1}^{16}{q^2_{i}}}{\sum_{i=1}^{16}{Q^4_{i}}}[/math] (absolute error)
For all ADC channel we have [math]ADC_{avg}[/math]+-[math]ADC_{err}[/math]

Here are results:
Fc poisson error.png

Faraday Cup 3D plot

Below is the plot of the charge in Faraday cup (pC) as a function of magnet current (vertical axis, A) (basically magnetic field) and ADC (horizontal axis).

Far.jpg

What about moving beam? splitting beam?

in time charge distribution plots (Run 1477 -5A)

ADC6.pngADC7.pngADC8.png
ADC9.pngADC10.pngADC11.pngADC12.png

As we can see from plots above:
1. There are no beam on blocks 6 and 12.
2. There are beam on blocks (7, 8, 9, 10)
3. The total charge for blocks (7, 8, 9, 10, 11) is strongly fluctuates with time.

Also it is interesting to note the switching in mean value of charge distribution from "up" (blocks 7,8) to "down" (blocks (9,10,11)
This switching in mean value position is easy to see if we plot the histogram of charge distribution. (don't mind about that, wrong observation)

ADC6h.pngADC7h.pngADC8h.png
ADC9h.pngADC10h.pngADC11h.pngADC12h.png

Question:
1. Is it moving beam?
2. Or is it maybe splitting of beam into two beams with very unstable current fluctuation?

correlation plots (Run 1477 -5A)

ADC6-ADC8.pngADC6-ADC9.pngADC6-ADC10.png
ADC7-ADC8.pngADC7-ADC9.pngADC7-ADC10.pngADC7-ADC11.pngADC7-ADC12.png
Empty.pngADC8-ADC9.pngADC8-ADC10.pngADC8-ADC11.pngADC8-ADC12.png
Empty.pngEmpty.pngADC9-ADC10.pngADC9-ADC11.pngADC9-ADC12.png
Empty.pngEmpty.pngEmpty.pngADC10-ADC11.pngADC10-ADC12.png

As we can see from plots above, there are strong correlations between blocks (7, 8, 10, 11)
There are weak correlation between blocks (9) and (7, 8, 10, 11)
And there are no correlation between blocks (6, 12) and all other blocks.

For moving beam I would expect the negative correlation between blocks (7, 8, 9, 10), if somewhere decrease, somewhere increase
So because the positive correlation between blocks (7, 8, 10, 11)
we have beam splitting into two beams with very unstable current?

This structure can be seen if I plot histogram for charge distribution for all ADC blocks
ADC channel distribution.png
Here is interesting to note that the total charge for first peak (blocks 7, 8) is about 3 times more then for second peak (blocks 10, 11)
It would be interesting to see this ratio for all other kicker magnet current values. Is it constant?
Also it would be interesting to see the structure of the splitting for different current values.
Let see how it looks

Error Analysis

The ADC measures the charge deposited on each of the 16 Aluminum blocks. The ADC is a 12 bit ADC with a max input of 400 pc. THe means that the charge per channel is:

[math]\frac{400 \mbox{pC}}{2^{12}\mbox{Channels}} = \frac{400 \mbox{pC}}{4096\mbox{Channels}} =\frac{1 \mbox{pC}}{10.24 \mbox{Channels}} [/math]

Run1477 -5A

Raw ADC result for channel 8 :

Run1477 ADC8 Raw.png

ADC 8 using channel -> Coul conversion:

Run1477 ADC8 Raw pC.png

The above histogram shows an RMS of 66.74 pC.

Comment

The ADC basically counts the number of electrons collected by the aluminum FC bricks and transfered through the cables to the ADC. This is a poisson process with a large number of trials leading to a large mean value [math](\mu)[/math]. One would expect a gaussian parent population with
[math]\sigma = \sqrt{\mu}[/math]

The histogram shows

Mean = [math]172.7 =\frac{\sum_i^N x_i}{N}[/math]
RMS = [math]66.74 = \frac{\sqrt{\sum_i^N(x_i - 172.7)^2}}{N}[/math]
[math] \sqrt{172.7} = 13.1 \ll 66.74[/math]
What does this mean?
The theoretical distribution would be Guassian with [math]\sigma = \sqrt{\mu}[/math]. The above suggests that the beam charge delivered to the FC is not following this statistical parent distribution. Most likely the beam current is changing while you are measuring the charge with the FC.
Is the charge lost or is it just moving to different FC channels

Run1477 TotalCharge pC.png


Run1477 TotalCharge vsTime.png

Run1477 TotalCharge vsRFpulse.png

Within 100 RF pulses the total charge on the FC drops from 600 to 100 pC. It seems the beam current is very unstable.


There also appears to be a gaussian distribution centered around 600 pC and a wider one centered around 100 pC.


Could it be that I can put a cut on total beam current and see if the charge moves between FC elements?


Let see what the total charge looks like for 3 cases.

Go Up Go Back