Determining Momentum and Energy Corrections for g1c Using Kinematic Fitting

Size: px
Start display at page:

Download "Determining Momentum and Energy Corrections for g1c Using Kinematic Fitting"

Transcription

1 CLAS-NOTE 4-17 Determining Momentum and Energy Corrections for g1c Using Kinematic Fitting Mike Williams, Doug Alegate and Curtis A. Meyer Carnegie Mellon University June 7, 24 Abstract We have used the CLAS Kinematic Fit and the exclusive reaction γ π + π to derive energy and momentum corrections for the g1c run eriod. By selectively leaving articles out of the fit, and then comaring fit and measured quantities, we have been able to derive corrections for the hoton energy and the final state articles. Contents 1 Introduction 2 2 The Procedure The π Sectrum The π + sectrum The sectrum The Photon Beam Energy Comarison of Results Mass eaks for known resonances Summary and Conclusions 18 A π Corrections 19 B π + Corrections 27 C Alying the π + corrections to the roton 3 1

2 1 Introduction In late 23, a Kinematic Fitting ackage for the CLAS exeriments was released by the authors. The details of this ackage can be found in an earlier CLAS-Note [1], and background information on the mathematics of the fitting rocedure are described in a book by Brandt [2]. The urose of this study was to understand what if any momentum and energy corrections should be alied to the data from g1c. 2 The Procedure The γ π + π data set is an excellent environment in which to study the CLAS detector. The reaction is fully reconstructed, which in turn allows one to carry out a 4-constraint (4-C) fit to the data. The ull [3] and confidence-level distributions from this data can then oint one in the correct direction to small roblems in the data set. These small systematic errors can then be studied in detail by selectively leaving out articles in a series of 1-C fits erformed on the data. This study started by examining the momentum resolution of the π by using the reaction: γ π + (π ) missing. Treating the π as an undetected article in a 1-C fit yields a momentum, π (fit), that is not biased by the measured value. The difference, π = π (fit) π (meas) (1) can then be studied as a function of ion angles, momentum and sector number. Following this work, the π is corrected, and one can then study the correction to the π + momentum. Because of the asymmetric bending of ositive and negative articles, these two corrections are not exected to be the same. With this comlete, the roton is examined. While we nominally exect the roton and the π + corrections to be the same, there may be residual issues in the energy loss corrections of rotons versus ions that introduce small differences. All energy loss corrections are erformed using the ELOSS [6] ackage and are erformed before any fits or comarisons. Finally, we can look at the deviation in the tagged hoton energy. This rocess can then be iterated tbtain the best set of corrections. We note here that the changes to the ELOSS corrections along with the corrections to the tagged hoton energy were found to be significant enough to warrant restarting the rocess with first aroximations for these already inserted. The results resented in the following sections were obtained from the final iteration. As a reference oint, we resent the confidence-level and ull distributions obtained by erforming a 4-C fit on the data rior to making any of the corrections resented in this aer (standard ELOSS corrections have been made). Figure 1 shows the confidence level distribution, while Figure 2 shows the ull distributions, the means and widths of the ull distributions are given in Table 1. Both the confidence level and ull distributions look quite reasonable, while it could be argued that some of the ulls are skewed, they are within accetable limits for carrying out kinematic fits. 2

3 C Fit No Corrections Figure 1: The confidence level distribution for the 4-C fit to γ π + π before any momentum corrections have been made(standard ELOSS corrections have been made). Quantity Mean Sigma λ φ π λ π φ π π λ π φ π E γ Table 1: The means and widths of the ull distributions, shown in Figure 2, before momentum corrections have been alied to the data. 3

4 2 x1 2 x1 2 x λ 12 1 φ x x x π + 1 λ π + 1 φ π x π - 2 x λ π - 2 x φ π x E γ Figure 2: The ull distributions for the 4-C fit to γ π + π before any momentum corrections have been made. The ulls have been fit with a Gaussian function which is overlaid as the solid line. The means and σ s are given in Table 1. 4

5 2.1 The π Sectrum When the π is treated as a missing article in a 1-C fit, we can form the quantity π (see Eq. 1). To get an idea of where roblems might be, we made 2-D histograms of vs., θ and φ, then fit each bin to a gaussian and lotted the means. Figure 3 shows the mean of this difference as a function of the π momentum, while Figures 4 and show it as a function of θ and φ resectively. One striking feature seen in all of these lots is the large offset in sector 2 relative to the other sectors in CLAS. Figure 4 shows that this aears to be concentrated at forward angles (θ), and that the deviation can get as large as 3 MeV/c. All other sectors show a significantly smaller deviation. This effect was actually first observed by Kim [4] and is attributed to an alignment issue in the Region-3 drift chambers. Clearly any corrections that we make are going to have to be sector deendent. It is also clear that the corrections have an aarent comlicated deendence on θ, φ and, something noted by Taylor []. Rather than trying to correct the momentum in terms of its magnitude,, we chose to use the, where Bdl B d l is the total amount of magnetic field the article travels through in CLAS. This deendence was seen in earlier work by Kim [4]. The quantity Bdl is comuted by crudely swimming the articles through the CLAS magnetic field. Figure 6 shows a lot of the result of several of these tracks. quantity which is actually measured, Bdl We can then correct Bdl by dividing each sector u into 1 φ bins giving us a correction of the form, Bdl Bdl = α ij + β ij Bdl = a i + b i θ + c i θ 2 + d i θ 3 + e i θ 4 (3) (2) where i = 1, 2,.., 6 is the sector number and j = 1, 2,...1 is the φ bin. The arameters α, β, a, b, c, d and e are constants which are fit to the available data. Terms where Bdl deends only on θ, φ and sector number are due to drift chamber misalignments. Thus, Equation (3) can be interreted as the θ-deendence of the misalignments, while α ij in Equation (2) can be interreted as the φ-deendence Bdl of the misalignments averaged over all θ. The term β ij in Equation (2) accounts for roblems in calculating B d l that could be caused either by inaccuracies in the CLAS magnetic field ma or by drift chamber misalignments leading to systematic shifts in the articles trajectories with resect to the field ma. Figures in Aendix A show the results of lotting Bdl, dividing each sector into 1 φ bins, for the π. Each bin is, to first order, linear in Bdl Bdl vs. as exected. We fit each bin to Equation(2), then aly this correction and lot Bdl vs. θ for each sector (see Figure 24 in Aendix A). We do not correct the measured angles θ and φ. We found the systematic shifts in these angles to be less than 1mrad, which is not enough to ose any roblems during the kinematic fitting rocess. It is imortant to kee in mind that EXCEPT for some forward angle regions in sector 2, the tyical size of these momentum corrections are a few MeV. These are NOT large corrections. Figure 17 shows the corrections to the π momentum as a function of, θ and φ resectively summed over all sectors.

6 Sector 1 Sector Sector 3 Sector Sector Sector Figure 3: vs. : The change in momentum of the π as a function of the π momentum for each of the six sectors in the CLAS detector. 6

