Problem 562

Summary: G4ParticleChange::CheckIt
Product: Geant4 Reporter: galante
Component: processes/hadronicAssignee: Hans-Peter.Wellisch
Status: RESOLVED FIXED    
Severity: critical    
Priority: P1    
Version: 5.0   
Hardware: PC   
OS: Linux   
URL: http://hep.fi.infn.it/PAMELA/galante.tar

Description galante 2004-01-05 09:33:03 CET
Hello,
I'm simulating the cosmic ray flux in a magnetic spectrometer and i obtain the
following error:
G4ParticleChange::CheckIt  : the Momentum Change is not unit vector !!
  Difference:  1
 G4ParticleChange::CheckIt
      -----------------------------------------------
        G4ParticleChange Information
      -----------------------------------------------
        # of 2ndaries       :                    0
      -----------------------------------------------
        Energy Deposit (MeV):                    0
        Track Status        :                Alive
        True Path Length (mm) :                0.929
        Stepping Control     :                    0

        Mass (GeV)   :                 1.12
        Charge (eplus)   :                   -1
        Position - x (mm)   :                  105
        Position - y (mm)   :                -53.1
        Position - z (mm)   :                  114
        Time (ns)           :             6.69e+06
        Proper Time (ns)    :             6.69e+06
        Momentum Direct - x :                    0
        Momentum Direct - y :                    0
        Momentum Direct - z :                    0
        Kinetic Energy (MeV):             6.04e-18
        Polarization - x    :                    0
        Polarization - y    :                    0
        Polarization - z    :                    0
        Track Weight      :                    1


The pid of the particle is 3122(lambda)
I really don't understand what's wrong...
Can you help me please? It's really urgent because this is my last job for the
degree and i have to finish before 26th of january...
U can find my code at the following address:

http://hep.fi.infn.it/PAMELA/galante.tar

Thanks a lot
Annamaria Galante
Comment 1 galante 2004-01-07 12:20:59 CET
Hello,
as i'm intersted only in primary protons that i generate i deleted the physics
about lambda and now it seems to work..
Is it correct?
Anyway i simulated 500000 protons withouth energy cut :
in trackerPrimaryGenerationAction.cc i put

 while(abs(x+sin(theta)*cos(phi)*l)>x_max ||
abs(y+sin(theta)*sin(phi)*l)>y_max);

instead of

while(abs(x+sin(theta)*cos(phi)*l)>x_max || abs(y+sin(theta)*sin(phi)*l)>y_max
|| en_k>=2000*GeV);

and i obteined the following error:

GetHadronicInteraction: No model found for this energy range

*** G4Exception: Aborting execution ***
Aborted


What's wrong now??
Please can you help?
Thanks a lot,
Annamaria Galante
Comment 2 Hans-Peter.Wellisch 2004-01-23 10:06:59 CET
Hi,

   the problem with the lambda should be gone in geant4 6.0.
Otherwise, when you get the message you see, it means you use an energy
for which no model was registered. This can happen beyone 10 TeV.

Many greetings,

Hans-Peter.