42 nCutMax(7),ThresholdParameter(0.000*
GeV),
43 QGSMThreshold(3*
GeV),theNucleonRadius(1.5*
fermi)
50 ThresholdParameter(right.ThresholdParameter), QGSMThreshold(right.QGSMThreshold),
51 theNucleonRadius(right.theNucleonRadius)
116 #ifdef debug_G4QGSParticipants
124 std::pair<G4double, G4double> theImpactParameter;
126 G4double impactX = theImpactParameter.first;
127 G4double impactY = theImpactParameter.second;
131 G4int nucleonCount = 0;
133 #ifdef debug_G4QGSParticipants
134 intNuclMom=aPrimaryMomentum;
146 G4double s_nucleon = (aPrimaryMomentum + nucleonMomentum).mag2();
149 G4double Probability = theProbability.GetInelasticProbability(s_nucleon, Distance2);
153 if (Probability > rndNumber)
155 #ifdef debug_G4QGSParticipants
156 G4cout <<
"DEBUG p="<< Probability<<
" r="<<rndNumber<<
" d="<<std::sqrt(Distance2)<<
G4endl;
157 G4cout <<
" qgspart+ " << aPrimaryMomentum <<
" # " << aPrimaryMomentum.mag()
159 intNuclMom += nucleonMomentum;
165 pNucleon->
Hit(aTarget);
166 if ((theProbability.GetDiffractiveProbability(s_nucleon, Distance2)/Probability >
G4UniformRand()
191 for(nCut = 0; nCut <
nCutMax; nCut++)
193 running[nCut] = theProbability.GetCutPomeronProbability(s_nucleon, Distance2, nCut + 1);
194 if(nCut!=0) running[nCut] += running[nCut-1];
197 for(nCut = 0; nCut <
nCutMax; nCut++)
199 if(running[nCut] > random)
break;
210 #ifdef debug_G4QGSParticipants
211 impactUsed=Distance2;
216 #ifdef debug_G4QGSParticipants
221 #ifdef debug_G4QGSParticipants
223 G4cout <<
"Impact Parameter used = "<<impactUsed<<
G4endl;
244 #ifdef debug_G4QGSPart_PDiffColl
262 #ifdef debug_G4QGSPart_PDiffColl
277 std::vector<G4InteractionContent*>::iterator i;
292 #ifdef debug_G4QGSPart_PSoftColl
300 #ifdef debug_G4QGSParticipants
307 #ifdef debug_G4QGSPart_PSoftColl
315 #ifdef debug_G4QGSParticipants
327 #ifdef debug_G4QGSPart_PSoftColl
G4bool ExciteParticipants(G4VSplitableHadron *aPartner, G4VSplitableHadron *bPartner) const
G4Parton * GetParton2(void)
void SetMomentum(G4LorentzVector &aMomentum)
virtual G4bool StartLoop()=0
const G4LorentzVector & Get4Momentum() const
const G4double QGSMThreshold
virtual const G4LorentzVector & Get4Momentum() const
virtual G4Parton * GetNextAntiParton()=0
virtual const G4ThreeVector & GetPosition() const
G4bool IsSingleDiffractive()
G4ParticleDefinition * GetDefinition() const
G4Parton * GetParton1(void)
virtual G4double GetOuterRadius()=0
virtual ~G4QGSParticipants()
G4GLOB_DLL std::ostream G4cout
virtual G4VSplitableHadron * SelectInteractions(const G4ReactionProduct &thePrimary)
virtual G4Parton * GetNextParton()=0
void PerformDiffractiveCollisions()
void IncrementCollisionCount(G4int aCount)
G4QGSDiffractiveExcitation theDiffExcitaton
virtual G4bool ExciteParticipants(G4VSplitableHadron *aPartner, G4VSplitableHadron *bPartner) const
G4V3DNucleus * theNucleus
G4double GetKineticEnergy() const
std::pair< G4double, G4double > ChooseImpactXandY(G4double maxImpact)
G4double GetTotalEnergy() const
void SetNumberOfDiffractiveCollisions(int)
std::vector< G4PartonPair * > thePartonPairs
G4double GetPDGMass() const
void SetNumberOfSoftCollisions(int)
std::vector< G4VSplitableHadron * > theTargets
G4VSplitableHadron * GetTarget() const
G4VSplitableHadron * GetProjectile() const
void BuildInteractions(const G4ReactionProduct &thePrimary)
G4ThreeVector GetMomentum() const
virtual G4ParticleDefinition * GetDefinition() const
virtual G4Nucleon * GetNextNucleon()=0
void PerformSoftCollisions()
void Hit(G4VSplitableHadron *aHit)
G4SingleDiffractiveExcitation theSingleDiffExcitation
G4double GetBindingEnergy() const
std::vector< G4InteractionContent * > theInteractions
const G4double theNucleonRadius
const G4double ThresholdParameter
static const double fermi
G4ThreadLocal G4int G4QGSParticipants_NPart
G4int GetNumberOfSoftCollisions()
void SetTarget(G4VSplitableHadron *aTarget)
CLHEP::HepLorentzVector G4LorentzVector