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

#include <G4DNAWaterDissociationDisplacer.hh>

+ Inheritance diagram for G4DNAWaterDissociationDisplacer:

Public Member Functions

 G4DNAWaterDissociationDisplacer ()
 
 ~G4DNAWaterDissociationDisplacer () override
 
std::vector< G4ThreeVector > GetProductsDisplacement (const G4MolecularDissociationChannel *) const override
 
G4ThreeVector GetMotherMoleculeDisplacement (const G4MolecularDissociationChannel *) const override
 
G4ThreeVector radialDistributionOfElectron () const
 
G4ThreeVector radialDistributionOfProducts (G4double r_rms) const
 
- Public Member Functions inherited from G4VMolecularDissociationDisplacer
virtual ~G4VMolecularDissociationDisplacer ()=default
 
void SetVerbose (G4int verbose)
 

Static Public Member Functions

static G4double ElectronProbaDistribution (G4double r)
 

Additional Inherited Members

- Protected Member Functions inherited from G4VMolecularDissociationDisplacer
 G4VMolecularDissociationDisplacer ()=default
 
- Protected Attributes inherited from G4VMolecularDissociationDisplacer
G4int fVerbose {0}
 

Detailed Description

Definition at line 57 of file G4DNAWaterDissociationDisplacer.hh.

Constructor & Destructor Documentation

◆ G4DNAWaterDissociationDisplacer()

G4DNAWaterDissociationDisplacer::G4DNAWaterDissociationDisplacer ( )

◆ ~G4DNAWaterDissociationDisplacer()

G4DNAWaterDissociationDisplacer::~G4DNAWaterDissociationDisplacer ( )
override

Definition at line 176 of file G4DNAWaterDissociationDisplacer.cc.

177{
178}

Member Function Documentation

◆ ElectronProbaDistribution()

static G4double G4DNAWaterDissociationDisplacer::ElectronProbaDistribution ( G4double r)
static

◆ GetMotherMoleculeDisplacement()

G4ThreeVector G4DNAWaterDissociationDisplacer::GetMotherMoleculeDisplacement ( const G4MolecularDissociationChannel * theDecayChannel) const
overridevirtual

Implements G4VMolecularDissociationDisplacer.

Definition at line 182 of file G4DNAWaterDissociationDisplacer.cc.

184{
185 G4int decayType = theDecayChannel->GetDisplacementType();
186 G4double RMSMotherMoleculeDisplacement(0.0);
187
188 switch (decayType) {
189 case Ionisation_DissociationDecay:
190 RMSMotherMoleculeDisplacement = 2.0 * nanometer;
191 break;
192 case A1B1_DissociationDecay:
193 RMSMotherMoleculeDisplacement = 0.0 * nanometer;
194 break;
195 case B1A1_DissociationDecay:
196 RMSMotherMoleculeDisplacement = 0.0 * nanometer;
197 break;
198 case B1A1_DissociationDecay2:
199 RMSMotherMoleculeDisplacement = 0.0 * nanometer;
200 break;
201 case AutoIonisation:
202 RMSMotherMoleculeDisplacement = 2.0 * nanometer;
203 break;
204 case DissociativeAttachment:
205 RMSMotherMoleculeDisplacement = 0.0 * nanometer;
206 break;
207 case DoubleIonisation_DissociationDecay1:
208 RMSMotherMoleculeDisplacement = 2.0 * nanometer;
209 break;
210 case DoubleIonisation_DissociationDecay2:
211 RMSMotherMoleculeDisplacement = 2.0 * nanometer;
212 break;
213 case DoubleIonisation_DissociationDecay3:
214 RMSMotherMoleculeDisplacement = 2.0 * nanometer;
215 break;
216 case TripleIonisation_DissociationDecay:
217 RMSMotherMoleculeDisplacement = 2.0 * nanometer;
218 break;
219 case QuadrupleIonisation_DissociationDecay:
220 RMSMotherMoleculeDisplacement = 2.0 * nanometer;
221 break;
222 }
223
224 if (RMSMotherMoleculeDisplacement == 0) {
225 return G4ThreeVector(0, 0, 0);
226 }
227
228 auto RandDirection =
229 radialDistributionOfProducts(RMSMotherMoleculeDisplacement);
230
231 return RandDirection;
232}
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition G4ThreeVector.hh:36
G4double
double G4double
Definition G4Types.hh:83
G4int
int G4int
Definition G4Types.hh:85
G4DNAWaterDissociationDisplacer::radialDistributionOfProducts
G4ThreeVector radialDistributionOfProducts(G4double r_rms) const
Definition G4DNAWaterDissociationDisplacer.cc:683

