Computing a differential from plotted data

alex3
Messages
43
Reaction score
0
I have a differential that depends only on \cos{\theta}

\frac{\operatorname{d}\sigma}{\operatorname{d} \Omega} = f(\cos{\theta})

I am numerically solving this differential equation for \sigma, which physically is a cross section, for 0 \leq \theta \leq \pi. The differential contains a parameter s. I am solving the differential for a given value of s, then incrementing this parameter and solving the differential again.

I am solving the differential over a range of s (3 to 200), in increments of 0.1 (s is, for the curious, my particle accelerator energy in GeV).

My problem is that I'm asked to now plot the transverse momentum

p_{T} = \lvert p_{f} \rvert \sin{\theta}

I've been told that I am effectively trying to plot

\frac{\operatorname{d} \sigma}{\operatorname{d} p_{T}}

and that I'll have to do some sort of histogram/binning to make my plot.

I'm not sure what I should binning; the differentials, the final integrated cross section?

So, in essence, my questions is how can I plot \frac{\operatorname{d} \sigma}{\operatorname{d} p_{T}} given that I have numerical data for \frac{\operatorname{d}\sigma}{\operatorname{d} \Omega} (and \sigma)?
 
Last edited:
Physics news on Phys.org
Normally posts which don't use the homework template tend to be at low priority, so keep that in mind for the future... but anyway: \mathrm{d}\sigma/\mathrm{d}p_T can be physically interpreted as the "amount" of cross section which falls into a certain (small) range of p_T, divided by the size of that range. When you are plotting this differential cross section, it makes sense to plot it versus p_T.
 
diazona said:
Normally posts which don't use the homework template tend to be at low priority, so keep that in mind for the future
Apologies, I had trouble fitting the question into the format. I shall try harder next time :)

\mathrm{d}\sigma/\mathrm{d}p_T can be physically interpreted as the "amount" of cross section which falls into a certain (small) range of p_T, divided by the size of that range. When you are plotting this differential cross section, it makes sense to plot it versus p_T.

Do you mean I should plot \operatorname{d}\sigma/\operatorname{d} p_T against p_T? I know I need to do that, my problem is trying to find the differential values. i.e. How could I figure out "the 'amount' of cross section which falls into a certain (small) range of p_T"?

Would I calculate a small range in p_T and then see what values of \sigma (the integrated cross section) have a value in this range?
 
alex3 said:
Do you mean I should plot \operatorname{d}\sigma/\operatorname{d} p_T against p_T? I know I need to do that, my problem is trying to find the differential values. i.e. How could I figure out "the 'amount' of cross section which falls into a certain (small) range of p_T"?

Would I calculate a small range in p_T and then see what values of \sigma (the integrated cross section) have a value in this range?
Yeah, that sounds about right. The details vary somewhat depending on exactly what kind of data you have, so it's kind of hard for me to be specific without actually seeing your numbers (or inventing an example); unfortunately I don't have time for that right now but if you will still be working on this in a couple days, perhaps I can get back to it).
 
Hello everyone, I’m considering a point charge q that oscillates harmonically about the origin along the z-axis, e.g. $$z_{q}(t)= A\sin(wt)$$ In a strongly simplified / quasi-instantaneous approximation I ignore retardation and take the electric field at the position ##r=(x,y,z)## simply to be the “Coulomb field at the charge’s instantaneous position”: $$E(r,t)=\frac{q}{4\pi\varepsilon_{0}}\frac{r-r_{q}(t)}{||r-r_{q}(t)||^{3}}$$ with $$r_{q}(t)=(0,0,z_{q}(t))$$ (I’m aware this isn’t...
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
Back
Top