Problem 460 - Problem on G4LElastic
Summary: Problem on G4LElastic
Status: CLOSED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models/low_energy (show other problems)
Version: other
Hardware: PC Linux
: P1 normal
Assignee: Hans-Peter.Wellisch
URL:
Depends on:
Blocks:
 
Reported: 2003-02-24 07:57 CET by laurent.desorgher
Modified: 2012-02-15 09:56 CET (History)
1 user (show)

See Also:


Attachments

Note You need to log in before you can comment on or make changes to this problem.
Description laurent.desorgher 2003-02-24 07:57:15 CET
Dear Hans-Peter,

When simulating in geant4.5.0  Elastic scattering of particle on  material
heavier than hydrogen (and probably also on hydrogen) a bug occurs if the
incoming particle has exacltly the direction (0,0,-1.).  In this case the
particle scatters backward  while a  forward direction is expected.
The bug occurs in the file G4LElastic.cc, in the method  G4LElastic::Defs1 in
line 373-377.
As the code  is written (0,0,1) and (0,0,-1.) incident direction are consider as
equivalent to (0,0,1).

To solve the bug I propose to change the lines 373-377 by

else {
       G4double cost=pzinc/p;
       *pxnew = cost*px;
       *pynew = py;
       *pznew = cost*pz;
   }



Best Regards



Laurent Desorgher
Comment 1 Hans-Peter.Wellisch 2003-05-25 12:10:59 CEST
was included already some time ago. Closing the 'ticket' now.