LHC beam operation of the SLB boards



Collisions

First candidate collisions: here
SLBs in ECAL barrel registered them at bx 52

Beam splashes

Analysis

The particles produced after a beam splashed traverse the CMS detector resulting in an increase of the production of trigger primitives and in consequence in the increase in the occupancy of the SLB accumulators. The excess in the trigger primitive production rate can be searched using a method based on the computation of a significant deviation from background only counts or by computing the correlation with respect to a reference distribution.

Description (show) Description (hide)

Search for deviation from background only counts The instant of the beam splash can be found by summing up all the SLB accumulators and getting the most probable value of the obtained distribution. The accumulators are then divided into 2 regions:

  • A: [0,bx-50[ U [bx+50,1023[ - background region (dominated by cosmics and electronics noise counts)
  • B: [bx-50,bx+50[ - signal region (expected region for the splash)

The average number of counts per bunch recorded by the accumulators in these regions is computed as:

$\begin{equation*} N =\int_{b_{1}}^{b_{2}} \frac{H(b)}{(b_{2}-b_{1})}\end{equation*}$.

The significance of the increase of counts due to beam can thus be estimated from:

$\begin{equation*}S_{occupancy} = \frac{N_{B}-N_{A}}{\sqrt{N_{A}}} = \frac{signal}{\sqrt{noise+background}}\end{equation*}$

One can ask if the maximum observed in each accumulator is compatible with the noise counts expected, as measured in the background dominated region A. The answer is found based on Poisson statistics: the probability that the counts observed at the maximum of an accumulator histogram - $x_{obs}$ - are due to background effects (cosmics+electronics noise) is given by:

$\begin{equation*}p(x_{obs}=cts\; \max, N_{A}=background\; cts) = Poisson(x_{obs},N_{A}) = \frac{N_{A}^{x_{obs}}}{x_{obs}!}e^{-N_{A}}\end{equation*}$

From the expression above the following P-value is computed for the total probability that background effects yield $x$ or above trigger primitives above the threshold:

$\begin{equation*}P(x\geq x_{obs}) = 1-\sum_{y=0}^{x_{obs}-1} \frac{N_{A}^{y}}{y!}e^{-N_{A}}\end{equation*}$

Correlation method Details on the correlation method can be found in Test beam operation of the CMS calorimeter trigger synchronization boards.

The results obtained with these methods are shown below.

2009 results

Friday 20 Nov 2009
Comments
Threshold set to 4 (1 GeV)
Sum Maximum distribution Signal intervals defined Rate in Background Region (A) Rate in Signal Region (B) Significance
[293; 393[
Delays (individual maximum - maximum in the total sum) p-value for maximum p-value as function of the accumulators summed up
Delays (from correlation with maximum distribution) Correlation Significance Case studies for some channels
 
Delays (from correlation with ECAL sum) Correlation Significance Case studies for some channels
 
List of channels falling behind most probable bunch channel list
List of channels laying ahead most probable bunch channel list

From Sat Nov 7 2009 to Sun Nov 8 2009
Comments
Three sets of histograms have been stored and are shown. Each row in the plots correspond to a different acquisiton time. After each acquisiton the histograms are reseted.
  • First row aquired in Sat Nov 7 23:37:19 2009 and corresponding run is 120020
  • Second row aquired in Sun Nov 8 11:20:21 2009 and corresponding run is 120025
  • Third row aquired in Mon Nov 9 01:17:57 2009 and corresponding run is 120042
  • Threshold set to 3 (corresponds to approximately 750 MeV)
  • Endcap histograms were not logged
Sum Maximum distribution Signal intervals defined Rate in Background Region (A) Rate in Signal Region (B) Significance
EB
[54; 154[
[54; 154[
[54; 154[
Delays (individual maximum - maximum in the total sum) p-value for maximum p-value as function of the accumulators summed up
Delays (from correlation with maximum distribution) Correlation Significance Case studies for some channels (2nd dataset only)
Clear evidence (no delay)
Clear evidence (delay=-1)
Noisy channel (delay=-6)
Noise channel (delay=0)
Almost no counts (delay=+1)
List of channels falling behind most probable bunch channel list
List of channels laying ahead most probable bunch channel list
List of channels not ok (empty histograms mostly) channel list

First LHC beam splash (2008)

(show) (hide)

A consequence of incoming beam is the increase in the occupancy of the SLB accumulators. First beam was expected around bunch 680 ($17 \mu s$) so an increase of the occupancy in this region was expected. Dividing the accumulator into 2 regions:

  • A: 0 < bx < 500 - cosmics and noise dominated region of the accumulator;
  • B: 650 < bx < 750 - expected region for beam;

one can compute the average number of counts per bunch recorded by the accumulators in these regions as

$\begin{equation *} N =\int_{b_{1}}^{b_{2}} H(b) / (b_{2}-b_{1}) \end{equation*}$.

The significance of the increase of counts due to beam can be estimated from:

$\begin{equation*}S_{occupancy} = \frac{N_{B}-N_{A}}{\sqrt{N_{A}}} = \frac{signal}{\sqrt{noise+background}} \end{equation*}$

The plots in the table below show the results for $N_{A}$,$N_{B}$, and $S_{occupancy}$ and one can conclude that something was observed in the region 650 < bx < 750 but with relatively low significance ($\approx 1-2$) on a single accumulator basis.

Occupancy increase in 650<bx<750
Region A: 0 < bx < 500 Region B: 650 < bx < 750 Significance for occupancy in region B

Timing

Looking for the bunch corresponding to maximum number of counts in each accumulator we get the results shown in the table below.

Bunch for maximum in accumulators
Bunch $N_{trigger\; towers}=N_{accumulators}$ Total
681 45
682 2152
683 3
<681 or > 683 393

One can ask if the maximum observed is compatible with the noise counts expected, as measured in the region A (0<bx<500). We can give an answer based on Poisson statistics: the probability that the counts observed at the maximum of an accumulator histogram - $x_{obs}$ - are due to background effects (cosmics+electronics noise) is given by:

$\begin{equation*}p(x_{obs}=cts\; \max, N_{A}=background\; cts) = Poisson(x_{obs},N_{A}) = \frac{N_{A}^{x_{obs}}}{x_{obs}!}e^{-N_{A}}\end{equation*}$

From the expression above we compute the following P-value: the total probability that background effects yield $x$ or above trigger primitives above the threshold is:

$\begin{equation*}P(x\geq x_{obs}) = 1-\sum_{y=0}^{x_{obs}-1} \frac{N_{A}^{y}}{y!}e^{-N_{A}}\end{equation*}$

Below we show the result obtained for the distribution of the maximum of each SLB accumulator and the corresponding P-value.

Timing results
Maximum bunch position with respect to bx=680 $P-value$

From the distribution for P-value we see that the large part of the accumulators has a low P-value which indicates that the counts observed at the maximum on the histogram cannot be due to background effects only. There is, however a considerable ammount of accumulators with high P-value (counts are dominated by background effects).

Instead of looking to individual timing we can look to the full ECAL barrel integrated timing by summing up all the SLB accumulators. In this case we proceed as follows:

  • sum the histograms one by one
  • for each new sum re-compute $N_{A}$ and $x_{obs}=cts\; max$
  • compute the P-value corresponding to $x_{obs}$
  • plot the P-value as function of the number of trigger towers used in the sum

The plots below summarize the results obtained:

Timing results (barrel integrated)
Accumulator sum Zoom for the accumulator sum in region B (650<bx<750) $P-value$ evolution with the number of accumulators used

From the behavior of the P-value distribution we conclude that the first LHC beam was correctly identified by the ECAL SLBs. The delay with respect to the expected bunch (680) was measured to be 2 bunches (50ns).

-- PedroSilva - 12 Sep 2008

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng AllEcalZoomH_small.png r1 manage 11.6 K 2009-11-24 - 10:56 PedroSilva  
PNGpng AverageCounts_0to500.png r2 r1 manage 10.2 K 2008-09-16 - 02:01 PedroSilva  
PNGpng AverageCounts_600to700.png r1 manage 11.8 K 2008-09-12 - 16:47 PedroSilva  
PNGpng AverageCounts_650to750.png r1 manage 11.8 K 2008-09-16 - 01:49 PedroSilva  
PNGpng BunchForAccMax.png r2 r1 manage 5.9 K 2008-09-16 - 01:49 PedroSilva  
PNGpng Correlation.png r1 manage 10.8 K 2008-09-12 - 16:48 PedroSilva  
PNGpng Correlation3_20.11.09.png r1 manage 19.3 K 2009-11-20 - 23:09 PedroSilva  
PNGpng CorrelationSignificanceEvol.png r1 manage 8.4 K 2008-09-12 - 16:48 PedroSilva  
PNGpng Correlation_08.11.09.png r2 r1 manage 47.6 K 2009-11-19 - 11:59 PedroSilva  
PNGpng Correlation_20.11.09.png r1 manage 19.0 K 2009-11-20 - 22:38 PedroSilva  
PNGpng CountsSignificance.png r3 r2 r1 manage 15.2 K 2008-09-16 - 01:50 PedroSilva  
PNGpng DelayMeasured.png r2 r1 manage 8.0 K 2008-09-16 - 01:51 PedroSilva  
PNGpng DelaySignificance.png r1 manage 13.5 K 2008-09-12 - 16:49 PedroSilva  
PNGpng DelaysInECAL2_08.11.09.png r2 r1 manage 34.1 K 2009-11-19 - 12:00 PedroSilva  
PNGpng DelaysInECAL2_20.11.09.png r1 manage 16.7 K 2009-11-20 - 22:39 PedroSilva  
PNGpng DelaysInECAL3_20.11.09.png r1 manage 16.0 K 2009-11-20 - 23:09 PedroSilva  
PNGpng DelaysInECAL_08.11.09.png r3 r2 r1 manage 31.0 K 2009-11-17 - 23:00 PedroSilva  
PNGpng ECALMaxDistributions_08.11.09.png r3 r2 r1 manage 9.4 K 2009-11-17 - 23:04 PedroSilva  
PNGpng ECALMaxDistributions_20.11.09.png r1 manage 2.9 K 2009-11-20 - 22:39 PedroSilva  
PNGpng ECALSums_08.11.09.png r3 r2 r1 manage 10.3 K 2009-11-17 - 23:04 PedroSilva  
PNGpng ECALSums_20.11.09.png r1 manage 3.3 K 2009-11-20 - 22:40 PedroSilva  
PNGpng EtaSpectrum.png r1 manage 12.9 K 2009-11-19 - 23:19 PedroSilva  
PNGpng FirstBeam.png r1 manage 62.5 K 2009-11-24 - 10:56 PedroSilva  
PNGpng IndividualPvalue.png r1 manage 17.4 K 2008-09-16 - 02:02 PedroSilva  
Texttxt ListAhead_08.11.2009.txt r2 r1 manage 0.1 K 2009-11-19 - 13:54 PedroSilva  
Texttxt ListAhead_20.11.2009.txt r1 manage 0.2 K 2009-11-20 - 23:09 PedroSilva  
Texttxt ListBehind_08.11.2009.txt r2 r1 manage 11.2 K 2009-11-19 - 13:58 PedroSilva  
Texttxt ListBehind_20.11.2009.txt r1 manage 1.4 K 2009-11-20 - 23:10 PedroSilva  
Texttxt NotOkList_08.11.2009.txt r1 manage 0.5 K 2009-11-19 - 12:04 PedroSilva  
PNGpng PtSpectrum.png r1 manage 17.1 K 2009-11-19 - 23:18 PedroSilva  
PNGpng PvalueEvol.png r1 manage 15.5 K 2008-09-16 - 01:50 PedroSilva  
PNGpng PvaluesEvolution_08.11.09.png r4 r3 r2 r1 manage 9.6 K 2009-11-17 - 23:20 PedroSilva  
PNGpng PvaluesEvolution_20.11.09.png r1 manage 3.3 K 2009-11-20 - 22:40 PedroSilva  
PNGpng RateBackgroundRegion_08.11.09.png r3 r2 r1 manage 25.8 K 2009-11-17 - 23:01 PedroSilva  
PNGpng RateBackgroundRegion_20.11.09.png r1 manage 16.5 K 2009-11-20 - 22:40 PedroSilva  
PNGpng RateSignalRegion_08.11.09.png r3 r2 r1 manage 19.8 K 2009-11-17 - 23:02 PedroSilva  
PNGpng RateSignalRegion_20.11.09.png r1 manage 17.3 K 2009-11-20 - 22:40 PedroSilva  
PNGpng SignificancesInECAL.png r1 manage 67.9 K 2009-11-16 - 13:22 PedroSilva  
PNGpng SignificancesInECAL2_08.11.09.png r2 r1 manage 62.6 K 2009-11-19 - 12:00 PedroSilva  
PNGpng SignificancesInECAL2_20.11.09.png r1 manage 22.1 K 2009-11-20 - 22:41 PedroSilva  
PNGpng SignificancesInECAL3_20.11.09.png r1 manage 22.6 K 2009-11-20 - 23:09 PedroSilva  
PNGpng SignificancesInECAL_08.11.09.png r3 r2 r1 manage 25.5 K 2009-11-17 - 23:03 PedroSilva  
PNGpng SignificancesInECAL_20.11.09.png r1 manage 22.3 K 2009-11-20 - 22:41 PedroSilva  
PNGpng Singlepvalues_08.11.09.png r3 r2 r1 manage 39.1 K 2009-11-17 - 23:01 PedroSilva  
PNGpng Singlepvalues_20.11.09.png r1 manage 16.6 K 2009-11-20 - 22:41 PedroSilva  
PNGpng SynchHisto.png r2 r1 manage 9.6 K 2008-09-16 - 01:52 PedroSilva  
PNGpng SynchHisto_650to750.png r1 manage 10.1 K 2008-09-16 - 01:52 PedroSilva  
PNGpng TCC37CH17.png r1 manage 13.6 K 2009-11-19 - 13:32 PedroSilva  
PNGpng TCC37CH65.png r1 manage 11.8 K 2009-11-19 - 13:30 PedroSilva  
PNGpng TCC43CH10.png r1 manage 13.2 K 2009-11-19 - 11:21 PedroSilva  
PNGpng TCC47CH7.png r1 manage 13.2 K 2009-11-19 - 11:42 PedroSilva  
PNGpng TCC53CH63.png r1 manage 11.4 K 2009-11-19 - 13:14 PedroSilva  
Edit | Attach | Watch | Print version | History: r13 < r12 < r11 < r10 < r9 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r13 - 2009-11-24 - PedroSilva
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright &© 2008-2024 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback