00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 // 00027 // $Id$ 00028 // 00029 // by V. Lara 00030 // 00031 // Modif (03 September 2008) by J. M. Quesada for external choice of inverse 00032 // cross section option 00033 // JMQ (06 September 2008) Also external choice has been added for: 00034 // - superimposed Coulomb barrier (if useSICB=true) 00035 // 23.08.2010 V.Ivanchenko general cleanup, move constructor and destructor 00036 // the source, make GetReactionProduct() and IsItPossible inlined 00037 00038 inline G4bool G4VPreCompoundFragment:: 00039 IsItPossible(const G4Fragment & aFragment) const 00040 { 00041 G4int pplus = aFragment.GetNumberOfCharged(); 00042 G4int pneut = aFragment.GetNumberOfParticles()-pplus; 00043 return ((pneut >= theA - theZ) && (pplus >= theZ) 00044 && (theMaximalKineticEnergy > 0)); 00045 } 00046 00047 00048 inline 00049 G4ReactionProduct * G4VPreCompoundFragment::GetReactionProduct() const 00050 { 00051 G4ReactionProduct * theReactionProduct = 00052 new G4ReactionProduct(const_cast<G4ParticleDefinition*>(particle)); 00053 theReactionProduct->SetMomentum(GetMomentum().vect()); 00054 theReactionProduct->SetTotalEnergy(GetMomentum().e()); 00055 return theReactionProduct; 00056 } 00057 00058 inline G4int G4VPreCompoundFragment::GetA() const 00059 { 00060 return theA; 00061 } 00062 00063 inline G4int G4VPreCompoundFragment::GetZ() const 00064 { 00065 return theZ; 00066 } 00067 00068 inline G4int G4VPreCompoundFragment::GetRestA() const 00069 { 00070 return theRestNucleusA; 00071 } 00072 00073 inline G4int G4VPreCompoundFragment::GetRestZ() const 00074 { 00075 return theRestNucleusZ; 00076 } 00077 00078 inline G4double G4VPreCompoundFragment::ResidualA13() const 00079 { 00080 return theRestNucleusA13; 00081 } 00082 00083 inline G4double G4VPreCompoundFragment::GetCoulombBarrier() const 00084 { 00085 return theCoulombBarrier; 00086 } 00087 00088 inline G4double G4VPreCompoundFragment::GetBindingEnergy() const 00089 { 00090 return theBindingEnergy; 00091 } 00092 00093 inline G4double G4VPreCompoundFragment::GetMaximalKineticEnergy() const 00094 { 00095 return theMaximalKineticEnergy; 00096 } 00097 00098 inline G4double G4VPreCompoundFragment::GetEnergyThreshold() const 00099 { 00100 return theMaximalKineticEnergy - theCoulombBarrier; 00101 } 00102 00103 inline G4double G4VPreCompoundFragment::GetEmissionProbability() const 00104 { 00105 return theEmissionProbability; 00106 } 00107 00108 inline G4double G4VPreCompoundFragment::GetNuclearMass(void) const 00109 { 00110 return theMass; 00111 } 00112 00113 inline G4double G4VPreCompoundFragment::GetRestNuclearMass() const 00114 { 00115 return theRestNucleusMass; 00116 } 00117 00118 inline G4double G4VPreCompoundFragment::GetReducedMass() const 00119 { 00120 return theReducedMass; 00121 } 00122 00123 inline 00124 const G4LorentzVector& G4VPreCompoundFragment::GetMomentum() const 00125 { 00126 return theMomentum; 00127 } 00128 00129 inline 00130 void G4VPreCompoundFragment::SetMomentum(const G4LorentzVector & value) 00131 { 00132 theMomentum = value; 00133 } 00134 00135 inline const G4String G4VPreCompoundFragment:: 00136 GetName() const 00137 { 00138 return particle->GetParticleName(); 00139 } 00140 00141 //for inverse cross section choice 00142 inline void G4VPreCompoundFragment::SetOPTxs(G4int opt) 00143 { 00144 OPTxs=opt; 00145 } 00146 00147 //for superimposed Coulomb Barrier for inverse cross sections 00148 inline void G4VPreCompoundFragment::UseSICB(G4bool use) 00149 { 00150 useSICB=use; 00151 }