A fifth analysis of 2D GEM position resolution

Dean Karlen / January 15, 2001

This document summarizes the analysis of the GEM data taken November 20-24, 2000, following the general approaches described in the earlier analyses. This analysis is intended to form the basis of a publication. The gas is P10 and the preamps were ALEPH. The scope time scale was decreased by a factor of 2 compared to the October data. The data recorded here therefore correspond to 250 MHz sampling. Sets of calibration data taken on December 14 and January 5 are also examined.
 

Index


 

GEM pad layout

The figure linked here shows the GEM layout and coordinate system used in the analysis. The coordinate system (shown in blue) is centred on pad 1. The ticks are shown with 0.1 mm spacing. The pads are numbered from 1 to 8, according to the readout channel. The central pad is read out by both oscilloscopes (channels 1 and 5), to provide a common trigger. The red dots show the locations of the x-ray collimator during the data taking.

GEM data

The data sets taken on November 20-24 were taken with P10 gas with Vdrift = 3475V, Vgem=3375, using the ALEPH preamplifier, and the xray tube set at 6 kV. At total of 750 triggers are taken for each run. Some 10% of triggers are not coincident between the scopes - therefore fewer than 750 events are analyzed in each event. The analysis also rejects events that fail in the fits of the pulses. In total, about 37000 events are analyzed.

The 60 data runs are summarized below. The collimator location is indicated using the coordinate system described above. The first two digits of the run number is the day of the month that the data was taken.

Before the run 2403, the xay collimator was brought back to the vertex position. A realignment was necessary at that time to bring the signal sizes on the 3 pads to agreement. The x and y position of the collimator was shifted in the negative direction by 30 microns. Unfortunately, pad-pad gain measurements were not made. The GEM system failed after the runs recorded here, requiring the GEM box to be opened for investigation of the problem.

As an example, the first event from run 2201 is shown here.

run number x_coll (mm) y_coll (mm)
2001 0. 1.443
2002 0. 1.043
2003 0. 0.843
2004 0. 0.643
2005 0. 0.443
2006 0. 0.243
2007 0. 0.043
2008 0. -0.157
2009 0.200 -0.157
2010 0.200 0.043
2011 0.200 0.243
2012 0.200 0.443
2013 0.200 0.643
2014 0.200 0.843
2015 0.200 1.043

run number x_coll (mm) y_coll (mm)
2101 0.200 1.243
2102 0.200 1.443
2103 0. 1.443
2104 0. 1.643
2105 0. 1.243
2106 0.400 1.243
2107 0.400 1.043
2108 0.400 0.843
2109 0.400 0.643
2110 0.400 0.443
2111 0.400 0.243
2112 0.400 0.043
2113 0.400 -0.157
2114 0.600 -0.157
2115 0.600 0.043
2116 0.600 0.243
2117 0.600 0.443
2118 0.600 0.643


run number x_coll (mm) y_coll (mm)
2201 0.600 0.843
2202 0.600 1.043
2203 0.800 0.843
2204 0.800 0.643
2205 0.800 0.443
2206 0.800 0.243
2207 0.800 0.043
2208 0.800 -0.157
2209 1.000 -0.157
2210 1.000 0.043

run number x_coll (mm) y_coll (mm)
2301 1.000 0.243
2302 1.000 0.443
2303 1.000 0.643
2304 1.200 0.443
2305 1.200 0.243
2306 1.200 0.043
2307 1.200 -0.157
2308 1.400 -0.157
2309 1.400 0.043
2310 1.400 0.243
2311 1.600 0.043


run number x_coll (mm) y_coll (mm)
2401 1.600 -0.157
2402 1.800 -0.157
2403 0.0 1.443
2404 0.0 1.243
2405 0.0 1.043
2406 0.0 0.843

 

 

GEM data analysis programs

The results shown below come from the gemanal program (version 0.9) located in the directory /home/karlen/gem. An associated paw kumac file, gem.kumac, is found in the same area.

