| Summary: | Kinematics error in Geant 4.4.1 simulation | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | clundst |
| Component: | processes/hadronic/processes | Assignee: | Hans-Peter.Wellisch |
| Status: | CLOSED FIXED | ||
| Severity: | major | ||
| Priority: | P3 | ||
| Version: | 4.1 | ||
| Hardware: | PC | ||
| OS: | Linux | ||
|
Description
clundst
2002-10-22 11:31:08 CEST
I did run a set of reactions with 10^-10MeV neutron into B. The result is below. The numbers after "SECONDARIES info" are total energy, and 3-momentum. As you see, the fragments come out essentially back-to-back. May be you can put your application someplace, where I can find it, and I'll have a closer look. Many greetings, Hans-Peter. Last chance before DoIt call: NUMBER OF SECONDARIES=3FINAL STATE = 1.000000e-10 (0.000000e+00,0.000000e+00,1.000000e+00) SECONDARIES info3.728866e+03(-2.570862e+01,6.699680e+01,-7.519707e+01)3 SECONDARIES info6.534667e+03(2.595904e+01,-6.691105e+01,7.560655e+01)0 SECONDARIES info4.775995e-01(-2.579513e-01,-8.036663e-02,-3.938320e-01)-3 QVALUE = 2.758158e+00 Last chance before DoIt call: NUMBER OF SECONDARIES=3FINAL STATE = 1.000000e-10 (0.000000e+00,0.000000e+00,1.000000e+00) SECONDARIES info3.728866e+03(3.045075e+01,5.231214e+01,8.450027e+01)3 SECONDARIES info6.534655e+03(-3.019973e+01,-5.219066e+01,-8.410274e+01)0 SECONDARIES info4.775996e-01(-2.405916e-01,-1.035754e-01,-3.993611e-01)-3 QVALUE = 2.746057e+00 Last chance before DoIt call: NUMBER OF SECONDARIES=3FINAL STATE = 1.000000e-10 (0.000000e+00,0.000000e+00,1.000000e+00) SECONDARIES info3.728866e+03(5.870308e+01,7.327475e+01,-4.459672e+01)3 SECONDARIES info6.534659e+03(-5.870543e+01,-7.284930e+01,4.486167e+01)0 SECONDARIES info4.775989e-01(8.899067e-03,-4.089418e-01,-2.465525e-01)-3 QVALUE = 2.750400e+00 Last chance before DoIt call: NUMBER OF SECONDARIES=3FINAL STATE = 1.000000e-10 (0.000000e+00,0.000000e+00,1.000000e+00) SECONDARIES info3.728866e+03(4.476266e+01,-9.281023e+01,-1.366277e+01)3 SECONDARIES info6.534658e+03(-4.460716e+01,9.269016e+01,1.320569e+01)0 SECONDARIES info4.775991e-01(-1.631135e-01,9.763042e-02,4.381360e-01)-3 QVALUE = 2.749615e+00 Last chance before DoIt call: NUMBER OF SECONDARIES=3FINAL STATE = 1.000000e-10 (0.000000e+00,0.000000e+00,1.000000e+00) SECONDARIES info3.728866e+03(-6.070129e+01,4.966587e+01,-6.821115e+01)3 SECONDARIES info6.534658e+03(6.072128e+01,-4.982315e+01,6.774538e+01)0 SECONDARIES info4.775993e-01(-1.344920e-02,1.558518e-01,4.512543e-01)-3 QVALUE = 2.749860e+00 Last chance before DoIt call: NUMBER OF SECONDARIES=3FINAL STATE = 1.000000e-10 (0.000000e+00,0.000000e+00,1.000000e+00) SECONDARIES info3.728866e+03(1.022882e+01,-1.013865e+02,-2.049985e+01)3 SECONDARIES info6.534655e+03(-1.002068e+01,1.009540e+02,2.046022e+01)0 SECONDARIES info4.775999e-01(-2.208317e-01,4.232616e-01,1.358680e-02)-3 QVALUE = 2.746179e+00 Last chance before DoIt call: NUMBER OF SECONDARIES=3FINAL STATE = 1.000000e-10 (0.000000e+00,0.000000e+00,1.000000e+00) SECONDARIES info3.728866e+03(-7.864378e+01,3.812353e+01,5.626590e+01)3 SECONDARIES info6.534656e+03(7.855544e+01,-3.793079e+01,-5.585317e+01)0 SECONDARIES info4.776006e-01(9.472475e-02,-1.921710e-01,-4.268488e-01)-3 QVALUE = 2.747597e+00 Hi Carl, I have your aoolication running, and used the new tool for checking the geometry. It is just a on-line command: /geometry/test/recursive_test which looks for overlaps in the geometrical description. It tells me that there are overlaps between the silicon and the BC layers. Most of the parameters of the geant4 solids are actually half-sizes, not sizes, so may be there was some confusion. Please use the geometry-checking, and sort out the overlaps. If you then still have problems, please let me know. Many greetings, Hans-Peter. /geometry/test/recursive_test |