cern.ch

Performance of the LHCb Vertex Locator

[to restricted-access page]

Abstract

The Vertex Locator (VELO) is a silicon microstrip detector that surrounds the proton-proton interaction region in the LHCb experiment. The performance of the detector during the first years of its physics operation is reviewed. The system is operated in vacuum, uses a bi-phase CO2 cooling system, and the sensors are moved to 7 mm from the LHC beam for physics data taking. The performance and stability of these characteristic features of the detector are described, and details of the material budget are given. The calibration of the timing and the data processing algorithms that are implemented in FPGAs are described. The system performance is fully characterised. The sensors have a signal to noise ratio of approximately 20 and a best hit resolution of 4 microns is achieved at the optimal track angle. The typical detector occupancy for minimum bias events in standard operating conditions in 2011 is around 0.5, and the detector has less than 1 of faulty strips. The proximity of the detector to the beam means that the inner regions of the n+-on-n sensors have undergone space-charge sign inversion due to radiation damage. The VELO performance parameters that drive the experiment's physics sensitivity are also given. The track finding efficiency of the VELO is typically above 98 and the modules have been aligned to a precision of 1 micron for translations in the plane transverse to the beam. A primary vertex resolution of 13 microns in the transverse plane and 71 microns along the beam axis is achieved for vertices with 25 tracks. An impact parameter resolution of less than 35 microns is achieved for particles with transverse momentum greater than 1 GeV/c.

Figures and captions

(top left) The LHCb VELO vacuum tank. The cut-away view allows the VELO sensors, hybrids and module support on the left-hand side to be seen. (top right) A photograph of one side of the VELO during assembly showing the silicon sensors and readout hybrids. (bottom) Cross-section in the $xz$ plane at $y=0$ of the sensors and a view of the sensors in the $xy$ plane. The detector is shown in its closed position. $R$ ($\Phi$ ) sensors are shown with solid blue (dashed red) lines. The modules at positive (negative) $x$ are known as the left or A-side (right or C-side).

VeloFu[..].pdf [424 KiB]
HiDef png [5 MiB]
Thumbnail [1 MiB]
VeloFull_Cut.pdf
VeloAs[..].pdf [248 KiB]
HiDef png [6 MiB]
Thumbnail [1 MiB]
VeloAssemblyPhoto.pdf
VELO_d[..].png [60 KiB]
HiDef png [61 KiB]
Thumbnail [23 KiB]
VELO_detector_layout_crop-reduce.png

Schematic representation of an $R$ and a $\Phi$ sensor. The $R$ sensor strips are arranged into four approximately $45 ^{\circ} $ segments and have routing lines perpendicular to the strips. The $\Phi$ sensor has two zones with inner and outer strips. The routing lines of the inner strips are orientated parallel to the outer strips.

randph[..].pdf [165 KiB]
HiDef png [3 MiB]
Thumbnail [1 MiB]
randphisensors.pdf

Illustration of the influence of the proton beams on the beam volume pressure. The readings of the beam volume (PE 412) and detector volume (PE 422) are shown, and the increase in pressure in the beam volume coinciding with the beam injections (each containing $2\times10^{14}$ protons) at $t = 6 \mathrm{hours}$ can be seen. The beams were dumped 24 hours later.

Velova[..].pdf [16 KiB]
HiDef png [147 KiB]
Thumbnail [117 KiB]
Velovacuum130130.pdf

The horizontal position of LHC collision vertices as reconstructed by the A-Side (left) and C-Side (right). These online screenshots of data illustrate the "stop-measure-move" cycle of the closing procedure, where each peak in the distribution corresponds to a "stop", as well as the degradation of vertex resolution with opening distance.

Malcol[..].png [47 KiB]
HiDef png [111 KiB]
Thumbnail [59 KiB]
Malcolm_Closing_ASide.png
Malcol[..].png [50 KiB]
HiDef png [116 KiB]
Thumbnail [61 KiB]
Malcolm_Closing_CSide.png

(left) Material passed through on trajectories originating from the interaction point and terminating at the radius of the first active strip in the detector expressed in terms of the fraction of a radiation length ( $ X_0$ ) and as a function of pseudorapidity ($\eta$) and azimuthal angle ($\phi$). The material is higher in the regions around $\pm90 ^{\circ} $ where the two halves of the detector overlap. (right) Breakdown of the total material budget by component of the VELO . The number given for each component is the percentage of the total average VELO material budget ($0.227 X_0 $).

radlen[..].pdf [165 KiB]
HiDef png [1 MiB]
Thumbnail [659 KiB]
radlength-fmp-simple-colour.pdf
velo_p[..].pdf [14 KiB]
HiDef png [346 KiB]
Thumbnail [216 KiB]
velo_pie_blacktext.pdf

Vertices of hadronic interactions in the LHCb VELO material. (top left) Full VELO system with entrance and exit regions visible in data. (top right) Zoom in to a group of sensors downstream from the interaction point in data. (bottom left) the same region reconstructed using simulation events. (bottom right) Zoom onto an pile-up module consisting of a single R-sensor to check the distance between the sensor and foil.

scan-f[..].png [65 KiB]
HiDef png [130 KiB]
Thumbnail [67 KiB]
scan-full-data.png
scan-z[..].png [70 KiB]
HiDef png [138 KiB]
Thumbnail [63 KiB]
scan-zoom-data.png
scan-z[..].png [52 KiB]
HiDef png [93 KiB]
Thumbnail [46 KiB]
scan-zoom-MC.png
PU02_C.png [48 KiB]
HiDef png [90 KiB]
Thumbnail [43 KiB]
PU02_C.png

An example of the scan over the digitisation phases (see text). Data for 16 analogue links are shown. The horizontal axis represents the sampling time, where the integer numbers indicate the location of the serially transmitted Beetle channels, and between two channels there are 16 bins. The solid red vertical line indicates the chosen sampling point and the dashed red lines the position of the previous and next samples.

delays[..].png [15 KiB]
HiDef png [117 KiB]
Thumbnail [72 KiB]
delayscan-edit-reduce.png

The FHS distribution before and after the first gain calibration. The data were taken two months apart in late 2009. One entry is made in the histogram for each link.

GainCa[..].pdf [37 KiB]
HiDef png [235 KiB]
Thumbnail [76 KiB]
GainCalibration_BeforeAfter-edit.pdf

Pulse-shape obtained by combining data from all sensors, with a 1 $ {\rm ns}$ time step and after applying for each sensor the time shift determined by the procedure. The plot on the left shows the ADC spectra versus time. The plot on the right shows the most probable value for each time bin obtained from a fit of a Landau convolved with a Gaussian function. A bifurcated Gaussian function is fitted to this distribution from which the rise and the fall widths are obtained. The fitted function also predicts the amount of pre-spill (25 $ {\rm ns}$ before the chosen point) and the spillover (25 ns after), which are marked with red vertical dotted lines. The chosen sampling time is also marked with a red vertical dotted line.

pulses[..].pdf [100 KiB]
HiDef png [695 KiB]
Thumbnail [293 KiB]
pulseshape-edit.pdf

Typical raw (Non-Zero Suppressed) data ADC values. Data from two ASICs, each with four analogue links, is shown. Variations in the pedestal values at the ASIC, analogue link, and channel level are seen.

rawDat[..].pdf [83 KiB]
HiDef png [624 KiB]
Thumbnail [280 KiB]
rawData-edit.pdf

Distribution of common mode noise per link. The blue histogram is for all links. The green histogram is for links in a ASIC that has large signals (see text).

common[..].pdf [39 KiB]
HiDef png [153 KiB]
Thumbnail [146 KiB]
commonModeDist-edit.pdf

(left) Event snapshot for an ASIC that has large signals (channels 1024-1152). The blue dots are the raw ADC values. The red lines are the pedestals. (right) ASIC baseline shift in response to injected charge in units of the most probable charge deposited by a MIP in 300 $ \upmu\rm m$ thick silicon.

largeE[..].pdf [45 KiB]
HiDef png [278 KiB]
Thumbnail [247 KiB]
largeEnergyDepEvt-edit.pdf
cmShif[..].pdf [71 KiB]
HiDef png [170 KiB]
Thumbnail [128 KiB]
cmShiftTestPulse-edit.pdf

