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
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
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.
