jump to navigation

Top production: gluon fusion or qq annihilation ? August 31, 2006

Posted by dorigo in news, physics, science.
add a comment

Ok, I vowed a more varied reporting of physics results here in the last post. So why not show one I was unaware of, from my own collaboration, which is also quite interesting and new ?

Toronto University physicists in CDF have studied observable quantities showing a dependence on the production mechanism of top quark pairs at the Tevatron. They have found that when top quark pairs are produced by gluon fusion, there is a higher number of soft tracks in the interaction, ultimately due (or so I guess) to the higher color charge of the initial state, and the resulting higher probability of soft gluon radiation.

They studied soft tracks in several high-Pt processes and concluded that they could indeed extract a measurement of the fraction of top quark pairs due to gluon fusion processes.

Their fit is shown in the plot below. They find a fraction sigma(gg–>tt)/sigma(tt) of 0.25+-0.24+-0.10, when QCD predicts this fraction to be 15% at the Tevatron energy (will be 85% at the LHC). Quite interesting, isn’t it ?

If you want more information on this measurement, please visit the web site of the analysis:

http://www-cdf.fnal.gov/physics/new/top/2005/topProp/ttProduction/public.html

H limits, Mt, H limits, Mt… It’s getting boring August 31, 2006

Posted by dorigo in news, personal, physics, science.
add a comment

I believe I may be starting to get the occasional readers of this blog bored as of late, since as far as particle physics topics are concerned, I have been talking mostly, if not exclusively, of top quark mass determinations and Higgs boson cross section limits.

My excuses. I will try to be more varied in the topics I pick in the future… But I thought today I would post here the summary slide I sent to the convener of my session at the QCHS 2006 conference. So here it is:

I think it is kind of hard to read… However, the top right plot has been posted here a week ago, and the bottom left also appeared here maybe a month ago. The most interesting one is the new view of the top mass / W mass plane, with the new world averages direct measurements (blue ellipse, which shows 68% limits).

The summary Top Search Plot (C) August 30, 2006

Posted by dorigo in news, physics, science.
add a comment

A mail message by Patrizia, who asks me about a plot I put together for a conference a while ago, stimulates me into reproducing it here. It represents a huge amount of information, so let’s try to discuss it bottoms up.

First of all, this plot has not been updated since 2004 - so the latest measurements of the top quark mass have not been inserted. They are not needed though, because by now the precision on the mass measurements make it impossible to see the error bars on a plot such as this.

The plot represent, indeed, 10 years of searches for the top quark and of theoretical estimates for its mass, as well as ten years of direct measurements. The horizontal scale goes from year 1987 to 2005, and the vertical scale from 0 to 300 GeV of mass.

In dark blue, the area shows the region excluded by direct searches at electron-positron colliders: a business ended in the early nineties, when LEP found no evidence, and set a limit of the top mass at half the Z mass value (if Mt<MZ/2, the Z could decay to two top quarks, and indeed they would have found it).

In light blue is the region excluded by proton-antiproton colliders - basically, the Tevatron, which sought it with more and more statistics and until 1994 could not find it, and only exclude its existence for low mass - a low mass would have implied a high production rate, and they would have not missed it.

The upper green region is instead picturing theoretical upper limits on the top quark mass, extracted from the combination of theoretical arguments and experimental evidence in other areas of particle physics (for instance, B mixing results).

The yellow area is the most intriguing: it shows theoretical determinations of the top quark mass as a function of time, based on dozens of different studies, which incorporated theory and the latest measurements of electroweak physics observables to infer the most probable mass of the top quark - in case it existed. Notice how the red line (the most probable value by these theoretical arguments) bounces around until the CDF evidence in 1994, when all at once the indirect evidence by LEP started to be a whole lot more compelling… Hmmm.

Finally, the points with error bars show the direct measurements by CDF and D0 until 1995. Having been in CDF for 14 years, I cannot help pointing out here that our very first evidence had a mass measurement which was right on the money, from the very start! Instead, theoretical estimates seemed to follow the experimental lower limits in an ascending trend.

Quite an informative plot, if you ask me!

Lightning on a plane August 30, 2006

Posted by dorigo in Blogroll, internet, news, physics, science, travel.
3 comments

