55 default_evaluation =
"ENDF.B-VII.0";
66 for ( std::map< G4int , G4LENDUsedTarget* >::iterator
77 for ( std::map< G4int , G4LENDUsedTarget* >::iterator
98 for (
size_t i = 0 ; i < numberOfElements ; ++i )
101 const G4Element* anElement = (*theElementTable)[i];
104 if ( numberOfIsotope > 0 )
107 for (
G4int i_iso = 0 ; i_iso < numberOfIsotope ; i_iso++ )
114 if ( allow_nat ==
true ) aTarget->
AllowNat();
115 if ( allow_any ==
true ) aTarget->
AllowAny();
127 for (
G4int ii = 0 ; ii < numberOfNistIso ; ii++ )
137 if ( allow_nat ==
true ) aTarget->
AllowNat();
138 if ( allow_any ==
true ) aTarget->
AllowAny();
151 G4cout <<
"Requested Evaluation, Z , A -> Actual Evaluation, Z , A(0=Nat) , Pointer of Target" <<
G4endl;
152 for ( std::map< G4int , G4LENDUsedTarget* >::iterator
156 <<
" " << it->second->GetWantedEvaluation()
157 <<
", " << it->second->GetWantedZ()
158 <<
", " << it->second->GetWantedA()
159 <<
" -> " << it->second->GetActualEvaluation()
160 <<
", " << it->second->GetActualZ()
161 <<
", " << it->second->GetActualA()
162 <<
", " << it->second->GetTarget()
224 G4double cmsMom = std::sqrt(the3CMS*the3CMS);
225 G4double sqrts = std::sqrt((totE-cmsMom)*(totE+cmsMom));
229 theNeutron.
Lorentz(theNeutron, theCMS);
230 theTarget.
Lorentz(theTarget, theCMS);
236 tempVector.
setX(std::cos(theta)*std::sin(cms_theta)*std::cos(cms_phi)
237 +std::sin(theta)*std::cos(phi)*std::cos(cms_theta)*std::cos(cms_phi)
238 -std::sin(theta)*std::sin(phi)*std::sin(cms_phi) );
239 tempVector.
setY(std::cos(theta)*std::sin(cms_theta)*std::sin(cms_phi)
240 +std::sin(theta)*std::cos(phi)*std::cos(cms_theta)*std::sin(cms_phi)
241 +std::sin(theta)*std::sin(phi)*std::cos(cms_phi) );
242 tempVector.
setZ(std::cos(theta)*std::cos(cms_theta)
243 -std::sin(theta)*std::cos(phi)*std::sin(cms_theta) );
250 theNeutron.
Lorentz(theNeutron, -1.*theCMS);
251 theTarget.
Lorentz(theTarget, -1.*theCMS);
size_t GetNumberOfIsotopes() const
ThreeVector shoot(const G4int Ap, const G4int Af)
void SetMomentum(const G4ThreeVector &momentum)
G4double GetTotalMomentum() const
void Lorentz(const G4ReactionProduct &p1, const G4ReactionProduct &p2)
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
void SetKineticEnergy(const G4double en)
void SetMomentum(const G4double x, const G4double y, const G4double z)
G4NistElementBuilder * GetNistElementBuilder()
void recreate_used_target_map()
G4double GetIsotopeAbundance(G4int Z, G4int N) const
G4int GetNumberOfNistIsotopes(G4int Z) const
G4ParticleDefinition * proj
const G4String & GetModelName() const
G4LENDModel(G4String name="LENDModel")
G4int GetNucleusEncoding(G4int iZ, G4int iA, G4int iM)
void SetMinEnergy(G4double anEnergy)
void SetMass(const G4double mas)
G4GLOB_DLL std::ostream G4cout
static size_t GetNumberOfElements()
const G4ParticleDefinition * GetDefinition() const
double getElasticFinalState(double e_in, double temperature, double(*rng)(void *), void *rngState)
std::map< G4int, G4LENDUsedTarget * > usedTarget_map
G4double GetKineticEnergy() const
G4int GetNistFirstIsotopeN(G4int Z) const
void SetTotalEnergy(const G4double en)
G4ErrorTarget * theTarget
G4ParticleDefinition * FindIon(G4int atomicNumber, G4int atomicMass, G4double excitationEnergy)
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
const G4LorentzVector & Get4Momentum() const
G4double GetKineticEnergy() const
void SetEnergyChange(G4double anEnergy)
G4double GetTotalEnergy() const
static G4LENDManager * GetInstance()
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
void create_used_target_map()
const G4Isotope * GetIsotope()
void SetMaxEnergy(const G4double anEnergy)
G4ThreeVector GetMomentum() const
const G4Isotope * GetIsotope(G4int iso) const
G4double GetTemperature() const
const G4Material * GetMaterial() const
std::vector< G4Element * > G4ElementTable
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
G4bool RequestChangeOfVerboseLevel(G4int)
G4LENDManager * lend_manager
static G4ElementTable * GetElementTable()
void SetMomentumChange(const G4ThreeVector &aV)
void AddSecondary(G4DynamicParticle *aP)