Separation of direct and induced components of signals

As in the previous analyses, the direct charge component of a signal is deduced from the amplitude measured a fixed time after the peak (referred to as the "late" amplitude). In this analysis the delay is chosen to be 600 ns.

The left plot from the figure linked here shows the mean ratio of the "late" to the peak amplitudes on pads 1,2, 7 and 8 as a function of the distance from the centre of the pad to the collimator position. The error bars indicate the standard deviations of the ratio. (The mean and standard deviations are found by fitting each ratio distribution to a Gaussian). The standard deviation is less than 1%; so the late amplitude is used instead of the peak amplitude for the charge fraction determination of both direct and mixed signals.

The ratio for pad 1 is independent of the location of the xray collimator provided all of the charge is collected by the pad. Surprisingly this ratio initially increases as the xray collimator gets beyond 0.7 mm from the pad centre. The scaling factor for the "late" amplitude for pad 1 is taken to be 0.681 from this data. For the other channels the scaling factor is somewhat larger. The right plot from the figure linked here shows the ratios after a multiplicative correction is applied to bring the channels into better agreement. The reason that channel 1 is different from the rest may be due to the fact that the scope is terminated with 1 MOhm for that channel (because the signal is shared with channel 5) whereas the rest of the channels are terminated with 50 Ohm.

The scaling factors deduced for channels 1,2,5,7,8 are 0.681, 0.686, 0.689, 0.691, 0.689 respectively.

Gain variation

The gain of the system was not constant over these runs, as can be seen in the figure linked here. The plots show the total charge collected by all pads, as deduced by the "late" amplitude and scaled by the factors from the previous section. There is a rapid rise in the gain during the runs on the first day. Gain fluctuations continued afterwards.

Pedestals

The data for each channel is corrected by using a pedestal defined by the average of measurements before the pulse (time bins 10-60).

Some plots have been produced to check if the baseline changes from run to run: The pulse shape is fit to a quadratic 300ns - 900 ns after the induced pulse. The value of the fit at the point 600 ns after the pulse defines the "baseline". The value of the baseline for the different runs are shown in the figure linked here. The baseline is seen to be independent of position of the x-ray collimator as long as it is more than about 3.1 mm from the centre of the pad for the for pads, 2,3,4, and 6. The other pads do not have measurements with the x-ray collimator far enough away. There is no evidence for cross talk, which in the past showed up as a dependence of the baseline from one channel on the amplitude of the pulse on another channel.

Position analysis from direct charge sharing

Observed charge fraction in pad 1 - determination of cloud size
The figure linked here shows the observed charge fraction in pad 1, as a function of the y-coordinate of the collimator. The left plot is the data for runs on November 20 and 21, the right plot is for runs on November 24. The solid curve shows the prediction of the model when sigma_x=0.60 mm (left figure) and sigma_x=0.56 (right figure). The dashed curves show the predictions for sigma_x  increased and decreased by 0.03 mm from these values. Apparently the transverse diffusion has reduced over the period of the 5 day run. For the analyses presented here, the average cloud size is assumed to be 0.58 mm.
Determining position from charge fractions
Events with the signal shared amongst exactly three pads are used to estimate the location of each x-ray event, according the the procedure applied in the earlier analyses. The criteria used in this analysis to define that a pad has direct charge is that the value of the "late" amplitude is below -0.020 V. In earlier analyses, the cut value was -0.010 V. The value was changed in this analysis to reduce the number of events designated as having 4 pads that share charge.

The efficiency for events to have its position determined in this way is summarized in the plot linked here. The lines indicate the pad boundaries and the efficiencies are printed at the micrometer readings of the x-ray collimator. Low efficiency is found in regions where the charge is shared amongst fewer or more than 3 pads. Some inefficiency would be recovered by considering events with 4 pads that share the charge.

