Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
G4PiMinusStopMaterial Class Referenceabstract

#include <G4PiMinusStopMaterial.hh>

+ Inheritance diagram for G4PiMinusStopMaterial:

Public Member Functions

 G4PiMinusStopMaterial ()
 
virtual ~G4PiMinusStopMaterial ()
 
virtual std::vector< G4ParticleDefinition * > * DefinitionVector ()
 
virtual std::vector< G4LorentzVector * > * P4Vector (const G4double binding, const G4double mass)
 
virtual G4double FinalNucleons ()=0
 

Protected Member Functions

G4double GenerateAngle (G4double range)
 
G4LorentzVector MakeP4 (G4double p, G4double theta, G4double phi, G4double e)
 
G4double RecoilEnergy (const G4double mass)
 

Protected Attributes

std::vector< G4ParticleDefinition * > * _definitions
 
std::vector< G4LorentzVector * > * _momenta
 
G4DistributionGenerator * _distributionE
 
G4DistributionGenerator * _distributionAngle
 
G4double theR
 

Detailed Description

Definition at line 49 of file G4PiMinusStopMaterial.hh.

Constructor & Destructor Documentation

◆ G4PiMinusStopMaterial()

G4PiMinusStopMaterial::G4PiMinusStopMaterial ( )

Definition at line 55 of file G4PiMinusStopMaterial.cc.

56{
57 _definitions = 0;
58 _momenta = 0;
59 _distributionE = 0;
60 _distributionAngle = 0;
61 theR = 0.5;
62}
G4PiMinusStopMaterial::_definitions
std::vector< G4ParticleDefinition * > * _definitions
Definition: G4PiMinusStopMaterial.hh:80
G4PiMinusStopMaterial::_distributionE
G4DistributionGenerator * _distributionE
Definition: G4PiMinusStopMaterial.hh:82
G4PiMinusStopMaterial::_momenta
std::vector< G4LorentzVector * > * _momenta
Definition: G4PiMinusStopMaterial.hh:81
G4PiMinusStopMaterial::_distributionAngle
G4DistributionGenerator * _distributionAngle
Definition: G4PiMinusStopMaterial.hh:83

◆ ~G4PiMinusStopMaterial()

G4PiMinusStopMaterial::~G4PiMinusStopMaterial ( )
virtual

Definition at line 67 of file G4PiMinusStopMaterial.cc.

68{
69 if (_definitions != 0) delete _definitions;
70 _definitions = 0;
71
72 //A.R. 26-Jul-2012 Coverity fix
73 if (_momenta != 0) {
74 for (unsigned int i=0; i<_momenta->size(); i++) delete(*_momenta)[i];
75 delete _momenta;
76 }
77
78 delete _distributionE;
79 delete _distributionAngle;
80}

Member Function Documentation

◆ DefinitionVector()

std::vector< G4ParticleDefinition * > * G4PiMinusStopMaterial::DefinitionVector ( )
virtual

Definition at line 82 of file G4PiMinusStopMaterial.cc.

83{
84
85 _definitions->push_back(G4Neutron::Neutron());
86
87 G4double ranflat = G4UniformRand();
88 if (ranflat < theR)
89 { _definitions->push_back(G4Proton::Proton()); }
90 else
91 { _definitions->push_back(G4Neutron::Neutron()); }
92
93 return _definitions;
94
95}
G4double
double G4double
Definition: G4Types.hh:64
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:53
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:93

Referenced by G4PiMinusStopAbsorption::DoAbsorption().

◆ FinalNucleons()

virtual G4double G4PiMinusStopMaterial::FinalNucleons ( )
pure virtual

Implemented in G4PiMinusStopAl, G4PiMinusStopC, G4PiMinusStopCo, G4PiMinusStopCu, G4PiMinusStopLi, G4PiMinusStopN, G4PiMinusStopO, G4PiMinusStopPb, and G4PiMinusStopTa.

Referenced by G4PiMinusStopAbsorption::DoAbsorption().

◆ GenerateAngle()

G4double G4PiMinusStopMaterial::GenerateAngle ( G4double  range)
protected

Definition at line 165 of file G4PiMinusStopMaterial.cc.

166{
167 G4double ranflat = G4UniformRand();
168 G4double value = ranflat * x;
169 return value;
170}

Referenced by P4Vector().

◆ MakeP4()

G4LorentzVector G4PiMinusStopMaterial::MakeP4 ( G4double  p,
G4double  theta,
G4double  phi,
G4double  e 
)
protected

Definition at line 172 of file G4PiMinusStopMaterial.cc.

