Geant4-11
Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
G4CollisionMesonBaryonToResonance Class Reference

#include <G4CollisionMesonBaryonToResonance.hh>

Inheritance diagram for G4CollisionMesonBaryonToResonance:
G4CollisionComposite G4VCollision

Public Member Functions

void AddComponent (G4VCollision *aC)
 
void establish_G4MT_TLS_G4VCollision ()
 
virtual G4KineticTrackVector * FinalState (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
 G4CollisionMesonBaryonToResonance ()
 
virtual const G4VAngularDistribution * GetAngularDistribution () const
 
virtual const G4CollisionVector * GetComponents () const
 
virtual const G4VCrossSectionSource * GetCrossSectionSource () const
 
virtual G4bool IsInCharge (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
G4bool operator!= (const G4VCollision &right) const
 
G4bool operator== (const G4VCollision &right) const
 
virtual void Print () const
 
virtual void Print (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
virtual ~G4CollisionMesonBaryonToResonance ()
 

Protected Member Functions

virtual const std::vector< G4String > & GetListOfColliders (G4int) const
 
G4int GetNumberOfPartons (const G4ParticleDefinition *aP) const
 

Private Member Functions

void BufferCrossSection (const G4ParticleDefinition *aP, const G4ParticleDefinition *bP)
 
G4double BufferedCrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
virtual G4double CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
 G4CollisionMesonBaryonToResonance (const G4CollisionMesonBaryonToResonance &)
 
virtual G4String GetName () const
 
G4CollisionMesonBaryonToResonance & operator= (const G4CollisionMesonBaryonToResonance &)
 

Private Attributes

G4Mutex bufferMutex
 
G4CollisionVector components
 
std::vector< G4CrossSectionBuffer > theBuffer
 
G4XpimNTotal thepimp
 
G4XpipNTotal thepipp
 

Static Private Attributes

static const G4int nPoints = 32
 
static const G4double theT []
 

Detailed Description

Definition at line 44 of file G4CollisionMesonBaryonToResonance.hh.

Constructor & Destructor Documentation

◆ G4CollisionMesonBaryonToResonance() [1/2]

G4CollisionMesonBaryonToResonance::G4CollisionMesonBaryonToResonance ( )

Definition at line 47 of file G4CollisionMesonBaryonToResonance.cc.

48{
49 // the particle types used are indicative for the particle class;
50 // i.e. pip stands for pions, proton for nucleon, etc..
51
52 G4ParticleDefinition * aProton = G4Proton::ProtonDefinition();
53 G4ParticleDefinition * aPionp = G4PionPlus::PionPlusDefinition();
54
55 G4ParticleDefinition * aDeltap = G4ParticleTable::GetParticleTable()->FindParticle(2214); // D+
56 G4ParticleDefinition * aD1600 = G4ParticleTable::GetParticleTable()->FindParticle(32214); // D+
57 G4ParticleDefinition * aD1620 = G4ParticleTable::GetParticleTable()->FindParticle(2122); // D+
58 G4ParticleDefinition * aD1700 = G4ParticleTable::GetParticleTable()->FindParticle(12214); // D+
59 G4ParticleDefinition * aD1900 = G4ParticleTable::GetParticleTable()->FindParticle(12122); // D+
60 G4ParticleDefinition * aD1905 = G4ParticleTable::GetParticleTable()->FindParticle(2126); // D+
61 G4ParticleDefinition * aD1910 = G4ParticleTable::GetParticleTable()->FindParticle(22122); // D+
62 G4ParticleDefinition * aD1920 = G4ParticleTable::GetParticleTable()->FindParticle(22214); // D+
63 G4ParticleDefinition * aD1930 = G4ParticleTable::GetParticleTable()->FindParticle(12126); // D+
64 G4ParticleDefinition * aD1950 = G4ParticleTable::GetParticleTable()->FindParticle(2218); // D+
65
66 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aDeltap, "D1232_Npi"));
67 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aD1600, "D1600_Npi"));
68 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aD1620, "D1620_Npi"));
69 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aD1700, "D1700_Npi"));
70 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aD1900, "D1900_Npi"));
71 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aD1905, "D1905_Npi"));
72 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aD1910, "D1910_Npi"));
73 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aD1920, "D1920_Npi"));
74 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aD1930, "D1930_Npi"));
75 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aD1950, "D1950_Npi"));
76
77
78 G4ParticleDefinition * aN1440 = G4ParticleTable::GetParticleTable()->FindParticle(12112); // N+
79 G4ParticleDefinition * aN1520 = G4ParticleTable::GetParticleTable()->FindParticle(2124); // N+
80 G4ParticleDefinition * aN1535 = G4ParticleTable::GetParticleTable()->FindParticle(22212); // N+
81 G4ParticleDefinition * aN1650 = G4ParticleTable::GetParticleTable()->FindParticle(32212); // N+
82 G4ParticleDefinition * aN1675 = G4ParticleTable::GetParticleTable()->FindParticle(2216); // N+
83 G4ParticleDefinition * aN1680 = G4ParticleTable::GetParticleTable()->FindParticle(12216); // N+
84 G4ParticleDefinition * aN1700 = G4ParticleTable::GetParticleTable()->FindParticle(22124); // N+
85 G4ParticleDefinition * aN1710 = G4ParticleTable::GetParticleTable()->FindParticle(42212); // N+
86 G4ParticleDefinition * aN1720 = G4ParticleTable::GetParticleTable()->FindParticle(32124); // N+
87 G4ParticleDefinition * aN1900 = G4ParticleTable::GetParticleTable()->FindParticle(42124); // N+
88 G4ParticleDefinition * aN1990 = G4ParticleTable::GetParticleTable()->FindParticle(12218); // N+
89 G4ParticleDefinition * aN2090 = G4ParticleTable::GetParticleTable()->FindParticle(52214); // N+
90 G4ParticleDefinition * aN2190 = G4ParticleTable::GetParticleTable()->FindParticle(2128); // N+
91 G4ParticleDefinition * aN2220 = G4ParticleTable::GetParticleTable()->FindParticle(100002210); // N+
92 G4ParticleDefinition * aN2250 = G4ParticleTable::GetParticleTable()->FindParticle(100012210); // N+
93
94 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1440, "N1440_Npi"));
95 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1520, "N1520_Npi"));
96 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1535, "N1535_Npi"));
97 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1650, "N1650_Npi"));
98 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1675, "N1675_Npi"));
99 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1680, "N1680_Npi"));
100 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1700, "N1700_Npi"));
101 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1710, "N1710_Npi"));
102 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1720, "N1720_Npi"));
103 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1900, "N1900_Npi"));
104 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN1990, "N1990_Npi"));
105 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN2090, "N2090_Npi"));
106 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN2190, "N2190_Npi"));
107 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN2220, "N2220_Npi"));
108 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aPionp, aN2250, "N2250_Npi"));
109
110/*
111 G4ParticleDefinition * aKaon = G4KaonPlus::KaonPlus();
112
113// sig=0 G4ParticleDefinition * aL1405 = G4ParticleTable::GetParticleTable()->FindParticle(13122);
114 G4ParticleDefinition * aL1520 = G4ParticleTable::GetParticleTable()->FindParticle(3124);
115 G4ParticleDefinition * aL1600 = G4ParticleTable::GetParticleTable()->FindParticle(23122);
116 G4ParticleDefinition * aL1670 = G4ParticleTable::GetParticleTable()->FindParticle(33122);
117 G4ParticleDefinition * aL1690 = G4ParticleTable::GetParticleTable()->FindParticle(13124);
118 G4ParticleDefinition * aL1800 = G4ParticleTable::GetParticleTable()->FindParticle(43122);
119 G4ParticleDefinition * aL1810 = G4ParticleTable::GetParticleTable()->FindParticle(53122);
120 G4ParticleDefinition * aL1820 = G4ParticleTable::GetParticleTable()->FindParticle(3126);
121 G4ParticleDefinition * aL1830 = G4ParticleTable::GetParticleTable()->FindParticle(13126);
122 G4ParticleDefinition * aL1890 = G4ParticleTable::GetParticleTable()->FindParticle(23124);
123 G4ParticleDefinition * aL2100 = G4ParticleTable::GetParticleTable()->FindParticle(3128);
124 G4ParticleDefinition * aL2110 = G4ParticleTable::GetParticleTable()->FindParticle(23126);
125 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL1520, "L1520_NKbar"));
126 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL1600, "L1600_NKbar"));
127 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL1690, "L1690_NKbar"));
128 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL1670, "L1670_NKbar"));
129 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL1800, "L1800_NKbar"));
130 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL1810, "L1810_NKbar"));
131 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL1820, "L1820_NKbar"));
132 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL1830, "L1830_NKbar"));
133 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL1890, "L1890_NKbar"));
134 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL2100, "L2100_NKbar"));
135 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aL2110, "L2110_NKbar"));
136
137//@ G4ParticleDefinition * aSigma = G4ParticleTable::GetParticleTable()->FindParticle(3222);
138//@ G4ParticleDefinition * aS1385 = G4ParticleTable::GetParticleTable()->FindParticle(3224);
139 G4ParticleDefinition * aS1660 = G4ParticleTable::GetParticleTable()->FindParticle(13222);
140 G4ParticleDefinition * aS1670 = G4ParticleTable::GetParticleTable()->FindParticle(13224);
141 G4ParticleDefinition * aS1750 = G4ParticleTable::GetParticleTable()->FindParticle(23222);
142 G4ParticleDefinition * aS1775 = G4ParticleTable::GetParticleTable()->FindParticle(3226);
143 G4ParticleDefinition * aS1915 = G4ParticleTable::GetParticleTable()->FindParticle(13226);
144 G4ParticleDefinition * aS1940 = G4ParticleTable::GetParticleTable()->FindParticle(23224);
145 G4ParticleDefinition * aS2030 = G4ParticleTable::GetParticleTable()->FindParticle(3228);
146 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aS1660, "S1660_NKbar"));
147 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aS1670, "S1670_NKbar"));
148 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aS1750, "S1750_NKbar"));
149 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aS1775, "S1775_NKbar"));
150 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aS1915, "S1915_NKbar"));
151 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aS1940, "S1940_NKbar"));
152 G4CollisionComposite::AddComponent(new G4ConcreteMesonBaryonToResonance(aProton, aKaon, aS2030, "S2030_NKbar"));
153*/
154}
G4CollisionComposite::AddComponent
void AddComponent(G4VCollision *aC)
Definition: G4CollisionComposite.hh:72
G4ParticleTable::FindParticle
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
Definition: G4ParticleTable.cc:586
G4ParticleTable::GetParticleTable
static G4ParticleTable * GetParticleTable()
Definition: G4ParticleTable.cc:87
G4PionPlus::PionPlusDefinition
static G4PionPlus * PionPlusDefinition()
Definition: G4PionPlus.cc:92
G4Proton::ProtonDefinition
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:87

References G4CollisionComposite::AddComponent(), G4ParticleTable::FindParticle(), G4ParticleTable::GetParticleTable(), G4PionPlus::PionPlusDefinition(), and G4Proton::ProtonDefinition().

◆ ~G4CollisionMesonBaryonToResonance()

virtual G4CollisionMesonBaryonToResonance::~G4CollisionMesonBaryonToResonance ( )
inlinevirtual

Definition at line 50 of file G4CollisionMesonBaryonToResonance.hh.

50{};

◆ G4CollisionMesonBaryonToResonance() [2/2]

G4CollisionMesonBaryonToResonance::G4CollisionMesonBaryonToResonance ( const G4CollisionMesonBaryonToResonance &  )
private

Member Function Documentation

◆ AddComponent()

void G4CollisionComposite::AddComponent ( G4VCollision *  aC)
inlineinherited

Definition at line 72 of file G4CollisionComposite.hh.

72{components.push_back(aC);}
G4CollisionComposite::components
G4CollisionVector components
Definition: G4CollisionComposite.hh:120

References G4CollisionComposite::components.

Referenced by G4CollisionMesonBaryon::G4CollisionMesonBaryon(), G4CollisionMesonBaryonToResonance(), G4CollisionComposite::Register::operator()(), and G4CollisionComposite::Resolve::operator()().

◆ BufferCrossSection()

void G4CollisionComposite::BufferCrossSection ( const G4ParticleDefinition *  aP,
const G4ParticleDefinition *  bP 
)
privateinherited

Definition at line 141 of file G4CollisionComposite.cc.

143{
144 // check if already buffered
145 size_t i;
146 for(i=0; i<theBuffer.size(); i++)
147 {
148 if(theBuffer[i].InCharge(aP, bP)) return;
149 }
150// G4cerr << "Buffering for "<<aP->GetParticleName()<<" "<<bP->GetParticleName()<<G4endl;
151
152 // buffer the new one.
153 G4CrossSectionBuffer aNewBuff(aP, bP);
154 size_t maxE=nPoints;
155 for(size_t tt=0; tt<maxE; tt++)
156 {
157 G4double aT = theT[tt]*GeV;
158 G4double crossSect = 0;
159 // The total cross-section is summed over all the component channels
160
161 //A.R. 28-Sep-2012 Fix reproducibility problem
162 // Assign the kinetic energy to the lightest of the
163 // two particles, instead to the first one always.
164 G4double atime = 0;
165 G4double btime = 0;
166 G4ThreeVector aPosition(0,0,0);
167 G4ThreeVector bPosition(0,0,0);
168 G4double aM = aP->GetPDGMass();
169 G4double bM = bP->GetPDGMass();
170 G4double aE = aM;
171 G4double bE = bM;
172 G4ThreeVector aMom(0,0,0);
173 G4ThreeVector bMom(0,0,0);
174 if ( aM <= bM ) {
175 aE += aT;
176 aMom = G4ThreeVector(0,0,std::sqrt(aE*aE-aM*aM));
177 } else {
178 bE += aT;
179 bMom = G4ThreeVector(0,0,std::sqrt(bE*bE-bM*bM));
180 }
181 G4LorentzVector a4Momentum(aE, aMom);
182 G4LorentzVector b4Momentum(bE, bMom);
183 G4KineticTrack a(aP, atime, aPosition, a4Momentum);
184 G4KineticTrack b(bP, btime, bPosition, b4Momentum);
185
186 for (i=0; i<components.size(); i++)
187 {
188 if(components[i]->IsInCharge(a,b))
189 {
190 crossSect += components[i]->CrossSection(a,b);
191 }
192 }
193 G4double sqrts = (a4Momentum+b4Momentum).mag();
194 aNewBuff.push_back(sqrts, crossSect);
195 }
196 theBuffer.push_back(aNewBuff);
197// theBuffer.back().Print();
198}
GeV
static constexpr double GeV
Definition: G4SIunits.hh:203
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:36
G4double
double G4double
Definition: G4Types.hh:83
G4CollisionComposite::theBuffer
std::vector< G4CrossSectionBuffer > theBuffer
Definition: G4CollisionComposite.hh:121
G4CollisionComposite::IsInCharge
virtual G4bool IsInCharge(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
Definition: G4CollisionComposite.cc:121
G4CollisionComposite::nPoints
static const G4int nPoints
Definition: G4CollisionComposite.hh:123
G4CollisionComposite::theT
static const G4double theT[]
Definition: G4CollisionComposite.hh:124

References G4CollisionComposite::components, G4ParticleDefinition::GetPDGMass(), GeV, G4CollisionComposite::IsInCharge(), G4CollisionComposite::nPoints, G4CrossSectionBuffer::push_back(), G4CollisionComposite::theBuffer, and G4CollisionComposite::theT.

Referenced by G4CollisionComposite::CrossSection().

◆ BufferedCrossSection()

G4double G4CollisionComposite::BufferedCrossSection ( const G4KineticTrack &  trk1,
const G4KineticTrack &  trk2 
) const
privateinherited

Definition at line 201 of file G4CollisionComposite.cc.

203{
204 for(size_t i=0; i<theBuffer.size(); i++)
205 {
206 if(theBuffer[i].InCharge(trk1.GetDefinition(), trk2.GetDefinition()))
207 {
208 return theBuffer[i].CrossSection(trk1, trk2);
209 }
210 }
211 throw G4HadronicException(__FILE__, __LINE__, "G4CollisionComposite::BufferedCrossSection - Blitz !!");
212 return 0;
213}
G4KineticTrack::GetDefinition
const G4ParticleDefinition * GetDefinition() const
Definition: G4KineticTrack.hh:213

References G4KineticTrack::GetDefinition(), and G4CollisionComposite::theBuffer.

Referenced by G4CollisionComposite::CrossSection().

◆ CrossSection()

virtual G4double G4CollisionMesonBaryonToResonance::CrossSection ( const G4KineticTrack &  trk1,
const G4KineticTrack &  trk2 
) const
inlineprivatevirtual

Reimplemented from G4CollisionComposite.

Definition at line 58 of file G4CollisionMesonBaryonToResonance.hh.

60 {
61 G4double result=0;
62// if(trk1.GetDefinition() == G4PionPlus::PionPlus() && trk2.GetDefinition()==G4Proton::Proton())
63// {
64// result=thepipp.CrossSection(trk1, trk2);
65// }
66// else if(trk2.GetDefinition() == G4PionPlus::PionPlus() && trk1.GetDefinition()==G4Proton::Proton())
67// {
68// result=thepipp.CrossSection(trk2, trk1);
69// }
70// else if(trk1.GetDefinition() == G4PionMinus::PionMinus() && trk2.GetDefinition()==G4Proton::Proton())
71// {
72// result=thepimp.CrossSection(trk1, trk2);
73// }
74// else if(trk2.GetDefinition() == G4PionMinus::PionMinus() && trk1.GetDefinition()==G4Proton::Proton())
75// {
76// result=thepimp.CrossSection(trk2, trk1);
77// }
78// else
79 {
80 result = G4CollisionComposite::CrossSection(trk1, trk2);
81 }
82 return result;
83 }
G4CollisionComposite::CrossSection
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
Definition: G4CollisionComposite.cc:57

References G4CollisionComposite::CrossSection().

◆ establish_G4MT_TLS_G4VCollision()

void G4VCollision::establish_G4MT_TLS_G4VCollision ( )
inherited

Definition at line 143 of file G4VCollision.cc.

143{ }

Referenced by G4VScatteringCollision::establish_G4MT_TLS_G4VScatteringCollision().

◆ FinalState()

G4KineticTrackVector * G4CollisionComposite::FinalState ( const G4KineticTrack &  trk1,
const G4KineticTrack &  trk2 
) const
virtualinherited

Implements G4VCollision.

Definition at line 79 of file G4CollisionComposite.cc.

81{
82 std::vector<G4double> cxCache;
83 G4double partialCxSum = 0.0;
84
85 size_t i;
86 for (i=0; i<components.size(); i++)
87 {
88 G4double partialCx;
89// cout << "comp" << i << " " << components[i]()->GetName();
90 if (components[i]->IsInCharge(trk1,trk2))
91 {
92 partialCx = components[i]->CrossSection(trk1,trk2);
93 }
94 else
95 {
96 partialCx = 0.0;
97 }
98// cout << " cx=" << partialCx << endl;
99 partialCxSum += partialCx;
100 cxCache.push_back(partialCx);
101 }
102
103 G4double random = G4UniformRand()*partialCxSum;
104 G4double running = 0;
105 for (i=0; i<cxCache.size(); i++)
106 {
107 running += cxCache[i];
108 if (running > random)
109 {
110 return components[i]->FinalState(trk1, trk2);
111 }
112 }
113// G4cerr <<"in charge = "<<IsInCharge(trk1, trk2)<<G4endl;
114// G4cerr <<"Cross-section = "<<CrossSection(trk1, trk2)/millibarn<<" "<<running<<" "<<cxCache.size()<<G4endl;
115// G4cerr <<"Names = "<<trk1.GetDefinition()->GetParticleName()<<", "<<trk2.GetDefinition()->GetParticleName()<<G4endl;
116// throw G4HadronicException(__FILE__, __LINE__, "G4CollisionComposite: no final state found!");
117 return NULL;
118}
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:52

References G4CollisionComposite::components, G4UniformRand, and G4CollisionComposite::IsInCharge().

◆ GetAngularDistribution()

virtual const G4VAngularDistribution * G4CollisionComposite::GetAngularDistribution ( ) const
inlinevirtualinherited

Implements G4VCollision.

Reimplemented in G4CollisionNN, and G4CollisionPN.

Definition at line 76 of file G4CollisionComposite.hh.

76{ return 0; }

◆ GetComponents()

virtual const G4CollisionVector * G4CollisionComposite::GetComponents ( ) const
inlinevirtualinherited

Reimplemented from G4VCollision.

Reimplemented in G4CollisionNN, G4CollisionNNToDeltaDelta1600, G4CollisionNNToDeltaDelta1620, G4CollisionNNToDeltaDelta1700, G4CollisionNNToDeltaDelta1900, G4CollisionNNToDeltaDelta1905, G4CollisionNNToDeltaDelta1910, G4CollisionNNToDeltaDelta1920, G4CollisionNNToDeltaDelta1930, G4CollisionNNToDeltaDelta1950, G4CollisionNNToNDelta, G4CollisionNNToNDelta1600, G4CollisionNNToNDelta1620, G4CollisionNNToNDelta1700, G4CollisionNNToNDelta1900, G4CollisionNNToNDelta1905, G4CollisionNNToNDelta1910, G4CollisionNNToNDelta1920, G4CollisionNNToNDelta1930, and G4CollisionNNToNDelta1950.

Definition at line 78 of file G4CollisionComposite.hh.

78{ return &components;}

References G4CollisionComposite::components.

Referenced by G4CollisionComposite::IsInCharge().

◆ GetCrossSectionSource()

virtual const G4VCrossSectionSource * G4CollisionComposite::GetCrossSectionSource ( ) const
inlinevirtualinherited

Implements G4VCollision.

Reimplemented in G4CollisionNN, and G4CollisionPN.

Definition at line 75 of file G4CollisionComposite.hh.

75{ return 0; }

Referenced by G4CollisionComposite::CrossSection().

◆ GetListOfColliders()

virtual const std::vector< G4String > & G4CollisionMesonBaryonToResonance::GetListOfColliders ( G4int  ) const
inlineprotectedvirtual

Implements G4VCollision.

Definition at line 87 of file G4CollisionMesonBaryonToResonance.hh.

88 {
89 throw G4HadronicException(__FILE__, __LINE__, "Tried to call G4CollisionMesonBaryonToResonance::GetListOfColliders. Please find out why!");
90 std::vector<G4String> * aList = new std::vector<G4String>;
91 return *aList;
92 }

◆ GetName()

virtual G4String G4CollisionMesonBaryonToResonance::GetName ( ) const
inlineprivatevirtual

Implements G4VCollision.

Definition at line 57 of file G4CollisionMesonBaryonToResonance.hh.

57{ return "mN -> baryon Resonance Collision"; }

◆ GetNumberOfPartons()

G4int G4VCollision::GetNumberOfPartons ( const G4ParticleDefinition *  aP) const
inlineprotectedinherited

Definition at line 68 of file G4VCollision.hh.

69 {
70 G4int result = 0;
71 for(G4int i=0; i<6; i++)
72 {
73 result += aP->GetQuarkContent(i+1);
74 result += aP->GetAntiQuarkContent(i+1);
75 }
76 return result;
77 }
G4int
int G4int
Definition: G4Types.hh:85
G4ParticleDefinition::GetQuarkContent
G4int GetQuarkContent(G4int flavor) const
G4ParticleDefinition::GetAntiQuarkContent
G4int GetAntiQuarkContent(G4int flavor) const

References G4ParticleDefinition::GetAntiQuarkContent(), and G4ParticleDefinition::GetQuarkContent().

Referenced by G4CollisionMesonBaryonElastic::IsInCharge().

◆ IsInCharge()

G4bool G4CollisionComposite::IsInCharge ( const G4KineticTrack &  trk1,
const G4KineticTrack &  trk2 
) const
virtualinherited

Implements G4VCollision.

Reimplemented in G4GeneralNNCollision.

Definition at line 121 of file G4CollisionComposite.cc.

123{
124 G4bool isInCharge = false;
125
126 // The composite is in charge if any of its components is in charge
127
128 const G4CollisionVector* comps = GetComponents();
129 if (comps)
130 {
131 G4CollisionVector::const_iterator iter;
132 for (iter = comps->begin(); iter != comps->end(); ++iter)
133 {
134 if ( ((*iter))->IsInCharge(trk1,trk2) ) isInCharge = true;
135 }
136 }
137
138 return isInCharge;
139}
G4CollisionVector
std::vector< G4VCollision * > G4CollisionVector
Definition: G4CollisionVector.hh:48
G4bool
bool G4bool
Definition: G4Types.hh:86
G4CollisionComposite::GetComponents
virtual const G4CollisionVector * GetComponents() const
Definition: G4CollisionComposite.hh:78

References G4CollisionComposite::GetComponents(), and G4CollisionComposite::IsInCharge().

Referenced by G4CollisionComposite::BufferCrossSection(), G4CollisionComposite::FinalState(), and G4CollisionComposite::IsInCharge().

◆ operator!=()

G4bool G4VCollision::operator!= ( const G4VCollision &  right) const
inherited

Definition at line 48 of file G4VCollision.cc.

49{
50 return (this != (G4VCollision *) &right);
51}

◆ operator=()

G4CollisionMesonBaryonToResonance & G4CollisionMesonBaryonToResonance::operator= ( const G4CollisionMesonBaryonToResonance &  )
private

◆ operator==()

G4bool G4VCollision::operator== ( const G4VCollision &  right) const
inherited

Definition at line 42 of file G4VCollision.cc.

43{
44 return (this == (G4VCollision *) &right);
45}

◆ Print() [1/2]

void G4VCollision::Print ( ) const
virtualinherited

Definition at line 70 of file G4VCollision.cc.

71{
72 G4String name = GetName();
73
74 G4cout << "---- " << name << "---- Cross section" << G4endl;
75
76 const G4VCrossSectionSource* xSource = GetCrossSectionSource();
77 if (xSource) xSource->Print();
78
79 G4int nComponents = 0;
80 const G4CollisionVector* components = GetComponents();
81 if (components)
82 {
83 nComponents = components->size();
84 }
85 G4cout << "---- " << name << "---- has " << nComponents << " components" <<G4endl;
86 G4int i = 0;
87 G4CollisionVector::const_iterator iter;
88 if (components)
89 {
90 for (iter = components->begin(); iter != components->end(); ++iter)
91 {
92 G4cout << "---- " << name << " ---- Component " << i << G4endl;
93 ((*iter))->Print();
94 i++;
95 }
96 }
97
98}
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4VCollision::GetComponents
virtual const G4CollisionVector * GetComponents() const
Definition: G4VCollision.hh:79
G4VCollision::Print
virtual void Print() const
Definition: G4VCollision.cc:70
G4VCollision::GetCrossSectionSource
virtual const G4VCrossSectionSource * GetCrossSectionSource() const =0
G4VCollision::GetName
virtual G4String GetName() const =0
G4VCrossSectionSource::Print
virtual void Print() const
Definition: G4VCrossSectionSource.cc:86
G4InuclParticleNames::name
const char * name(G4int ptype)
Definition: G4InuclParticleNames.hh:76

References G4cout, G4endl, G4VCollision::GetComponents(), G4VCollision::GetCrossSectionSource(), G4VCollision::GetName(), G4InuclParticleNames::name(), G4VCollision::Print(), and G4VCrossSectionSource::Print().

Referenced by G4VCollision::Print().

◆ Print() [2/2]

void G4VCollision::Print ( const G4KineticTrack &  trk1,
const G4KineticTrack &  trk2 
) const
virtualinherited

Definition at line 101 of file G4VCollision.cc.

103{
104 G4String name = GetName();
105
106 if (IsInCharge(trk1,trk2))
107 {
108 G4cout << "---- " << name << "is in charge ---- " << G4endl;
109 }
110 else
111 {
112 G4cout << "---- " << name << "is not in charge ---- " << G4endl;
113 }
114
115 G4cout << "---- " << name << "---- Cross section" << G4endl;
116
117 const G4VCrossSectionSource* xSource = GetCrossSectionSource();
118 if (xSource) xSource->Print();
119 G4cout << "Cross section = " << CrossSection(trk1,trk2) << G4endl;
120
121 G4int nComponents = 0;
122 const G4CollisionVector* components = GetComponents();
123 if (components)
124 {
125 nComponents = components->size();
126 }
127 G4cout << "---- " << name << "has " << nComponents << " components" <<G4endl;
128
129 G4int i = 0;
130 G4CollisionVector::const_iterator iter;
131 if (components)
132 {
133 for (iter = components->begin(); iter != components->end(); ++iter)
134 {
135 G4cout << "Component " << i << G4endl;
136 ((*iter))->Print();
137 i++;
138 }
139 }
140}
G4VCollision::CrossSection
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
Definition: G4VCollision.cc:54
G4VCollision::IsInCharge
virtual G4bool IsInCharge(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0

References G4VCollision::CrossSection(), G4cout, G4endl, G4VCollision::GetComponents(), G4VCollision::GetCrossSectionSource(), G4VCollision::GetName(), G4VCollision::IsInCharge(), G4InuclParticleNames::name(), G4VCollision::Print(), and G4VCrossSectionSource::Print().

Field Documentation

◆ bufferMutex

G4Mutex G4CollisionComposite::bufferMutex
mutableprivateinherited

Definition at line 126 of file G4CollisionComposite.hh.

Referenced by G4CollisionComposite::CrossSection(), G4CollisionComposite::G4CollisionComposite(), and G4CollisionComposite::~G4CollisionComposite().

◆ components

G4CollisionVector G4CollisionComposite::components
privateinherited

Definition at line 120 of file G4CollisionComposite.hh.

Referenced by G4CollisionComposite::AddComponent(), G4CollisionComposite::BufferCrossSection(), G4CollisionComposite::FinalState(), G4CollisionComposite::GetComponents(), and G4CollisionComposite::~G4CollisionComposite().

◆ nPoints

const G4int G4CollisionComposite::nPoints = 32
staticprivateinherited

Definition at line 123 of file G4CollisionComposite.hh.

Referenced by G4CollisionComposite::BufferCrossSection().

◆ theBuffer

std::vector<G4CrossSectionBuffer> G4CollisionComposite::theBuffer
privateinherited

Definition at line 121 of file G4CollisionComposite.hh.

Referenced by G4CollisionComposite::BufferCrossSection(), and G4CollisionComposite::BufferedCrossSection().

◆ thepimp

G4XpimNTotal G4CollisionMesonBaryonToResonance::thepimp
private

Definition at line 95 of file G4CollisionMesonBaryonToResonance.hh.

◆ thepipp

G4XpipNTotal G4CollisionMesonBaryonToResonance::thepipp
private

Definition at line 94 of file G4CollisionMesonBaryonToResonance.hh.

◆ theT

const G4double G4CollisionComposite::theT
staticprivateinherited
Initial value:
=
{.01, .03, .05, .1, .15, .2, .3, .4, .5, .6, .7, .8, .9, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.5, 3.0, 3.5, 4.0, 5.0, 6.0, 8.0, 10., 15, 20, 50, 100}

Definition at line 124 of file G4CollisionComposite.hh.

Referenced by G4CollisionComposite::BufferCrossSection().

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.3