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

#include <G4LorentzConvertor.hh>

Public Member Functions

 G4LorentzConvertor ()
 
 G4LorentzConvertor (const G4LorentzVector &bmom, G4double bmass, const G4LorentzVector &tmom, G4double tmass)
 
 G4LorentzConvertor (const G4InuclParticle *bullet, const G4InuclParticle *target)
 
void setVerbose (G4int vb=0)
 
void setBullet (const G4InuclParticle *bullet)
 
void setTarget (const G4InuclParticle *target)
 
void setBullet (const G4InuclParticle &bullet)
 
void setTarget (const G4InuclParticle &target)
 
void setBullet (const G4LorentzVector &bmom)
 
void setTarget (const G4LorentzVector &bmom)
 
void setBullet (const G4LorentzVector &bmom, G4double bmass)
 
void setTarget (const G4LorentzVector &tmom, G4double tmass)
 
void toTheCenterOfMass ()
 
void toTheTargetRestFrame ()
 
void fillKinematics ()
 
G4LorentzVector backToTheLab (const G4LorentzVector &mom) const
 
const G4LorentzVector & getBullet () const
 
const G4LorentzVector & getTarget () const
 
G4double getKinEnergyInTheTRS () const
 
G4double getTotalSCMEnergy () const
 
G4double getSCMMomentum () const
 
G4double getTRSMomentum () const
 
G4LorentzVector rotate (const G4LorentzVector &mom) const
 
G4LorentzVector rotate (const G4LorentzVector &mom1, const G4LorentzVector &mom) const
 
G4bool reflectionNeeded () const
 
G4bool trivial () const
 
void printBullet () const
 
void printTarget () const
 

Detailed Description

Definition at line 45 of file G4LorentzConvertor.hh.

Constructor & Destructor Documentation

◆ G4LorentzConvertor() [1/3]

G4LorentzConvertor::G4LorentzConvertor ( )

Definition at line 47 of file G4LorentzConvertor.cc.

48 : verboseLevel(0), v2(0.), ecm_tot(0.), valong(0.), degenerated(false) {}

◆ G4LorentzConvertor() [2/3]

G4LorentzConvertor::G4LorentzConvertor ( const G4LorentzVector &  bmom,
G4double  bmass,
const G4LorentzVector &  tmom,
G4double  tmass 
)

Definition at line 50 of file G4LorentzConvertor.cc.

53 : verboseLevel(0), v2(0.), ecm_tot(0.), valong(0.), degenerated(false) {
54 setBullet(bmom, bmass);
55 setTarget(tmom, tmass);
56}
G4LorentzConvertor::setBullet
void setBullet(const G4InuclParticle *bullet)
Definition: G4LorentzConvertor.cc:67
G4LorentzConvertor::setTarget
void setTarget(const G4InuclParticle *target)
Definition: G4LorentzConvertor.cc:71

◆ G4LorentzConvertor() [3/3]

G4LorentzConvertor::G4LorentzConvertor ( const G4InuclParticle *  bullet,
const G4InuclParticle *  target 
)

Definition at line 58 of file G4LorentzConvertor.cc.

61 : verboseLevel(0), v2(0.), ecm_tot(0.), valong(0.), degenerated(false) {
62 setBullet(bullet);
63 setTarget(target);
64}

Member Function Documentation

◆ backToTheLab()

G4LorentzVector G4LorentzConvertor::backToTheLab ( const G4LorentzVector &  mom) const

Definition at line 135 of file G4LorentzConvertor.cc.

