115 G4cout <<
"G4ChargeExchange::DoIt: Incident particle plab="
116 << plab/
GeV <<
" GeV/c "
117 <<
" ekin(MeV) = " << ekin/
MeV <<
" "
127 <<
" PDGcode= " << projPDG <<
" on nucleus Z= " << Z
128 <<
" A= " << A <<
" N= " << N
143 G4bool theHyperon =
false;
159 }
else if(theParticle ==
theKPlus) {
161 else theSecondary =
theK0L;
165 else theSecondary =
theK0L;
167 }
else if(theParticle ==
theK0S || theParticle ==
theK0L) {
181 }
else if(theParticle ==
theL) {
185 theSecondary =
theS0;
186 }
else if (x < 0.4) {
189 }
else if (x < 0.6) {
194 }
else if (x < 0.8) {
207 theSecondary =
theS0;
208 }
else if (x < 0.4) {
211 }
else if (x < 0.6) {
216 }
else if (x < 0.8) {
230 if (Z == 1 && A == 2) theDef =
theD;
231 else if (Z == 1 && A == 3) theDef =
theT;
232 else if (Z == 2 && A == 3) theDef =
theHe3;
233 else if (Z == 2 && A == 4) theDef =
theA;
242 if(theRecoil) { m21 += theRecoil->
GetPDGMass(); }
243 else { theRecoil = theDef; }
248 if(etot < m11 + m21) {
255 G4double e1 = 0.5*etot*(1.0 - (m21*m21 - m11*m11)/(etot*etot));
257 G4double ptot = std::sqrt(e1*e1 - m11*m11);
265 G4cout <<
"## G4ChargeExchange t= " << t <<
" tmax= " << tmax
266 <<
" ptot= " << ptot <<
G4endl;
271 if(std::abs(cost) > 1.0) cost = 1.0;
272 G4double sint = std::sqrt((1.0-cost)*(1.0+cost));
277 G4ThreeVector v1(sint*std::cos(phi),sint*std::sin(phi),cost);
303 if(erec < 0.0) erec = 0.0;
330 const G4int maxNumberOfLoops = 10000;
331 G4int loopCounter = 0;
334 }
while ( (t > tmax) &&
335 ++loopCounter < maxNumberOfLoops );
336 if ( loopCounter >= maxNumberOfLoops ) {
static G4Pow * GetInstance()
G4ParticleDefinition * theSMinus
G4double powA(G4double A, G4double y) const
static G4double GetNuclearMass(const G4double A, const G4double Z)
G4ParticleDefinition * theL
G4ParticleDefinition * theOmega
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
G4ParticleDefinition * theASPlus
static G4AntiOmegaMinus * AntiOmegaMinus()
G4ParticleDefinition * theXiMinus
CLHEP::Hep3Vector G4ThreeVector
G4double lowEnergyRecoilLimit
G4ParticleDefinition * theNeutron
G4ParticleDefinition * theAOmega
G4ParticleDefinition * theS0
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
static G4OmegaMinus * OmegaMinus()
G4ParticleDefinition * theAProton
static G4KaonZeroLong * KaonZeroLong()
G4ParticleDefinition * theAS0
G4ParticleDefinition * theT
G4ParticleDefinition * theA
G4int GetPDGEncoding() const
G4ParticleDefinition * theKMinus
static G4AntiSigmaPlus * AntiSigmaPlus()
G4ParticleDefinition * theK0S
const G4String & GetParticleName() const
static G4SigmaZero * SigmaZero()
G4ParticleDefinition * theAXi0
static G4KaonMinus * KaonMinus()
void SetStatusChange(G4HadFinalStateStatus aS)
static G4AntiSigmaMinus * AntiSigmaMinus()
static G4XiZero * XiZero()
void SetMinEnergy(G4double anEnergy)
G4IonTable * GetIonTable() const
G4GLOB_DLL std::ostream G4cout
G4ParticleDefinition * theK0L
double A(double temperature)
const G4ParticleDefinition * GetDefinition() const
static G4KaonZeroShort * KaonZeroShort()
static G4AntiProton * AntiProton()
G4double GetKineticEnergy() const
static const double twopi
static G4XiMinus * XiMinus()
static G4AntiXiMinus * AntiXiMinus()
static G4Triton * Triton()
static G4Proton * Proton()
static G4PionPlus * PionPlus()
G4ParticleDefinition * theXi0
static G4Neutron * Neutron()
G4ParticleDefinition * theKPlus
G4ParticleDefinition * theASMinus
const G4LorentzVector & Get4Momentum() const
static G4PionZero * PionZero()
static G4Deuteron * Deuteron()
void SetKineticEnergy(G4double aEnergy)
G4ParticleDefinition * theANeutron
static G4SigmaMinus * SigmaMinus()
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
void SetEnergyChange(G4double anEnergy)
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
G4double SampleT(G4double p, G4double A)
static G4AntiLambda * AntiLambda()
G4ParticleDefinition * theD
static G4PionMinus * PionMinus()
virtual ~G4ChargeExchange()
static G4AntiSigmaZero * AntiSigmaZero()
const G4double x[NPOINTSGL]
G4ParticleDefinition * thePiPlus
void SetMaxEnergy(const G4double anEnergy)
G4HadFinalState theParticleChange
void SetLocalEnergyDeposit(G4double aE)
static G4AntiXiZero * AntiXiZero()
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
G4ParticleDefinition * thePiZero
static G4SigmaPlus * SigmaPlus()
G4ParticleDefinition * theAXiMinus
static G4Lambda * Lambda()
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
static G4KaonPlus * KaonPlus()
G4ParticleDefinition * theSPlus
void SetMomentumChange(const G4ThreeVector &aV)
G4double lowestEnergyLimit
static G4AntiNeutron * AntiNeutron()
G4double GetTotalMomentum() const
G4ParticleDefinition * theAntiL
CLHEP::HepLorentzVector G4LorentzVector
G4ParticleDefinition * theProton
G4ParticleDefinition * thePiMinus
G4ParticleDefinition * theHe3