Problem 728

Summary: G4HadronicProcess failed in ApplyYourself call for - Particle energy[GeV] = 0.00613 - Material = Lsim - Particle type = pi+
Product: Geant4 Reporter: timmyf
Component: processes/hadronic/modelsAssignee: Hans-Peter.Wellisch
Status: CLOSED REMIND    
Severity: major    
Priority: P1    
Version: 6.0   
Hardware: PC   
OS: Linux   

Description timmyf 2005-02-27 21:08:27 CET
Dear Sir,
  I am Timmy Francis, Project Associate in the Indian Moon Mission Program,
Chandrayaan -1. My group is trying to simulate the lunar gamma ray back ground
continuum   in the Energy region of our interest (20-250 keV ), using Geant4.6.0.
  Our detector consists of a cylinder filled with material having average lunar
composition. We are firing protons with a given input flux from random positions
from a hemisphere such that protons are hitting uniformly on a 1 sq.cm area on
the top face of the cylinder which is filled with the material having lunar
composition.
  In the physics processes simulated, we have included hadronic  processes also
and since we found that the hadronic processes simulated in the LISA example
which is given in geant4.6.2/examples, are satisfying our requirement we have
simply copied that hadronic processes to our PhysicsList.
  But while the energy of the protons fired reaches about .460GeV, the program
is aborting its execution, with the following lines:
  Event: 434272
In
/root/geant4.6.0.p01/source/processes/hadronic/models/util/include/G4Fragment.hh,
line 261:
===> G4Fragment::SetNumberOfCharged: Number of charged particles can't be
greater than number of particles
 G4HadronicProcess failed in ApplyYourself call for
 - Particle energy[GeV] = 0.00613
 - Material = Lsim
 - Particle type = pi+

*** G4Exception : 007
      issued by : G4HadronicProcess
GeneralPostStepDoIt failed.
*** Fatal Exception *** core dump ***

*** G4Exception: Aborting execution ***
Aborted

   Sir, there is another problem which we are facing.
 Once we go to higer energies, (a few GeVs), the program is aborting many events
reporting that the momentum direction was illegal.
   Sir, In our DetectorConstruction.cc, our fractional mass is not added up to
1, it is less than 1 by some .00000001. Could this be a reason for  the above
mentioned problems.
  Kindly send me a reply at your earliest convenience.

Thanking You,
Yours faithfully
 TIMMY FRANCIS.
 Project Associate, PLANEX
 Physical Research Laboratory
 India.

 E-mail:timmyf@prl.ernet.in
       :timmyf@rediffmail.com
Comment 1 Hans-Peter.Wellisch 2005-07-05 06:56:59 CEST
At the first glance this looks like there is an inconsistancy between cross-
sections and final state generators used in the LISA physics list. Can you find
a physics list within the educated guess physics lists, that would fulfill your
requirements? It should resolve the problem.

With kind regards,

Hans-Peter.