135 {
136 if (verboseLevel > 2)
137 G4cout << " >>> G4LorentzConvertor::backToTheLab" << G4endl;
138
139 if (verboseLevel > 3)
140 G4cout << " at rest: px " << mom.x() << " py " << mom.y() << " pz "
141 << mom.z() << " e " << mom.e() << G4endl
142 << " v2 " << v2 << G4endl;
143
144 G4LorentzVector mom1 = mom;
145 if (v2 > small) mom1.boost(velocity);
146
147 if (verboseLevel > 3)
148 G4cout << " at lab: px " << mom1.x() << " py " << mom1.y() << " pz "
149 << mom1.z() << G4endl;
150
151 return mom1;
152}
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout
CLHEP::HepLorentzVector::boost
HepLorentzVector & boost(double, double, double)
Definition: LorentzVector.cc:59

Referenced by G4CollisionOutput::boostToLabFrame(), G4ElementaryParticleCollider::collide(), G4EquilibriumEvaporator::collide(), G4NonEquilibriumEvaporator::collide(), and G4NucleiModel::initializeCascad().

◆ fillKinematics()

void G4LorentzConvertor::fillKinematics ( )

Definition at line 113 of file G4LorentzConvertor.cc.

113 {
114 ecm_tot = (target_mom+bullet_mom).m();
115
116 scm_direction = scm_momentum.vect().unit();
117 valong = velocity.dot(scm_direction);
118
119 v2 = velocity.mag2();
120
121 G4double pvsq = v2 - valong*valong; // velocity perp to scm_momentum
122 if (verboseLevel > 3) G4cout << " pvsq " << pvsq << G4endl;
123
124 degenerated = (pvsq < small);
125 if (degenerated && verboseLevel > 2)
126 G4cout << " degenerated case (already along Z) " << G4endl;
127
128 if (verboseLevel > 3) {
129 G4cout << " v2 " << v2 << " valong " << valong
130 << " valong*valong " << valong*valong << G4endl;
131 }
132}
G4double
double G4double
Definition: G4Types.hh:64
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CLHEP::Hep3Vector::mag2
double mag2() const
CLHEP::Hep3Vector::dot
double dot(const Hep3Vector &) const
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const

Referenced by toTheCenterOfMass(), and toTheTargetRestFrame().

◆ getBullet()

const G4LorentzVector & G4LorentzConvertor::getBullet ( ) const
inline

Definition at line 93 of file G4LorentzConvertor.hh.

93{ return bullet_mom; }

◆ getKinEnergyInTheTRS()

G4double G4LorentzConvertor::getKinEnergyInTheTRS ( ) const

Definition at line 157 of file G4LorentzConvertor.cc.

157 {
158 if (verboseLevel > 2)
159 G4cout << " >>> G4LorentzConvertor::getKinEnergyInTheTRS" << G4endl;
160
161 G4LorentzVector bmom = bullet_mom;
162 bmom.boost(-target_mom.boostVector());
163 return bmom.e()-bmom.m();
164}
CLHEP::HepLorentzVector::boostVector
Hep3Vector boostVector() const
Definition: LorentzVector.cc:177

Referenced by G4InuclCollider::collide(), and G4ElementaryParticleCollider::collide().

◆ getSCMMomentum()

G4double G4LorentzConvertor::getSCMMomentum ( ) const
inline

Definition at line 98 of file G4LorentzConvertor.hh.

98{ return scm_momentum.rho(); }

Referenced by G4ElementaryParticleCollider::collide().

◆ getTarget()

const G4LorentzVector & G4LorentzConvertor::getTarget ( ) const
inline

Definition at line 94 of file G4LorentzConvertor.hh.

94{ return target_mom; }

◆ getTotalSCMEnergy()

G4double G4LorentzConvertor::getTotalSCMEnergy ( ) const
inline

Definition at line 97 of file G4LorentzConvertor.hh.

97{ return ecm_tot; }

Referenced by G4ElementaryParticleCollider::collide().

◆ getTRSMomentum()

G4double G4LorentzConvertor::getTRSMomentum ( ) const

Definition at line 166 of file G4LorentzConvertor.cc.

166 {
167 if (verboseLevel > 2)
168 G4cout << " >>> G4LorentzConvertor::getTRSMomentum" << G4endl;
169
170 G4LorentzVector bmom = bullet_mom;
171 bmom.boost(-target_mom.boostVector());
172 return bmom.rho();
173}