173{
174 // G4LorentzVector p4;
175 G4double px = p * std::sin(theta) * std::cos(phi);
176 G4double py = p * std::sin(theta) * std::sin(phi);
177 G4double pz = p * std::cos(theta);
178 G4LorentzVector p4(px,py,pz,e);
179 return p4;
180}

Referenced by P4Vector().

◆ P4Vector()

std::vector< G4LorentzVector * > * G4PiMinusStopMaterial::P4Vector ( const G4double  binding,
const G4double  mass 
)
virtual

Definition at line 98 of file G4PiMinusStopMaterial.cc.

100{
101 // Generate energy of direct absorption products according to experimental
102 // data. The energy distribution of the two nucleons is assumed to be the
103 // same for protons and neutrons.
104
105 G4double eKin1;
106 G4double eKin2;
107 G4double eRecoil;
108
109 // Assume absorption on two nucleons
110 G4int nNucleons = 2;
111 G4double availableE = G4PionMinus::PionMinus()->GetPDGMass() - nNucleons * binding;
112 G4LorentzVector p1;
113 G4LorentzVector p2;
114
115 do
116 {
117 G4double ranflat;
118 G4double p;
119 G4double energy;
120 G4double mass;
121
122 ranflat = G4UniformRand();
123 eKin1 = _distributionE->Generate(ranflat);
124 mass = (*_definitions)[0]->GetPDGMass();
125 energy = eKin1 + mass;
126 p = std::sqrt(energy*energy - mass*mass);
127 G4double theta1 = pi*G4UniformRand();
128 G4double phi1 = GenerateAngle(2.*pi);
129 p1 = MakeP4(p,theta1,phi1,energy);
130
131 ranflat = G4UniformRand();
132 eKin2 = _distributionE->Generate(ranflat);
133 mass = (*_definitions)[1]->GetPDGMass();
134 energy = eKin2 + mass;
135 p = std::sqrt(energy*energy - mass*mass);
136 ranflat = G4UniformRand();
137 G4double opAngle = _distributionAngle->Generate(ranflat);
138 G4double theta2 = theta1 + opAngle;
139 G4double phi2 = phi1 + opAngle;
140
141 p2 = MakeP4(p,theta2,phi2,energy);
142
143 G4double pNucleus = (p1.vect() + p2.vect()).mag();
144 eRecoil = std::sqrt(pNucleus*pNucleus + massNucleus*massNucleus) - massNucleus;
145
146 // ---- Debug
147 // G4cout << " ---- binding = " << binding << ", nucleus mass = " << massNucleus
148 // << ", p nucleus = " << pNucleus << G4endl;
149 // G4cout << "eKin1,2 " << eKin1 << " " << eKin2 << " eRecoil " << eRecoil
150 // << " availableE " << availableE << G4endl;
151 // ----
152
153 } while ((eKin1 + eKin2 + eRecoil) > availableE);
154
155 //A.R. 26-Jul-2012 Coverity fix
156 if (_momenta != 0) {
157 _momenta->push_back(new G4LorentzVector(p1));
158 _momenta->push_back(new G4LorentzVector(p2));
159 }
160
161 return _momenta;
162
163}
G4LorentzVector
CLHEP::HepLorentzVector G4LorentzVector
Definition: G4LorentzVector.hh:35
G4int
int G4int
Definition: G4Types.hh:66
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
G4DistributionGenerator::Generate
G4double Generate(G4double ranflat)
Definition: G4DistributionGenerator.cc:91
G4PiMinusStopMaterial::GenerateAngle
G4double GenerateAngle(G4double range)
Definition: G4PiMinusStopMaterial.cc:165
G4PiMinusStopMaterial::MakeP4
G4LorentzVector MakeP4(G4double p, G4double theta, G4double phi, G4double e)
Definition: G4PiMinusStopMaterial.cc:172
G4PionMinus::PionMinus
static G4PionMinus * PionMinus()
Definition: G4PionMinus.cc:98
G4INCL::Math::pi
const G4double pi
Definition: G4INCLGlobals.hh:70

Referenced by G4PiMinusStopAbsorption::DoAbsorption().

◆ RecoilEnergy()

G4double G4PiMinusStopMaterial::RecoilEnergy ( const G4double  mass)
protected

Definition at line 182 of file G4PiMinusStopMaterial.cc.

