Autocalibration of high precision drift tubes

UCLEAR PHYSIC~

ELSEVIER

Nuclear Physics B (Proc. Suppl.) 54B (1997) 311-316

PROCEEDINGS SUPPLEMENTS

Autocalibration of high precision drift tubes C.Bacci, C.Bini, G.Ciapetti, G.De Zorzi, P.Gauzzi, F.Lacava, A.Nisati, L.Pontecorvo*, S.Rosati and S.Veneziano Universita' di R o m a "La Sapienza" and I N F N R o m a M.Cambiaghi, G.Casellotti, C.Conta, M.Fraternali, A.Lanza, M.Livan, G.Polesello, A.Rimoldi and V.Vercesi Universita' di Pavia and INFN Pavia We present the results on MDT (monitored drift tubes) autocalibration studies obtained from the analisys of the data collected in Summer 1995 on the H8B Muon Test Beam. In particular we studied the possibility of autocalibration of the MDT using four or three layers of tubes, and we compared the calibration obtained using a precise external tracker with the output of the autocalibration procedure. Results show the feasibility of autocalibration with four and three tubes and the good accuracy of the autocalibration procedure.

1. I n t r o d u c t i o n In the framework of the Atlas muon activities we have performed in the H8B test b e a m at CERN, measurements using an M D T [1] multilayer of four layers, operated with different gas mixtures at three Bar absolute pressure. The main aim of the measurements was to show the feasibility and the accuracy of the autocalibration [2] procedure, and the possibility of autocalibration with only three layers of tubes per multilayer; both problems are of relevance for the design of the Atlas Muon Spectrometer. In the following we present the experimental set up, the autocalibration method, the result of autocalibration obtained using both three and four layers of tubes and the comparison between the M D T Radius-Time (R-T) relation obtained with autocalibration and the one measured with the external tracker. 2. E x p e r i m e n t a l

set-up

In the H8 b e a m line at CERN, we have set up a muon test facility equipped with a trigger hodoscope made out of four scintillators (2 x 2, 4 × 4, 10 x 10, 10 × 10 cm ~ ), two b e a m chambers and a reference system (external tracker) made out of two matrices of four by four drift tubes, *Presented by L.Pontecorvo INFN Roma, c/o Dip. Fisica Univ. Roma "La Sapienza" P.le A. Moro 4 00100 Rome (Italy), fax# +39 6 4957697, e-mail [email protected] .infn.it 0920-5632/97/$17.00 © 1997 Elsevier Science B.V All rights reserved. Pll: S0920-5632(97)00128-X

operated in streamer mode. The muon b e a m mom e n t u m varied between 100 and 180 GeV/c. The d a t a aquisition was based on Spider under the OS9 operating system. We have performed measurements on an M D T made out of four layers containing fifteen drift tubes each. The chamber was made of single alluminum tubes 400 p m thick, 3.010 cm diameter and 80 cm in lenght. We used 50 p m gold plated W-Re wires. We used the following gas mixtures : Ar 94.6 % C2H6 0.4 % COs 5% at a pressure of 3 Bar absolute. The tubes were operated in proportional mode at 3.1 KV. The front end electronics was based on the VTX chip [3], the amplification was 1 m V / f C and the peaking time 8 ns. After the preamplification board there was a voltage amplifier and discriminator card in wieh the signal was amplified by a factor of ten and then discriminated. The discrimination threshold was set quite high (300 mV) due to noise and cross talk problems between channels, corresponding, at a tipicM gas gain of 4 x 104, to about 100 electrons. Finally the discriminated ECL signals were fed in a LeCroy 2277 T D C with a LSB of 1 ns and a range up to 64 ps. Some signals were split before the discriminator and were sent to a LRS 2249W ADC to measure the charge spectra and the streamer fraction.

312

C. Bacci et al./Nuclear Physics B (Proc. Suppl.) 54B (1997) 311 316

3. A u t o c a l i b r a t i o n m e t h o d The aim of autocalibration is to obtain the RT relation, the position of the wires and the zero time (T0s) from the data of the tubes themselves without the help of external detectors. To obtain this we used an iterative procedure in which we first get the 0-approximation of the R-T relation using the well known technique of the integration of the time spectrum, assuming a uniform illumination of the tubes [4]; then using it we fit the tracks in the detector and we measure the systematic error of the current R-T relation from the value of the residuals vs the drift time in the tubes. The 1-approximation of the R-T relation is then obtained correcting the 0-approximation with the value of the systematic error measured at each drift time. This procedure is iterated until the difference between the n-approximation and the n-1 approximation of the R-T relation is negligible. The T0s of each wires is set at the beginning of the procedure to the value of the rise of the time spectrum, known with an accuracy of few ns. The first approximation position of the wires is given by their nominal geometrical position; a better determination of the T0s and of the wire positions is obtained during the iterative procedure studing the residuals of the track fit separately for tracks crossing the tubes on the upper and the lower side of the wires. An error on the TO is reflected in a separation of the up and down residuals by the quantity:

the time spectrum and assuming that this correspond to the physical dimension of the detector. This approach is not very precise because of several reasons: 1) it is always quite difficult to precisely determine the end point of the time spectrum; 2) the tube geometrical dimensions could be different from the dimension where the tube is really efficient and 3) the time spectrum could be distorted close to the end wall of the detector by the presence of ~ electrons. To overcame this difficulties we have used a different approach. From our geometry and from the fact that the tracks are almost perpendicular to the detector it is easy to see that tracks crossing the first tube at R/2, where R is the geometrical tube radius, will cross all the other tubes almost in the same position. In this condition it is possible to select all the tracks passing at the center of the tubes asking that the drift time of the tubes being nearly equal (+ 5 ns). In this way we obtain a gaussian distribution of drift time correspondig to a drift distance equal to R/2 Fig (1), which sets the absolute scale with

R/2 DistanceScale Setting

tD

20

2O

I@

10 170

130

140

150

16~

17,o 13o

14o

15o

16o

120

130

140

150

160

12o

l~to

14o

15o

16o

A T O - ~ (ResUP -- Resd°wn)

(2 ×

vdr /t)

(1) tD

where V dri/t is the average drift velocity in the tube. After this correction the accuracy on the T0s is about 300 ps. The same technique is used to determine the error on the position of the wires with an accuracy of about 10 micron (assuming no systematic error introduced by the fitting procedure): A Y O = (Res~P + Resd°wn)

2

1o 5

5

n$

Figure 1. Time distribution corresponding to R/2 drift distance.

(2)

One of the problem in the autocalibration procedure is the setting of the absolute scale of distances. Usually this is done using the end point of

high precision. Another common problem of the autocalibration method is due to the least square fitting procedure in presence of systematic errors coming from the R-T relation. It is possible to

313

C. Bacci et al./Nuclear Physics B (Proc. Suppl.) 54B (1997) 311 316

demonstrate that the reconstruction of a track, obtained using four measured points all affected by systematic errors of the same magnitude, results in a wrong reconstructed direction of the track while the intercept of the track in the mid plane of the detector is not affected. To correct for this we have used different weights in the fit for the different measured points, namelly the two outer layers' weight is set to 1 while the internal layers' weight is set to 3. In the case of three tube autocalibration the situation is different, in fact the direction of the reconstructed track is not affected by the systematic errors on the measured point while the intercept at the central plane of the detector is wrong. In this case the weights used are 1 for the external layers and 2 for the internal one. This weighting technique is used only in the autocalibration procedure and not during the track fitting after calibration. 4.

Result

on

Residuols vs Drift Time LI

..O.~S .0.1 0

1O0

200

300

0.1

O75 O.OS 0.025 .,0.02S .,0.8~ .4.0?5 -0.1

-4L! R$

Figure 2. 0-iteration residuals.

time is flat and is centered to zero. Moreover if we study the variable

external trackers

The external tracker is made of two identical detectors each of which consist in a matrix of 4 × 4,500 p m thick alluminum drift tubes, 3.030 cm diameter, 30 cm long, wired with a 100 p m alluminum wire. The gas mixture used was 600£ I-Butane and 40% Argon at normal pressure, the high voltage was set to 5200 V, the detector was operated in streamer mode. The readout chain for each channel consisted simply in a discriminator whose threshold was set to 30 mV, and a LRS 2277 TDC. We have autocalibrated the two trackers using the procedure previously explained on groups of four aligned tubes (quadruplet). Using the 0-iteration R-T relation we measure the residuals Fig (2) with respect to the fitted tracks for each tube defined as: R e s -- y / i t _ y m e a ,

0.1 ~ 0.075 0.0~ O.iJ~ 0 .,,O.4J~S

(3)

From this figure we see that the systematic error on the 0-iteration R-T relation is as high as 500 pm, and is rapidly varing with the drift time; at the distance of R/2, the error is zero by definition because that point was used to set the space absolute scale. In Fig (3) we present the situation of the residuals of the four tubes after the first iteration. The residual distribution vs the drift

R = (y~e~, + y~ea,)/2 + y~eo,

(4)