Referenced by G4InuclCollider::collide().

◆ printBullet()

void G4LorentzConvertor::printBullet ( ) const

Definition at line 289 of file G4LorentzConvertor.cc.

289 {
290 G4cout << " G4LC bullet: px " << bullet_mom.px() << " py " << bullet_mom.py()
291 << " pz " << bullet_mom.pz() << " e " << bullet_mom.e()
292 << " mass " << bullet_mom.m() << G4endl;
293 }

Referenced by setBullet().

◆ printTarget()

void G4LorentzConvertor::printTarget ( ) const

Definition at line 295 of file G4LorentzConvertor.cc.

295 {
296 G4cout << " G4LC target: px " << target_mom.px() << " py " << target_mom.py()
297 << " pz " << target_mom.pz() << " e " << target_mom.e()
298 << " mass " << target_mom.m() << G4endl;
299}

Referenced by setTarget().

◆ reflectionNeeded()

G4bool G4LorentzConvertor::reflectionNeeded ( ) const

Definition at line 265 of file G4LorentzConvertor.cc.

265 {
266 if (verboseLevel > 2)
267 G4cout << " >>> G4LorentzConvertor::reflectionNeeded (query)" << G4endl;
268
269 if (verboseLevel > 3) {
270 G4cout << " v2 = " << v2 << " SCM z = " << scm_momentum.z()
271 << " degenerated? " << degenerated << G4endl;
272 }
273
274 if (v2 < small && !degenerated)
275 throw G4HadronicException(__FILE__, __LINE__, "G4LorentzConvertor::reflectionNeeded - return value undefined");
276
277 if (verboseLevel > 2) {
278 G4cout << " reflection across XY is"
279 << ((v2>=small && (!degenerated || scm_momentum.z()<0.0))?"":" NOT")
280 << " needed" << G4endl;
281 }
282
283 return (v2>=small && (!degenerated || scm_momentum.z()<0.0));
284}

Referenced by G4CollisionOutput::boostToLabFrame().

◆ rotate() [1/2]

G4LorentzVector G4LorentzConvertor::rotate ( const G4LorentzVector &  mom) const

Definition at line 175 of file G4LorentzConvertor.cc.

175 {
176 if (verboseLevel > 2)
177 G4cout << " >>> G4LorentzConvertor::rotate(G4LorentzVector)" << G4endl;
178
179 if (verboseLevel > 3) {
180 G4cout << " valong " << valong << " degenerated " << degenerated << G4endl
181 << " before rotation: px " << mom.x() << " py " << mom.y()
182 << " pz " << mom.z() << G4endl;
183 }
184
185 G4LorentzVector mom_rot = mom;
186 if (!degenerated) {
187 if (verboseLevel > 2)
188 G4cout << " rotating to align with reference z axis " << G4endl;
189
190 G4ThreeVector vscm = velocity - valong*scm_direction;
191 G4ThreeVector vxcm = scm_direction.cross(velocity);
192
193 if (vscm.mag() > small && vxcm.mag() > small) { // Double check
194 if (verboseLevel > 3) {
195 G4cout << " reference z axis " << scm_direction
196 << " vscm " << vscm << " vxcm " << vxcm << G4endl;
197 }
198
199 mom_rot.setVect(mom.x()*vscm.unit() + mom.y()*vxcm.unit() +
200 mom.z()*scm_direction);
201 } else {
202 if (verboseLevel)
203 G4cerr << ">>> G4LorentzVector::rotate zero with !degenerated" << G4endl;
204 }
205 }
206
207 if (verboseLevel > 3) {
208 G4cout << " after rotation: px " << mom_rot.x() << " py " << mom_rot.y()
209 << " pz " << mom_rot.z() << G4endl;
210 }
211
212 return mom_rot;
213}
G4cerr
G4DLLIMPORT std::ostream G4cerr
CLHEP::Hep3Vector::cross
Hep3Vector cross(const Hep3Vector &) const
CLHEP::Hep3Vector::mag
double mag() const
CLHEP::HepLorentzVector::setVect
void setVect(const Hep3Vector &)

Referenced by G4CollisionOutput::boostToLabFrame().

◆ rotate() [2/2]

G4LorentzVector G4LorentzConvertor::rotate ( const G4LorentzVector &  mom1,
const G4LorentzVector &  mom 
) const

Definition at line 215 of file G4LorentzConvertor.cc.

216 {
217 if (verboseLevel > 2)
218 G4cout << " >>> G4LorentzConvertor::rotate(G4LorentzVector,G4LorentzVector)"
219 << G4endl;
220
221 if (verboseLevel > 3) {
222 G4cout << " before rotation: px " << mom.x() << " py " << mom.y()
223 << " pz " << mom.z() << G4endl;
224 }
225
226 G4ThreeVector mom1_dir = mom1.vect().unit();
227 G4double pv = velocity.dot(mom1_dir);
228
229 G4double vperp = v2 - pv*pv; // velocity perpendicular to mom1
230 if (verboseLevel > 3) {
231 G4cout << " vperp " << vperp << " small? " << (vperp <= small) << G4endl;
232 }
233
234 G4LorentzVector mom_rot = mom;
235
236 if (vperp > small) {
237 if (verboseLevel > 2)
238 G4cout << " rotating to align with first z axis " << G4endl;
239
240 G4ThreeVector vmom1 = velocity - pv*mom1_dir;
241 G4ThreeVector vxm1 = mom1_dir.cross(velocity);
242
243 if (vmom1.mag() > small && vxm1.mag() > small) { // Double check
244 if (verboseLevel > 3) {
245 G4cout << " first z axis " << mom1_dir << G4endl
246 << " vmom1 " << vmom1 << " vxm1 " << vxm1 << G4endl;
247 }
248
249 mom_rot.setVect(mom.x()*vmom1.unit() + mom.y()*vxm1.unit() +
250 mom.z()*mom1_dir );
251 } else {
252 if (verboseLevel)
253 G4cerr << ">>> G4LorentzVector::rotate zero with !degenerated" << G4endl;
254 }
255 }
256
257 if (verboseLevel > 3) {
258 G4cout << " after rotation: px " << mom_rot.x() << " py " << mom_rot.y()
259 << " pz " << mom_rot.z() << G4endl;
260 }
261
262 return mom_rot;
263}

◆ setBullet() [1/4]

void G4LorentzConvertor::setBullet ( const G4InuclParticle &  bullet)
inline

Definition at line 60 of file G4LorentzConvertor.hh.

60{ setBullet(&bullet); }

Referenced by setBullet().

◆ setBullet() [2/4]

void G4LorentzConvertor::setBullet ( const G4InuclParticle *  bullet)

Definition at line 67 of file G4LorentzConvertor.cc.

67 {
68 setBullet(bullet->getMomentum());
69}
G4InuclParticle::getMomentum
G4LorentzVector getMomentum() const
Definition: G4InuclParticle.hh:126

Referenced by G4ElementaryParticleCollider::collide(), G4EquilibriumEvaporator::collide(), G4NonEquilibriumEvaporator::collide(), G4LorentzConvertor(), and setBullet().

◆ setBullet() [3/4]

void G4LorentzConvertor::setBullet ( const G4LorentzVector &  bmom)
inline

Definition at line 64 of file G4LorentzConvertor.hh.

64 {
65 bullet_mom = bmom;
66 if (verboseLevel > 3) printBullet();
67 }

◆ setBullet() [4/4]

void G4LorentzConvertor::setBullet ( const G4LorentzVector &  bmom,
G4double  bmass 
)
inline

Definition at line 75 of file G4LorentzConvertor.hh.

75 {
76 bullet_mom.setVectM(bmom.vect(), bmass);
77 if (verboseLevel > 3) printBullet();
78 }
CLHEP::HepLorentzVector::setVectM
void setVectM(const Hep3Vector &spatial, double mass)

◆ setTarget() [1/4]

void G4LorentzConvertor::setTarget ( const G4InuclParticle &  target)
inline

Definition at line 61 of file G4LorentzConvertor.hh.

61{ setTarget(&target); }

Referenced by setTarget().

◆ setTarget() [2/4]

void G4LorentzConvertor::setTarget ( const G4InuclParticle *  target)

Definition at line 71 of file G4LorentzConvertor.cc.

71 {
72 setTarget(target->getMomentum());
73}

Referenced by G4ElementaryParticleCollider::collide(), G4EquilibriumEvaporator::collide(), G4NonEquilibriumEvaporator::collide(), G4LorentzConvertor(), and setTarget().

◆ setTarget() [3/4]

void G4LorentzConvertor::setTarget ( const G4LorentzVector &  bmom)
inline

Definition at line 69 of file G4LorentzConvertor.hh.

69 {
70 target_mom = bmom;
71 if (verboseLevel > 3) printTarget();
72 }

◆ setTarget() [4/4]

void G4LorentzConvertor::setTarget ( const G4LorentzVector &  tmom,
G4double  tmass 
)
inline

Definition at line 80 of file G4LorentzConvertor.hh.

80 {
81 target_mom.setVectM(tmom.vect(), tmass);
82 if (verboseLevel > 3) printTarget();
83 }

◆ setVerbose()

void G4LorentzConvertor::setVerbose ( G4int  vb = 0)
inline

Definition at line 55 of file G4LorentzConvertor.hh.

55{ verboseLevel = vb; }

Referenced by G4ElementaryParticleCollider::collide().

◆ toTheCenterOfMass()

void G4LorentzConvertor::toTheCenterOfMass ( )

Definition at line 78 of file G4LorentzConvertor.cc.

78 {
79 if (verboseLevel > 2)
80 G4cout << " >>> G4LorentzConvertor::toTheCenterOfMass" << G4endl;
81
82 velocity = (target_mom+bullet_mom).boostVector();
83 if (verboseLevel > 3) G4cout << " boost " << velocity << G4endl;
84
85 // "SCM" is reverse target momentum in the CM frame
86 scm_momentum = target_mom;
87 scm_momentum.boost(-velocity);
88 scm_momentum.setVect(-scm_momentum.vect());
89
90 if (verboseLevel > 3) G4cout << " pscm " << scm_momentum.vect() << G4endl;
91
92 fillKinematics();
93}

Referenced by G4ElementaryParticleCollider::collide().

◆ toTheTargetRestFrame()

void G4LorentzConvertor::toTheTargetRestFrame ( )

Definition at line 95 of file G4LorentzConvertor.cc.

95 {
96 if (verboseLevel > 2)
97 G4cout << " >>> G4LorentzConvertor::toTheTargetRestFrame" << G4endl;
98
99 velocity = target_mom.boostVector();
100 if (verboseLevel > 3) G4cout << " boost " << velocity << G4endl;
101
102 // "SCM" is bullet momentum in the target's frame
103 scm_momentum = bullet_mom;
104 scm_momentum.boost(-velocity);
105
106 if (verboseLevel > 3) G4cout << " pseudo-pscm " << scm_momentum.vect() << G4endl;
107
108 fillKinematics();
109}

Referenced by G4InuclCollider::collide(), G4EquilibriumEvaporator::collide(), G4NonEquilibriumEvaporator::collide(), and G4NucleiModel::initializeCascad().

◆ trivial()

G4bool G4LorentzConvertor::trivial ( ) const
inline

Definition at line 108 of file G4LorentzConvertor.hh.

108{ return degenerated; }

Referenced by G4InuclCollider::collide().

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