59    for ( 
int i = 0 ; i < nSecondaries ; i++ )
    93    if( getenv(
"G4ParticleHPDebug")) 
G4cout  << 
"G4ParticleHPFinalState::adjust_final_stat  BaseZ " << baseZNew << 
" BaseA " << baseANew << 
" sum_Z " << sum_Z << 
" sum_A " << sum_A << 
G4endl;
    96    G4bool needOneMoreSec = 
false;
    98    if ( (
int)(baseZNew - sum_Z) == 0 && (
int)(baseANew - sum_A) == 0 )
   105       if ( max_SecA > 
int(baseANew - sum_A) ) 
   109          needOneMoreSec = 
true;
   117       if ( needOneMoreSec )
   119          if ( 
int(baseZNew - sum_Z) == 0 && (
int)(baseANew - sum_A) > 0 )
   121             if ( 
int(baseANew - sum_A) > 1 ) 
G4cout << 
"More than one neutron is required for the balance of baryon number!" << 
G4endl;
   127        if( getenv(
"G4ParticleHPDebug")) 
G4cout << 
this << 
"G4ParticleHPFinalState oneMoreSec_pd Z " << baseZNew << 
" - " << sum_Z << 
" A " << baseANew << 
" - " << sum_A << 
" projectile " << 
theProjectile->
GetParticleName() << 
G4endl;
   130        if( !oneMoreSec_pd ) {
   132          G4Exception(
"G4ParticleHPFinalState:adjust_final_state",
   135              "No adjustment will be done!");
   141       if ( resi_pd == NULL )
   162         if ( (
int)(ndlZNew - sum_Z) == 0 && (
int)(ndlANew - sum_A) == 0 )
   169          G4Exception(
"G4ParticleHPFinalState:adjust_final_state",
   172              "No adjustment will be done!");
   176            for ( 
int i = 0 ; i < nSecondaries ; i++ )
   201    for ( 
G4int i = 0 ; i < n_sec ; i++ )
   224       if ( slow > beta && beta != 0 ) 
   268       if ( slow > p4.beta() ) 
   273       dif_4p = init_4p_lab - ( secs_4p_lab + p4 );  
   276    if ( needOneMoreSec && oneMoreSec_pd)
   288       if ( slow > p4.beta() ) 
   291          islow = nSecondaries-1; 
   293       dif_4p = init_4p_lab - ( secs_4p_lab + p4 );  
   298    if ( dif_4p.
v().
mag() < std::abs( dif_4p.
e() ) )
   303       if ( minimum_energy < dif_4p.
v().
mag() && dif_4p.
v().
mag() < 1*
MeV )
   331    G4double dif_e = dif_4p.
e() - ( dif_4p.
v() ).mag();
   342       if ( minimum_energy < e1 )  
   347                             std::sin(std::acos(costh))*std::sin(phi),
   368       if ( ke0 + dif_e > 0 )
 static G4ParticleHPManager * GetInstance()
 
G4Cache< G4HadFinalState *> theResult
 
void SetMomentum(const G4ThreeVector &momentum)
 
G4int GetNumberOfSecondaries() const
 
G4HadSecondary * GetSecondary(size_t i)
 
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
 
G4ParticleDefinition * theProjectile
 
G4double GetKineticEnergy() const
 
G4ThreeVector GetMomentum() const
 
const G4String & GetParticleName() const
 
G4GLOB_DLL std::ostream G4cout
 
static const double twopi
 
static G4Triton * Triton()
 
static G4Proton * Proton()
 
static G4Neutron * Neutron()
 
static G4Deuteron * Deuteron()
 
void SetKineticEnergy(G4double aEnergy)
 
static G4IonTable * GetIonTable()
 
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
 
G4int GetAtomicNumber() const
 
G4double GetPDGMass() const
 
G4DynamicParticle * GetParticle()
 
G4LorentzVector Get4Momentum() const
 
G4ParticleDefinition * GetDefinition() const
 
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
 
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
 
G4int GetAtomicMass() const
 
G4GLOB_DLL std::ostream G4cerr