7 Sector 1 Sector Sector 3 Sector Sector Sector Figure 4: vs. θ: The change in momentum of the π as a function of the π θ angle for each of the six sectors in the CLAS detector. 7

8 Sector 1 Sector Sector 3 Sector Sector Sector Figure : vs. φ: The change in momentum of the π as a function of the π φ angle for each of the six sectors in the CLAS detector. 8

9 3 2 q = 1 P = 1 GeV/c o Θ = q = 1 P = MeV/c o Θ = q = -1 P = 1 GeV/c o Θ = q = -1 P = MeV/c o Θ = Figure 6: Outut of the mini tracking to comute the B d l along a article trajectory through the magnetic field. 9

10 2.2 The π + sectrum We followed the same rocedure for the π + as we did for the π. We treated the π + as the missing article in a 1-C fit and used the quantity Bdl to determine the corrections to the magnitude of the π + momentum. Figures 26-31, in Aendix B show the results of lotting Bdl Bdl vs., dividing each sector into 1 φ bins, for the π +. Each bin is, to first order, linear in Bdl just as for the π. We fit each bin to Equation(2), then aly this correction and lot Bdl vs. θ for each sector (see Figure 32 in Aendix B). Again, excet for some forward angle regions in sector 2, the tyical size of these momentum corrections are a few MeV. Figure 2 shows the corrections to the π + momentum as a function of, θ and φ resectively summed over all 6 sectors. 2.3 The sectrum If all of the roblems in the measured momentum are due to either magnetic field or alignment issues in CLAS, then one would exect that the corrections for all ositive articles (π +, K + and ) would be the same, and all corrections for negative articles, (π and K ) will be the same. Initially, we were seeing differences between and π +. However, we soon realized that it aeared the ELOSS corrections were still somewhat too small. In Figure 7 we show the ELOSS correction ( ) as a function of the article momentum for the roton. Figure 8 shows the correction to ELOSS needed on the roton as a function of the size of the alied correction. The corrections needed are linear with sector deendent sloes corresonding to about 18%. We attribute this additional correction to energy loss in traversing the tracking region of CLAS. Note that this is only really relevant for low momentum (< 3 MeV/c) rotons. If we aly these sector deendent scale factors to the ELOSS ackage corrections, then the corrections that we derived for the π + in the revious section also work for the rotons. Thus, our corrections deend only on which sector in CLAS the article traveled through and the article s charge. This is consistent with our icture of how these corrections should work. (GeV) P eloss P measured (GeV) Figure 7: Eloss Correction vs. : Correction alied to the magnitude of the roton momenta by the ELOSS ackage. 1

11 ..4.3 Sector Sector Sector Sector Sector..4.3 Sector Figure 8: P vs. Eloss Correction : Correction needed for the roton momentum, after the π + momentum corrections and standard ELOSS corrections have been made, versus the energy loss correction obtained from the ELOSS ackage. 11

12 Sector ELOSS Scale Factor Table 2:. Scale Factors this analysis shows need to be alied to the energy loss corrections obtained from the ELOSS ackage. 2.4 The Photon Beam Energy In the final round of analysis, all final state articles in the detector have had their momentum corrected, and we then want to look at the measured hoton energy. Because it is not trivial to modify the kinematic fit to leave the hoton energy free, we took a somewhat different aroach here. We have used the target mass and the detected final state articles to redict what the hoton beam energy, E γ(fit), should be and then comared it to the measured energy, E γ(measured). So E γ(fit) is defined as, E γ(fit) = E + E π + + E π m. (4) We can then form the difference: E γ = E γ(fit) E γ(measured). () Tbtain a correction for the hoton energy, we first histogrammed E γ vs. E γ. We then fit each bin of this histogram to a Gaussian then lotted the mean versus the hoton energy. Figure 9 shows this histogram fitted to a iecewise quadratic function in E γ. It is interesting that this function is not constant. The double hum structure is similar in shae to the results of Steanyan [7]. (As an aside, there is an absolute offset between the values in this study and those of Steanyan). Since this observation, the shift can be exlained by a calculable sage in the tagger [9] E γ E γ (GeV) Figure 9: Correction to the tagged hoton beam energy as a function of the beam energy. 12

13 3 Comarison of Results After all of the corrections described in this note have been erformed, the resulting confidence-level is shown in Figure 1. Clearly, some of the errors going into the fit were slightly too large. A result that should be exected from the fact that we have corrected 4 of the 1 measured quantities. In our earlier work [1], the errors on the energy-loss were fit to the data. After this study, we refit those and come u with numbers that are about 6% of the earlier error. After we have adjusted the errors to comensate, the resulting confidence level distribution is shown in Figure 11 and the ull distributions are shown in Figure 12. The means and σ s of the ull distributions are shown in Table 3. It should be noted that the error that we adjusted is only one of several comonents that feed into the total error. The effect on the total error is much smaller than imlied by the 6% number, but it is the one degree of freedom that is not a-riori totally constrained by hysics C Fit With Corrections Figure 1: The confidence level distribution for the 4-C fit to γ π + π after all momentum corrections have been made. The rise near 1 demonstrates a need to scale down our errors slightly now that we ve imroved the data being inut to the fit C Fit With Corrections With Adjusted Errors Figure 11: The confidence level distribution for the 4-C fit to γ π + π after all momentum corrections have been made and the errors have been adjusted. 13

14 λ 12 1 φ π λ π φ π π λ π φ π E γ Figure 12: The ull distributions for the 4-C fit to γ π + π after all momentum corrections have been made. 14

15 Quantity Mean(initial) Mean (final) Sigma(initial) Sigma(final) λ φ π λ π φ π π λ π φ π E γ Table 3: The means and widths of the ull distributions before and after all momentum corrections have been alied to the data. 1

16 3.1 Mass eaks for known resonances In order to demonstrate that the kinematic fit is robust against these small errors in the momentum and energy, we have examined the reaction: γ π + π X miss. (6) where the missing article can be either a π or an η. In the case of the π, there are two three ion resonances of interest, η π + π π and ω π + π π. In the case of the missing η, the main resonance of interest is the η ηπ + π. We select the data for dislay by lacing a 1% confidence level cut on the 1C fit to the two channels. In Figure 13, we show the three mass eaks as roduced using unfit and uncorrected momentum. This is effectively the missing mass off the γ X miss sectrum. The eaks have been fit using a Voigtian function [8] with the natural width of the resonance as an inut arameter. Figure 14 shows the same sectrum after all momentum corrections have been made, but again without fitting. There is significant imrovement in the accuracy of masses and a slight inrovement in the widths. In Figure 1 we show the kinematically fit sectrum starting from the uncorrected momentum, while in Figure 16 we show the fits starting from the corrected momentum. In Table 4 we summarize the results of all of these fits. It should be ointed out that the value of the fit mean can vary by about ±.2 MeV/c 2 due to assumtions about background shae and variations in fit ranges. Excet for the case of uncorrected measured data, all of the masses are extremely accurate. The slight imrovements in the width from the momentum corrections translates into a narrower resolution in the fit data as well. Particle Mass [1] Uncor., Measured Cor., Measured Uncor., Fit Cor., Fit Mean σ Mean σ Mean σ Mean σ η ω η Table 4: The mass and widths of known resonances fit using a Voigtian function to derive the mass and resolution width. The four blocks of data corresond to uncorrected momentum using measured values, corrected measured momentum, uncorrected momentum that has been kinematically fit and corrected momentum that has been kinematically fit. All quantities are given in MeV/c Measured M = 1.9 MeV σ = MeV M(π + π - π ) GeV Measured M = MeV σ = MeV M(π + π - π ) GeV Measured M = 96.2 MeV σ =.81 MeV M(π + π - η) GeV Figure 13: The invariant mass sectrum for the η (left), ω (center) and η (right) made using the unfit and uncorrected momentum. 16