where Y n 'nea" indicates the measured distance in tube n, and R, in our geometrical situation should be equal for each event, to the distance between the wires of two staggered planes. We see that the mean value of this distribution 1.514 cm is in very good agreement with what expected from the tracker mechanical construction, proving the good accuracy of the absolute distance scale. We have repeated the procedure for other two iterations and the measured R-T relation shows no relevant difference with the previous ones, indicating the very fast convergence of the autocalibration procedure. We have applied the autocalibration procedure to different quadruplets in the first and the second tracker obtaining always consistent results, so that we can use a single RT relation for each detector (there is a small systematic difference between the R-T relation of the two trackers). We have also studied the variations with time of the R-T relation finding it very stable. After the calibration procedure we have fit tracks in the two detectors separately, using an average single point resolution of 65/~m, the resuiting distribution of the X2 probability Fig (4) is flat showing the good space accuracy of the

314

C. Bacci et al./Nuclear Physics B (Proc. Suppl.) 54B (1997) 311-316 Residuols vs Drift Time 1 Iterotion

0.1

~ O.075

•";';2 ;':: - " . " " : " •

Chi Square Probability Tracker 1

.0"I$

~:;:~...:....:- .

O.OS O.O75

,.>, 40O

e, -O.O25

350

-O.OS ..O.O~S

3~

-0.1 O

lOO

200

300

0

100

200

300 ~e

n~,

0.1

0.1 lee

O.OS o.ff25

.i.,

:i.... :~:.

~

,':~.:-',~':i~:.

IS0

e 100