The figure linked here shows histograms of the x and y coordinate estimates, with respect to the x and y collimator position, for run 2106. Fitting the distributions to Gaussians, gives central values of (0.007 mm,0.022 mm) and standard deviations of 66 and 64 microns in x and y, respectively.

There are regions where the biases in x and y for this measurement are large, as shown in the figures linked here. The standard deviations of the measurements are between 50 and 80 microns throughout the region. The observed biases can result from gain inequalities from pad to pad. It appears, for example, that the gain in pad 7 is smaller than for other pads, resulting in biases towards -x in regions where pad 7 is used to estimate the position.

In the section on position analysis from induced pulses below, the relative gains of the outer 6 pads is deduced by matching their induced pulse height functions. Unfortunately, this does not tie the gain relative to the central pad, leaving that parameter to be determined. That analysis does not support the hypothesis that pad 7 has a smaller gain than pad 2 and 8. The calibration data analysis, however, does support the hypothesis that pad 7 has a smaller gain than pads 2 and 8.

A map of all reconstructed events appears in the figure linked here. The gaps at 30 degrees and 90 degrees are due to a problem with the process of inverting the mapping between cloud centroid and pad charge fractions.

Position analysis from induced pulses

To determine the coordinate from induced pulses, the same approach is used as in the previous analyses. The amplitude of the induced pulse is assumed to be proportional to the total charge of the event and a function of the distance from the cloud centroid to the pad centre. An algorithm combines the radial information from all pads that have an induced pulse to determine the x-ray location.

The ratio of the peak amplitude of the induced pulse to the total charge of the event as a function of distance to the pad centre is shown in the figure linked here. Some of the scatter seen for the pads (for example pad 7) is reduced if the first 9 runs and the last 6 runs are ignored. The response functions are overlayed for all the pads in the figure linked here. In order to match the curves, the signals in pads 2, 7, and 8 were scaled by factors 0.93, 0.95, and 0.95 respectively. The curve is the result of a fit to a 4th order polynomial that is used in the analysis to describe the response function.

The efficiency for events to have the x-ray position successfully determined with the induced pulse method is quite high, as summarized in the figure linked here. The bias and standard deviations of the measurements are shown in the figure linked here. The region centered on pad 1, but near pad 7 has a bias in the +x direction, whereas in the direct charge analysis a bias was seen in the opposite direction.

A map of all reconstructed events from induced pulses appears in the figure linked here. An unusual artifact appears near y=1. The grid nature of the data taking is evident.
 

Calibration Data - Part 1

On December 14, data was taken with the collimator centred on the outer pads to study the relative gain of the various pads. Unfortunately, the data was taken with the scope at 125MHz sampling (instead of 250 MHz sampling that was used for the data recorded in November). The first run (1401) was recorded with x-ray source at 6.8 kV by mistake. The rest had the x-ray source at 6.0 kV. A summary of the runs are shown in the table below.
 
 
run
centre pad
events
R1
R5
< signal>
1401
1
624
0.495
0.509
0.762 +/- 0.009
1402
3
755
0.507
0.514
0.688 +/- 0.008
1403
7
763
0.504
0.511
0.665 +/- 0.008
1404
6
797
0.503
0.510
0.595 +/- 0.007
1405
4
755
0.499
0.506
0.674 +/- 0.009
1406
2
749
0.496
0.502
0.734 +/- 0.007
1407
8
734
0.503
0.509
0.703 +/- 0.008
1408
1
51
 
 
0.654 +/- 0.029

The "late amplitude" is defined as 1200 ns after the peak pulse. The values R1 and R5 are the ratios of the late amplitude to the peak amplitude for channels 1 and 5 (both connected to the centre pad). The values are within about a percent of each other. The average signal (which includes the corrections R1,R5) shows much larger variation, over 10%. The variation could arise from the following sources: (1) changes in operating conditions - gas mixture and charging of GEM; (2) variation in the pre-amplifier gain from channel to channel; and (3) variation in the GEM amplification as a function of position over the GEM.