Example of the graphical user interface used in the VELO monitoring. Shown are distributions of the number of hits per VELO track segment (top left), the track angle with respect to the beam axis (top right), and the ADC counts of $R$ and $\Phi$ sensor clusters associated to a track (bottom left and right, respectively). The results of the current run are shown in red filled symbols and the reference data for comparison are shown as black open symbols.

VeloMo[..].pdf [100 KiB]
HiDef png [862 KiB]
Thumbnail [341 KiB]
VeloMoniGUI.pdf

(left) Fit of the signal distribution for Sensor 104 (a $\Phi$ sensor) for clusters on long tracks with a track intercept point on the sensor $10 < r (\rm mm ) <15$. (right) The most probable value of the fits to sensors 40 and 104, the $R$ and $\Phi$ sensor of the same module respectively, as a function of the track intercept point radius. The cluster ADC counts were normalised to a track crossing 300 $ \upmu\rm m$ of silicon.

Sensor[..].pdf [59 KiB]
HiDef png [175 KiB]
Thumbnail [147 KiB]
Sensor104R10-15_ADCFit-edit.pdf
data_s[..].pdf [60 KiB]
HiDef png [103 KiB]
Thumbnail [92 KiB]
data_sens_MPV_40-edit.pdf

Noise in ADC counts averaged across the 42 installed R (left) and $\Phi$ (right) sensors, with the error bars indicating the RMS of the distribution.

RNoise[..].pdf [107 KiB]
HiDef png [135 KiB]
Thumbnail [67 KiB]
RNoise_StripAve-edit.pdf
PhiNoi[..].pdf [421 KiB]
HiDef png [422 KiB]
Thumbnail [211 KiB]
PhiNoise_StripAve-edit.pdf

Signal to noise (S/N) ratio from the MPV of the signal for single strip clusters on tracks divided by the noise of that strip. Shown are the S/N values for sensor 40 (R) and sensor 104, the $\Phi$ sensor of the same module, as a function of impact point radius.

dataSN[..].pdf [57 KiB]
HiDef png [118 KiB]
Thumbnail [105 KiB]
dataSN_sens_40-edit.pdf

(left) The VELO resolution for two projected angle bins for the $R$ sensors as a function of the readout pitch compared with binary resolution. (right) Resolution divided by pitch as function of the track projected angle for four different strip pitches.

VeloHi[..].pdf [42 KiB]
HiDef png [209 KiB]
Thumbnail [192 KiB]
VeloHitReso2010_withbinary-edit.pdf
ResoRs[..].pdf [38 KiB]
HiDef png [232 KiB]
Thumbnail [217 KiB]
ResoRsensor2010_projectedangle-edit.pdf

The percentage of one (left) and two (right) strip clusters as a function of the track projected angle for four different strip pitches.

percen[..].pdf [38 KiB]
HiDef png [228 KiB]
Thumbnail [207 KiB]
percentageClsize1_projectedangle-edit.pdf
percen[..].pdf [38 KiB]
HiDef png [229 KiB]
Thumbnail [210 KiB]
percentageClsize2_projectedangle-edit.pdf

Cluster occupancy in the VELO silicon strip detectors. (left plots) Average occupancy in the sensor as a function of the position of the sensor along the beam-line. (right plots) Occupancy as a function of the local radius of strips on the $R$ sensors, with a negative sign applied to sensors on one half. The upper plots show the occupancy when fully closed using 2011 data with a $\mu$ of 1.7 and for events selected using a random trigger or with events passing the high level trigger. The lower plots show the occupancy as a function of closing distance with fully closed labelled as 0 $\rm mm$ , the points below then follow in order of the retraction distance indicated in the key.

occ_vs[..].pdf [17 KiB]
HiDef png [109 KiB]
Thumbnail [81 KiB]
occ_vs_zr_calib_hlt.pdf
occ_vs[..].pdf [24 KiB]
HiDef png [137 KiB]
Thumbnail [95 KiB]
occ_vs_r_calib_hlt.pdf
occ_vs[..].pdf [24 KiB]
HiDef png [193 KiB]
Thumbnail [133 KiB]
occ_vs_zr_closing.pdf
occ_vs[..].pdf [27 KiB]
HiDef png [266 KiB]
Thumbnail [161 KiB]
occ_vs_r_closing.pdf