-0.025 -O.05 -O.~/5 -0.1 (I

~.:~:-*'.i'/~.;.~'~:U *~l I~ 200 300

I *.1

(';"i', ~ "'"'.1"::'" ", I 1~ :200 300

0 O.I

0.2

0.3

0,4

05

06

07

0.8

09

n$

Figure 3. 1-iteration residuals.

single tubes. After aligning the two detectors we have fit tracks using both trackers. The two detectors were set at a distance of about 9 metres on the beam line to assure a good determination of the direction of the beam tracks. Unfortunatelly, due to the presence on the beam line of material (radioactive source shielding) introducing a large multiple scattering, the precision of the interpolation of the tracks on the MDT in between the two trackers is not very accurate ( 2 300 #m) because of the big distance between the two trackers.

Figure 4. Distribution of the X 2 probability assuming 6 5 / l m resolution per point.

lute space scale. After the calibration procedure we fit tracks in the MDT, the resulting distribution of the X 2 probability is flat using an average single point resolution of 110 /tm. The resolution is mainly limited by the high threshold used, infact our resulution figure is in good agreement with the resolution calculated by a detailed simulation of the drift tube response [5] in which the threshold was set to 100 electrons.

5. R e s u l t s o n M D T

6. T h r e e vs f o u r t u b e s a u t o c a l i b r a t i o n

A typical time spectrum for an MDT tube is shown in Fig (5). We notice that the rise of the spectrum is rather slow because of the high threshold used. It gives problems in the first determination of the T0s and introduces a dead region close to the wire (~_ 1 mm). The residuals vs the drift time using the 0-iteration R-T relation are shown in Fig (6), also in this case we see the big systematic errors due to the wrong 0-iteration R-T relation. After few iterations the autocalibration procedure converges and the residuals vs drift time, shown in Fig (8), are flat and centered to zero. The distribution of the variable (Y1 m*a~ + Y 3 meas )/2+Y2 m ~ is centered to 1.504 cm showing again the precision of the abso-

We have repeated the calibration procedure on both the external tracker and the MDT using only three contiguous tubes instead of four tubes and we have compared the resulting R-T relations bin by bin. The R-T relations obtained with the two procedures give very similar results on the whole drift distance (about 10 p m difference) except close to the wire and to the end wall of the tube where we notice a peak structure of about 40 p m in the R-T relations difference. We still have to investigate the nature of this systematic difference, but from the practical point of view we think that it is not relevant because it influences a very small part of the tube where anyway the space resolution is quite degraded ( > 100/tm).

315

C. Bacci et al./Nuclear Physics B (Proc. Suppl.) 54B (1997) 311 316

MDT Time Distribution

Residuals vs Drift Time and R-T Relation

0.1

(0 C

~

:';:

0.ff75

::- ":.

0.025 0

~b

-0.025 80

-0.05 -0.075 -O.!

n~ EE 1.4

~' "'?

1.2

40

" "~,~:." .

.

~.~'~

: ~'¢C-'" ': ; :~ .. :~i.~ f : .,, : • "

,',.

. . d : ,::.~.'~ ..... ,,I ; ~ :-..~,,'..~ ~' ..."~?::~"

*'~':'~ i7.:~ :~-::;~,'3~.~'

. .."-

20 0.4

1~

200

3~

4~

~

6~

h 70~

~.::

-600

-400

.200

~~*,~.'

0

7-

200

4~

600 fl$

n$

Figure 5. M D T time spectrum.

Figure 6. 0-iteration residuals vs drift time.

7. C o m p a r i s o n between autocalibration and external tracker calibration

the wire with a constant velocity V dri'ft the mini m u m time detectable will be the time needed by N electrons to reach the wire:

The easiest way to measure the R - T relation for a drift tube is to use a precise tracking detector to determine, event by event, the distance of closest approach of a track to the wire of the tube, and to correlate it with the corresponding drift time. We used our external tracker to define tracks inpinging on the M D T and we measured the R - T relation for the M D T tubes in an independent way with respect to the autocalibration procedure. The distribution of the difference between the distance of closest approach defined by the external tracker and the one measured using the autocalibration is shown in Fig(8). The width of this distribution (2300/~m) is due, as previously stated (cfr 4), to multiple scattering on the source shielding. We see that there is a good agreement between the two R - T relations in the whole tube apart from a region of +1.5 m m from the wire and the corresponding region close to the tube wall. This systematic error is due to the high threshold used for the MDT; this is explained in the following. Consider a track crossing the tube right on the wire, then ideally the measured time and distance from the wire would be zero. Now if we set a threshold of N electrons and a crossing muon produces M electrons/cm wich drift toward

N

T..i. = (2 × M x Vd~il t)

(5)

the factor 2 being due to the fact that the electrons are produced on both sides of the wire. Correspondigly all the tracks passing very close to the wire will be reconstructed at a m i n i m u m distance from the wire given by: N P ~ i , - (2 × M-------~

(6)

In our case, with a threshold of 100 electrons, 300 electrons/cm produced and an average drift velocity V driyt of about 30 # m / n s the m i n i m u m distance from the wire is 1.6 m m and the mini m u m time to collect the 100 electrons is about 50 ns. Those values gives us respectively the size of the zone interested by the sistematic error in the R - T relation obtained with the autocalibration procedure and the slowness of the rise of the time spectrum for the MDT.

8. C o n c l u s i o n s During the s u m m e r 1995 we have set up and succesfully operated a muon test facility in the H8B zone, where we made measurements on one

316

C. Bacci et al./Nuclear Physics B (Proc. Suppl.) 54B (1997) 311-316 Finol Residuole ve Drift Time and R - T Relotion

E

0.1

E

o.o"/$

O

t .;,.>~,;.

:. ~

o.~J

• .

,

.

o.~$

Meosured - Autocalibroted R - T

.: ....:: ,0.1

o.o~q

• -o.o25

I...:~ 14

,.0~$ -o, o 7 5

4

T

/

,

b

"

;.?

: : - : : ,:.:!3

0.05

. .

,

, : i;~,

,.~1 0

I00

~

,t~

400

500

600

n~ 1-2

-0.05

I

~ ~"..' ~i7.'.~,.r:~:~.. ::~:~4:"./"-' ".:: -~c~:;~:,:

::~

oI 0.6

":?

0.4

i~:~(!!:

0.2 0

I

-0.15

'

.500

~

.7,00

0

200

400

-1-$

~00

-I

I O-q

1

l,.,q CrTI

n$

Figure 7. 3-iteration residuals vs drift time.

-0.$

Figure 8. Difference between R-T relation obtained with trackers and autocalibration R-T relation for MDT.

REFERENCES

MDT prototype. The focus of our studies was centered on the autocalibration procedure used to calibrate both the MDT prototype and a set of drift tubes used as a precision external tracker. The results on the tracker shows the very good spacial resolution of this device, 65 # m each single point , and the fast convergence of the autocalibration procedure. The MDT results shows a resolution of about 110 m m and also in this case the convergence of the calibration procedure is very fast. We demonstrated the possibility of autocalibrating the tubes using three layers of tubes showing that there is a very small difference between the R-T relations obtained using three or four layers. Eventually we checked the autocalibration procedure comparing the R-T relation for the MDT measured using the externM tracker with the one obtained by autocalibration. For the whole drift distance, except for a small region close to the wire and the end wall, the two R-T relation are in good agreement, showing the correctness of the autocalibration procedure• A c k n o w l e d g e m e n t s We wish to acknowledge the highly qualified contribution of our technical staff, in particular we would like to thank G.Capradossi, A.Mattei and A.Zullo of the Rome mechanical workshop, and G.Iuvino, A.Freddi and F.Vercellati of the Pavia mechanical workshop.

1 2

3 4 5

Atlas technical proposal, Cern\lhcc\94 -43 lhcc\P2 15 Dec 1994 S.Baranov et al.,Atlas internal note Muon-No040, Z.Y.Feng et all., Atlas internal note MuonNo-056, Z.Y.Feng et al., Atlas internal note Muon-No-095. R.J.Yarema et al., IEEE Trans on Nucl. Science Vol 39 No 4 1992 F.Sauli Cern Report 77-09(1977) V.Tcherniatine private comunication