◆ GetProductsDisplacement()

vector< G4ThreeVector > G4DNAWaterDissociationDisplacer::GetProductsDisplacement ( const G4MolecularDissociationChannel * pDecayChannel) const
overridevirtual

Implements G4VMolecularDissociationDisplacer.

Definition at line 236 of file G4DNAWaterDissociationDisplacer.cc.

238{
239 G4int nbProducts = pDecayChannel->GetNbProducts();
240 vector<G4ThreeVector> theProductDisplacementVector(nbProducts);
241
242 typedef map<const G4MoleculeDefinition*, G4double> RMSmap;
243 RMSmap theRMSmap;
244
245 G4int decayType = pDecayChannel->GetDisplacementType();
246
247 switch (decayType) {
248 case Ionisation_DissociationDecay:
249 {
250 if (fVerbose != 0) {
251 G4cout << "Ionisation_DissociationDecay" << G4endl;
252 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
253 }
254 G4double RdmValue = G4UniformRand();
255
256 if (RdmValue < 0.5) {
257 // H3O
258 theRMSmap[G4H3O::Definition()] = 0.0 * nanometer;
259 // OH
260 theRMSmap[G4OH::Definition()] = 0.8 * nanometer;
261 } else {
262 // H3O
263 theRMSmap[G4H3O::Definition()] = 0.8 * nanometer;
264 // OH
265 theRMSmap[G4OH::Definition()] = 0.0 * nanometer;
266 }
267
268 for (G4int i = 0; i < nbProducts; i++) {
269 auto pProduct = pDecayChannel->GetProduct(i);
270 G4double theRMSDisplacement = theRMSmap[pProduct->GetDefinition()];
271
272 if (theRMSDisplacement == 0.0) {
273 theProductDisplacementVector[i] = G4ThreeVector();
274 } else {
275 auto RandDirection =
276 radialDistributionOfProducts(theRMSDisplacement);
277 theProductDisplacementVector[i] = RandDirection;
278 }
279 }
280 break;
281 }
282 case A1B1_DissociationDecay:
283 {
284 if (fVerbose != 0) {
285 G4cout << "A1B1_DissociationDecay" << G4endl;
286 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
287 }
288
289 constexpr G4double theRMSDisplacement = 2.4 * nanometer;
290 auto RandDirection = radialDistributionOfProducts(theRMSDisplacement);
291
292 for (G4int i = 0; i < nbProducts; i++) {
293 auto pProduct = pDecayChannel->GetProduct(i);
294
295 if (pProduct->GetDefinition() == G4OH::Definition()) {
296 theProductDisplacementVector[i] = -1.0 / 18.0 * RandDirection;
297 } else if (pProduct->GetDefinition() == G4Hydrogen::Definition()) {
298 theProductDisplacementVector[i] = +17.0 / 18.0 * RandDirection;
299 }
300 }
301 break;
302 }
303 case B1A1_DissociationDecay:
304 {
305 if (fVerbose != 0) {
306 G4cout << "B1A1_DissociationDecay" << G4endl;
307 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
308 }
309
310 constexpr G4double theRMSDisplacement = 0.8 * nanometer;
311 auto RandDirection = radialDistributionOfProducts(theRMSDisplacement);
312
313 G4int NbOfOH = 0;
314 for (G4int i = 0; i < nbProducts; ++i) {
315 auto pProduct = pDecayChannel->GetProduct(i);
316 if (pProduct->GetDefinition() == G4H2::Definition()) {
317 // In the paper of Kreipl (2009)
318 // theProductDisplacementVector[i] = -2.0 / 18.0 * RandDirection;
319
320 // Based on momentum conservation
321 theProductDisplacementVector[i] = -16.0 / 18.0 * RandDirection;
322 } else if (pProduct->GetDefinition() == G4OH::Definition()) {
323 // In the paper of Kreipl (2009)
324 // G4ThreeVector OxygenDisplacement = +16.0 / 18.0 * RandDirection;
325
326 // Based on momentum conservation
327 G4ThreeVector OxygenDisplacement = +2.0 / 18.0 * RandDirection;
328 constexpr G4double OHRMSDisplacement = 1.1 * nanometer;
329
330 auto OHDisplacement =
331 radialDistributionOfProducts(OHRMSDisplacement);
332
333 if (NbOfOH == 0) {
334 OHDisplacement = 0.5 * OHDisplacement;
335 } else {
336 OHDisplacement = -0.5 * OHDisplacement;
337 }
338
339 theProductDisplacementVector[i] =
340 OHDisplacement + OxygenDisplacement;
341
342 ++NbOfOH;
343 }
344 }
345 break;
346 }
347 case B1A1_DissociationDecay2:
348 {
349 if (fVerbose != 0){
350 G4cout << "B1A1_DissociationDecay2" << G4endl;
351 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
352 }
353
354 G4int NbOfH = 0;
355 for (G4int i = 0; i < nbProducts; ++i) {
356 auto pProduct = pDecayChannel->GetProduct(i);
357 if (pProduct->GetDefinition() == G4Oxygen::Definition()) {
358 // O(3p)
359 theProductDisplacementVector[i] = G4ThreeVector(0, 0, 0);
360 } else if (pProduct->GetDefinition() == G4Hydrogen::Definition()) {
361 // H
362 constexpr G4double HRMSDisplacement = 1.6 * nanometer;
363
364 auto HDisplacement =
365 radialDistributionOfProducts(HRMSDisplacement);
366
367 if (NbOfH == 0) {
368 HDisplacement = 0.5 * HDisplacement;
369 } else {
370 HDisplacement = -0.5 * HDisplacement;
371 }
372 theProductDisplacementVector[i] = HDisplacement;
373
374 ++NbOfH;
375 }
376 }
377 break;
378 }
379 case AutoIonisation:
380 {
381 if (fVerbose != 0) {
382 G4cout << "AutoIonisation" << G4endl;
383 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
384 }
385
386 G4double RdmValue = G4UniformRand();
387
388 if (RdmValue < 0.5) {
389 // H3O
390 theRMSmap[G4H3O::Definition()] = 0.0 * nanometer;
391 // OH
392 theRMSmap[G4OH::Definition()] = 0.8 * nanometer;
393 } else {
394 // H3O
395 theRMSmap[G4H3O::Definition()] = 0.8 * nanometer;
396 // OH
397 theRMSmap[G4OH::Definition()] = 0.0 * nanometer;
398 }
399
400 for (G4int i = 0; i < nbProducts; i++) {
401 auto pProduct = pDecayChannel->GetProduct(i);
402 auto theRMSDisplacement = theRMSmap[pProduct->GetDefinition()];
403
404 if (theRMSDisplacement == 0) {
405 theProductDisplacementVector[i] = G4ThreeVector();
406 } else {
407 auto RandDirection =
408 radialDistributionOfProducts(theRMSDisplacement);
409 theProductDisplacementVector[i] = RandDirection;
410 }
411 if (pProduct->GetDefinition() == G4Electron_aq::Definition()) {
412 theProductDisplacementVector[i] = radialDistributionOfElectron();
413 }
414 }
415 break;
416 }
417 case DissociativeAttachment:
418 {
419 if (fVerbose != 0) {
420 G4cout << "DissociativeAttachment" << G4endl;
421 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
422 }
423 constexpr G4double theRMSDisplacement = 0.8 * nanometer;
424 auto RandDirection = radialDistributionOfProducts(theRMSDisplacement);
425
426 G4int NbOfOH = 0;
427 for (G4int i = 0; i < nbProducts; ++i) {
428 auto pProduct = pDecayChannel->GetProduct(i);
429 if (pProduct->GetDefinition() == G4H2::Definition()) {
430 // In the paper of Kreipl (2009)
431 // theProductDisplacementVector[i] = -2.0 / 18.0 * RandDirection;
432
433 // Based on momentum conservation
434 theProductDisplacementVector[i] = -16.0 / 18.0 * RandDirection;
435 } else if (pProduct->GetDefinition() == G4OH::Definition()) {
436 // In the paper of Kreipl (2009)
437 // G4ThreeVector OxygenDisplacement = +16.0 / 18.0 * RandDirection;
438
439 // Based on momentum conservation
440 G4ThreeVector OxygenDisplacement = +2.0 / 18.0 * RandDirection;
441 constexpr G4double OHRMSDisplacement = 1.1 * nanometer;
442
443 auto OHDisplacement =
444 radialDistributionOfProducts(OHRMSDisplacement);
445
446 if (NbOfOH == 0) {
447 OHDisplacement = 0.5 * OHDisplacement;
448 } else {
449 OHDisplacement = -0.5 * OHDisplacement;
450 }
451 theProductDisplacementVector[i] = OHDisplacement +
452 OxygenDisplacement;
453 ++NbOfOH;
454 }
455 }
456 break;
457 }
458 case DoubleIonisation_DissociationDecay1:
459 {
460 if (fVerbose != 0) {
461 G4cout << "DoubleIonisation_DissociationDecay1" << G4endl;
462 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
463 }
464
465 // Ref.) B. Gervais, et al., DOI: 10.1016/j.radphyschem.2005.09.01
466 // Decay Channel #1: H2O^2+ -> 2H+ + O(3P) -> 2H3O+ + O(3P)
467
468 theRMSmap[G4H3O::Definition()] = 1.2 * nanometer;
469 theRMSmap[G4Oxygen::Definition()] = 0.0 * nanometer;
470
471 for (G4int i = 0, num_H3O = 0; i < nbProducts; i++) {
472
473 const auto prod_def = pDecayChannel->GetProduct(i)->GetDefinition();
474
475 if (prod_def == G4H3O::Definition()) {
476 num_H3O++;
477 if (num_H3O == 2) {
478 constexpr G4double H3Op_rms = 0.3 * nanometer;
479 theRMSmap[G4H3O::Definition()] = H3Op_rms;
480 }
481 }
482
483 theProductDisplacementVector[i]
484 = radialDistributionOfProducts(theRMSmap[prod_def]);
485
486 }
487
488 break;
489 }
490 case DoubleIonisation_DissociationDecay2:
491 {
492 if (fVerbose != 0) {
493 G4cout << "DoubleIonisation_DissociationDecay2" << G4endl;
494 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
495 }
496
497 // Ref.) B. Gervais, et al., DOI: 10.1016/j.radphyschem.2005.09.01
498 // Decay Channel #2: H2O^2+ -> H+ + H* + O+ -> 2H3O+ + H* + *OH + O(3P)
499
500 theRMSmap[G4H3O::Definition()] = 1.2 * nanometer;
501 theRMSmap[G4Oxygen::Definition()] = 0.0 * nanometer;
502 theRMSmap[G4Hydrogen::Definition()] = 0.8 * nanometer;
503 theRMSmap[G4OH::Definition()] = 0.3 * nanometer;
504
505 const auto OH_disp = radialDistributionOfProducts(
506 theRMSmap[G4OH::Definition()]);
507
508 for (G4int i = 0, num_H3O = 0; i < nbProducts; i++) {
509
510 const auto prod_def = pDecayChannel->GetProduct(i)->GetDefinition();
511
512 if (prod_def == G4H3O::Definition()) {
513 num_H3O++;
514 if (num_H3O == 2) {
515 constexpr G4double OH_rms = 0.3 * nanometer;
516 theProductDisplacementVector[i]
517 = radialDistributionOfProducts(OH_rms);
518 theProductDisplacementVector[i] += OH_disp;
519 continue;
520 }
521 } else if (prod_def == G4OH::Definition()) {
522 theProductDisplacementVector[i] = OH_disp;
523 continue;
524 }
525
526 theProductDisplacementVector[i]
527 = radialDistributionOfProducts(theRMSmap[prod_def]);
528
529 }
530
531 break;
532 }
533 case DoubleIonisation_DissociationDecay3:
534 {
535 if (fVerbose != 0) {
536 G4cout << "DoubleIonisation_DissociationDecay3" << G4endl;
537 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
538 }
539
540 // Ref.) B. Gervais, et al., DOI: 10.1016/j.radphyschem.2005.09.01
541 // Decay Channel #3: H2O^2+ -> H+ + OH+ -> 2H3O+ + O(3P)
542
543 theRMSmap[G4H3O::Definition()] = 1.2 * nanometer;
544 theRMSmap[G4Oxygen::Definition()] = 0.0 * nanometer;
545
546 for (G4int i = 0; i < nbProducts; i++) {
547
548 const auto prod_def = pDecayChannel->GetProduct(i)->GetDefinition();
549
550 theProductDisplacementVector[i]
551 = radialDistributionOfProducts(theRMSmap[prod_def]);
552
553 }
554
555 break;
556 }
557 case TripleIonisation_DissociationDecay:
558 {
559 if (fVerbose != 0) {
560 G4cout << "TripleIonisation_DissociationDecay" << G4endl;
561 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
562 }
563
564 // Ref.) B. Gervais, et al., DOI: 10.1016/j.radphyschem.2005.09.01
565 // Decay Channel: H2O^3+ -> 3H3O+ + *OH + O(3P)
566 //
567 // ** Detaied decay chain **
568 // H2O^3+ -> O+ + H3O+ + H3O+
569 // |
570 // --> H2O+ + O(3P)
571 // |
572 // --> *OH + H3O+
573
574 theRMSmap[G4H3O::Definition()] = 1.2 * nanometer;
575 theRMSmap[G4Oxygen::Definition()] = 0.0 * nanometer;
576 theRMSmap[G4OH::Definition()] = 0.3 * nanometer;
577
578 const auto OH_disp = radialDistributionOfProducts(
579 theRMSmap[G4OH::Definition()]);
580
581 for (G4int i = 0, num_H3O = 0; i < nbProducts; i++) {
582
583 const auto prod_def = pDecayChannel->GetProduct(i)->GetDefinition();
584
585 if (prod_def == G4H3O::Definition()) {
586 num_H3O++;
587 if (num_H3O == 3) {
588 constexpr G4double H3Op_rms = 0.3 * nanometer;
589 theProductDisplacementVector[i]
590 = radialDistributionOfProducts(H3Op_rms);
591 theProductDisplacementVector[i] += OH_disp;
592 continue;
593 }
594 } else if (prod_def == G4OH::Definition()) {
595 theProductDisplacementVector[i] = OH_disp;
596 continue;
597 }
598
599 theProductDisplacementVector[i]
600 = radialDistributionOfProducts(theRMSmap[prod_def]);
601
602 }
603
604 break;
605 }
606 case QuadrupleIonisation_DissociationDecay:
607 {
608 if (fVerbose != 0) {
609 G4cout << "QuadrupleIonisation_DissociationDecay" << G4endl;
610 G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
611 }
612
613 // Ref.) B. Gervais, et al., DOI: 10.1016/j.radphyschem.2005.09.01
614 // Decay Channel: H2O^4+ -> 4H3O+ + 2*OH + O(3P)
615 //
616 // ** Detaied decay chain **
617 // H2O^4+ -> O^2+ + H3O+ + H3O+
618 // |
619 // --> H2O+ + H2O+ + O(3P)
620 // | |
621 // | --> *OH + H3O+
622 // |
623 // --> *OH + H3O+
624
625 theRMSmap[G4H3O::Definition()] = 1.2 * nanometer;
626 theRMSmap[G4Oxygen::Definition()] = 0.0 * nanometer;
627 theRMSmap[G4OH::Definition()] = 0.3 * nanometer;
628
629 const auto OH_disp1 = radialDistributionOfProducts(
630 theRMSmap[G4OH::Definition()]);
631 const auto OH_disp2 = radialDistributionOfProducts(
632 theRMSmap[G4OH::Definition()]);
633
634 for (G4int i = 0, num_H3O = 0, num_OH = 0; i < nbProducts; i++) {
635
636 const auto prod_def = pDecayChannel->GetProduct(i)->GetDefinition();
637
638 if (prod_def == G4H3O::Definition()) {
639
640 num_H3O++;
641 constexpr G4double H3Op_rms = 0.3 * nanometer;
642
643 if (num_H3O == 3) {
644 theProductDisplacementVector[i]
645 = radialDistributionOfProducts(H3Op_rms);
646 theProductDisplacementVector[i] += OH_disp1;
647 continue;
648 } else if (num_H3O == 4) {
649 theProductDisplacementVector[i]
650 = radialDistributionOfProducts(H3Op_rms);
651 theProductDisplacementVector[i] += OH_disp2;
652 continue;
653 }
654
655 } else if (prod_def == G4OH::Definition()) {
656
657 num_OH++;
658 if (num_OH == 1) {
659 theProductDisplacementVector[i] = OH_disp1;
660 } else {
661 theProductDisplacementVector[i] = OH_disp2;
662 }
663
664 continue;
665 }
666
667 theProductDisplacementVector[i]
668 = radialDistributionOfProducts(theRMSmap[prod_def]);
669
670 }
671
672 break;
673 }
674
675 }
676 return theProductDisplacementVector;
677}
G4endl
#define G4endl
Definition G4ios.hh:67
G4cout
G4GLOB_DLL std::ostream G4cout
G4UniformRand
#define G4UniformRand()
Definition Randomize.hh:52
G4Electron_aq::Definition
static G4Electron_aq * Definition()
Definition G4Electron_aq.cc:45
G4H3O::Definition
static G4H3O * Definition()
Definition G4H3O.cc:46
G4Hydrogen::Definition
static G4Hydrogen * Definition()
Definition G4Hydrogen.cc:45
G4MolecularConfiguration::GetDefinition
const G4MoleculeDefinition * GetDefinition() const
Definition G4MolecularConfiguration.hh:475
G4OH::Definition
static G4OH * Definition()
Definition G4OH.cc:45
G4Oxygen::Definition
static G4Oxygen * Definition()
Definition G4Oxygen.cc:44

◆ radialDistributionOfElectron()

G4ThreeVector G4DNAWaterDissociationDisplacer::radialDistributionOfElectron ( ) const

Definition at line 696 of file G4DNAWaterDissociationDisplacer.cc.

697{/*
698 G4double rand_value = G4UniformRand();
699 size_t nBins = fElectronThermalization.size();
700 size_t bin = size_t(floor(rand_value * nBins));
701 size_t bin_p1 = min(bin + 1, nBins - 1);
702
703 return (fElectronThermalization[bin] * (1. - rand_value)
704 + fElectronThermalization[bin_p1] * rand_value) *
705 G4RandomDirection();*/
706
707 G4ThreeVector pdf = G4ThreeVector(0, 0, 0);
708
709 if (dnaSubType == fRitchie1994eSolvation) {
710 DNA::Penetration::Ritchie1994::GetPenetration(ke, pdf);
711 } else if (dnaSubType == fTerrisol1990eSolvation) {
712 DNA::Penetration::Terrisol1990::GetPenetration(ke, pdf);
713 } else if (dnaSubType == fMeesungnoensolid2002eSolvation) {
714 DNA::Penetration::Meesungnoen2002_amorphous::GetPenetration(ke, pdf);
715 } else if (dnaSubType == fKreipl2009eSolvation) {
716 DNA::Penetration::Kreipl2009::GetPenetration(ke, pdf);
717 } else {
718 DNA::Penetration::Meesungnoen2002::GetPenetration(ke, pdf);
719 }
720 return pdf;
721}
fKreipl2009eSolvation
@ fKreipl2009eSolvation
Definition G4DNAModelSubType.hh:48
fMeesungnoensolid2002eSolvation
@ fMeesungnoensolid2002eSolvation
Definition G4DNAModelSubType.hh:49
fRitchie1994eSolvation
@ fRitchie1994eSolvation
Definition G4DNAModelSubType.hh:45
fTerrisol1990eSolvation
@ fTerrisol1990eSolvation
Definition G4DNAModelSubType.hh:46
DNA::Penetration::Kreipl2009::GetPenetration
static void GetPenetration(G4double energy, G4ThreeVector &displacement)
Definition G4DNAOneStepThermalizationModel.cc:157
DNA::Penetration::Meesungnoen2002_amorphous::GetPenetration
static void GetPenetration(G4double energy, G4ThreeVector &displacement)
Definition G4DNAOneStepThermalizationModel.cc:150
DNA::Penetration::Meesungnoen2002::GetPenetration
static void GetPenetration(G4double energy, G4ThreeVector &displacement)
Definition G4DNAOneStepThermalizationModel.cc:144
DNA::Penetration::Ritchie1994::GetPenetration
static void GetPenetration(G4double energy, G4ThreeVector &displacement)
Definition G4DNAOneStepThermalizationModel.cc:176
DNA::Penetration::Terrisol1990::GetPenetration
static void GetPenetration(G4double energy, G4ThreeVector &displacement)
Definition G4DNAOneStepThermalizationModel.cc:238

Referenced by GetProductsDisplacement().

◆ radialDistributionOfProducts()

G4ThreeVector G4DNAWaterDissociationDisplacer::radialDistributionOfProducts ( G4double r_rms) const

Definition at line 682 of file G4DNAWaterDissociationDisplacer.cc.

684{
685 static const double inverse_sqrt_3 = 1.0 / sqrt(3.0);
686 double sigma = Rrms * inverse_sqrt_3;
687 double x = G4RandGauss::shoot(0.0, sigma);
688 double y = G4RandGauss::shoot(0.0, sigma);
689 double z = G4RandGauss::shoot(0.0, sigma);
690 return G4ThreeVector(x, y, z);
691}

Referenced by GetMotherMoleculeDisplacement(), and GetProductsDisplacement().

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