Re: Changing t distribution slopes [message #12533 is a reply to message #12532] |
Sun, 21 August 2011 10:35 |
Ingo Froehlich
Messages: 167 Registered: March 2004 Location: IKF - Frankfurt
|
first-grade participant |
From: *dip.t-dialin.net
|
|
Hi,
the method is based on sampling (but with phase space only) and rejection (pure Monte-Carlo). For this one needs a ceiling value. The rejected events are those, for which a random number between 0 and the ceiling value is higher as the function. If the value is larger as the ceiling, the ceiling is adjusted. If this happens during the first events, it does not play a big role. But in your case this is too much.
Maybe it would be helpful if you could attach your macro. It could be that it depends on the beam energy. I tested it at T_kin=2.2 GeV only. There the slope was OK over some kEvents.
--
Ingo Froehlich
IKF - University of Frankfurt
069-798-47027, FAX: -47024
|
|
|