Home » Hades » Pluto » [SOLVED] Beam Smear and Cross Sections
Re: Bear Smear and Cross Sections [message #13893 is a reply to message #13892] |
Mon, 27 August 2012 22:59 |
Ingo Froehlich
Messages: 167 Registered: March 2004 Location: IKF - Frankfurt
|
first-grade participant |
From: *dip.t-dialin.net
|
|
Michael Kunkel wrote on Mon, 27 August 2012 21:07 | Is _f the density function? If so, wouldn't using Input : _x s cos(theta), _y is differential cross section
Output : cross section suffice?
|
No, it's a very simple implementation: there are 2 input variables, mapped on the 2 axis of TF2, and one output variable, named _f, and this is the return value of the TF2. How you fill the _f is up to you, but it must be defined.
Michael Kunkel wrote on Mon, 27 August 2012 21:07 |
I ask this, because I have 64 different models I want to use, each a segment of c.m. energy covering the range 1.68 - 2.84 GeV in the c.m. frame.
|
Than maybe the simplest method is that you use a 2-dimensional histogram, if the binning is the same for all the models, it should be no problem to merge them. Keep in mind that there is no linear interpolation implemented yet.
Michael Kunkel wrote on Mon, 27 August 2012 21:07 |
Would you please elaborate more on how to add more histograms into the model to ensure I get the right topology with the scaleing of the cross sections, the beam smearing, etc. I have having a hard time understanding everything so far.
|
You can concat the calculation:
model->AddHistogram(distribution1,"_f = Eval(_x. _y);");
model->AddHistogram(distribution2,"_f = _f * Eval(_x. _y);");
or
model->AddHistogram(distribution3,"_f = _f * Eval(_y);");
as indicated in the example, you can also make transformations:
model->AddEquation("t_lab = (_y*_y - g.mass*g.mass - p.mass*p.mass)/(2*p.mass*p.mass) - g.mass;");
model->AddHistogram(profile,"_f = _f * Eval(t_lab);");
...just in case that your histogram profile is a function of t_lab.
In the same way you can replace a beam profile by an analytical function.
--
Ingo Froehlich
IKF - University of Frankfurt
069-798-47027, FAX: -47024
[Updated on: Mon, 27 August 2012 23:01] Report message to a moderator
|
|
|
Goto Forum:
Current Time: Thu Dec 12 22:38:05 CET 2024
Total time taken to generate the page: 0.00800 seconds
|