Geant4-11
Public Member Functions | Protected Member Functions
G4VElasticCollision Class Referenceabstract

#include <G4VElasticCollision.hh>

Inheritance diagram for G4VElasticCollision:
G4VCollision G4CollisionMesonBaryonElastic G4CollisionNNElastic G4CollisionnpElastic

Public Member Functions

virtual G4double CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
void establish_G4MT_TLS_G4VCollision ()
 
virtual G4KineticTrackVector * FinalState (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
 G4VElasticCollision ()
 
virtual G4String GetName () const =0
 
virtual G4bool IsInCharge (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0
 
G4bool operator!= (const G4VCollision &right) const
 
G4bool operator!= (const G4VElasticCollision &right) const
 
G4bool operator== (const G4VCollision &right) const
 
G4bool operator== (const G4VElasticCollision &right) const
 
virtual void Print () const
 
virtual void Print (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
virtual ~G4VElasticCollision ()
 

Protected Member Functions

virtual const G4VAngularDistribution * GetAngularDistribution () const =0
 
virtual const G4CollisionVector * GetComponents () const
 
virtual const G4VCrossSectionSource * GetCrossSectionSource () const =0
 
virtual const std::vector< G4String > & GetListOfColliders (G4int whichOne) const =0
 
G4int GetNumberOfPartons (const G4ParticleDefinition *aP) const
 

Detailed Description

Definition at line 39 of file G4VElasticCollision.hh.

Constructor & Destructor Documentation

◆ G4VElasticCollision()

G4VElasticCollision::G4VElasticCollision ( )

Definition at line 44 of file G4VElasticCollision.cc.

45{
46}

◆ ~G4VElasticCollision()

G4VElasticCollision::~G4VElasticCollision ( )
virtual

Definition at line 49 of file G4VElasticCollision.cc.

50{ }

Member Function Documentation

◆ CrossSection()

G4double G4VCollision::CrossSection ( const G4KineticTrack &  trk1,
const G4KineticTrack &  trk2 
) const
virtualinherited

Reimplemented in G4CollisionComposite, G4CollisionMesonBaryonToResonance, and G4CollisionNN.

Definition at line 54 of file G4VCollision.cc.

56{
57 G4double sigma = 0.;
58
59 const G4VCrossSectionSource* xSource = GetCrossSectionSource();
60
61 if (xSource != 0)
62 {
63 // There is a cross section for this Collision
64 sigma = xSource->CrossSection(aTrk1,aTrk2);
65 }
66 return sigma;
67}
G4double
double G4double
Definition: G4Types.hh:83
G4VCollision::GetCrossSectionSource
virtual const G4VCrossSectionSource * GetCrossSectionSource() const =0
G4VCrossSectionSource::CrossSection
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0

References G4VCrossSectionSource::CrossSection(), and G4VCollision::GetCrossSectionSource().

Referenced by G4Scatterer::GetCrossSection(), G4Scatterer::GetTimeToInteraction(), G4VCollision::Print(), and G4Scatterer::Scatter().

◆ 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 * G4VElasticCollision::FinalState ( const G4KineticTrack &  trk1,
const G4KineticTrack &  trk2 
) const
virtual

Implements G4VCollision.

Definition at line 53 of file G4VElasticCollision.cc.

55{
56 const G4VAngularDistribution* angDistribution;
57
58 angDistribution = GetAngularDistribution();
59
60
61 G4LorentzVector pCM=trk1.Get4Momentum() + trk2.Get4Momentum();
62
63 G4LorentzRotation toLabFrame(pCM.boostVector());
64 G4LorentzVector Ptmp=toLabFrame.inverse() * trk1.Get4Momentum(); //trk1 in CMS
65 G4LorentzRotation toZ;
66 toZ.rotateZ(-Ptmp.phi());
67 toZ.rotateY(-Ptmp.theta());
68 toLabFrame *= toZ.inverse();
69
70 G4double S = pCM.mag2();
71 G4double m10 = trk1.GetDefinition()->GetPDGMass();
72 G4double m20 = trk2.GetDefinition()->GetPDGMass();
73 if(S-(m10+m20)*(m10+m20) < 0) return new G4KineticTrackVector;
74
75 G4double m_1 = trk1.GetActualMass();
76 G4double m_2 = trk2.GetActualMass();
77
78 // Angles of outgoing particles
79 G4double cosTheta = angDistribution->CosTheta(S,m_1,m_2);
80
81 if ( (trk1.GetDefinition() == G4Proton::Proton() || trk1.GetDefinition() == G4Neutron::Neutron() )
82 &&(trk2.GetDefinition() == G4Proton::Proton() || trk2.GetDefinition() == G4Neutron::Neutron() ) )
83 {
84 if ( trk1.GetDefinition() == trk2.GetDefinition() )
85 {
86 if ( trk1.GetDefinition() == G4Proton::Proton() )
87 {
88// G4cout << "scatterangle pp " << cosTheta
89// << " " << typeid(*angDistribution).name() << G4endl;
90 } else {
91// G4cout << "scatterangle nn " << cosTheta
92// << " " << typeid(*angDistribution).name() << G4endl;
93 }
94 } else {
95// G4cout << "scatterangle pn " << cosTheta
96// << " " << typeid(*angDistribution).name() << G4endl;
97 }
98 } else {
99// G4cout << "scatterangle other " << cosTheta
100// << " " << typeid(*angDistribution).name() << G4endl;
101 }
102
103 G4double phi = angDistribution->Phi();
104 G4double Theta = std::acos(cosTheta);
105
106 // Unit vector of three-momentum
107 G4ThreeVector pFinal1(std::sin(Theta)*std::cos(phi), std::sin(Theta)*std::sin(phi), cosTheta);
108 // Three momentum in cm system
109 G4double pInCM = std::sqrt((S-(m10+m20)*(m10+m20))*(S-(m10-m20)*(m10-m20))/(4.*S));
110 pFinal1 = pFinal1 * pInCM;
111 G4ThreeVector pFinal2 = -pFinal1;
112
113 G4double eFinal1 = std::sqrt(pFinal1.mag2() + m10*m10);
114 G4double eFinal2 = std::sqrt(pFinal2.mag2() + m20*m20);
115
116 G4LorentzVector p4Final1(pFinal1, eFinal1);
117 G4LorentzVector p4Final2(pFinal2, eFinal2);
118
119 // Lorentz transformation
120 p4Final1 *= toLabFrame;
121 p4Final2 *= toLabFrame;
122
123 // Final tracks are copies of incoming ones, with modified 4-momenta
124 G4KineticTrack* final1 = new G4KineticTrack(trk1);
125 final1->Set4Momentum(p4Final1);
126 G4KineticTrack* final2 = new G4KineticTrack(trk2);
127 final2->Set4Momentum(p4Final2);
128
129 G4KineticTrackVector* finalTracks = new G4KineticTrackVector;
130 finalTracks->push_back(final1);
131 finalTracks->push_back(final2);
132
133 return finalTracks;
134}
S
G4double S(G4double temp)
Definition: G4DNAElectronHoleRecombination.cc:76
CLHEP::Hep3Vector::mag2
double mag2() const
CLHEP::HepLorentzRotation::rotateY
HepLorentzRotation & rotateY(double delta)
Definition: LorentzRotation.cc:158
CLHEP::HepLorentzRotation::rotateZ
HepLorentzRotation & rotateZ(double delta)
Definition: LorentzRotation.cc:170
CLHEP::HepLorentzRotation::inverse
HepLorentzRotation inverse() const
CLHEP::HepLorentzVector::theta
double theta() const
CLHEP::HepLorentzVector::boostVector
Hep3Vector boostVector() const
Definition: LorentzVector.cc:173
G4KineticTrack::Set4Momentum
void Set4Momentum(const G4LorentzVector &a4Momentum)
Definition: G4KineticTrack.hh:253
G4KineticTrack::GetDefinition
const G4ParticleDefinition * GetDefinition() const
Definition: G4KineticTrack.hh:213
G4KineticTrack::Get4Momentum
const G4LorentzVector & Get4Momentum() const
Definition: G4KineticTrack.hh:243
G4KineticTrack::GetActualMass
G4double GetActualMass() const
Definition: G4KineticTrack.hh:322
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:103
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:92
G4VAngularDistribution::Phi
virtual G4double Phi() const
Definition: G4VAngularDistribution.hh:66
G4VAngularDistribution::CosTheta
virtual G4double CosTheta(G4double s, G4double m1, G4double m2) const =0
G4VCollision::GetAngularDistribution
virtual const G4VAngularDistribution * GetAngularDistribution() const =0

References CLHEP::HepLorentzVector::boostVector(), G4VAngularDistribution::CosTheta(), G4KineticTrack::Get4Momentum(), G4KineticTrack::GetActualMass(), G4VCollision::GetAngularDistribution(), G4KineticTrack::GetDefinition(), G4ParticleDefinition::GetPDGMass(), CLHEP::HepLorentzRotation::inverse(), CLHEP::HepLorentzVector::mag2(), CLHEP::Hep3Vector::mag2(), G4Neutron::Neutron(), CLHEP::HepLorentzVector::phi(), G4VAngularDistribution::Phi(), G4Proton::Proton(), CLHEP::HepLorentzRotation::rotateY(), CLHEP::HepLorentzRotation::rotateZ(), S(), G4KineticTrack::Set4Momentum(), and CLHEP::HepLorentzVector::theta().

◆ GetAngularDistribution()

virtual const G4VAngularDistribution * G4VCollision::GetAngularDistribution ( ) const
protectedpure virtualinherited

Implemented in G4CollisionComposite, G4CollisionMesonBaryonElastic, G4CollisionNN, G4CollisionNNElastic, G4CollisionnpElastic, G4CollisionPN, G4VAnnihilationCollision, and G4VScatteringCollision.

Referenced by FinalState().

◆ GetComponents()

virtual const G4CollisionVector * G4VCollision::GetComponents ( ) const
inlineprotectedvirtualinherited

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

Definition at line 79 of file G4VCollision.hh.

79{ return 0;}

Referenced by G4VCollision::Print().

◆ GetCrossSectionSource()

virtual const G4VCrossSectionSource * G4VCollision::GetCrossSectionSource ( ) const
protectedpure virtualinherited

Implemented in G4CollisionComposite, G4CollisionMesonBaryonElastic, G4CollisionNN, G4CollisionNNElastic, G4CollisionnpElastic, G4CollisionPN, G4ConcreteMesonBaryonToResonance, and G4ConcreteNNTwoBodyResonance.

Referenced by G4VCollision::CrossSection(), and G4VCollision::Print().

◆ GetListOfColliders()

virtual const std::vector< G4String > & G4VCollision::GetListOfColliders ( G4int  whichOne) const
protectedpure virtualinherited

Implemented in G4CollisionNN, G4CollisionNNElastic, G4CollisionnpElastic, G4CollisionPN, G4CollisionMesonBaryon, G4CollisionMesonBaryonElastic, G4CollisionMesonBaryonToResonance, G4CollisionNNToDeltaDelta, G4CollisionNNToDeltaDelta1600, G4CollisionNNToDeltaDelta1620, G4CollisionNNToDeltaDelta1700, G4CollisionNNToDeltaDelta1900, G4CollisionNNToDeltaDelta1905, G4CollisionNNToDeltaDelta1910, G4CollisionNNToDeltaDelta1920, G4CollisionNNToDeltaDelta1930, G4CollisionNNToDeltaDelta1950, G4CollisionNNToDeltaDeltastar, G4CollisionNNToDeltaNstar, G4CollisionNNToNDelta, G4CollisionNNToNDelta1600, G4CollisionNNToNDelta1620, G4CollisionNNToNDelta1700, G4CollisionNNToNDelta1900, G4CollisionNNToNDelta1905, G4CollisionNNToNDelta1910, G4CollisionNNToNDelta1920, G4CollisionNNToNDelta1930, G4CollisionNNToNDelta1950, G4CollisionNNToNDeltastar, G4CollisionNNToNNstar, G4CollisionNStarNToNN, G4ConcreteMesonBaryonToResonance, and G4ConcreteNNTwoBodyResonance.

◆ GetName()

virtual G4String G4VCollision::GetName ( ) const
pure virtualinherited

Implemented in G4CollisionMesonBaryon, G4CollisionMesonBaryonElastic, G4CollisionMesonBaryonToResonance, G4CollisionNN, G4CollisionNNElastic, G4CollisionNNToDeltaDelta, G4CollisionNNToDeltaDelta1600, G4CollisionNNToDeltaDelta1620, G4CollisionNNToDeltaDelta1700, G4CollisionNNToDeltaDelta1900, G4CollisionNNToDeltaDelta1905, G4CollisionNNToDeltaDelta1910, G4CollisionNNToDeltaDelta1920, G4CollisionNNToDeltaDelta1930, G4CollisionNNToDeltaDelta1950, G4CollisionNNToDeltaDeltastar, G4CollisionNNToDeltaNstar, G4CollisionNNToNDelta, G4CollisionNNToNDelta1600, G4CollisionNNToNDelta1620, G4CollisionNNToNDelta1700, G4CollisionNNToNDelta1900, G4CollisionNNToNDelta1905, G4CollisionNNToNDelta1910, G4CollisionNNToNDelta1920, G4CollisionNNToNDelta1930, G4CollisionNNToNDelta1950, G4CollisionNNToNDeltastar, G4CollisionNNToNNstar, G4CollisionnpElastic, G4CollisionNStarNToNN, G4CollisionPN, G4ConcreteMesonBaryonToResonance, G4ConcreteNNToDeltaDelta, G4ConcreteNNToDeltaDeltastar, G4ConcreteNNToDeltaNstar, G4ConcreteNNToNDelta, G4ConcreteNNToNDeltaStar, G4ConcreteNNToNNStar, G4ConcreteNNTwoBodyResonance, and G4ConcreteNStarNToNN.

Referenced by G4VCollision::Print().

◆ 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()

virtual G4bool G4VCollision::IsInCharge ( const G4KineticTrack &  trk1,
const G4KineticTrack &  trk2 
) const
pure virtualinherited

Implemented in G4CollisionComposite, G4CollisionMesonBaryonElastic, G4CollisionNNElastic, G4CollisionnpElastic, G4ConcreteMesonBaryonToResonance, G4ConcreteNNTwoBodyResonance, and G4GeneralNNCollision.

Referenced by G4Scatterer::FindCollision(), and G4VCollision::Print().

◆ operator!=() [1/2]

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

Definition at line 48 of file G4VCollision.cc.

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

◆ operator!=() [2/2]

G4bool G4VElasticCollision::operator!= ( const G4VElasticCollision &  right) const

◆ operator==() [1/2]

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

Definition at line 42 of file G4VCollision.cc.

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

◆ operator==() [2/2]

G4bool G4VElasticCollision::operator== ( const G4VElasticCollision &  right) const

◆ 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}
G4CollisionVector
std::vector< G4VCollision * > G4CollisionVector
Definition: G4CollisionVector.hh:48
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::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().

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