183{
184 G4ThreeVector p(0.,0.,0.);
185
186 //A.R. 26-Jul-2012 Coverity fix
187 if (_momenta != 0) {
188 for (unsigned int i = 0; i< _momenta->size(); i++)
189 {
190 p = p + (*_momenta)[i]->vect();
191 }
192 }
193
194 G4double pNucleus = p.mag();
195 G4double eNucleus = std::sqrt(pNucleus*pNucleus + mass*mass);
196
197 return eNucleus;
198}

Member Data Documentation

◆ _definitions

std::vector<G4ParticleDefinition* >* G4PiMinusStopMaterial::_definitions
protected

Definition at line 80 of file G4PiMinusStopMaterial.hh.

Referenced by DefinitionVector(), G4PiMinusStopAl::G4PiMinusStopAl(), G4PiMinusStopC::G4PiMinusStopC(), G4PiMinusStopCo::G4PiMinusStopCo(), G4PiMinusStopCu::G4PiMinusStopCu(), G4PiMinusStopLi::G4PiMinusStopLi(), G4PiMinusStopMaterial(), G4PiMinusStopN::G4PiMinusStopN(), G4PiMinusStopO::G4PiMinusStopO(), G4PiMinusStopPb::G4PiMinusStopPb(), G4PiMinusStopTa::G4PiMinusStopTa(), and ~G4PiMinusStopMaterial().

◆ _distributionAngle

G4DistributionGenerator* G4PiMinusStopMaterial::_distributionAngle
protected

Definition at line 83 of file G4PiMinusStopMaterial.hh.

Referenced by G4PiMinusStopAl::G4PiMinusStopAl(), G4PiMinusStopC::G4PiMinusStopC(), G4PiMinusStopCo::G4PiMinusStopCo(), G4PiMinusStopCu::G4PiMinusStopCu(), G4PiMinusStopLi::G4PiMinusStopLi(), G4PiMinusStopMaterial(), G4PiMinusStopN::G4PiMinusStopN(), G4PiMinusStopO::G4PiMinusStopO(), G4PiMinusStopPb::G4PiMinusStopPb(), G4PiMinusStopTa::G4PiMinusStopTa(), P4Vector(), and ~G4PiMinusStopMaterial().

◆ _distributionE

G4DistributionGenerator* G4PiMinusStopMaterial::_distributionE
protected

Definition at line 82 of file G4PiMinusStopMaterial.hh.

Referenced by G4PiMinusStopAl::G4PiMinusStopAl(), G4PiMinusStopC::G4PiMinusStopC(), G4PiMinusStopCo::G4PiMinusStopCo(), G4PiMinusStopCu::G4PiMinusStopCu(), G4PiMinusStopLi::G4PiMinusStopLi(), G4PiMinusStopMaterial(), G4PiMinusStopN::G4PiMinusStopN(), G4PiMinusStopO::G4PiMinusStopO(), G4PiMinusStopPb::G4PiMinusStopPb(), G4PiMinusStopTa::G4PiMinusStopTa(), P4Vector(), and ~G4PiMinusStopMaterial().

◆ _momenta

std::vector<G4LorentzVector* >* G4PiMinusStopMaterial::_momenta
protected

Definition at line 81 of file G4PiMinusStopMaterial.hh.

Referenced by G4PiMinusStopAl::G4PiMinusStopAl(), G4PiMinusStopC::G4PiMinusStopC(), G4PiMinusStopCo::G4PiMinusStopCo(), G4PiMinusStopCu::G4PiMinusStopCu(), G4PiMinusStopLi::G4PiMinusStopLi(), G4PiMinusStopMaterial(), G4PiMinusStopN::G4PiMinusStopN(), G4PiMinusStopO::G4PiMinusStopO(), G4PiMinusStopPb::G4PiMinusStopPb(), G4PiMinusStopTa::G4PiMinusStopTa(), P4Vector(), RecoilEnergy(), and ~G4PiMinusStopMaterial().

◆ theR

G4double G4PiMinusStopMaterial::theR
protected

Definition at line 84 of file G4PiMinusStopMaterial.hh.

Referenced by DefinitionVector(), G4PiMinusStopAl::G4PiMinusStopAl(), G4PiMinusStopC::G4PiMinusStopC(), G4PiMinusStopCo::G4PiMinusStopCo(), G4PiMinusStopCu::G4PiMinusStopCu(), G4PiMinusStopLi::G4PiMinusStopLi(), G4PiMinusStopMaterial(), G4PiMinusStopN::G4PiMinusStopN(), G4PiMinusStopO::G4PiMinusStopO(), G4PiMinusStopPb::G4PiMinusStopPb(), and G4PiMinusStopTa::G4PiMinusStopTa().

The documentation for this class was generated from the following files: