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

#include <G4CollisionnpElastic.hh>

Inheritance diagram for G4CollisionnpElastic:
G4VElasticCollision G4VCollision

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
 
 G4CollisionnpElastic ()
 
virtual G4String GetName () const
 
virtual G4bool IsInCharge (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
G4bool operator!= (const G4CollisionnpElastic &right) const
 
G4bool operator!= (const G4VCollision &right) const
 
G4bool operator!= (const G4VElasticCollision &right) const
 
G4bool operator== (const G4CollisionnpElastic &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 ~G4CollisionnpElastic ()
 

Protected Member Functions

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

Private Member Functions

 G4CollisionnpElastic (const G4CollisionnpElastic &)
 
G4CollisionnpElastic & operator= (const G4CollisionnpElastic &)
 

Private Attributes

G4VAngularDistribution * angularDistribution
 
std::vector< G4String > colliders1
 
std::vector< G4String > colliders2
 
G4VCrossSectionSource * crossSectionSource
 

Detailed Description

Definition at line 41 of file G4CollisionnpElastic.hh.

Constructor & Destructor Documentation

◆ G4CollisionnpElastic() [1/2]

G4CollisionnpElastic::G4CollisionnpElastic ( )

Definition at line 42 of file G4CollisionnpElastic.cc.

43{
44 // Subtype of interacting particles
45 G4String subType1 = G4Proton::ProtonDefinition()->GetParticleSubType();
46 G4String subType2 = G4Neutron::NeutronDefinition()->GetParticleSubType();
47
48 colliders1.push_back(subType1);
49 colliders2.push_back(subType2);
50
51// angularDistribution = new G4AngularDistribution(true);
52 angularDistribution = new G4AngularDistributionNP();
53 crossSectionSource = new G4XnpElastic();
54}
G4CollisionnpElastic::crossSectionSource
G4VCrossSectionSource * crossSectionSource
Definition: G4CollisionnpElastic.hh:72
G4CollisionnpElastic::colliders1
std::vector< G4String > colliders1
Definition: G4CollisionnpElastic.hh:75
G4CollisionnpElastic::angularDistribution
G4VAngularDistribution * angularDistribution
Definition: G4CollisionnpElastic.hh:73
G4CollisionnpElastic::colliders2
std::vector< G4String > colliders2
Definition: G4CollisionnpElastic.hh:76
G4Neutron::NeutronDefinition
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:98
G4ParticleDefinition::GetParticleSubType
const G4String & GetParticleSubType() const
Definition: G4ParticleDefinition.hh:131
G4Proton::ProtonDefinition
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:87

References angularDistribution, colliders1, colliders2, crossSectionSource, G4ParticleDefinition::GetParticleSubType(), G4Neutron::NeutronDefinition(), and G4Proton::ProtonDefinition().

◆ ~G4CollisionnpElastic()

G4CollisionnpElastic::~G4CollisionnpElastic ( )
virtual

Definition at line 57 of file G4CollisionnpElastic.cc.

58{
59 delete angularDistribution;
60 delete crossSectionSource;
61}

References angularDistribution, and crossSectionSource.

◆ G4CollisionnpElastic() [2/2]

G4CollisionnpElastic::G4CollisionnpElastic ( const G4CollisionnpElastic &  )
private

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
virtualinherited

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 * G4CollisionnpElastic::GetAngularDistribution ( ) const
inlineprotectedvirtual

Implements G4VCollision.

Definition at line 65 of file G4CollisionnpElastic.hh.

65{ return angularDistribution; }

References angularDistribution.

◆ 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 * G4CollisionnpElastic::GetCrossSectionSource ( ) const
inlineprotectedvirtual

Implements G4VCollision.

Definition at line 64 of file G4CollisionnpElastic.hh.

64{ return crossSectionSource; }

References crossSectionSource.

◆ GetListOfColliders()

const std::vector< G4String > & G4CollisionnpElastic::GetListOfColliders ( G4int  whichOne) const
protectedvirtual

Implements G4VCollision.

Definition at line 84 of file G4CollisionnpElastic.cc.

85{
86 if (whichOne == 1) {
87 return colliders1;
88 } else if (whichOne == 2) {
89 return colliders2;
90 }
91
92 throw G4HadronicException(__FILE__, __LINE__, "G4CollisionnpElastic::GetListOfColliders - Argument outside valid range");
93}

References colliders1, and colliders2.

◆ GetName()

virtual G4String G4CollisionnpElastic::GetName ( ) const
inlinevirtual

Implements G4VCollision.

Definition at line 56 of file G4CollisionnpElastic.hh.

56{ return "np Elastic 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 G4CollisionnpElastic::IsInCharge ( const G4KineticTrack &  trk1,
const G4KineticTrack &  trk2 
) const
virtual

Implements G4VCollision.

Definition at line 64 of file G4CollisionnpElastic.cc.

66{
67 G4bool isInCharge = false;
68
69 const G4ParticleDefinition* def1 = trk1.GetDefinition();
70 const G4ParticleDefinition* def2 = trk2.GetDefinition();
71
72 if ( (def1 == G4Neutron::NeutronDefinition() &&
73 def2 == G4Proton::ProtonDefinition() )
74 ||
75 (def1 == G4Proton::ProtonDefinition() &&
76 def2 == G4Neutron::NeutronDefinition() ) )
77 {
78 isInCharge = true;
79 }
80 return isInCharge;
81}
G4bool
bool G4bool
Definition: G4Types.hh:86

References G4KineticTrack::GetDefinition(), G4Neutron::NeutronDefinition(), and G4Proton::ProtonDefinition().

◆ operator!=() [1/3]

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

◆ operator!=() [2/3]

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

Definition at line 48 of file G4VCollision.cc.

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

◆ operator!=() [3/3]

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

◆ operator=()

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

◆ operator==() [1/3]

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

◆ operator==() [2/3]

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

Definition at line 42 of file G4VCollision.cc.

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

◆ operator==() [3/3]

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

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

Field Documentation

◆ angularDistribution

G4VAngularDistribution* G4CollisionnpElastic::angularDistribution
private

Definition at line 73 of file G4CollisionnpElastic.hh.

Referenced by G4CollisionnpElastic(), GetAngularDistribution(), and ~G4CollisionnpElastic().

◆ colliders1

std::vector<G4String> G4CollisionnpElastic::colliders1
private

Definition at line 75 of file G4CollisionnpElastic.hh.

Referenced by G4CollisionnpElastic(), and GetListOfColliders().

◆ colliders2

std::vector<G4String> G4CollisionnpElastic::colliders2
private

Definition at line 76 of file G4CollisionnpElastic.hh.

Referenced by G4CollisionnpElastic(), and GetListOfColliders().

◆ crossSectionSource

G4VCrossSectionSource* G4CollisionnpElastic::crossSectionSource
private

Definition at line 72 of file G4CollisionnpElastic.hh.

Referenced by G4CollisionnpElastic(), GetCrossSectionSource(), and ~G4CollisionnpElastic().

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