17 Measured with Corrections M = 48.4 MeV σ = MeV M(π + π - π ) GeV Measured with Corrections M = MeV σ = 6.32 MeV M(π + π - π ) GeV Measured with Corrections M = 97.2 MeV σ =.88 MeV M(π + π - η) GeV Figure 14: The invariant mass sectrum for the η (left), ω (center) and η (right) made using the unfit but corrected momentum. 22 Kinematically Fit 2 M = 48.7 MeV 18 σ = 4.9 MeV M(π + π - π ) GeV 2 Kinematically Fit 18 M = MeV 16 σ =.869 MeV M(π + π - π ) GeV Kinematically Fit M = 96.6 MeV σ = 4.18 MeV M(π + π - η) GeV Figure 1: The invariant mass sectrum for the η (left), ω (center) and η (right) made using the kinematically-fit uncorrected momentum Kinematically Fit with Corrections M = 47.8 MeV σ = 3.21 MeV M(π + π - π ) GeV 2 Kinematically Fit with 18 Corrections 16 M = MeV 14 σ =.62 MeV M(π + π - π ) GeV Kinematically Fit with Corrections M = 97.3 MeV σ = 4.7 MeV M(π + π - η) GeV Figure 16: The invariant mass sectrum for the η (left), ω (center) and η (right) made using the kinematically-fit and corrected momentum. 17

18 4 Summary and Conclusions We have carried out a detailed study of momentum and energy corrections for the CLAS g1c data. The result of this work is three-fold. First, the hoton beam energy needs to be corrected by the interesting function shown in Figure 9. If the sloe and offset are removed from this function, the wiggles are almost certainly within the design goals of the tagger. However, using this function, we can actually imrove uon that slightly. Second, we believe that the energy loss of the articles as they travel through the tracking region is not accurately handled in CLAS. We have identified a sector deendent effect that is on average about 18% of the ELOSS correction that needs to be added. This correction is needed if we want the residual momentum corrections for all articles of the same sign to be identical. Finally, we have develoed residual momentum corrections for articles that deend only on the sector in CLAS and the charge of the articles. All of these functions will be checked into the CLAS CVS reository. Alying these corrections will rovide imroved reconstructed quantities in CLAS data. We also oint out that even before these corrections, the error functions worked out in our kinematic fit study [1] were of high enough quality that a good kinematic fit could be made. While slight imrovements are achieved by alying the corrections before using a fit, they are NOT necessary tbtain good results kinematically fitting the data. References [1] Mike Williams and Curtis A. Meyer, Kinematic Fitting in CLAS, CLAS-Note 3-17, November 23. [2] Sigmund Brandt, Statistical and Comutational Methods in Data Analysis, North Holland Publishing Comany, (197), (see Chater 9). [3] A. G. Frodesen and O. Skjeggestad, Probability and Statistics in Particle Physics, Columbia University Press, (1979), (see Chater 1 on stretch functions). [4] K. Y. Kim, H. Denizli, J. A. Mueller and S. A. Dytman, General momentum correction form for CLAS: A study of momentum corrections for the Prod 1-9 cooking of e1c data, CLAS-Note 1-18, October 21. [] S. Taylor and G. Mutchler, Analysis udate for E89-24 (radiative decays of excitedstate hyerons), CLAS Analysis Note 23-18, February 24. [6] Eugene Pasyuk, The CLAS ELOSS Package. [7] Stean Steanyan, Evidence for an Exotic Pentaquark Baryon State Θ + (13), in the γd K + K n reaction, CLAS Analysis Note 23-1, July 23. [8] Curtis A. Meyer, Mike Williams and Robert Bradford, Comuting Invariant Masses and Missing Masses, CLAS-Note 4-8, February 23. [9] Dan Sober, rivate communication. [1] The 22 version of the PDG. 18

19 A π Corrections This aendix shows detailed lots of the π momentum corrections. Figures 18, 19,2,21,22 and 23 show the quantity Bdl versus Bdl in sectors 1-6. The lots are broken into ten bins in φ in each sector. Figure 17 shows the corections alied to data from one run in g1c summed over all 6 sectors. These corrections are small, tyically only a few MeV. One interesting feature of these lots is that the sloe decreases as φ increases for each sector. We interret this as a rotation of the drift chambers with resect to the torus coils. Thus, one side of each sector is closer to the coils while the other side is farther away than what is inut to the tracking code. This results in a reorted Bdl that is too large on the side of the sector rotated away from the coils and too small on the side of the sector rotated towards the coils. While the decreasing sloe with resect to φ is common to all sectors, the actual values of the sloes are different in each sector. One interretation for this effect is that the current isn t the same in all six torus coils. In fact, to cause the affect that we see in the data the differences in the current would have to be.2%, or a few ams (GeV) θ φ Figure 17: Corrections to the π momentum as a function of, θ and φ summed over all 6 sectors. 19

20 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 18: Bdl Bdl vs. : The change in Bdl, for the π Bdl, in Sector 1. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 2

21 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 19: Bdl Bdl vs. : The change in Bdl, for the π Bdl, in Sector 2. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 21

22 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 2: Bdl Bdl vs. : The change in Bdl, for the π Bdl, in Sector 3. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 22

23 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 21: Bdl Bdl vs. : The change in Bdl, for the π Bdl, in Sector 4. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 23

24 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 22: Bdl Bdl vs. : The change in Bdl, for the π Bdl, in Sector. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 24

25 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 23: Bdl Bdl vs. : The change in Bdl, for the π Bdl, in Sector 6. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 2

26 Sector Sector Sector Sector Sector Sector Figure 24: Bdl units of kg cm GeV. vs. θ : The change in Bdl, for the π Bdl, in each sector as a function of θ. has 26

27 B π + Corrections This aendix shows detailed lots of the π + momentum corrections. Figures 34, 3,36,37,38 and 39 show the quantity Bdl versus Bdl in sectors 1-6. The lots are broken into ten bins in φ in each sector. Figure 2 shows the corrections made tne run from the g1c data summed over all 6 sectors. The corrections are small, tyically only a few MeV. The π + corrections have a decreasing sloe, in the Bdl Bdl vs. lots, as φ increases. Recall that in Aendix A we also saw this affect in the π corrections and attributed it to a rotation of the drift chambers with resect to the torus coils. A rotation of the drift chambers would be seen in both ostive and negative articles (GeV) θ φ Figure 2: Corrections to the π + momentum as a function of, θ and φ summed over all 6 sectors. 27

