86 G4double cmsMom = std::sqrt(the3CMS*the3CMS);
87 G4double sqrts = std::sqrt((totE-cmsMom)*(totE+cmsMom));
91 theNeutron.
Lorentz(theNeutron, theCMS);
92 theTarget.
Lorentz(theTarget, theCMS);
98 tempVector.
setX( std::cos(theta)*std::sin(cms_theta)*std::cos(cms_phi)
99 +std::sin(theta)*std::cos(phi)*std::cos(cms_theta)*std::cos(cms_phi)
100 -std::sin(theta)*std::sin(phi)*std::sin(cms_phi) );
101 tempVector.
setY( std::cos(theta)*std::sin(cms_theta)*std::sin(cms_phi)
102 +std::sin(theta)*std::cos(phi)*std::cos(cms_theta)*std::sin(cms_phi)
103 +std::sin(theta)*std::sin(phi)*std::cos(cms_phi) );
104 tempVector.
setZ( std::cos(theta)*std::cos(cms_theta)
105 -std::sin(theta)*std::cos(phi)*std::sin(cms_theta) );
114 theNeutron.
Lorentz(theNeutron, -1.*theCMS);
122 theTarget.
Lorentz(theTarget, -1.*theCMS);
129 theTarget.
Lorentz(theTarget, -1.*theCMS);
ThreeVector shoot(const G4int Ap, const G4int Af)
void SetMomentum(const G4ThreeVector &momentum)
G4double GetTotalMomentum() const
void Lorentz(const G4ReactionProduct &p1, const G4ReactionProduct &p2)
void SetKineticEnergy(const G4double en)
void SetMomentum(const G4double x, const G4double y, const G4double z)
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
G4int GetNucleusEncoding(G4int iZ, G4int iA, G4int iM)
void SetMass(const G4double mas)
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
void SetTotalEnergy(const G4double en)
G4ErrorTarget * theTarget
G4ParticleDefinition * FindIon(G4int atomicNumber, G4int atomicMass, G4double excitationEnergy)
const G4LorentzVector & Get4Momentum() const
G4double GetKineticEnergy() const
void SetEnergyChange(G4double anEnergy)
G4double GetTotalEnergy() const
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
const G4Isotope * GetIsotope()
G4ThreeVector GetMomentum() const
G4HadFinalState theParticleChange
G4double GetTemperature() const
const G4Material * GetMaterial() const
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
G4LENDManager * lend_manager
void SetMomentumChange(const G4ThreeVector &aV)
void AddSecondary(G4DynamicParticle *aP)