Message: Re: Energy of secondary nuclei produced in G4.9.4 inelastic hadronic processes Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

More Re: Energy of secondary nuclei produced in G4.9.4 inelastic hadronic processes 

Forum: Hadronic Processes
Re: Question Energy of secondary nuclei produced in G4.9.4 inelastic hadronic processes (Daniel Matthiń)
Re: Idea Re: Energy of secondary nuclei produced in G4.9.4 inelastic hadronic processes (Vladimir Ivanchenko)
Date: 07 Jan, 2011
From: Daniel Matthiń <Daniel Matthiń>

Dear Vladimir,
the output is created by the application that I am using (PLANETOCOSMICS) for the transport of cosmic ray particles through the atmosphere.
I found out that it is the Bertini model which is responsible for the discrepancies between the total energy of a particle and the sum of its kinetic energy and rest mass.
If I use the QGSP_BIC_HP list instead of QGSP_BERT_HP or if I use GEANT4.9.3.p02 the warnings disappear.
Primary particles are protons with 1GeV energy.
Here is the part of the code which produces the message:

 const G4DynamicParticle* thePart =aTrack->GetDynamicParticle();
 double E =thePart->GetTotalEnergy();
 double p =thePart->GetMomentum().mag();
 double Ekin = thePart->GetKineticEnergy();
 double m0= thePart->GetDefinition()->GetPDGMass();
 if (std::abs(E-m0-Ekin) >100.*eV ) {
  if (bug_verbose >0){
    G4cout<<"Warning a particle has been produced with E!=Ekin+m0";
    G4cout<<"particle: "
         << thePart->GetDefinition()->GetParticleName()<<std::endl;
    G4cout<<"Creator process:"
         << aTrack->GetCreatorProcess()->GetProcessName()<<std::endl;
      
    G4cout<<"E[MeV] Ekin P[MeV/c] m0 E-Ekin-E0: "
         <<E/MeV<<'\t'
                                 <<Ekin/MeV<<'\t'
   <<p/MeV<<'\t'
         <<m0/MeV<<'\t'
         <<(E-Ekin-m0)/MeV<<std::endl;
            
  }
 }
 
 

Regards,
Daniel

 Add Message Add Message
to: "Re: Energy of secondary nuclei produced in G4.9.4 inelastic hadronic processes"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]