The mean signal values in the table above were used to define scale factors for the direct charge signals from the November data, and the position analysis from charge sharing was repeated. The relative scale factor for pad 1 could not be determined from the calibration data. A few values were considered, and a value corresponding to mean signal size of 0.71 was selected, as it gave the smallest variance of the position residuals. The average cloud size was also varied to reduce the variance of the position residuals, and the best value was found to be 0.55 mm, instead of the value 0.58 mm that was used above. An offset of (-6,26) microns (x,y) is also applied to bring the central values of the residuals to zero.

The calibrations described above reduce the scatter in the biases seen across the GEM pad, as shown in the figure linked here. The distributions of residuals for the entire data set is shown in the figure linked here, with the standard deviations in x and y found to be about 60 and 70 microns respectively. A more careful calibration procedure could possibly reduce the biases further. It should be kept in mind that after moving the micrometer stage and returning to its origin there are shifts of order 20 microns observed.
 

Calibration Data - Part 2

On January 5, further calibration data was taken. The detector had been opened since the previous calibration run, to fix an internal breakdown. Once again the pinhole was centred on the different pads:
 
 
 
run
centre pad
sampling (MHz)
events
R1
R5
< signal>
pad1 correction
501
1
125
393
0.497
0.506
0.731 +/- 0.010
502
1
125
374
0.498
0.504
0.752 +/- 0.010
0.752
503
3
250
306
0.688
0.696
0.774 +/- 0.012
0.757 0.978
504
7
250
358
0.680
0.686
0.807 +/- 0.010
0.761 0.942
505
6
250
308
0.681
0.687
0.736 +/- 0.009
0.766 1.041
506
4
250
351
0.679
0.687
0.777 +/- 0.011
0.770 0.991
507
2
250
213
0.692
0.698
0.786 +/- 0.013
0.775 0.986
508
8
250
243
0.682
0.689
0.800 +/- 0.011
0.779 0.974
509
1
250
223
0.678
0.686
0.784 +/- 0.012
0.784

The second and last calibration run shows an increase in gain by 4%. The assumed pad 1 gain is taken to increase linearly during the calibration runs, with the value shown in the column labelled "pad1".  Pad 6 is again seen to have the lowest gain. However, the two sets of calibrations show greater variance in the relative mean signals than would be expected from the errors in the means. Comparing the pads that were measured for both sets of runs:
 
 

centre pad
< sig>/av <sig> calib 1
< sig>/av <sig> calib 2
pull 1 vs 2 (sigma)
3
1.017 +/- 0.012
0.992 +/- 0.015
+1.3
7
0.983 +/- 0.012
1.035 +/- 0.013
-2.9
6
0.880 +/- 0.010
0.944 +/- 0.012
-4.1
4
0.996 +/- 0.013
0.996 +/- 0.014
  0.0
2
1.085 +/- 0.010
1.008 +/- 0.017
+3.9
8
1.039 +/- 0.012
1.026 +/- 0.014
+0.7

The gain variation from pad to pad is not reproducible between the two sets of calibration runs. Further studies of the gain stability of the gem would be useful. By taking the second set of pad to pad gain corrections, the overall residual deviations from charge sharing (61,73) microns (x,y) is achieved when the cloud size is reduced to 0.52 mm. The width of the residuals is not quite as good as was found for the previous set of calibration constants. An offset of (25,-33) microns is applied to bring the residuals to a mean of zero.

The induced pulse analysis is repeated using these calibrations as the starting point. No improvement in residual widths is observed.

Conclusion

The data taken in November 20-24, 2000 shows promise. Biases seen in the measurements are reduced by a calibration procedure. The calibration procedure itself is not well reproducible.
© 2006 Carleton University 1125 Colonel By Drive, Ottawa, Ontario, K1S 5B6 Canada (613) 520-7400
| Contacts |
Canada's Capital University