28 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 26: Bdl Bdl vs. : The change in Bdl, for the π+ Bdl, in Sector 1. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 28

29 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 27: Bdl Bdl vs. : The change in Bdl, for the π+ Bdl, in Sector 2. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 29

30 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 28: Bdl Bdl vs. : The change in Bdl, for the π+ Bdl, in Sector 3. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 3

31 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 29: Bdl Bdl vs. : The change in Bdl, for the π+ Bdl, in Sector 4. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 31

32 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 3: Bdl Bdl vs. : The change in Bdl, for the π+ Bdl, in Sector. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 32

33 2 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 31: Bdl Bdl vs. : The change in Bdl, for the π+ Bdl, in Sector 6. has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. 33

34 Sector Sector Sector Sector Sector Sector Figure 32: Bdl units of kg cm GeV. vs. θ : The change in Bdl, for the π+ Bdl, in each sector as a function of θ. has 34

35 C Alying the π + corrections to the roton The lots in this aendix show the affects of alying the π + corrections along with the extra eloss correction to the roton. Figures 34, 3,36,37,38 and 39 show the quantity Bdl versus Bdl in sectors 1-6. The lots are broken into ten bins in φ in each sector. The red circles were calculated using the measured momentum with the standard eloss correction alied. The blue triangles were calculated after alying the π + corrections to the roton momentum and using the increased eloss correction. These corrections work quite well, in almost every bin Bdl is consistent with zero. Figure 33 shows the results of alying the π + momentum corrections along with the increase to the ELOSS correction (not the correction itself, just the extra 18%) to the roton momentum for one run in the g1c data. Excet for rotons with mometum less than about 3 MeV/c 2, these corrections are small, tyically less than MeV/c (GeV) θ φ Figure 33: Corrections, including our increase to the eloss correction (but not the standard ELOSS correction), to the roton momentum as a function of, θ and φ summed over all 6 sectors. 3

36 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 34: Bdl Bdl vs. : The change in Bdl, for the roton, in Sector 1. Bdl has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. The red circles were calculated using the measured momentum with the standard eloss correction. The blue triangles were calculated after alying the π + momentum corrections to the roton and using the increased eloss correction. 36

37 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 3: Bdl Bdl vs. : The change in Bdl, for the roton, in Sector 2. Bdl has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. The red circles were calculated using the measured momentum with the standard eloss correction. The blue triangles were calculated after alying the π + momentum corrections to the roton and using the increased eloss correction. 37

38 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 36: Bdl Bdl vs. : The change in Bdl, for the roton, in Sector 3. Bdl has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. The red circles were calculated using the measured momentum with the standard eloss correction. The blue triangles were calculated after alying the π + momentum corrections to the roton and using the increased eloss correction. 38

39 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 37: Bdl Bdl vs. : The change in Bdl, for the roton, in Sector 4. Bdl has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. The red circles were calculated using the measured momentum with the standard eloss correction. The blue triangles were calculated after alying the π + momentum corrections to the roton and using the increased eloss correction. 39

40 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 38: Bdl Bdl vs. : The change in Bdl, for the roton, in Sector. Bdl has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. The red circles were calculated using the measured momentum with the standard eloss correction. The blue triangles were calculated after alying the π + momentum corrections to the roton and using the increased eloss correction. 4

41 φ (-2,-2 ) φ (-2,-1 ) φ (-1,-1 ) φ (-1,- ) φ (-, ) φ (, ) φ (,1 ) φ (1,1 ) φ (1,2 ) φ (2,2 ) Figure 39: Bdl Bdl vs. : The change in Bdl, for the roton, in Sector 6. Bdl has units of kg cm GeV and φ is the sector deendent azimuthal angle, with φ = corresonding to the center of the sector. The red circles were calculated using the measured momentum with the standard eloss correction. The blue triangles were calculated after alying the π + momentum corrections to the roton and using the increased eloss correction. 41

42 Sector Sector Sector Sector Sector Sector Figure 4: Bdl units of kg cm GeV vs. θ : The change in Bdl, for the roton, in each sector as a function of θ. Bdl has. The red circles were calculated using the measured momentum with the standard eloss correction. The blue triangles were calculated after alying the π + momentum corrections to the roton and using the increased eloss correction. 42

Invariant yield calculation

Invariant yield calculation Chater 6 Invariant yield calculation he invariant yield of the neutral ions and η mesons er one minimum bias collision as a function of the transverse momentum is given by E d3 N d 3 = d 3 N d dydφ = d

More information

pp physics, RWTH, WS 2003/04, T.Hebbeker

pp physics, RWTH, WS 2003/04, T.Hebbeker 1. PP TH 03/04 Accelerators and Detectors 1 hysics, RWTH, WS 2003/04, T.Hebbeker 2003-12-03 1. Accelerators and Detectors In the following, we concentrate on the three machines SPS, Tevatron and LHC with

More information

Feedback-error control

Feedback-error control Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller

More information

The Quark-Parton Model

The Quark-Parton Model The Quark-Parton Model Before uarks and gluons were generally acceted Feynman roosed that the roton was made u of oint-like constituents artons Both Bjorken Scaling and the Callan-Gross relationshi can

More information

The Noise Power Ratio - Theory and ADC Testing

The Noise Power Ratio - Theory and ADC Testing The Noise Power Ratio - Theory and ADC Testing FH Irons, KJ Riley, and DM Hummels Abstract This aer develos theory behind the noise ower ratio (NPR) testing of ADCs. A mid-riser formulation is used for

More information

An Improved Calibration Method for a Chopped Pyrgeometer

An Improved Calibration Method for a Chopped Pyrgeometer 96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum

More information

Single and double coincidence nucleon spectra in the weak decay of Λ hypernuclei

Single and double coincidence nucleon spectra in the weak decay of Λ hypernuclei Single and double coincidence nucleon sectra in the weak decay of hyernuclei E. Bauer 1, G. Garbarino 2, A. Parreño 3 and A. Ramos 3 1 Deartamento de Física, Universidad Nacional de La Plata, C. C. 67

More information

One-way ANOVA Inference for one-way ANOVA

One-way ANOVA Inference for one-way ANOVA One-way ANOVA Inference for one-way ANOVA IPS Chater 12.1 2009 W.H. Freeman and Comany Objectives (IPS Chater 12.1) Inference for one-way ANOVA Comaring means The two-samle t statistic An overview of ANOVA

More information

Solved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.

Solved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points. Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the

More information

On Using FASTEM2 for the Special Sensor Microwave Imager (SSM/I) March 15, Godelieve Deblonde Meteorological Service of Canada

On Using FASTEM2 for the Special Sensor Microwave Imager (SSM/I) March 15, Godelieve Deblonde Meteorological Service of Canada On Using FASTEM2 for the Secial Sensor Microwave Imager (SSM/I) March 15, 2001 Godelieve Deblonde Meteorological Service of Canada 1 1. Introduction Fastem2 is a fast model (multile-linear regression model)

More information

arxiv: v1 [hep-ph] 19 Nov 2012

arxiv: v1 [hep-ph] 19 Nov 2012 Evidence for a narrow D 0 state in K η near threshold Bo-Chao Liu 1,, and Ju-Jun Xie,, 1 Deartment of Alied Physics, Xi an Jiaotong University, Xi an, Shanxi 710049, China Theoretical Physics Center for

More information

Higgs Modeling using EXPER and Weak Fusion. by Woody Stanford (c) 2016 Stanford Systems.

Higgs Modeling using EXPER and Weak Fusion. by Woody Stanford (c) 2016 Stanford Systems. iggs Modeling using EXPER and Weak Fusion by Woody Stanford (c) 2016 Stanford Systems. Introduction The EXPER roject, even though its original findings were inconclusive has lead to various ideas as to

More information

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional

More information

4. Score normalization technical details We now discuss the technical details of the score normalization method.

4. Score normalization technical details We now discuss the technical details of the score normalization method. SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules

More information

arxiv: v1 [physics.data-an] 26 Oct 2012

arxiv: v1 [physics.data-an] 26 Oct 2012 Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch

More information

Possibilities for studying few nucleon correlations and Δ isobars in processes with several final state baryons.

Possibilities for studying few nucleon correlations and Δ isobars in processes with several final state baryons. Possibilities for studying few nucleon correlations and Δ isobars in rocesses with several final state baryons. Short-Range Structure of Nuclei at 12 GeV October 26-27, 2007 Jefferson Lab, Newort News,

More information

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment

More information

Performance of a First-Level Muon Trigger with High Momentum Resolution Based on the ATLAS MDT Chambers for HL-LHC

Performance of a First-Level Muon Trigger with High Momentum Resolution Based on the ATLAS MDT Chambers for HL-LHC Performance of a First-Level Muon rigger with High Momentum Resolution Based on the ALAS MD Chambers for HL-LHC P. Gadow, O. Kortner, S. Kortner, H. Kroha, F. Müller, R. Richter Max-Planck-Institut für

More information

Production of Strange Particles in Jets in Heavy-ion Collisions in ALICE

Production of Strange Particles in Jets in Heavy-ion Collisions in ALICE WD'3 Proceedings of Contributed Paers, Part III, 5 57, 23. IBN 978-8-737852-8 MAFYZPRE Production of trange Particles in Jets in Heavy-ion Collisions in ALICE V. učera Charles University in Prague, Faculty

More information

The Forward Detectors of CDF and DØ

The Forward Detectors of CDF and DØ The Forward Detectors of CDF and DØ Konstantin Goulianos The Rockefeller University, 123 York Avenue, New York, NY 165-9965, USA DOI: will be assigned The forward detectors of CDF II are resented with

More information

arxiv: v3 [physics.data-an] 23 May 2011

arxiv: v3 [physics.data-an] 23 May 2011 Date: October, 8 arxiv:.7v [hysics.data-an] May -values for Model Evaluation F. Beaujean, A. Caldwell, D. Kollár, K. Kröninger Max-Planck-Institut für Physik, München, Germany CERN, Geneva, Switzerland

More information

arxiv:hep-ex/ v1 18 Jan 2007

arxiv:hep-ex/ v1 18 Jan 2007 Jet roerties from di-hadron correlations in + collisions at s= GeV Jan Rak for the PHENIX collaboration arxiv:he-ex/77v 8 Jan 7 Abstract Deartment of Physics P.O.Box 5 (YFL) Jyväskylä FI-44 University

More information

Improved search of COSY-TOF on Ө +

Improved search of COSY-TOF on Ө + Pentaquark 05 International Worksho, October 20-22, 2005, Jefferson Lab, U.S. Imroved search of COSY-TOF on Ө + Wolfgang Eyrich Physikalisches Institut, University of Erlangen-Nuremberg, Germany for the

More information

STA 250: Statistics. Notes 7. Bayesian Approach to Statistics. Book chapters: 7.2

STA 250: Statistics. Notes 7. Bayesian Approach to Statistics. Book chapters: 7.2 STA 25: Statistics Notes 7. Bayesian Aroach to Statistics Book chaters: 7.2 1 From calibrating a rocedure to quantifying uncertainty We saw that the central idea of classical testing is to rovide a rigorous

More information

A Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition

A Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino

More information

arxiv: v1 [hep-ex] 1 Feb 2018

arxiv: v1 [hep-ex] 1 Feb 2018 arxiv:8.6v [he-ex] Feb 8 MA Wigner RCP E-mail: varga-kofarago.monika@wigner.mta.hu In heavy-ion collisions, the quark gluon lasma is exected to be roduced, which is an almost erfect liquid that made u

More information

Diffractive results from CDF

Diffractive results from CDF International Journal of Modern Physics A Vol. 3, No. 8 (5) 543 ( ages) c World Scientific Publishing Comany DOI:.4/S775X5438 Diffractive results from CDF Konstantin Goulianos Exerimental High Energy Physics,

More information

University of North Carolina-Charlotte Department of Electrical and Computer Engineering ECGR 4143/5195 Electrical Machinery Fall 2009

University of North Carolina-Charlotte Department of Electrical and Computer Engineering ECGR 4143/5195 Electrical Machinery Fall 2009 University of North Carolina-Charlotte Deartment of Electrical and Comuter Engineering ECG 4143/5195 Electrical Machinery Fall 9 Problem Set 5 Part Due: Friday October 3 Problem 3: Modeling the exerimental

More information

State Estimation with ARMarkov Models

State Estimation with ARMarkov Models Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,

More information

Outline Solutions to Particle Physics Problem Sheet 1

Outline Solutions to Particle Physics Problem Sheet 1 2010 Subatomic: Particle Physics 1 Outline Solutions to Particle Physics Problem Sheet 1 1. List all fundamental fermions in the Standard Model There are six letons and six quarks. Letons: e, ν e, µ, ν

More information

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO) Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment

More information

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College

More information

Maximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming

Maximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming Maximum Entroy and the Stress Distribution in Soft Disk Packings Above Jamming Yegang Wu and S. Teitel Deartment of Physics and Astronomy, University of ochester, ochester, New York 467, USA (Dated: August

More information

MEASUREMENT OF THE INCLUSIVE ELECTRON (POSITRON) +PROTON SCATTERING CROSS SECTION AT HIGH INELASTICITY y USING H1 DATA *

MEASUREMENT OF THE INCLUSIVE ELECTRON (POSITRON) +PROTON SCATTERING CROSS SECTION AT HIGH INELASTICITY y USING H1 DATA * Romanian Reorts in Physics, Vol. 65, No. 2, P. 420 426, 2013 MEASUREMENT OF THE INCLUSIVE ELECTRON (POSITRON) +PROTON SCATTERING CROSS SECTION AT HIGH INELASTICITY y USING H1 DATA * IVANA PICURIC, ON BEHALF

More information

MODEL-BASED MULTIPLE FAULT DETECTION AND ISOLATION FOR NONLINEAR SYSTEMS

MODEL-BASED MULTIPLE FAULT DETECTION AND ISOLATION FOR NONLINEAR SYSTEMS MODEL-BASED MULIPLE FAUL DEECION AND ISOLAION FOR NONLINEAR SYSEMS Ivan Castillo, and homas F. Edgar he University of exas at Austin Austin, X 78712 David Hill Chemstations Houston, X 77009 Abstract A

More information

The non-stochastic multi-armed bandit problem

The non-stochastic multi-armed bandit problem Submitted for journal ublication. The non-stochastic multi-armed bandit roblem Peter Auer Institute for Theoretical Comuter Science Graz University of Technology A-8010 Graz (Austria) auer@igi.tu-graz.ac.at

More information

A Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression

A Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Journal of Modern Alied Statistical Methods Volume Issue Article 7 --03 A Comarison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Ghadban Khalaf King Khalid University, Saudi

More information

Session 5: Review of Classical Astrodynamics

Session 5: Review of Classical Astrodynamics Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve

More information

PHYSICAL REVIEW LETTERS

PHYSICAL REVIEW LETTERS PHYSICAL REVIEW LETTERS VOLUME 81 20 JULY 1998 NUMBER 3 Searated-Path Ramsey Atom Interferometer P. D. Featonby, G. S. Summy, C. L. Webb, R. M. Godun, M. K. Oberthaler, A. C. Wilson, C. J. Foot, and K.

More information

arxiv: v1 [hep-ex] 8 Jun 2017

arxiv: v1 [hep-ex] 8 Jun 2017 UCHEP 17 05 6 Aril 017 Prosects for time-deendent mixing and CP-violation measurements at Belle II arxiv:1706.0363v1 [he-ex] 8 Jun 017 Physics Deartment, University of Cincinnati, Cincinnati, Ohio 51 E-mail:

More information

Computer arithmetic. Intensive Computation. Annalisa Massini 2017/2018

Computer arithmetic. Intensive Computation. Annalisa Massini 2017/2018 Comuter arithmetic Intensive Comutation Annalisa Massini 7/8 Intensive Comutation - 7/8 References Comuter Architecture - A Quantitative Aroach Hennessy Patterson Aendix J Intensive Comutation - 7/8 3

More information

Bayesian Spatially Varying Coefficient Models in the Presence of Collinearity

Bayesian Spatially Varying Coefficient Models in the Presence of Collinearity Bayesian Satially Varying Coefficient Models in the Presence of Collinearity David C. Wheeler 1, Catherine A. Calder 1 he Ohio State University 1 Abstract he belief that relationshis between exlanatory

More information

FE FORMULATIONS FOR PLASTICITY

FE FORMULATIONS FOR PLASTICITY G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND

More information

arxiv: v2 [hep-ph] 13 Dec 2018

arxiv: v2 [hep-ph] 13 Dec 2018 Predictions for azimuthal anisotroy in Xe+Xe collisions at = 5.44 ev using a multihase transort model Sushanta riathy, Sudian De, Mohammed Younus, and Raghunath Sahoo Disciline of Physics, School of Basic

More information

Radial Basis Function Networks: Algorithms

Radial Basis Function Networks: Algorithms Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.

More information

A New Asymmetric Interaction Ridge (AIR) Regression Method

A New Asymmetric Interaction Ridge (AIR) Regression Method A New Asymmetric Interaction Ridge (AIR) Regression Method by Kristofer Månsson, Ghazi Shukur, and Pär Sölander The Swedish Retail Institute, HUI Research, Stockholm, Sweden. Deartment of Economics and

More information

Lower bound solutions for bearing capacity of jointed rock

Lower bound solutions for bearing capacity of jointed rock Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying

More information

Hard Diffraction Results and Prospects at the Tevatron

Hard Diffraction Results and Prospects at the Tevatron Hard Diffraction Results and Prosects at the Tevatron University of Manchester On behalf of the CDF & DØ Collaboration ISMD 2005 14 th August, Kromeriz Tevatron Tevatron at Fermilab -collisions Run II:

More information

The Poisson Regression Model

The Poisson Regression Model The Poisson Regression Model The Poisson regression model aims at modeling a counting variable Y, counting the number of times that a certain event occurs during a given time eriod. We observe a samle

More information

A mgh ENERGY NEUTRON DETECTOR USING PROPORTIONAL WIRE CHAMBERS (Ii 727) Presented by M. Atac National Accelerator Laboratory* Batavia, Illinois

A mgh ENERGY NEUTRON DETECTOR USING PROPORTIONAL WIRE CHAMBERS (Ii 727) Presented by M. Atac National Accelerator Laboratory* Batavia, Illinois A mgh ENERGY NEUTRON DETECTOR USNG PROPORTONAL WRE CHAMBERS (i 727) Presented by M. Atac National Accelerator Laboratory* Batavia, llinois An exeriment to study the roerties of negative hyerons roduced

More information

Measuring center and spread for density curves. Calculating probabilities using the standard Normal Table (CIS Chapter 8, p 105 mainly p114)

Measuring center and spread for density curves. Calculating probabilities using the standard Normal Table (CIS Chapter 8, p 105 mainly p114) Objectives Density curves Measuring center and sread for density curves Normal distributions The 68-95-99.7 (Emirical) rule Standardizing observations Calculating robabilities using the standard Normal

More information

Developing A Deterioration Probabilistic Model for Rail Wear

Developing A Deterioration Probabilistic Model for Rail Wear International Journal of Traffic and Transortation Engineering 2012, 1(2): 13-18 DOI: 10.5923/j.ijtte.20120102.02 Develoing A Deterioration Probabilistic Model for Rail Wear Jabbar-Ali Zakeri *, Shahrbanoo

More information

Sound speed profile structure and variability measured over flat terrain

Sound speed profile structure and variability measured over flat terrain Sound seed rofile structure and variability measured over flat terrain Bradley, SG, Waddington, DC and Von Hunerbein, S Title Authors Tye URL Published Date 6 Sound seed rofile structure and variability

More information

Meshless Methods for Scientific Computing Final Project

Meshless Methods for Scientific Computing Final Project Meshless Methods for Scientific Comuting Final Project D0051008 洪啟耀 Introduction Floating island becomes an imortant study in recent years, because the lands we can use are limit, so eole start thinking

More information

u-channel Omega Meson Production from the Fpi-2 Experiment

u-channel Omega Meson Production from the Fpi-2 Experiment u-channel Omega Meson Production from the Fi- Exeriment Bill (Wenliang) Li Hall C Worksho. January, 015. Outline Where the data come from Theoretical justification Plan for data analysis Wenliang Li, Det.

More information

Recent flow results at RHIC

Recent flow results at RHIC Recent flow results at RHIC Hiroshi Masui / University of sukuba Flow and heavy flavour worksho in high energy heavy ion collisions: GRN worksho, Inchon, Feb./24-26, 215 H. Masui / Univ. of sukuba 1 /3

More information

arxiv:cond-mat/ v2 25 Sep 2002

arxiv:cond-mat/ v2 25 Sep 2002 Energy fluctuations at the multicritical oint in two-dimensional sin glasses arxiv:cond-mat/0207694 v2 25 Se 2002 1. Introduction Hidetoshi Nishimori, Cyril Falvo and Yukiyasu Ozeki Deartment of Physics,

More information

Monopolist s mark-up and the elasticity of substitution

Monopolist s mark-up and the elasticity of substitution Croatian Oerational Research Review 377 CRORR 8(7), 377 39 Monoolist s mark-u and the elasticity of substitution Ilko Vrankić, Mira Kran, and Tomislav Herceg Deartment of Economic Theory, Faculty of Economics

More information

Positive decomposition of transfer functions with multiple poles

Positive decomposition of transfer functions with multiple poles Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.

More information

Chapter 4 Randomized Blocks, Latin Squares, and Related Designs Solutions

Chapter 4 Randomized Blocks, Latin Squares, and Related Designs Solutions Solutions from Montgomery, D. C. (008) Design and Analysis of Exeriments, Wiley, NY Chater 4 Randomized Blocks, Latin Squares, and Related Designs Solutions 4.. The ANOVA from a randomized comlete block

More information

Pressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids

Pressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids American Journal of Alied Sciences (9): 19-195, 5 ISSN 1546-939 5 Science Publications Pressure-sensitivity Effects on Toughness Measurements of Comact Tension Secimens for Strain-hardening Solids Abdulhamid

More information

New technique and results of cosmic ray investigations in the energy interval ev

New technique and results of cosmic ray investigations in the energy interval ev EPJ Web of Conferences 53, 08001 (2013) DOI: 10.1051/ejconf/20135308001 C Owned by the authors, ublished by EDP Sciences, 2013 New technique and results of cosmic ray investigations in the energy interval

More information

Pulse Propagation in Optical Fibers using the Moment Method

Pulse Propagation in Optical Fibers using the Moment Method Pulse Proagation in Otical Fibers using the Moment Method Bruno Miguel Viçoso Gonçalves das Mercês, Instituto Suerior Técnico Abstract The scoe of this aer is to use the semianalytic technique of the Moment

More information

NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS

NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372

More information

Period-two cycles in a feedforward layered neural network model with symmetric sequence processing

Period-two cycles in a feedforward layered neural network model with symmetric sequence processing PHYSICAL REVIEW E 75, 4197 27 Period-two cycles in a feedforward layered neural network model with symmetric sequence rocessing F. L. Metz and W. K. Theumann Instituto de Física, Universidade Federal do

More information

Mukesh Saini. Florida State University, Tallahassee, FL. February 26, FSU Nuclear Physics Seminar. February 26

Mukesh Saini. Florida State University, Tallahassee, FL. February 26, FSU Nuclear Physics Seminar. February 26 C o l u m n 3 Mukesh Saini Florida State University, Tallahassee, FL, 2010 1 OUTLINE Introduction Meson Spectroscopy Strangeonia Experiment CEBAF & CLAS HyCLAS & g12 Calibrations Analysis Summary 2 QCD

More information

OPTIMIZATION OF EARTH FLIGHT TEST TRAJECTORIES TO QUALIFY PARACHUTES FOR USE ON MARS

OPTIMIZATION OF EARTH FLIGHT TEST TRAJECTORIES TO QUALIFY PARACHUTES FOR USE ON MARS OPTIMIZATION OF EARTH FLIGHT TEST TRAJECTORIES TO QUALIFY PARACHUTES FOR USE ON MARS Christoher L. Tanner (1) (1) Sace Systems Design Laboratory, Daniel Guggenheim School of Aerosace Engineering Georgia

More information

Submicrometer Position Control of Single Trapped Neutral Atoms

Submicrometer Position Control of Single Trapped Neutral Atoms Dotsenko, I and Alt, W and Khudaverdyan, M and Kuhr, S and Meschede, D and Miroshnychenko, Y and Schrader, D and Rauschenbeutel, A (25) Submicrometer osition control of single traed neutral atoms. Physical

More information

Structure of 11 Be studied in β-delayed neutron- and γ- decay from polarized 11 Li

Structure of 11 Be studied in β-delayed neutron- and γ- decay from polarized 11 Li Nuclear Physics A 46 (4) c c Structure of Be studied in β-delayed neutron- and γ- decay from olarized Li Y. Hirayama a, T. Shimoda a,h.izumi a,h.yano a,m.yagi a, A. Hatakeyama b, C.D.P. Levy c,k.p.jackson

More information

1 University of Edinburgh, 2 British Geological Survey, 3 China University of Petroleum

1 University of Edinburgh, 2 British Geological Survey, 3 China University of Petroleum Estimation of fluid mobility from frequency deendent azimuthal AVO a synthetic model study Yingrui Ren 1*, Xiaoyang Wu 2, Mark Chaman 1 and Xiangyang Li 2,3 1 University of Edinburgh, 2 British Geological

More information

from low Q (' 0:45GeV ) inclusive hotoroduction. Finally, we discuss the contribution that olarized HERA could make to the measurement of these high s

from low Q (' 0:45GeV ) inclusive hotoroduction. Finally, we discuss the contribution that olarized HERA could make to the measurement of these high s The Drell-Hearn-Gerasimov Sum-Rule at Polarized HERA S. D. Bass a, M. M. Brisudova b and A. De Roeck c a Institut fur Theoretische Kernhysik, Universitat Bonn, Nussallee 4{6, D-535 Bonn, Germany b Theoretical

More information

Plotting the Wilson distribution

Plotting the Wilson distribution , Survey of English Usage, University College London Setember 018 1 1. Introduction We have discussed the Wilson score interval at length elsewhere (Wallis 013a, b). Given an observed Binomial roortion

More information

Approximating min-max k-clustering

Approximating min-max k-clustering Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost

More information

Highlights from the ATLAS experiment

Highlights from the ATLAS experiment Nuclear Physics A Nuclear Physics A (28) 7 www.elsevier.com/locate/rocedia XXVIIth International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 28) Highlights from the ALAS exeriment

More information

Measuring center and spread for density curves. Calculating probabilities using the standard Normal Table (CIS Chapter 8, p 105 mainly p114)

Measuring center and spread for density curves. Calculating probabilities using the standard Normal Table (CIS Chapter 8, p 105 mainly p114) Objectives 1.3 Density curves and Normal distributions Density curves Measuring center and sread for density curves Normal distributions The 68-95-99.7 (Emirical) rule Standardizing observations Calculating

More information

Recursive Estimation of the Preisach Density function for a Smart Actuator

Recursive Estimation of the Preisach Density function for a Smart Actuator Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator

More information

Jet energy measurement and its systematic uncertainty in proton proton collisions at s = 7 TeV with the ATLAS detector

Jet energy measurement and its systematic uncertainty in proton proton collisions at s = 7 TeV with the ATLAS detector Eur. Phys. J. C (5) 75:7 DOI.4/ejc/s5-4-9-y Regular Article - Exerimental Physics Jet energy measurement and its systematic uncertainty in roton roton collisions at s = 7 ev with the ALAS detector ALAS

More information

The influence of cluster emission and the symmetry energy on neutron-proton spectral double ratios

The influence of cluster emission and the symmetry energy on neutron-proton spectral double ratios The influence of cluster emission and the symmetry energy on neutron-roton sectral double ratios Yingxun Zhang 1,,3, P. Danielewicz 1,,4, M. Famiano 5, Zhuxia Li 3, W.G. Lynch 1,,4, M.B. Tsang 1,,4* 1

More information

Temperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures

Temperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures Indian Journal of Pure & Alied Physics Vol. 44, December 2006,. 953-958 Temerature, current and doing deendence of non-ideality factor for n and nn unch-through structures Khurshed Ahmad Shah & S S Islam

More information

KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS

KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS 4 th International Conference on Earthquake Geotechnical Engineering June 2-28, 27 KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS Misko CUBRINOVSKI 1, Hayden BOWEN 1 ABSTRACT Two methods for analysis

More information

Morten Frydenberg Section for Biostatistics Version :Friday, 05 September 2014

Morten Frydenberg Section for Biostatistics Version :Friday, 05 September 2014 Morten Frydenberg Section for Biostatistics Version :Friday, 05 Setember 204 All models are aroximations! The best model does not exist! Comlicated models needs a lot of data. lower your ambitions or get

More information

LUMINOSITY DETERMINATION AT THE TEVATRON*

LUMINOSITY DETERMINATION AT THE TEVATRON* LUMINOSITY DETERMINATION AT THE TEVATRON* V. Paadimitriou #, Fermilab, Batavia, IL 60510, U.S.A. Abstract In this aer we discuss the luminosity determination at the Tevatron. We discuss luminosity measurements

More information

Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier

Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier Australian Journal of Basic and Alied Sciences, 5(12): 2010-2020, 2011 ISSN 1991-8178 Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doed Fiber Amlifier

More information

In this lesson you will use the Cloud Chamber applet to investigate the alpha decay process.

In this lesson you will use the Cloud Chamber applet to investigate the alpha decay process. Alha Decay In this lesson you will use the Cloud Chamber alet to investigate the alha decay rocess. Prerequisites: You should be familiar with the hysics of charged articles moving in magnetic fields and

More information

Notes on Instrumental Variables Methods

Notes on Instrumental Variables Methods Notes on Instrumental Variables Methods Michele Pellizzari IGIER-Bocconi, IZA and frdb 1 The Instrumental Variable Estimator Instrumental variable estimation is the classical solution to the roblem of

More information

Lower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data

Lower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data Quality Technology & Quantitative Management Vol. 1, No.,. 51-65, 15 QTQM IAQM 15 Lower onfidence Bound for Process-Yield Index with Autocorrelated Process Data Fu-Kwun Wang * and Yeneneh Tamirat Deartment

More information

Churilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics

Churilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated

More information

All-fiber Optical Parametric Oscillator

All-fiber Optical Parametric Oscillator All-fiber Otical Parametric Oscillator Chengao Wang Otical Science and Engineering, Deartment of Physics & Astronomy, University of New Mexico Albuquerque, NM 87131-0001, USA Abstract All-fiber otical

More information

where x i is the ith coordinate of x R N. 1. Show that the following upper bound holds for the growth function of H:

where x i is the ith coordinate of x R N. 1. Show that the following upper bound holds for the growth function of H: Mehryar Mohri Foundations of Machine Learning Courant Institute of Mathematical Sciences Homework assignment 2 October 25, 2017 Due: November 08, 2017 A. Growth function Growth function of stum functions.

More information

Ordered and disordered dynamics in random networks

Ordered and disordered dynamics in random networks EUROPHYSICS LETTERS 5 March 998 Eurohys. Lett., 4 (6),. 599-604 (998) Ordered and disordered dynamics in random networks L. Glass ( )andc. Hill 2,3 Deartment of Physiology, McGill University 3655 Drummond

More information

Characteristics of Beam-Based Flexure Modules

Characteristics of Beam-Based Flexure Modules Shorya Awtar e-mail: shorya@mit.edu Alexander H. Slocum e-mail: slocum@mit.edu Precision Engineering Research Grou, Massachusetts Institute of Technology, Cambridge, MA 039 Edi Sevincer Omega Advanced

More information

VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES

VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES Journal of Sound and Vibration (998) 22(5), 78 85 VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES Acoustics and Dynamics Laboratory, Deartment of Mechanical Engineering, The

More information

ECON Answers Homework #2

ECON Answers Homework #2 ECON 33 - Answers Homework #2 Exercise : Denote by x the number of containers of tye H roduced, y the number of containers of tye T and z the number of containers of tye I. There are 3 inut equations that

More information

The search for the η -mesic nuclei in the LEPS2/BGOegg experiment

The search for the η -mesic nuclei in the LEPS2/BGOegg experiment The search for the η -mesic nuclei in the LEPS2/BGOegg exeriment N. Tomida (RCNP, Osaka Univ.) 2018/Nov/13 QNP2018, Tsukuba Theory η -nucleus otical otential U(r) = ( V 0 + iw 0 ) x NJL model : V 0 = 150

More information

CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules

CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 13: PF Design II The Coil Solver

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 13: PF Design II The Coil Solver .615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 13: PF Design II The Coil Solver Introduction 1. Let us assume that we have successfully solved the Grad Shafranov equation for a fixed boundary

More information

ute measures of uncertainty called standard errors for these b j estimates and the resulting forecasts if certain conditions are satis- ed. Note the e

ute measures of uncertainty called standard errors for these b j estimates and the resulting forecasts if certain conditions are satis- ed. Note the e Regression with Time Series Errors David A. Dickey, North Carolina State University Abstract: The basic assumtions of regression are reviewed. Grahical and statistical methods for checking the assumtions

More information

Measurements of forward neutron and neutral pion productions with the LHCf detector

Measurements of forward neutron and neutral pion productions with the LHCf detector Measurements of forward neutron and neutral ion roductions with the LHCf detector Gaku Mitsuka (University of Florence, INFN Firenze, JSPS fellow) MPI4-7 Nov. 4, Krakow Outline Introduction and hysics

More information

Soft QCD Results from ATLAS and CMS

Soft QCD Results from ATLAS and CMS Soft QCD Results from ALAS and Moriond, March oics Proerties of minimum bias events - transverse momentum, seudoraidity and event-by-event multilicity distributions of charged articles Underlying event

More information

Christina Mesropian The Rockefeller University

Christina Mesropian The Rockefeller University Christina Mesroian The Rockefeller University Introduction Definition of diffractive rocesses Single Diffraction Hard Single Diffraction Soft Single diffraction Double Diffraction Hard Double Diffraction

More information