Quite an interesting post over at the Cocktail Party Physics site ( http://twistedphysics.typepad.com/cocktail_party_physics/ ). And pretty much in theme with the recent fashion of scaring people with flight fear by adding more scary ingredients (see http://www.snakesonaplane.com/ for a commercial example using snakes, or http://www.cnn.com/2006/WORLD/europe/08/11/terror.plot/index.html for liquid explosives)

The not-so-rare occurrence of lightning striking airplanes is discussed in detail. It appears that on average any commercial airplane is struck once a year, usually without any consequence to the safety of the flight. There is even a dramatic picture… Check it out.

Weather forecasts August 29, 2006

Posted by dorigo in astronomy, computers, mathematics, personal, science.
add a comment

Tonight two different computer models for weather forecasts over northern Italy were not predicting any rain at all… But my instinct was.

We went out for a pizza dinner with friends, and I decided I would trust computer science for once, leaving umbrellas behind. Of course, in the matter of an hour, rain was pouring. We did not get too wet anyways, because the restaurant is very close to our house. But the poor level of reliability of weather forecasts in Italy annoys me. We are talking about predictions for the next six hours, not six days ahead!

Isn’t weather prediction a problem that can be solved by just throwing kilospecints at it ? Evidently not. And this is for precipitations - cloud cover is even worse to predict, and as an astronomer that does not make me too happy.

Chisquared envelopes August 29, 2006

Posted by dorigo in computers, mathematics, news, personal, physics, science.
2 comments

I spent the afternoon discussing how to estimate a systematic uncertainty on our b-Jet energy scale factor from the Z–>bb signal with Julien.

The question is very technical, but it boils down to deciding how to deal with variations of a fit chisquared (actually a log-likelihood value, which we trivially convert to a chisquared) versus a fit parameter x to determine the error on x.

Normally, one would simply fit the chisquared variations around the minimum (where x is x_best) with a quadratic or cubic function chi2(x), determine the parameter values x_left and x_right when the chisquared has grown by one unit:

chi2(x_left)=chi2(x_right)=chi2(x_best)+1

and thus determine the parameter range as

x = x_best +- (x_right-x_left)/2.

The problem is that the values of chi2(x) are coming from different fits to the same base of data, using different functional forms to fit them, but the forms are correlated among each other.

Well, we concluded that correlations do not matter, and that a minimum of a chisquared is a minimum regardless. So, since he has hundreds of different fits (and so hundreds of chi2(x) values), rather than fitting all the points with a quadratic or a cubic, he should just find the lower envelope of the points, in the chi2 vs x plane.

The envelope could have any shape, but is by definition convex. I am not sure of the mathematical description of such a curve, but I know how to draw one: just imagine you throw straight lines at the set of points from the bottom of the plot, in all directions. The lines stop when they hit a point, or better, they rotate until they hit a second point and then stop.

Once all possible lines have been considered, the polygonal which goes from point to point is the sought envelope function.

Am I out of my mind, or does this have a reasonable mathematical formulation ? Did anybody out there used a similar method to find a parameter range ?

Search engine terms yesterday August 29, 2006

Posted by dorigo in games, humor, internet, language, personal.
2 comments

Just an update of my periodic survey of what people googled to end up visiting my page:

  • “easy physicists”: bingo. I sure am one.
  • “leaning to do cpr”: you may save lives, but you first have to learn how to spell.
  • “fortran program for gaussian smearing”: this guy will go a long way to avoid doing his homework. Come on! Programming is fun!
  • “universe BRAIN”: I beg your pardon ?
  • “8-27-2006 and mars conjuncture”: no idea what this guy was searching for. The funny thing is, my google does not produce *any* search results on this one, let alone my blog…
  • “beverage mia”: well, Mia is a student of mine, and she sure enjoys a drink or four every now and then. But I doubt she invented any.
  • “Hadronization GOD”: wow. Mystical theoreticians ? Skilled phenomenologists ?
  • “dorigo LA Marathon”: nope, you won’t see me there any time soon.

Quite fun to read what people search… I wonder why google does not have a “search of the day” feature… Or maybe they do ? 

The Say of the Week August 28, 2006

Posted by dorigo in games, humor.
1 comment so far

Sex without love is an empty experience, but, as empty experiences go, it’s one of the best

(W.Allen)

Luminosity and Cross Section August 28, 2006

Posted by dorigo in news, physics, science.
2 comments

In a comment to a post of a few days ago (http://dorigo.wordpress.com/2006/08/24/new-z-bb-signal-from-d0/), Markk writes:

“A kind of dumb question from a non-particle physics person. I see Pb-1 all over and understand it as an area, I even understand using it as luminosity but what is the top unit? I always see, say, L = 333 Pb-1 which makes sense as a luminosity if say energy is the “top” value. Or as you say above something like “700 pb-1 of data”, I don’t get it. What is the default value assumed that I am missing to make more sense of this stuff.”

I think I wrote about this in the past, but I will try to do my homework again here. Let me start by saying that particle fluxes are naturally computed as a number divided by area and time: N = [cm^(-2)][s^(-1)] is the way we write such a thing. These are the units of instantaneous luminosity. If you integrate over time (i.e., sum over a period of time) an instantaneous luminosity, you get what is called integrated luminosity: L = [cm^(-2)].

For instance, you could have a flux of a billion protons per square micron per second, and if you wanted to count the integrated flux per unit area in a minute you would do:

L = 1,000,000,000 * (100,000,000/cm^(2)) / s * (60 s) = 6*10^(18) cm^(-2).

That is beacuse there are 100 million squared microns in a squared cm, and you have multiplied the instantaneous luminosity by the total time (60s). Physicists readily use cm^(-2) s^(-1) as units for instantaneous luminosity, but when they time-integrate, they find out it is much less cumbersome to use a tiny fraction of the square cm as a reference area.

A picobarn is equal to 10^(-36) cm^2: a square of 10^(-18) cm across. By using picobarns, integrated luminosities take on much more manageable forms, as the following example shows.

A collider provides a luminosity of 1.055*10^(31) cm^(-2) s^(-1) and runs for a year. What is the integrated luminosity ? How many events of a kind that has a cross section of 10 pb were produced ?

A year is 86400×365.26 = 31,558,464 seconds. A cm^(-2) is 10^(36) picobarns. So,

L = 31,558,464 s * 1.055*10^(32) / (10^(36) pb s) = 333 pb^(-1).

To find how many events are produced of a process that has a cross section of 10 pb, just multiply:

N = sigma * L = 10 pb * 333 pb^(-1) = 3330.

What is, then, a ”cross section” ? It is something with area units. Picture a proton and imagine another particle (say, as at the Tevatron, an antiproton) coming towards it. They will collide if the axis of motion of the antiproton lies within a certain area, centered on the proton. That is the total proton cross section: about a tenth of a barn, or 10^(-25) cm^2.

When a proton and an antiproton collide, the energy released can produce many different results. Most of the times the two bounce off one another, but very rarely (with an increasing likelihood as the collision energy increases) they give rise to rare processes when massive particles are produced. For instance, only once every 10 million times at the Tevatron energy the collision gives rise to a Z boson. A simple measure of the probability of producing a Z is therefore the total proton cross section divided by this 10 million factor: 0.1 barns divided by 10 million makes 10 nanobarns (or 10,000 picobarns if you prefer). So we say that Z production at the Tevatron has a cross section of 10 nb.

Then, specifically for the case of Z->bb production in CDF (the plot this discussion originated from), the Z decays to a bb pair only 15% of the times, and then our selection cuts accept only one event every 150: that is why the “effective” cross section turns out to be 10 pb: 3330 events in 333 pb^(-1).

CPR, Heimlich, and Tracheotomy August 26, 2006

Posted by dorigo in food, internet, news, science.
3 comments

Anybody should know how to perform these. The more so if you are a parent.

So for your convenience, here are three sites to visit, print out, and keep in your kitchen - possibly attached to your bulletin board if you have one. These save lives!

The first is a site which shows a video on how to perform CPR on a unresponsive, non-breathing child (e.g. after an electric choc). You can even download the video and keep it on your computer: http://depts.washington.edu/learncpr/videodemo/child-cpr-video.html

The second site is a good explanation of how to perform the Heimlich maneuver, with several drawings. You can even learn how to do it on yourself with your fists or by leaning on a chair: http://www.heimlichinstitute.org/howtodo.html

The third site explains very clearly and concisely how to do a tracheotomy to a patient who is choking and to whom the Heimlich maneuver has not produced any result. After reading this, I think I would do it if conditions arose: http://www.tracheostomy.com/surgery/emergency.htm

I hope I will never need to do any of these. But choking is by no means a rare occurrence in children!