(left) Imaging of the LHC beams through the reconstruction of the production vertex of tracks in 2010. The vertices from beam-gas interactions with the beam travelling in the $+z$ ($-z$) direction in LHCb are shown in blue (red). Vertices arising from beam-beam interactions are in green. The vertical axis represents the horizontal direction ($x$) and the horizontal axis the beam direction ($z$). (right) An event display image showing clusters in $R$ and $\Phi$ sensors with a high occupancy in a localised region assumed to be due to a beam background splash event.

BeamGas.pdf [153 KiB]
HiDef png [341 KiB]
Thumbnail [147 KiB]
BeamGas.pdf

The cluster finding efficiency for each sensor with all strips (left) and when identified bad strips are excluded (right).

CFE_Al[..].pdf [16 KiB]
HiDef png [248 KiB]
Thumbnail [235 KiB]
CFE_AllStrips_0pb.pdf
CFE_BS[..].pdf [17 KiB]
HiDef png [266 KiB]
Thumbnail [256 KiB]
CFE_BSR_0pb.pdf

The intercept position of an extrapolated track when a cluster is not found for two example $R$ (left) and $\Phi$ (right) sensors. The bad channels are clearly visible.

BadStr[..].pdf [21 KiB]
HiDef png [368 KiB]
Thumbnail [284 KiB]
BadStrips_S36_2D.pdf
BadStr[..].pdf [24 KiB]
HiDef png [454 KiB]
Thumbnail [339 KiB]
BadStrips_S100_2D.pdf

Currents measured for each sensor as a function of time (bottom). The integrated luminosity delivered to LHCb and the average sensor temperature is shown over the same time scale (middle and top). Increases in the delivered luminosity are matched by increases in the sensor currents. The evolution of the mean measured current agrees well with the prediction from simulation. The mean measured value excludes sensors that are surface-current dominated.

Leakag[..].pdf [171 KiB]
Thumbnail [279 KiB]
LeakageCurrentWithData.pdf

The effective depletion voltage versus fluence for all VELO sensors up to 3.4 $ fb^{-1}$ delivered integrated luminosity.

Fluenc[..].pdf [66 KiB]
HiDef png [117 KiB]
Thumbnail [60 KiB]
Fluence_2013_bw-edit.pdf

Tracking efficiency for the 2011 data and simulation for the VELO as a function of the momentum, p (top left), the pseudorapidity, $\eta$ (top right), the azimuthal angle $\phi$ (bottom left) and the total number of tracks in the event, $N_{track}$ (bottom right). The simulation has been reweighted to the number of tracks observed in data for the p,$\eta$ and $\phi$ plots. The error bars indicate the statistical uncertainty.

Track_[..].pdf [13 KiB]
HiDef png [102 KiB]
Thumbnail [83 KiB]
Track_Velo_effVsP2011.pdf
Track_[..].pdf [13 KiB]
HiDef png [81 KiB]
Thumbnail [67 KiB]
Track_Velo_effVsEta2011.pdf
Track_[..].pdf [14 KiB]
HiDef png [97 KiB]
Thumbnail [79 KiB]
Track_Velo_effVsmyPhi2011.pdf
Track_[..].pdf [13 KiB]
HiDef png [95 KiB]
Thumbnail [77 KiB]
Track_Velo_effVstrackmult2011.pdf

(left) Fraction of ghost tracks versus number of VELO clusters in simulation. (right) VELO pattern recognition timing versus number of clusters, the pattern recognition alone (magenta triangles), the raw data decoding (blue squares) and the combined (red circles) timings are shown. The times are scaled to a 2.8 $ {\rm GHz}$ Xeon processor.

fracGh[..].pdf [8 KiB]
HiDef png [116 KiB]
Thumbnail [98 KiB]
fracGhostAccept-edit.pdf
time_clus.pdf [11 KiB]
HiDef png [195 KiB]
Thumbnail [125 KiB]
time_clus.pdf

Example unbiased sensor residuals as a function of the $\phi$ coordinate using only the survey information (left) and using the track-based software alignment (right). Results are given for two different example sensors. (top) A significant improvement in the residuals is seen in this sensor with the track-based alignment. (bottom) In this sensor the alignment quality using the survey information is already good.

Plotco[..].png [63 KiB]
HiDef png [249 KiB]
Thumbnail [102 KiB]
PlotcompareMeanResMetrology_TrackAlign_v30_col-reduce.png

Misalignment between the two VELO halves in each run, evaluated by fitting the PV separately with tracks in the two halves of the VELO . The run numbers shown here span the period of the last four months of operations in 2010.

PVXmis[..].pdf [71 KiB]
HiDef png [102 KiB]
Thumbnail [49 KiB]
PVXmisalignVsRun-edit.pdf

PV Resolution of events with exactly one PV in 2011 data as a function of track multiplicity. (left) $x$ (red) and $y$ (blue) resolution and (right) $z$ resolution. The fit parameters $A$, $B$ and $C$ for each coordinate are given.

ResXY_[..].pdf [68 KiB]
HiDef png [205 KiB]
Thumbnail [156 KiB]
ResXY_1PV_2011Data.pdf
ResZ_1[..].pdf [64 KiB]
HiDef png [99 KiB]
Thumbnail [57 KiB]
ResZ_1PV_2011Data.pdf

$IP_x$ and $IP_y$ resolution as a function of momentum (left) and $IP_x$ as a function of $1/p_T$ and compared with simulation (right). Determined with 2012 data.

IPRes-[..].pdf [21 KiB]
HiDef png [177 KiB]
Thumbnail [159 KiB]
IPRes-Vs-P-CompareIPxIPy-2012.pdf
IPXRes[..].pdf [17 KiB]
HiDef png [174 KiB]
Thumbnail [162 KiB]
IPXRes-Vs-InversePT-Compare2012DataToMC.pdf

$IP_x$ resolution as a function of azimuthal angle $\phi$, measured on 2011 data and compared to simulation.

IPXRes[..].pdf [16 KiB]
HiDef png [138 KiB]
Thumbnail [148 KiB]
IPXRes-Vs-Phi-Compare2011DataToMC.pdf

Decay time (left) and decay length (right) distribution for fake, prompt $ B ^0_s \rightarrow { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu} \phi\rightarrow \mu ^+\mu ^- K ^+ K ^- $ candidates in 2011 data (black points) and simulation (solid red histogram). Only events with a single PV are used. The simulated data are generated inclusive $ { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu} $ events from which signal $ B ^0_s \rightarrow { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu} \phi$ are removed. In the data contributions from non- $ { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu}$ {} di-muon combinations and from true $ B ^0_s \rightarrow { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu} \phi$ are subtracted using the sPlot technique \cite{Pivk:2004ty}.

decayt[..].pdf [63 KiB]
HiDef png [135 KiB]
Thumbnail [118 KiB]
decaytimepromptJPsiPhi-edit.pdf
decayl[..].pdf [47 KiB]
HiDef png [142 KiB]
Thumbnail [127 KiB]
decaylengthpromptJPsiPhi-edit.pdf

Decay time resolution (points) as a function of momentum (left) and as a function of the estimated decay time uncertainty (right) of fake, prompt $ B ^0_s \rightarrow { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu} \phi\rightarrow \mu ^+\mu ^- K ^+ K ^- $ candidates in data. Only events with a single PV are used. Contributions from non- $ { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu}$ {} di-muon combinations and from true $ B ^0_s \rightarrow { J \mskip -3mu/\mskip -2mu\psi \mskip 2mu} \phi$ are subtracted using the sPlot\ technique. The superimposed histogram shows the distribution of momentum (left) and estimated decay time uncertainty (right) on an arbitrary scale.

decayt[..].pdf [15 KiB]
HiDef png [133 KiB]
Thumbnail [124 KiB]
decaytimeresoVsMomentumJPsiPhi.pdf
decayt[..].pdf [15 KiB]
HiDef png [124 KiB]
Thumbnail [113 KiB]
decaytimeresoVsErrorJPsiPhi.pdf

Animated gif made out of all figures.

DP-2014-001.gif
Thumbnail
thumbnail_DP-2014-001.gif

Tables and captions

Fraction of faulty strips classified as dead or noisy. The results are obtained from the methods used at production, occupancy spectrum studies, and a cluster finding efficiency analysis. Values are given at the time of production or start of operations in 2010, at the end of 2011, and at the end of 2012.

Table_1.pdf [32 KiB]
HiDef png [63 KiB]
Thumbnail [31 KiB]
tex code
Table_1.pdf

Created on 18 October 2023.