Likelihood Template Fitting
|
|
- Lenard Miles
- 5 years ago
- Views:
Transcription
1 Momentum Scale Estimation Using Maximum Likelihood Template Fitting by Yu Zeng Department of Statistics Duke University Date: Approved: David L. Banks, Chair David B. Dunson Ashutosh V. Kotwal A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in the Department of Statistics in the Graduate School of Duke University 2010
2 Abstract (Statistics) Momentum Scale Estimation Using Maximum Likelihood Template Fitting by Yu Zeng Department of Statistics Duke University Date: Approved: David L. Banks, Chair David B. Dunson Ashutosh V. Kotwal An abstract of a thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in the Department of Statistics in the Graduate School of Duke University 2010
3 Copyright c 2010 by Yu Zeng All rights reserved except the rights granted by the Creative Commons Attribution-Noncommercial Licence
4 Abstract A maximum likelihood template fitting procedure is performed by using Υ µ + µ events to extract the momentum scale, a scale factor applied to measured momentum, of the CDF detector at Fermilab. The constructed invariant mass spectrum from data events is compared with the invariant mass spectrum from Monte Carlo simulated events, with the momentum scale varying as a free parameter in the simulation. The invariant mass spectrum from simulation which best matches the data spectrum gives the maximum likelihood estimation of the momentum scale. We find the momentum scale is p/p = ( ± (stat) ± (syst) ) iv
5 Contents Abstract List of Tables List of Figures List of Abbreviations and Symbols iv viii ix xii 1 Introduction Motivation of W boson mass measurement Relationship between momentum scale measurement and W mass measurement The computing Thesis structure Overview of Momentum Scale Measurement General procedure Some general issues Statistical error and systematic error Binning issue Choice of step size in p/p χ 2 and goodness of fit Parameter tuning The Data and Monte Carlo Samples 15 v
6 3.1 The Data sample Invariant mass spectrum from data The Monte Carlo samples Invariant mass spectrum from simulation Effect on χ 2 due to finite MC sample size and background description Template Fitting with Maximum Likelihood Template construction Construction of the inclusive template Construction of non-inclusive templates in parameter space R Maximum Likelihood method Variance of Maximum Likelihood Estimator Bias in the fitter Necessary steps before the final momentum scale measurement Background parameterization Resolution tuning in detector simulation Data bias removal Material simulation check Results Momentum scale and the statistical error Momentum scale systematic errors PYTHIA event generation Recycling technique Binning Step size in p/p Fit range vi
7 5.2.6 Background parameterization Resolution model Applied data correction Material simulation Selection cuts Combine momentum scale statistical and systematic errors Error propagation Summary of systematic and statistical errors Summary 65 A Invariant Mass 66 B Selection Cuts 67 C The Poisson Distribution 68 C.1 The Poisson postulates D Other Systematic Errors 71 D.1 Quantum electrodynamics and energy loss model D.2 Magnetic field nonuniformities D.3 Trigger efficiency (p T cut) Bibliography 76 vii
8 List of Tables 1.1 Comparison of required precisions in M t and M W to constrain M H precision to 100 GeV Summary of systematic errors and statistical error on p/p and on M W. 64 B.1 Requirements on muons from Υ µµ viii
9 List of Figures 1.1 m W and m t contours [4] The CDF detector and the coordinate system Transverse mass m T spectrum from simulation. Blue histogram is for input M W = 80 GeV, and red histogram is for input M W = 81 GeV. The Jacobian edge of the m T distribution is quite sensitive to the input M W value Flow chart of momentum scale (or p/p) measurement The final invariant mass spectrum from data Distribution of p Υ T before and after tuning. Blue points represent data and the red histogram represents events from PYTHIA Comparison of invariant mass spectrum m µµ at different simulation stages from PYTHIA and from PYTHIA plus detector simulation Comparison of not recycling and recycling PYTHIA-level generated events 10 times. Plots (a) and (b) show the Data-MC fit without and with the recycling, respectively. Plots (c) and (d) show the negative log-likelihood curve (ln L ln L max ) as a function of parameter p/p without and with the recycling, respectively (see discussion of likelihood in Sec. 4.2) Standard error ellipse for the estimators ˆθ i and ˆθ j when correlation coefficient ρ ij < 0 [6] m µµ spectrum from MC for p/p = and p/p = Inclusive Monte Carlo templates constructed with Υ µ + µ simulated events. The p/p axis is binned into 800 bins from to 0 with equal width The m µµ axis is binned into 320 bins from 8.8 GeV to 11.2 GeV with equal bin width GeV (7.5 MeV). 28 ix
10 4.3 Validation plots show that the fitter is unbiased Inclusive Υ µµ fit for resolution scale factor h 0 = 0.968, and The χ 2 as a function of hit resolution scale factor h 0. The red line is for χ 2 = χ 2 min + 1. Blues lines indicate the locations of the statistical error and the central value p/p vs. cot θ before correction where cot θ = cot θµ + cot θ µ. The black curve is a parabolic fit to the blue points with function p 0 + p 1 ( cot θ) + p 2 ( cot θ) p/p vs. cot θ after correction where cot θ = cot θµ + cot θ µ. Correction parameters: α = , β = The black curve is a parabolic fit to the blue points with function p 0 + p 1 ( cot θ) + p 2 ( cot θ) Vary best tuned α = by ±0.001 with β = fixed Vary β = by ± with α = fixed p/p vs. 1/p T (µ) Final central value of p/p with its statistical error Log-likelihood curve and signed χ Distribution of p T (Υ) when p 1 varies ±0.001 from ˆp 1 = χ 2 as a function of p Υ µµ fit when ˆp 1 varies by 1σ Binning effect of m µµ on p/p Vary fit range by 20% Fit range ( , ) GeV for B level ±3σ band tuning Back to normal fit range ( , ) GeV Fit range ( , ) for B slope ±3σ band tuning Fit range ( , ) for B slope ±3σ band tuning Vary h 0 = by ±3σ Vary α by ±3σ x
11 5.14 Vary β by ±3σ Vary M scale by ±3σ Vary the cut on N hit by ±5 from Vary the cut on d 0 by ±0.1 from 0.3 cm Vary the cut on z 0 by ±1.0 from 3.0 cm D.1 χ 2 as a function of Qscale factor q D.2 Inclusive Υ µµ fit with Qscale factor q 0 = D.3 Inclusive Υ µµ fit without magnetic field correction D.4 M W with and without applying the magnetic field correction D.5 Inclusive Υ µµ when p T cut is raised by 0.2 GeV xi
12 List of Abbreviations and Symbols Abbreviations SM CDF COT MC KeV MeV GeV The Standard Model. The Collider Detector at Fermilab. It is a 100-ton complex detector which measures particles coming out of the proton antiproton collison. Central Outer Tracker, a component of CDF detector. It is a 3-meter long cylinder for measuring tracks. Monte Carlo simulation. A technique to make inference by performing a large number of trail runs and analyzing their colletive results. Kiloelectron Volt. Unit of energy, momentum and mass when c = 1, = 1. Megaelectron Volt. Unit of energy, momentum and mass when c = 1, = 1. Gigaelectron Volt. Unit of energy, momentum and mass when c = 1, = 1 (1 GeV = 10 3 MeV = 10 6 KeV = 10 9 ev). Symbols W M W ν Υ M Υ W particle. It has a positive or negative charge. The mass of the W particle. Neutrino. It has no charge. Υ particle. It can decay to a µ + µ pair. The mass of the Υ particle. It is known to be 9.46 GeV with extremely high precision. xii
13 µ Muon particle. It has a positive charge or negative charge with its mass precisely known to be MeV. µ will leave a track in COT. µ + Positively charged muon particle. µ Negatively charged muon particle. m µµ Invariant mass of µ + µ pairs. See Appendix A. m T Transverse mass of µν pairs which are decay products of W. Transverse mass is constructed in a similar way as invariant mass. The only difference is that it only uses information on the plane perpendicular to beam axis z at z = 0. W µν A W particle decays to a charged µ and a neutral ν. Υ µµ A Υ particle decays to a µ + and a µ. p T Amount of momentum perpendicular to beam axis z. p T = p 2 x + p 2 y. p Total amount of momentum. p = p 2 x + p 2 y + p 2 z. p T (Υ) p T carried by Υ. p T (µ ± ) p T carried by µ + or µ. 1/p T (µ) Average of 1/p T (µ + ) and 1/p T (µ ). p/p Correction of momentum normalized to total momentum. 1 + p/p Momentum scale. θ The angle of µ track relative to the proton beam direction. See Fig cot θ The cot θ difference between µ + and µ, i.e., cot θ µ + cot θ µ. cot θ Track parameter. Cotangent of the angle of µ track w.r.t. beam axis z. z 0 Track parameter. Initial position of µ track along beam axis z. c d 0 φ 0 Track parameter. Curvature of µ track. Track parameter. Impact parameter of µ track. Track parameter. Angle of µ track in x y plane. xiii
14 1 Introduction 1.1 Motivation of W boson mass measurement In high energy physics, also known as particle physics, there is a theory called the Standard Model [1], which is a quantum field theory describing the physical properties of all fundamental particles and their mutual interactions. The Standard Model is quite successful at explaining experimentally observed phenomena, and more importantly, it gave successful predictions of some undiscovered particles which were later confirmed to exist. Though the Standard Model predictions agree well with almost all experimental observations until now, a fundamental particle, also known as the Higgs boson, is still missing from the observed particle spectrum. However, the existence of the Higgs boson is required to make the Standard Model theory self-consistent. If the Standard Model is a self-consistent theory, then the mass of the undiscovered Higgs boson is related to the W boson mass (M W ) and the top quark mass (M t ). By measuring M W and M t precisely, we can infer the possible Higgs boson mass (M H ) within the Standard Model framework. Their mass relationship (in units 1
15 of GeV, see pg. xii) is best illustrated by Eqn. (1.1) [3]: M W = [ ( ) 2 ( ) ( ) Mt MH 1] ln ln 2 MH (1.1) where other terms with small contributions are neglected. From Eqn. (1.1) we can see higher values of M W suggests smaller values of the Higgs mass for a fixed M t, and vice versa. Assume the central values M W = GeV, M t = GeV, and M H = 100 GeV. Table 1.1 shows the required precision of M W compared with the required precision of M t for the same constraining power on the precision of M H. The M W, M t and M H relationship is shown in Fig. 1.1, where the diagonal lines are for different values of M H and the blue contour corresponds to a 68% confidence level contour of M W and M t measurements. We can see the constraint on M H from the current measurement precision of M W is not as powerful as that from the measurement precision of M t. It is thus important to measure the W boson mass more precisely. Table 1.1: Comparison of required precisions in M t and M W to constrain M H precision to 100 GeV. Parameter Shift (GeV) M t (GeV) M W (GeV) M H = Relationship between momentum scale measurement and W mass measurement As one of the most important, yet difficult, analyses in experimental high energy physics, the W boson mass is measured at hadron colliders by constructing transverse 2
16 August LEP2 and Tevatron (prel.) LEP1 and SLD 68% CL m W [GeV] α m H [GeV] m t [GeV] Figure 1.1: m W and m t contours [4]. mass m T, the shape of which is a function of M W as well as other factors: m T = 2p µ T pν T (1 cos φ) (1.2) where µ (muon) and ν (neutrino) are decaying particles of W boson via W µν; p µ T (pν T ) is the momentum of µ (ν) in the transverse plane, a plane perpendicular to beam axis z; φ is the angle between µ and ν in the transverse plane. The coordinate system and the CDF detector are shown in Fig The m T spectrum constructed using Eqn. (1.2) is shown in Fig.1.3. The transverse mass is used because it is a kinematic quantity which is invariant under longitudinal lorentz boosts (boosts along the beam axis z). The m T distribution is not like any of the statistical distributions we already know. The shape of m T arises from a combined effect of the true m W value, dynamics of W boson production (p T (W )), W boson decay angular distribution, and 3
17 Figure 1.2: The CDF detector and the coordinate system. Figure 1.3: Transverse mass m T spectrum from simulation. Blue histogram is for input M W = 80 GeV, and red histogram is for input M W = 81 GeV. The Jacobian edge of the m T distribution is quite sensitive to the input M W value. 4
18 the detector effects which the W boson decay products will experience. The m T Jacobian edge position near 80 GeV is most sensitive to M W. Once the detector response and the W boson production and decay mechanisms are well understood and simulated, we can fit the m T distribution from Monte Carlo simulation against m T distribution obtained from real data to extract M W. By inspecting Eqn. (1.2), we can see that if the measured p µ T is bigger than the true p µ T, then the mmeasure T distribution will shift to the right compared with the m true T distribution. This will lead to a bigger M measure W smaller than the true p µ measure T, then we will get a smaller MW than MW true. If the measured pµ T is than M true. It is thus W important to measure p µ T (or pµ ) precisely. A proposed way to get the most precise p µ is using a particle, whose mass is precisely known, as a reference to calibrate the momentum p µ measured in the CDF detector. The process of determining how much the measured momentum should be corrected to match the true momentum, thus to match the mass of the particle, is called momentum scale (1 + p/p) measurement. Once this scale factor (1 + p/p) is found, it will be applied to the measured momentum to get the best estimate of true momentum p true = (1 + p/p) p measure. (1.3) One way to measure (1 + p/p) (or equivalently p/p) is to use Υ(1S) µ + µ events. The Υ(1S) particle is a flavorless b b meson which was first discovered at Fermilab in 1977 [5] with its mass precisely measured to be at 9.46 GeV by other experiments (see Particle Data Group [6]). Its well-known mass value makes it an ideal candidate to precisely calibrate CDF-measured momenta. (1S) is a label to distinguish it from other similar resonances with higher masses, which are known as Υ(2S) and Υ(3S). For simplicity we will just use Υ to denote Υ(1S) throughout this thesis when no confusion could arise. Using Υ µ + µ data events, we can construct an invariant mass spectrum m Data µµ (See Appendix A for the definition of 5
19 invariant mass). Similary, using Υ µ + µ simulated events, we can construct a m MC µµ spectrum. If we repeat the construction of the m MC µµ spectrum by applying a different scale factor (1 + p/p) to the measured momentum, we can construct a series of such m MC µµ spectra, with each of the spectra i generated at a pre-defined value (1 + p/p) i. The m Data µµ spectrum is then compared with each of the m MC µµ spectra to figure out which MC spectrum matches the data best. The m MC µµ spectrum k that matches the m Data µµ spectrum best then tells us the momentum scale should be (1 + p/p) k. The way to find this best match is by using a maximum-likelihood method for each m MC µµ template. See discussions in Chapter 4 for details. 1.3 The computing All analysis codes are written in C++ except the PYTHIA [7] event generation, which is written in FORTRAN. Some shell scripts are used for automatic job submission and automatic background tuning. Monte Carlo simulation jobs are submitted to the CDF Analysis Farm [8], a network of computer farms, to run hundreds of the jobs at the same time instead of running them one by one. All outputs are directed to the same location upon completion and are then recombined. The typical time to run 100 parallel jobs of the Υ µ + µ simulation is hours. 1.4 Thesis structure This thesis consists of 6 chapters. Chapter 1 gives the general introduction; Chapter 2 to Chapter 5 discusses the momentum scale estimation using maximum likelihood template fitting, as well as the statistical and systematic errors; Chapter 6 summarizes the obtained results. Additional discussions relevant to this thesis are presented in Appendices. Terminologies used in this thesis are collected and explained in the section of Abbreviations and Symbols on pg. xii and pg. xiii. 6
20 Overview of Momentum Scale Measurement General procedure As illustrated in Fig. 2.1, the momentum scale measurement contains 3 major steps: STEP 1: On the data side, we use Υ µ + µ data events collected by the CDF detector to construct a binned invariant mass spectrum m Data µµ. The location and shape of the m Data µµ spectrum is a combined effect of Υ µ + µ production and the interaction of decaying µ ± with the CDF detector (see the definition of invariant mass in Appendix A). STEP 2: On the simulation side, we take three sub-steps to construct a binned invariant mass spectrum m MC µµ. STEP 2.1: The first step is to use the Monte Carlo event generator PYTHIA [7] to simulate the production of Υ µ + µ events. At this stage, no detector effects are simulated. The invariant mass m P ythia µµ constructed at this stage is sharply peaked at m Υ and approximates a δ-function compared with the spread of m Data µµ. STEP 2.2: As a next step, these simulated events are propagated through a fast detector simulation to simulate how µ ± interacts with materials inside the detector. 7
21 Raw Data MC generator ( interested events) Reconstruction Track Parameters(TP) Condensed Data (Root trees) Detector simulation Reconstruction Event Selection Data TP Updated TP Invariant mass Invariant mass Data Histogram ML fitting Simulation Histograms Correct Data Yes Data bias? Simulation OK? No No Yes Final result! Figure 2.1: Flow chart of momentum scale (or p/p) measurement. 8
22 When this step is done, the mass spectrum m MC µµ will become much wider than m P ythia µµ due to detector resolution effects and it will look similar to m Data µµ if the detector simulation is correct. STEP 2.3: The background observed in the m Data µµ spectrum is parameterized and added to the m MC µµ spectrum to form an updated m MC µµ spectrum (see Sec for the discussion of background tuning). STEP 3: We then fit the m MC µµ spectrum against the data m Data µµ spectrum and can get a negative log-likelihood value ln L using the Poisson distribution as the probability function (see Sec. 4.2 for detailed discussion of likelihood function). If we scale the momentum in simulation by a scale factor (1 + p/p), and repeat STEP 2.2, we can get a new m MC µµ spectrum in STEP 2.3. Thus we can get a new ln L in STEP 3. If we repeat the above processes by varying p/p, we can get a series of ( p/p, ln L) pairs. If we plot them on x-y plane with the x axis as the trial value of p/p and corresponding ln L values on the y axis, we can get a ln L curve as a function of p/p. The x-coordinate of the minimum of ln L (equivalently the maximum of ln L) is the Maximum Likelihood Estimator (MLE) p/p. Thus we know momentum measured in our data should be corrected as p true = (1 + p/p) p measure. (2.1) Besides extracting the MLE p/p, the above three steps can be used to explore possible bias in data. They can also be used to validate the parameterization in our detector simulation. Actually, before getting a reliable p/p, we have to remove data bias first as well as make sure our detector simulation agrees with data. 9
23 2.2 Some general issues Statistical error and systematic error Statistical error is a random deviation between an observed value and the true value caused by the sampling or measurement noise. Statistical error can be reduced by repeated measurements. Systematic error is the bias in measurement caused by incomplete knowledge of nuisance parameters [9]. It cannot be reduced by repeating measurements. In our analysis, the m µµ spectrum is a function of p/p as well as other nuisance parameters listed in Table 5.1. We vary those parameters in the simulation to quantify the corresponding systematic errors Binning issue Our whole analysis is constructed based on histograms (binned spectrum) because the binning method is advantageous for fast calculations (see Fan and Marron [10], Wand [11] for detailed discussion). Generally speaking, we need to consider the following facts in choosing the number of bins N bin : 1) N bin should not be too small, otherwise we will lose the resolution in original measurements; 2) For presentation reasons and to make the shape of spectrum clear, N bin should not be too large [12]. We choose N bin to satisfy the above two conditions and confirm the bias caused by our choice of bins is negligible compared with the precision we want to reach. In our analysis, we deal with binning in m µµ and R with R to be cot θ or 1/p µ T. To get a high precision measurement, we need fine binning in m µµ. As to the binning in R (R = cot θ or 1/p µ T ), we just need N bin to be sufficient to capture the overall bias in data [13] [14] [15] (we use 16 cot θ bins in range (-1.6, 1.6)), or to check the material description in simulation [13] [14] [15] (we use 12 1/p µ T bins in range 10
24 (0, 0.3)). Binning in R space is discussed in the construction of the non-inclusive template in Sec As to the binning in m µµ, we force the binning of m MC µµ spectrum and m Data µµ spectrum to be the same, with N bin of m Data µµ spectrum chosen such that we can see a clear Υ peak with clear tails. The bias on fitted p/p due to the binning in m µµ, which turns out to be very small and can also be neglected, is studied in detail in Chapter Choice of step size in p/p We vary trial values of p/p from to 0 with a fixed step, which directly relates to the precision of p/p we want to measure. The range ( 0.004, 0) is chosen because other analysis [13] has shown p/p lies within a negative region close to 0. We choose the step size in p/p in such a way that the bias on p/p, when converted into M W, is negligible when compared with the goal of measuring M W within 25 MeV. In principle, if the size of step goes to 0 (thus infinite number of p/p trial values), we can remove this bias. In our analysis, we divide the p/p range ( 0.004, 0) into 800 regions with step size ɛ p/p = , which will introduce a maximum bias of 0.4 MeV on M W. Since the bias is so small compared with 25 MeV, its effect can be neglected and further reducing the step size in p/p would not help in improving the precision of M W measurement χ 2 and goodness of fit One commonly used goodness-of-fit test is based on Pearson s χ 2 statistic [16], which can tell the level of agreement between an observed data histogram and a MC simulated histogram. The χ 2 goodness-of-fit test approximates χ 2 in distribution under the null hypothesis when the following assumptions are satisfied [17]: 1) Random sampling of data; 2) A sufficiently large sample size; 3) Adequate number of entries 11
25 in a given category; 4) Independent observations. For m Data µµ and m MC µµ with the same binning and assuming the MC histogram bin has no associated statistical error, the Pearson s χ 2 is given by: χ 2 = n i=1 ( ) N Data i Ni MC 2 Ni MC (2.2) with n i=1 N Data i = n i=1 N MC i (2.3) where n is the total number of bins within a pre-selected fit range along m µµ axis; N Data i is the observed number of data events in bin i; N MC i is the number of simulated events in bin i with the background and the MC-to-Data normalization also taken into account. Details of how to choose the fit range, how to parameterize the background and how to calculate N MC i are discussed in Sec , Sec and Sec. 4.2, respectively. The change from the calculated χ 2 using Eqn. (2.2) due to our large but finite MC sample and due to uncertainties in background model is discussed in Sec Since χ 2 statistic gives the sum of squares of the differences between data and MC in units of the standard deviation, a large χ 2 usually means a large discrepancy between data and MC. It can also be shown that the χ 2 statistic asymptotically approaches a χ 2 distribution with degrees of freedom df = n p, where n is the number of bins and p is the reduction in degrees of freedom [18]. In case of extracting p/p via Maximum Likelihood method (see Sec. 4.2), the calculated χ 2 statistic will approximate a χ 2 n 2 distribution if the size of data sample is sufficiently large. The degree of freedom df equals to n 2 rather than n because of the normalization and the need to estimate p/p. Since the expected value of χ 2 n 2 distributed variable is n 2, the expectation value of the variable χ 2 /(n 2) should be 1 and thus χ 2 /(n 2) is often quoted as a measure of agreement between 12
26 data and MC. If χ 2 /(n 2) is much greater than 1, then the theoretical model used to generate MC samples is deficient; if χ 2 /(n 2) is much less than 1, we may need to check whether the variance used in calculating the χ 2 statistic (N MC i in our case) is overestimated. A P -value [18] [19] can also be calculated by comparing the value of the χ 2 statistic with a χ 2 n 2 distribution. Since the P -value is defined to be the probability of obtaining a test statistic, under the null hypothesis, at least as extreme as the one that was observed, we can thus calculate the P -value by integrating the χ 2 n 2 distribution from the observed χ 2 to infinity with P = χ 2 f(x; n 2)dx (2.4) f(x; n 2) χ 2 n 2. If one does not extract any parameter but just normalizes a MC histogram to the data histogram, as shown in Fig. 3.2 and Fig. 5.3, the χ 2 statistic will approximate a χ 2 n 1 distribution, with n as the total bin number. The P -value will then be with P = χ 2 f(x; n 1)dx (2.5) f(x; n 1) χ 2 n 1. The R [20] package supplies a command 1-pchisq(χ 2, df) to calculate this P -value. For presentation reasons, the calculated χ 2 is rounded to the nearest integer and added to each plot in this thesis. When we tune parameters and study systematic errors, the original χ 2 values are used Parameter tuning Throughout this thesis, our simulation parameters are tuned to their best values by minimizing the calculated χ 2 in a pre-selected fit range. This is based on the 13
27 least square method, which is equivalent to maximum likelihood method when the Poisson distributed N Data i approximates to a Gaussian [18]. In addition, it can be shown that one standard deviation from the best tuned parameter value will increase the χ 2 from χ 2 min to χ 2 min + 1 [21] [22]. We can thus vary a given parameter until χ 2 increases from χ 2 min to χ 2 min + 1 to find the corresponding one standard deviation of that parameter. 14
28 3 The Data and Monte Carlo Samples 3.1 The Data sample Data in High Energy Physics (HEP) take several forms [23] depending on the stages of data processing. In the time flow, it is called raw data, reconstructed data and skimmed (condensed) data in ROOT-based n-tuple formats, where ROOT is an object oriented framework for large scale data analysis in HEP [24]. We use the ROOT-based n-tuple in data analysis. A ROOT n-tuple is a collection of events categorized according to events properties and it is a derived data format to be used for the final analysis in high energy physics. For our analysis, we use n-tuples as the input data format and apply standard selection cuts (see Appendix B), then reconstruct m Data µµ. This approach is equivalent to first writing out µ ± track parameters for each event into a text file, then using those track parameters to construct m Data µµ spectrum. Without loss of generality, we can treat the data as stored in a text file with each line representing a Υ µ + µ event. For a given such event, there are 10 associated quantites with the first 5 quantities and last 5 quantities being the track parameters of µ + and µ, 15
29 respectively. Those track parameters of µ ± already contain the CDF detector effects Invariant mass spectrum from data The Υ µ + µ data used in this thesis were recorded by the CDF detector from Febuary 4, 2002 to August 4, 2007 [25]. We have a total of 464, 802 Υ µµ candidate events within the invariant mass range ( , ) GeV after our standard selection (see Appendix B). The range is chosen such that it can cover the observed spread of Υ peak and we use the same range to fit MC histograms against the data histogram. The number 464, 802 includes the background contribution. If we subtract the best tuned linear background, as discussed in Sec. 5.1, we get within the same fit range 238, 999 events, which can be viewed as the total number of pure Υ µµ events in the data. The final m µµ spectrum constructed from data is shown in Fig. 3.1, where the peak around 9.43 GeV is the Υ resonance and the peak around 10 GeV is Υ(2S) which is not of interest to us. Since Υ(2S) can also decay to µ + µ, they appear in the reconstructed m µ + µ spectrum. For our analysis, we will only consider the Υ peak. We choose the bin width such that there are 50 bins within the fitting window ( , ) GeV. The possible systematic bias caused by this binning method is studied by doubling the bin width and halving the bin width. A detailed discussion is shown in Chapter The Monte Carlo samples After the simulated production of Υ µ + µ and simulated interaction of µ ± with CDF detector, we write out the track parameters for each generated µ. Similar to the format of data events, each simulated event is stored in a 10-column text file with the first 5 columns being the track parameters associated with µ + and the last 5 columns being the track parameters associated with µ. 16
30 events m µµ (GeV) Figure 3.1: The final invariant mass spectrum from data Invariant mass spectrum from simulation The simulation of Υ µ + µ events within CDF detector has two stages. The first stage is event generation which only simulates the production of Υ µµ events without any detector effects. We generate 50 million Υ µ + µ events using the MC event generator PYTHIA [7]. Parameters which control the transverse momentum distribution of the generated Υ (p MC T (Υ)) are tuned to match p Data T (Υ) because PYTHIA does not describe it well. See Fig. 3.2 (a) and (b) for comparison before tuning and after tuning. The systematic uncertainty due to the tuning parameters are studied in Chapter 5. At this stage, if we construct the invariant mass m P ythia µµ using PYTHIA output, we will get a sharp spectrum with very narrow width, see Fig. 3.3 (a). The second stage of simulation is simulating the CDF detector effects on µ ± tracks produced by PYTHIA. At this stage, the same standard selection cuts (Appendix B), which are already applied to the data, are applied to the simulated events. After this second stage, track parameters of µ ± get updated from their original PYTHIAgenerated values. Using the updated track parameters, we can construct a new invariant mass spectrum m P ythia+detector µµ. As shown in Fig. 3.3 (b), the m P ythia+detector µµ 17
31 p (#ups) (GeV) T (a) Distribution of p Υ T before tuning ) χ 2 / DoF = / p (Υ (GeV) (b) Distribution of p Υ T after tuning T Figure 3.2: Distribution of p Υ T before and after tuning. Blue points represent data and the red histogram represents events from PYTHIA. spectrum is about 1200 times wider than the m P ythia µµ spectrum indicating most of the effects come from the detector simulation (the number 1200 is obtained by taking the ratio of the RMS on plot (a) and the RMS on plot (b)). This fact indicates that even if we feed the same PYTHIA event into the detector simulation, we will still get quite different behavior and quite different invariant mass m P ythia+detector µµ due to numerous random processes in the detector simulation [26]. The expected number of times we can use the same event in Fig. 3.3 (a) to simulate events in Fig. 3.3 (b) before the simulated events become identical is at least on the order of Thus, if we use the same PYTHIA event just a few times ( 10 3 ) through the detector simulation, we will get events with quite different behavior and they can be viewed, to a large extent, as uncorrelated. Considering the time and computing resources needed 1 to generate 10 times more PYTHIA events, and using the fact that almost all of the limitation comes from the detector simulation, we feed each of the 50 M PYTHIA generated events 10 times into the detector simulation to mimic the process of generating 500 M PYTHIA events and 1 (a). It is time consuming to generate 500 M events; (b). The CDF Analysis Farm (CAF) has limitations on the size of the input files that store the PYTHIA events. It also has limitations on the number of jobs to be submitted parallely at the same time. 18
32 events HepgMass Entries e+07 Mean 9.46 RMS events Upsmass Entries e+07 Mean RMS Pythia m µµ (a) m µµ after PYTHIA event generation Pythia+Sim m µµ (b) m µµ after PYTHIA event generation and detector simulation Figure 3.3: Comparison of invariant mass spectrum m µµ at different simulation stages from PYTHIA and from PYTHIA plus detector simulation. then feeding into the detector simulation. The possible bias caused by this technique is shown in Fig We see a change in p/p to be , which is a 0.8 MeV effect on M W. This change comes from a combined effect of the bias caused by recycling and the bias caused by the statistical fluctuation due to the size of the MC samples, which can be seen in the minimum region in Fig. 3.4 (c). This small change in p/p is negligible when compared with the 25 MeV precision goal of the M W measurement. We thus use the recycled MC samples as the default MC samples and quote 0.8 MeV as a systematic error along with other systematic errors in Chapter Effect on χ 2 due to finite MC sample size and background description For the MC histogram which matches with the data best (i.e., the 534th histogram), there are events within ( , ) GeV interval before normalizing to data entries and before adding backgrounds to MC. By comparing in the fit range ( , ) the total number of data events after background subtraction and the the total number of pure MC signal events, we can 19
33 events p/p = (-1.32 ± ) x 10 stat events p/p = (-1.33 ± ) x 10 stat χ 2 /dof = 61 / χ 2 /dof = 61 / m µµ (GeV) (a) Without recycling PYTHIA-level events m µµ (GeV) (b) Recycling PYTHIA-level events 10 times ln L max - ln L ln L max - ln L p/p (c) Without recycling PYTHIA-level events p/p (d) Recycling PYTHIA-level events 10 times Figure 3.4: Comparison of not recycling and recycling PYTHIA-level generated events 10 times. Plots (a) and (b) show the Data-MC fit without and with the recycling, respectively. Plots (c) and (d) show the negative log-likelihood curve (ln L ln L max ) as a function of parameter p/p without and with the recycling, respectively (see discussion of likelihood in Sec. 4.2). get the MC scaling factor S = (3.1) The effect of simulation statistical error on the calculated χ 2 (denoted as χ 2 to distinguish it from χ 2 in Sec ) can be illustrated by the following derivations: χ 2 = n i=1 ( N Data i N MC i ) Ni MC 2 + σ 2 Ni MC (3.2) 20
34 with N MC i = S Ni MC (0) + N Bkg i (3.3) σ 2 N MC i = σ 2 S Ni MC (0) + σ2 N bkg i = S 2 Ni MC (0) + σ 2 N bkg i (3.4) where m is the total number of bins within a pre-selected fit range, Ni MC (0) is the original number of pure MC events in bin i, N Bkg i i, σ 2 N bkg i is the number of bkg events in bin is the variance of background in bin i. If we plug Eqn. (3.3) and (3.4) into Eqn. (3.2), we will get χ 2 = = n ( i=1 n i=1 = χ 2 S N MC i [ ( Ni Data 1 + [ ( Ni Data (0) + N Bkg i S N MC i S Ni MC ) ( + S Ni MC (0) + N Bkg i 1 S2 N MC i S N MC i (0)+σ 2 N bkg i (0)+N Bkg i )] 2 (0) + N Bkg i S 2 N MC i )] 2 (0) + N Bkg i (0) + σ 2 N bkg i S2 N MC i S N MC i ) (0)+σ 2 N bkg i (0)+N Bkg i (3.5) with χ 2 the same as in equation as Eqn. (2.2), assuming the size of MC sample to be infinity and assuming no error on background parametrization: χ 2 = n i=1 [ ( Ni Data S N MC i S Ni MC (0) + N Bkg i )] 2 (0) + N Bkg i = n i=1 ( ) N Data i Ni MC 2 Ni MC. (3.6) As discussed in Sec , the background is parametrized as N bkg i = B level + B slope (m µµ 9.43) (3.7) 21
35 with ˆB level = 4517 (σ Blevel 7) and ˆB slope = 250 (σ Bslope 39), which are discussed in detail in Sec We thus have σ 2 N bkg i = σ 2 B level + σ 2 B slope (m µµ 9.43) (m µµ 9.43) Cov(B level, B slope ) (3.8) If there is no constraint, B level and B slope should be independent. When using χ 2 minimization as the criterion to tune parameters B level and B slope, they become weakly correlated. If we vary B level by ±1σ Blevel from its central value 4517, B slope needs to be changed by 2 to keep the χ 2 in the pre-selcted fit range minimized. As we can see from Fig. 3.5, the correlation coefficient between B level and B slope should satisfy ρ σ Bslope = 2. Thus we have ρ 0.05 and Cov(B level, B slope ) = ρ σ Blevel σ Bslope (3.9) Figure 3.5: Standard error ellipse for the estimators ˆθ i and ˆθ j when correlation coefficient ρ ij < 0 [6]. Thus Eqn. (3.8) will become σ 2 N bkg i (m µµ 9.43) (m µµ 9.43). (3.10) 22
36 Since we fit in the region < (m µµ 9.43) < 0.185, we have: σ 2 N bkg i (m µµ 9.43) (3.11) Using the value of S obtained in Eqn. (3.1) and with the m µµ peak estimation from Fig. 3.3 (b) to be , we have S 2 Ni MC (0) = 1.9, which is much smaller than σ 2. This means the variance from background dominates N bkg i the correction term in χ 2 calculation. When m µµ = 9.43 (the bin at the peak position), Eqn. (3.10) gives σ 2 N bkg i Thus Eqn. (3.5) can be approximated as 49. χ 2 χ χ 2 (3.12) where is read from the peak of Fig When m µµ = , which corresponds to the bin farthest away from the peak position in the left tail within the fit range ( , ), we get σ 2 N bkg i 106. Eqn. (3.5) then becomes χ 2 χ χ 2 (3.13) where 4563 is obtained by using Eqn. (3.7). When m µµ = , which corresponds to the bin farthest away from the peak position within the fit range ( , ), we get σ 2 N bkg i
37 Eqn. (3.5) then becomes χ 2 χ χ 2 (3.14) where 4471 is obtained by using Eqn. (3.7). The real χ 2 should thus be somewhere between the above extreme cases: 0.977χ 2 χ χ 2. (3.15) This means the finite MC sample and the background parameterization used in our analysis will reduce the χ 2, which is calculated by using Eqn. (2.2) with the assumption of infinite MC sample and accurate background theory, by at most 2.3%. Since the effect reduces χ 2 by 2% at most, it has a small effect (at most 1%) on the 1σ band of simulation paramters, which are obtained by increasing the χ 2 by 1 from χ 2 min. To cover this effect, we can scale those systematic errors (see Table 5.1) by a factor of The increase in those systematic errors turns out to be negligible when compared with the 25 MeV precision goal of M W measurement. In summary, we do not need to worry about the χ 2 calculation using Eqn. (2.2) since it is an approximation good enough for our purpose. A more general discussion of fitting using finite Monte Carlo samples can be found at [27]. 24
38 Template Fitting with Maximum Likelihood Template construction In our analysis, a template is a binned multi-dimensional histogram generated using MC simulated Υ µ + µ events [28]. It has discrete trial values of p/p with equal distances in the x axis and binned invariant mass spectrum m µµ in the y axis. The number of events is in the z axis. There are two kinds of templates used in this analysis. One is an inclusive template, which is generated using all simulated events; the other is a non-inclusive template in bins of physical quantity R, which can give valuable information on how the extracted p/p varies as a function of R. The sum of the non-inclusive templates over all bins in R will lead to the inclusive template and the splitting of the inclusive template according to the range of R values gives non-inclusive templates. We use non-inclusive templates in R = cot θ bins, where cot θ = cot θ µ + cot θ µ (see definitions of θ and cot θ on pg. xi), to study the possible bias in our data which has not been corrected by off-line calibration [29] and alignment [30]. Details of data correction in cot θ bins are shown in Sec After the data bias is removed, we 25
39 then use non-inclusive templates in R = 1/p µ T bins to check the material description in our detector simulation, since the linear dependence of p/p on 1/p µ T indicates the amount of material parameterized in our fast detector simulation does not match the real detector material (see reference [13] for the reasoning). The validation of material in our simulation is shown in Sec.4.4. Once the above two steps are successfully done, we can claim we understand the data and our simulation describes the data well. As a final step, we then use the inclusive template to fit against m Data µµ to extract the final momentum scale correction factor p/p Construction of the inclusive template For a given simulated event, we can calculate the momentum p µ and energy E µ using the track parameters associated with each µ. Once we have p µ and E µ, we can use the definition of invariant mass (see Appendix A) to calculate the invariant mass of µ + and µ : m µ + µ = (E µ + + E µ ) 2 ( p µ + + p µ ) 2 (4.1) with E µ ± = m 2 µ + (p x µ ±2 + p y µ ±2 + p z µ ±2 ). (4.2) If we repeat the calculation for each simulated event using Eqn. (4.1) and (4.2), we will then get a distribution of m MC µ + µ which is called the invariant mass spectrum of µ + µ. We bin this m MC µµ spectrum in such a way that it has the same binning as the m Data µµ spectrum. Using all simulated events, a m MC µµ inclusive template is formed by a series of m MC µµ binned histograms which are generated by multiplying momentum scale factor (1 + p/p) in front of p i µ ± (i = x, y, z) in Eqn. (4.1) and (4.2). As the p/p varies from to 0 with a step of , or equivalently as (1 + p/p) varies from to 1 with a step of , the 800 constructed m MC µ + µ 26
40 spectra will be slightly different from each other; see Fig. 4.1 for the comparison of m MC µµ spectrum generated at p/p = and m MC µµ spectrum generated at p/p = 0. All other m MC µµ spectra will lie between them events p/p=0 Entries e+08 Mean RMS p/p= Entries e+08 Mean RMS m µµ (GeV) Figure 4.1: m µµ spectrum from MC for p/p = and p/p = 0. For ease in accessing all p/p trial values in an efficient way, we re-arranged the inclusive template by turning the 1-to-1 correspondence between m MC µµ (i) and ( p/p) i into a binned 3-d histogram, with m MC µµ and p/p as its y and x axes, respectively, and N event the z axis. This is illustrated in Fig. 4.2, where a slice of 3-d histogram along p/p axis is the corresponding invariant mass spectrum m MC µµ (i) for a given value ( p/p) i. Since it is a 1-to-1 correspondence between p/p input value and the generated m MC µµ spectrum, we do not need to worry about the case where two or more spectra are generated with the same p/p input. Using this template, the fitter goes through all 800 m MC µµ histograms and fits against the m Data µµ spectrum by constructing the likelihood function (see Sec. 4.2). Thus, we will have 800 likelihood values. The m MC µµ spectrum that matches m Data µµ spectrum best will give the maximum likelihood. Once we find the best m MC µµ spectrum i, we will know the associated p/p value 27
41 Figure 4.2: Inclusive Monte Carlo templates constructed with Υ µ + µ simulated events. The p/p axis is binned into 800 bins from to 0 with equal width The m µµ axis is binned into 320 bins from 8.8 GeV to 11.2 GeV with equal bin width GeV (7.5 MeV). p/p = (i 1). This p/p is thus the Maximum Likelihood estimate of parameter p/p. The issues of step size choice in p/p and binning in m µµ are discussed in Chapter 2. Our step size choice in p/p and binning choice in m µµ turn out to be negligible when compared with the precision we want Construction of non-inclusive templates in parameter space R The same idea can be extended to the construction of the non-inclusive template in R (R = cot θ or R = 1/p µ T ) with different purposes. As mentioned before, templates in cot θ bins can help to figure out the bias in data by studying how the extracted p/p varies as a function of cot θ. Templates in 1/p T (µ) bins are used to determine whether the material description in our simulation matches the data. 28
42 We construct 16 non-inclusive templates in R = cot θ bins by selecting events passing requirement (i 1) < cot θ < i with i = 1,, 16 and construct 12 non-inclusive templates in R = 1/p T (µ) bins by selecting events passing requirement (j 1) < 1/p T (µ) < j with j = 1,, 12. There, each of the non-inclusive templates looks like Fig. 4.2, but with fewer N events since it is just a sub-sample of the inclusive template. Take R = cot θ as an example: we bin the cot θ space from 1.6 to 1.6 into 16 bins with equal bin width 0.2. The range (-1.6, 1.6) is chosen based on a previous analysis [13]. For a given cot θ bin i (i = 1,, 16), we first select simulated events passing the requirement (i 1) < cot θ < i, then use those selected events to construct m MC µµ spectrum as p/p varies in the same way as discussed in Sec For each cot θ bin, we can generate a corresponding template in forms of 3-d binned histogram and get a MLE of p/p associated with that cot θ bin. By repeating the same process over all 16 cot θ bins, we can get 16 MLEs of p/p. Using those 16 MLEs and 16 cot θ ranges, we can then study the dependence of MLEs on cot θ. Since different values of cot θ represent different decay topologies [13], the dependence of MLEs on cot θ means the MLEs depend on how Υ decays to µ + µ, which is contradictory to the physics principle that invariant mass is invariant in any reference frame. A conclusion is then drawn that there must exist a bias in the data in the measurement of mass, causing the measured mass to depend on cot θ µ ±. This bias needs to be corrected using a physicsbased hypothesis such that the dependence of MLE on cot θ is removed. See Sec for detailed discussion of removing the data bias. 4.2 Maximum Likelihood method The method of Maximum Likelihood (ML) is a statistical technique to estimate the parameters given a finite sample of data [18] [31] [32]. It is the core technique for 29
Measurement of Properties of Electroweak Bosons with the DØ Detector
Measurement of Properties of Electroweak Bosons with the DØ Detector Laboratoire de Physique Subatomique et de Cosmologie, 53, rue des Martyrs, 38026, Grenoble Cedex, France. E-mail: Hengne.Li@in2p3.fr
More informationMeasurement of the mass of the W boson at DØ
Measurement of the mass of the W boson at DØ 1 IPNL, Université de Lyon, Université Lyon 1, CNRS/IN2P3 4 Rue E. Fermi 69622 Villeurbanne, France E-mail: kurca@in2p3.fr We present a measurement of the mass
More informationCDF top quark " $ )(! # % & '
$% CDF quark 7 3 5 ( "#! Tevatron Run II Started Spring 1. proton-antiproton collider with (Run I :. antiproton recycler commissioning electron cooling operational by Summer 5. increase in luminosity.
More informationResults from the Tevatron: Standard Model Measurements and Searches for the Higgs. Ashutosh Kotwal Duke University
Results from the Tevatron: Standard Model Measurements and Searches for the Higgs Ashutosh Kotwal Duke University SLAC Summer Institute 31 July 2007 Why Build Accelerators? From Atoms to Quarks Scattering
More informationPoS(ICHEP2012)238. Search for B 0 s µ + µ and other exclusive B decays with the ATLAS detector. Paolo Iengo
Search for B s µ + µ and other exclusive B decays with the ATLAS detector. On behalf of the ATLAS Collaboration INFN Naples, Italy E-mail: paolo.iengo@cern.ch The ATLAS experiment, collecting data in pp
More informationTwo Early Exotic searches with dijet events at ATLAS
ATL-PHYS-PROC-2011-022 01/06/2011 Two Early Exotic searches with dijet events at ATLAS University of Toronto, Department of Physics E-mail: rrezvani@physics.utoronto.ca This document summarises two exotic
More informationElectroweak Physics at the Tevatron
Electroweak Physics at the Tevatron Adam Lyon / Fermilab for the DØ and CDF collaborations 15 th Topical Conference on Hadron Collider Physics June 2004 Outline Importance Methodology Single Boson Measurements
More informationPrecise measurements of the W mass at the Tevatron and indirect constraints on the Higgs mass. Rencontres de Moriond QCD and High Energy Interactions
Precise measurements of the W mass at the evatron and indirect constraints on the Higgs mass Rafael Lopes de Sá for the CDF and DØ Collaborations March 11, 212 Rencontres de Moriond QCD and High Energy
More informationOliver Stelzer-Chilton University of Oxford High Energy Physics Seminar Michigan State University
First Run II Measurement of the W Boson Mass by CDF Oliver Stelzer-Chilton University of Oxford High Energy Physics Seminar Michigan State University April 3 rd, 2007 1. Motivation Outline 2. W Production
More informationThe W-mass Measurement at CDF
2010-05 - 10 The W-mass Measurement at CDF Ilija Bizjak, University College London 1/33 Outline 1) Motivation for a W mass measurement Implications for the EW constraints on Higgs mass 2) Measurement of
More informationBoosted top quarks in the ttbar dilepton channel: optimization of the lepton selection
Boosted top quarks in the ttbar dilepton channel: optimization of the lepton selection DESY Summer School 24 9 September, 24 Author: Ibles Olcina Samblàs* Supervisor: Carmen Diez Pardos Abstract A study
More informationMeasurements of the W Boson Mass and Trilinear Gauge Boson Couplings at the Tevatron
Measurements of the Boson Mass and Trilinear Gauge Boson Couplings at the Tevatron John Ellison University of California, Riverside, USA Selection of and Z events Measurement of the mass Tests of the gauge
More informationTop production measurements using the ATLAS detector at the LHC
Top production measurements using the ATLAS detector at the LHC INFN, Sezione di Bologna and University of Bologna E-mail: romano@bo.infn.it This paper is an overview of recent results on top-quark production
More informationarxiv: v1 [hep-ex] 21 Aug 2011
arxiv:18.155v1 [hep-ex] 1 Aug 011 Early Searches with Jets with the ATLAS Detector at the LHC University of Chicago, Enrico Fermi Institute E-mail: georgios.choudalakis@cern.ch We summarize the analysis
More informationOverview of the Higgs boson property studies at the LHC
Overview of the Higgs boson property studies at the LHC Giada Mancini, a, Roberto Covarelli b a LNF-INFN and University of Roma Tor Vergata, b INFN and University of Torino E-mail: giada.mancini@lnf.infn.it,
More informationMuon reconstruction performance in ATLAS at Run-2
2 Muon reconstruction performance in ATLAS at Run-2 Hannah Herde on behalf of the ATLAS Collaboration Brandeis University (US) E-mail: hannah.herde@cern.ch ATL-PHYS-PROC-205-2 5 October 205 The ATLAS muon
More informationMeasurement of the W-mass with the ATLAS detector
Measurement of the W-mass with the ALAS detector Oleh Kivernyk CEA/Irfu/SPP July, 0 Oleh Kivernyk (CEA/Irfu/SPP) Measurement of the W-mass with the ALAS detector July, 0 / Motivation Relation between M
More informationMeasurement of the Inclusive Isolated Prompt Photon Cross Section at CDF
of the Inclusive Isolated Cross at IFAE Barcelona HEP Seminar University of Virginia Outline Theoretical introduction Prompt photon production The The Photon detection prediction The pqcd NLO prediction
More informationV cb. Determination of. and related results from BABAR. Masahiro Morii, Harvard University on behalf of the BABAR Collaboration
Determination of and related results from BABAR V cb Masahiro Morii, Harvard University on behalf of the BABAR Collaboration V cb from inclusive B semileptonic decays Lepton energy moments Hadron mass
More informationPoS(EPS-HEP2011)250. Search for Higgs to WW (lνlν, lνqq) with the ATLAS Detector. Jonas Strandberg
with the ATLAS Detector Royal Institute of Technology, KTH, Stockholm, Sweden E-mail: jostran@kth.se Higgs boson searches in the H WW ( ) lνlν (l = e, µ) and the H WW ( ) lνqq decay modes, using. fb of
More informationThe Compact Muon Solenoid Experiment. Conference Report. Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland. Rare B decays at CMS
Available on CMS information server CMS CR -2017/115 The Compact Muon Solenoid Experiment Conference Report Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland 17 April 2017 (v4, 10 May 2017) Rare
More informationNational Accelerator Laboratory
Fermi National Accelerator Laboratory FERMILAB-Conf-99/292-E CDF and D0 New W Mass Results from CDF and D0 Bill Carithers Lawrence Berkeley Laboratory Berkeley, California Fermi National Accelerator Laboratory
More informationHiggs and Z τ + τ in CMS
Higgs and Z τ + τ in CMS Christian Veelken for the CMS Collaboration Moriond EWK Conference, March 14 th 2011 Z τ + τ - Production @ 7 TeV τ + Z τ - CMS Measurement of Z/γ* l + l -, l = e/µ: σ BR(Z/γ*
More informationPhysics at Hadron Colliders
Physics at Hadron Colliders Part 2 Standard Model Physics Test of Quantum Chromodynamics - Jet production - W/Z production - Production of Top quarks Precision measurements -W mass - Top-quark mass QCD
More informationTHE FLORIDA STATE UNIVERSITY COLLEGE OF ARTS & SCIENCES W BOSON PRODUCTION CHARGE ASYMMETRY IN THE ELECTRON CHANNEL ASHLEY S HUFF
THE FLORIDA STATE UNIVERSITY COLLEGE OF ARTS & SCIENCES W BOSON PRODUCTION CHARGE ASYMMETRY IN THE ELECTRON CHANNEL By ASHLEY S HUFF A Thesis submitted to the Department of Physics In partial fulfillment
More informationInvariant Mass, Missing Mass, jet reconstruction and jet flavour tagging
1 Experimentelle Methods of Particle Physics HS 215 http://www.physik.uzh.ch/lectures/empp/ Wednesday 16.12.15 and Thursday 17.12.15 Invariant Mass, Missing Mass, jet reconstruction and jet flavour tagging
More informationDiscovery Physics at the Large Hadron Collider
+ / 2 GeV N evt 4 10 3 10 2 10 CMS 2010 Preliminary s=7 TeV -1 L dt = 35 pb R > 0.15 R > 0.20 R > 0.25 R > 0.30 R > 0.35 R > 0.40 R > 0.45 R > 0.50 10 1 100 150 200 250 300 350 400 [GeV] M R Discovery
More informationRecent QCD results from ATLAS
Recent QCD results from ATLAS PASCOS 2013 Vojtech Pleskot Charles University in Prague 21.11.2013 Introduction / Outline Soft QCD: Underlying event in jet events @7TeV (2010 data) Hard double parton interactions
More informationCombined Higgs Results
Chapter 2 Combined Higgs Results This chapter presents the combined ATLAS search for the Standard Model Higgs boson. The analysis has been performed using 4.7 fb of s = 7 TeV data collected in 2, and 5.8
More informationDirect measurement of the W boson production charge asymmetry at CDF
Direct measurement of the boson production charge asymmetry at CDF Eva Halkiadakis Rutgers University For the CDF collaboration Joint Experimental-Theoretical Physics Seminar Fermilab May 22 2009 Outline
More informationThe search for standard model Higgs boson in ZH l + l b b channel at DØ
The search for standard model Higgs boson in ZH l + l b b channel at DØ University of Science and Technology of China E-mail: pengj@fnal.gov on behalf of the DØ collaboration We present a search for the
More informationFERMI NATIONAL ACCELERATOR LABORATORY
FERMI NATIONAL ACCELERATOR LABORATORY arxiv:0908.1374v1 [hep-ex] 10 Aug 2009 TEVEWWG/WZ 2009/01 FERMILAB-TM-2439-E CDF Note 9859 D0 Note 5965 10 th August 2009 Updated Combination of CDF and D0 Results
More informationPoS(ICHEP2012)194. Measurements of the t t forward-backward asymmetry at CDF. Christopher Hays
Measurements of the t t forward-backward asymmetry at CDF On behalf of the CDF Collaboration Oxford University E-mail: hays@physics.ox.ac.uk At next-to-leading order in QCD, the t t forward-backward asymmetry
More informationMeasurement of the W boson mass at Tevatron
Cracow Warsaw Workshop on LHC Faculty of Physics, University of Warsaw 26 March 2010 Measurement of the W boson mass at Tevatron Mikołaj Ćwiok University of Warsaw University of Warsaw Motivation for precise
More informationarxiv: v1 [hep-ex] 8 Jan 2018
Top mass measurements in ATLAS and CMS arxiv:8.239v [hep-ex] 8 Jan 28 A. CASTRO University of Bologna and INFN, Bologna, Italy On behalf of the ATLAS and CMS Collaborations Abstract Top quarks are produced
More informationAn introduction to Bayesian reasoning in particle physics
An introduction to Bayesian reasoning in particle physics Graduiertenkolleg seminar, May 15th 2013 Overview Overview A concrete example Scientific reasoning Probability and the Bayesian interpretation
More informationSearch for H ± and H ±± to other states than τ had ν in ATLAS
Search for H ± and H to other states than τ had ν in On behalf of the collaboration Louisiana Tech University E-mail: Catrin.Bernius@cern.ch esults of searches for charged Higgs bosons based on data from
More informationDirect Measurement of the W Total Decay Width at DØ. Outline
Direct Measurement of the W Total Decay Width at DØ Introduction Junjie Zhu University of Maryland On behalf of the DØ Collaboration Monte Carlo Simulation Event Selection Outline Determination of the
More informationTop quark pair properties in the production and decays of t t events at ATLAS
ATL-PHYS-PROC-214-72 11 July 214 Top quark pair properties in the production and decays of t t events at DESY, Hamburg Universität Wuppertal E-mail: ralph.schaefer@cern.ch In proton-proton collisions at
More informationRelative branching ratio measurements of charmless B ± decays to three hadrons
LHCb-CONF-011-059 November 10, 011 Relative branching ratio measurements of charmless B ± decays to three hadrons The LHCb Collaboration 1 LHCb-CONF-011-059 10/11/011 Abstract With an integrated luminosity
More informationNicolas Berger SLAC, Menlo Park, California, U.S.A.
SLAC-PUB-11414 August 25 Frascati Physics Series Vol. VVVVVV (xxxx), pp. - DAΦNE 24: Physics at meson factories Frascati, June. 7-11, 24 Selected Contribution in Plenary Session INCLUSIVE HADRONIC RESULTS
More informationMeasurement of t-channel single top quark production in pp collisions
Measurement of t-channel single top quark production in pp collisions (on behalf of the CMS collaboration) INFN-Napoli & Università della Basilicata E-mail: Francesco.Fabozzi@cern.ch Measurements of t-channel
More informationLHCb Semileptonic Asymmetry
CERN E-mail: mika.vesterinen@cern.ch A recent measurement of the CP violating flavour specific asymmetry in B s decays, a s sl, is presented. This measurement is based on a data sample corresponding to
More informationSearches for Exotics in Upsilon Decays in BABAR
Searches for Exotics in Upsilon Decays in BABAR Yury Kolomensky UC Berkeley/LBNL For the BABAR Collaboration 14th Lomonosov Conference on Elementary Particle Physics August 24, 2009, Moscow Light Higgs
More informationDiscovery Potential for the Standard Model Higgs at ATLAS
IL NUOVO CIMENTO Vol.?, N.?? Discovery Potential for the Standard Model Higgs at Glen Cowan (on behalf of the Collaboration) Physics Department, Royal Holloway, University of London, Egham, Surrey TW EX,
More informationDiscovery potential of the SM Higgs with ATLAS
Discovery potential of the SM Higgs with P. Fleischmann On behalf of the Collaboration st October Abstract The discovery potential of the Standard Model Higgs boson with the experiment at the Large Hadron
More informationHighlights of top quark measurements in hadronic final states at ATLAS
Highlights of top quark measurements in hadronic final states at ATLAS Serena Palazzo 1,2,, on behalf of the ATLAS Collaboration 1 Università della Calabria 2 INFN Cosenza Abstract. Measurements of inclusive
More informationPrecision measurements of the top quark mass and width with the DØ detector
FERMILAB-CONF-10-588-E Precision measurements of the top quark mass and width with the DØ detector for the DØ collaboration Irfu/SPP/CEA Saclay E-mail: agrohsje@fnal.gov Since the discovery of the top
More informationResults on QCD and Heavy Flavors Production at the Tevatron
Results on QCD and Heavy Flavors Production at the Tevatron Donatella Lucchesi INFN and University of Padova October 21 Padova Outline: Machine and detector description Latest results on QCD Heavy Flavor:
More informationStudy of di-muon production with the ZEUS detector at HERA
Submitted to the International Conference on High Energy Physics 16-22 August 2004, Beijing, China Abstract: EW 4-0231, BSM 12-0230 Sessions: Electro Weak, Beyond the Standard Model Study of di-muon production
More informationCMS Event Simulation
CMS Event Simulation Nicole A. Larsen Department of Physics, Georgia Institute of Technology, Atlanta, GA, 33 (Dated: August 11, 6) The CMS Detector located at the Large Hadron Collider at CERN recently
More information14 Top Quark. Completing the Third Generation
14 Top Quark Completing the Third Generation No one could doubt that there would be a sixth quark, the top or t, but it was equally certain that initially no one knew where it would be found. With the
More informationStatistics and Data Analysis
Statistics and Data Analysis The Crash Course Physics 226, Fall 2013 "There are three kinds of lies: lies, damned lies, and statistics. Mark Twain, allegedly after Benjamin Disraeli Statistics and Data
More informationarxiv: v1 [hep-ex] 18 Jul 2016
Top quark mass measurements with the experiment at the LHC arxiv:67.4972v [hep-ex] 8 Jul 26 Deutsches Elektronen-Synchrotron DESY E-mail: simon.spannagel@desy.de Measurements of the top quark mass are
More informationHow to find a Higgs boson. Jonathan Hays QMUL 12 th October 2012
How to find a Higgs boson Jonathan Hays QMUL 12 th October 2012 Outline Introducing the scalar boson Experimental overview Where and how to search Higgs properties Prospects and summary 12/10/2012 2 The
More informationAnalysis of Top Quarks Using a Kinematic Likelihood Method
The Ohio State University Undergraduate Research Thesis Presented in Partial Fulfillment of the Requirements for Graduation with Research Distinction in Physics in the Undergraduate Colleges of The Ohio
More informationIntroduction. The Standard Model
Ph.D. Thesis in Engineering Physics Supervisor: Assoc. Prof. Dr. Ayda BEDDALL Co-supervisor: Assist. Prof. Dr. Andrew BEDDALL By Ahmet BNGÜL Page 1 Introduction Chapter 1-2 High energy collisions of sub-atomic
More informationSpectroscopy and Decay results from CDF
Quarkonium Spectroscopy and Decay results from CDF KIT Quarkonium Workshop December 3, 2008 Outline Tevatron and CDF Bc Mass Lifetime bc X(3872) Mass splitting and mass X page 2 Tevatron p s = 1.96 TeV
More informationLuminosity measurement and K-short production with first LHCb data. Sophie Redford University of Oxford for the LHCb collaboration
Luminosity measurement and K-short production with first LHCb data Sophie Redford University of Oxford for the LHCb collaboration 1 Introduction Measurement of the prompt Ks production Using data collected
More informationFirst V+jets results with CMS. Vitaliano Ciulli (Univ. & INFN Firenze) V+jets workshop, 8-10 Sep 2010, Durham
First V+jets results with CMS Vitaliano Ciulli (Univ. & INFN Firenze) V+jets workshop, 8- Sep 20, Durham Outline Jet and missing energy reconstruction W and Z selections and cross-section measurement First
More informationRecent results on search for new physics at BaBar
Recent results on search for new physics at BaBar B. Oberhof* *LNF-INFN, Italy on behalf of the BaBar Collaboration International Symposium Advances in Dark Matter and Particle Physics Messina, Italy 25th
More informationStatistics. Lent Term 2015 Prof. Mark Thomson. 2: The Gaussian Limit
Statistics Lent Term 2015 Prof. Mark Thomson Lecture 2 : The Gaussian Limit Prof. M.A. Thomson Lent Term 2015 29 Lecture Lecture Lecture Lecture 1: Back to basics Introduction, Probability distribution
More informationSingle top quark production at CDF
Available online at www.sciencedirect.com Nuclear and Particle Physics Proceedings 73 75 (16) 74 79 www.elsevier.com/locate/nppp Single top quark production at CDF Sandra Leone, on behalf of the CDF Collaboration
More informationPoS(Confinement X)171
in pp collisions at LHC energies Institute for Particle and Nuclear Physics, Wigner RCP, HAS 2 Budapest, XII. Konkoly Thege Miklos ut 29-33 E-mail: sona.pochybova@cern.ch With increasing luminosities at
More informationHiggs couplings and mass measurements with ATLAS. Krisztian Peters CERN On behalf of the ATLAS Collaboration
Higgs couplings and mass measurements with ATLAS CERN On behalf of the ATLAS Collaboration July observation: qualitative picture A single state observed around ~125 GeV Qualitatively all observations consistent
More informationMeasurement of Фs, ΔΓs and Lifetime in Bs J/ψ Φ at ATLAS and CMS
Measurement of Фs, ΔΓs and Lifetime in Bs J/ψ Φ at ATLAS and CMS Claudio Heller Excellence Cluster Universe LMU München For the ATLAS and CMS Collaborations Beauty 2013 14th International Conference on
More informationResults from D0: dijet angular distributions, dijet mass cross section and dijet azimuthal decorrelations
Results from D: dijet angular distributions, dijet mass cross section and dijet azimuthal decorrelations Zdenek Hubacek Czech Technical University in Prague E-mail: zdenek.hubacek@cern.ch on behalf of
More informationThe achievements of the CERN proton antiproton collider
The achievements of the CERN proton antiproton collider Luigi DiLella Scuola Normale Superiore, Pisa, Italy Motivation of the project The proton antiproton collider UA1 and UA2 detectors Discovery of the
More informationATLAS and CMS Higgs Mass and Spin/CP Measurement
ATLAS and CMS Higgs Mass and Spin/CP Measurement These proceedings summarize the ATLAS and CMS collaboration measurement of Higgs properties. The properties discussed include mass, spin, parity and the
More informationVBF SM Higgs boson searches with ATLAS
VBF SM Higgs boson searches with Stefania Xella (for the collaboration) Niels Bohr Institute, Copenhagen University, Denmark E-mail: xella@nbi.dk The observation of a Standard Model Higgs boson produced
More informationBeauty jet production at HERA with ZEUS
DESY and University of Hamburg E-mail: mlisovyi@mail.desy.de wo recent measurements of beauty production in deep inelastic ep scattering with the detector are presented. In the first analysis, beauty production
More informationSearch for Pair Production of Stop Quarks Mimicking Top Event Signatures
CDF/PUB/TOP/PUBLIC/9343 Search for Pair Production of Stop Quarks Mimicking Top Event Signatures The CDF Collaboration URL http://www-cdf.fnal.gov May 29, 2008 Abstract We present the search for the pair-produced
More informationMeasuring the Top Quark Mass using Kinematic Endpoints
Journal of Physics: Conference Series OPEN ACCESS Measuring the Top Quark Mass using Kinematic Endpoints To cite this article: Benjamin Nachman and the Cms Collaboration 13 J. Phys.: Conf. Ser. 45 156
More informationW mass measurement in the ATLAS experiment
W mass measurement in the ATLAS experiment DSM/Irfu/SPP, CEA/Saclay 99 Gif sur Yvette Cedex France E-mail: Nathalie.Besson@cern.ch A precise measurement of the mass of the W boson will be essential to
More informationMeasurement of charged particle spectra in pp collisions at CMS
Measurement of arged particle spectra in pp collisions at CMS for the CMS Collaboration Eötvös Loránd University, Budapest, Hungary E-mail: krisztian.krajczar@cern. We present the plans of the CMS collaboration
More informationNegative binomial distribution and multiplicities in p p( p) collisions
Negative binomial distribution and multiplicities in p p( p) collisions Institute of Theoretical Physics University of Wroc law Zakopane June 12, 2011 Summary s are performed for the hypothesis that charged-particle
More informationElectroweak results. Luca Lista. INFN - Napoli. LHC Physics
Electroweak results Luca Lista INFN - Napoli EWK processes at LHC p p W and Z production in pp collisions proceeds mainly form the scattering of a valence quark with a sea anti-quark The involved parton
More informationThesis. Wei Tang. 1 Abstract 3. 3 Experiment Background The Large Hadron Collider The ATLAS Detector... 4
Thesis Wei Tang Contents 1 Abstract 3 2 Introduction 3 3 Experiment Background 4 3.1 The Large Hadron Collider........................... 4 3.2 The ATLAS Detector.............................. 4 4 Search
More informationEvidence for Single Top Quark Production. Reinhard Schwienhorst
Evidence for Single Top Quark Production Reinhard Schwienhorst MSU high energy seminar, 1/9/2007 Outline Motivation Preparation Optimized Event Analysis Sensitivity Cross section measurement Vtb Conclusions/Outlook
More informationDiscovery and Significance. M. Witherell 5/10/12
Discovery and Significance M. Witherell 5/10/12 Discovering particles Much of what we know about physics at the most fundamental scale comes from discovering particles. We discovered these particles by
More informationQCD and jets physics at the LHC with CMS during the first year of data taking. Pavel Demin UCL/FYNU Louvain-la-Neuve
QCD and jets physics at the LHC with CMS during the first year of data taking Pavel Demin UCL/FYNU Louvain-la-Neuve February 8, 2006 Bon appétit! February 8, 2006 Pavel Demin UCL/FYNU 1 Why this seminar?
More informationStandard Model of Particle Physics SS 2012
Lecture: Standard Model of Particle Physics Heidelberg SS 2012 W- and Z-Bosons 1 2 Contents Discovery of real W- and Z-bosons Intermezzo: QCD at Hadron Colliders LEP + Detectors W- and Z- Physics at LEP
More informationDilepton Forward-Backward Asymmetry and electroweak mixing angle at ATLAS and CMS
Dilepton Forward-Backward Asymmetry and electroweak mixing angle at ATLAS and CMS University of Rochester E-mail: jyhan@fnal.gov We present measurements of the forward-backward asymmetry in Drell-Yan dilepton
More informationEvents with High P T Leptons and Missing P T and Anomalous Top at HERA
with High Leptons and Missing and Anomalous Top at HERA David South (DESY) On Behalf of the H Collaboration IIth International Workshop on Deep Inelastic Scattering (DIS 24) Štrbské Pleso, High Tatras,
More informationTutorial on Top-Quark Physics
Helmholtz Alliance at the Terascale Data Analysis Group Introductory School on Terascale Physics 21 25 February, 2011 Tutorial on Top-Quark Physics Introduction to the Tevatron, the CDF Detector and Top-Quark
More informationTop production at the Tevatron
Top production at the Tevatron Bob Hirosky University of Virginia for the CDF and D0 Collaborations Moriond Moriond QCD 2017 1 Top production at the Tevatron Top mass from differential cross section Polarization
More informationStudy of Higgs Boson Decaying to Four Muons at s =14 TeV
Study of Higgs Boson Decaying to Four Muons at s =14 TeV R.M. Aly 1, A.A. Abdelalim 1,2, M.N.El-Bakrey 1 and A. Mahrous 1 1 Department of physics, Faculty of science, Helwan University, Cairo, Egypt. 2
More informationV 0 production studies at LHCb. Mathias Knecht, EPFL , joint SPS-ÖPG-ÖGA 2 meeting, Innsbrück, Österreich, September 2-4, 2009
V 0 production studies at LHCb Mathias Knecht, EPFL 2 9 2009, joint SPS-ÖPG-ÖGA 2 meeting, Innsbrück, Österreich, September 2-4, 2009 Outline One of the first measurements to be done with the LHCb detector...
More informationThe Compact Muon Solenoid Experiment. Conference Report. Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland
Available on CMS information server CMS CR -2013/016 The Compact Muon Solenoid Experiment Conference Report Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland 18 January 2013 (v2, 21 January 2013)
More informationTop, electroweak and recent results from CDF and combinations from the Tevatron
IL NUOVO CIMENTO 40 C (2017) 181 DOI 10.1393/ncc/i20177181-7 Colloquia: LaThuile 2017 Top, electroweak and recent results from CDF and combinations from the Tevatron D. Lucchesi( 1 )( 2 ) on behalf of
More informationE. Santovetti lesson 4 Maximum likelihood Interval estimation
E. Santovetti lesson 4 Maximum likelihood Interval estimation 1 Extended Maximum Likelihood Sometimes the number of total events measurements of the experiment n is not fixed, but, for example, is a Poisson
More informationThe ATLAS Detector at the LHC
The ATLAS Detector at the LHC Results from the New Energy Frontier Cristina Oropeza Barrera Experimental Particle Physics University of Glasgow Somewhere near the Swiss Alps... A Toroidal LHC ApparatuS
More informationSearch for the very rare decays B s/d + - at LHCb
Search for the very rare decays B s/d + - at LHCb Justine Serrano Centre de Physique des Particules de Marseille On behalf of the LHCb collaboration EPS-HEP, July 22 th 2011, Grenoble Bs2MuMu @ LHCb Justine
More informationEarly SUSY searches at the LHC
on behalf of the ATLAS and CMS Collaborations Imperial College London E-mail: a.tapper@imperial.ac.uk Supersymmetry may give rise to striking events that could be discovered early in LHC running. Search
More information1 Measurements of sin 2 θ eff. (M z ), sin 2 θ w and indirect measurement of M w at the Tevatron (and Prospects for the LHC)
1 Measurements of sin 2 θ eff lept (M z ), sin 2 θ w and indirect measurement of M w at the Tevatron (and Prospects for the LHC) Arie Bodek, University of Rochester, Wed. March 15, 2017 11:00 AM BASIC
More informationPoS(CKM2016)117. Recent inclusive tt cross section measurements. Aruna Kumar Nayak
Recent inclusive tt cross section measurements Institute of Physics, Bhubaneswar, India E-mail: Aruna.Nayak@cern.ch Results of the recent measurements for the inclusive tt production cross section in the
More informationHiggs Searches at CMS
Higgs Searches at CMS Ashok Kumar Department of Physics and Astrophysics University of Delhi 110007 Delhi, India 1 Introduction A search for the Higgs boson in the Standard Model (SM) and the Beyond Standard
More informationDetermination of Electroweak Parameters
Determination of Electroweak Parameters Seminar Particle Physics at the LHC Proceedings R. Gugel Freiburg, 27.05.2014 Contents 1 Motivation 1 2 Template Method 4 3 Mass of the W Boson 5 3.1 Measurement
More informationWhat to do with Multijet Events?
DØ note 4465 2004-06-03 What to do with Multijet Events? T. Hebbeker RWTH Aachen ABSTRACT Multijet events are frequently produced in p p collisions. Events with several jets might indicate physics beyond
More informationMulti Lepton events at HERA
St. Petersburg, 25/04/2003 DIS03 Conference Multi Lepton events at HERA Andrea Parenti (on behalf of H and ZEUS Collabs.) Padova University and INFN A.Parenti - Multi Lepton Events at HERA p./?? Outline
More information