Loading...
Searching...
No Matches
G4LightIonQMDGroundStateNucleus Class Reference
#include <G4LightIonQMDGroundStateNucleus.hh>
Inheritance diagram for G4LightIonQMDGroundStateNucleus:
Additional Inherited Members | |
| Protected Attributes inherited from G4QMDSystem | |
| std::vector< G4QMDParticipant * > | participants |
Detailed Description
Definition at line 47 of file G4LightIonQMDGroundStateNucleus.hh.
Constructor & Destructor Documentation
◆ G4LightIonQMDGroundStateNucleus() [1/2]
| G4LightIonQMDGroundStateNucleus::G4LightIonQMDGroundStateNucleus | ( | ) |
◆ G4LightIonQMDGroundStateNucleus() [2/2]
Definition at line 43 of file G4LightIonQMDGroundStateNucleus.cc.
44: maxTrial ( 1000 )
45, r00 ( 1.124 ) // radius parameter for Woods-Saxon [fm]
46, r01 ( 0.5 ) // radius parameter for Woods-Saxon
47, saa ( 0.2 ) // diffuse parameter for initial Woods-Saxon shape
48, rada ( 0.9 ) // cutoff parameter
49, radb ( 0.3 ) // cutoff parameter
50, dsam ( 1.5 ) // minimum distance for same particle [fm]
51, ddif ( 1.0 ) // minimum distance for different particle
52, edepth ( 0.0 )
53, epse ( 0.000001 ) // torelance for energy in [GeV]
54, meanfield ( NULL )
55{
56
57 //std::cout << " G4LightIonQMDGroundStateNucleus( G4int z , G4int a ) Begin " << z << " " << a << std::endl;
58
59 dsam2 = dsam*dsam;
60 ddif2 = ddif*ddif;
61
62 G4LightIonQMDParameters* parameters = G4LightIonQMDParameters::GetInstance();
63
64 hbc = parameters->Get_hbc();
65 gamm = parameters->Get_gamm();
66 cpw = parameters->Get_cpw();
67 cph = parameters->Get_cph();
68 epsx = parameters->Get_epsx();
69 cpc = parameters->Get_cpc();
70
71 cdp = parameters->Get_cdp();
72 c0p = parameters->Get_c0p();
73 c3p = parameters->Get_c3p();
74 csp = parameters->Get_csp();
75 clp = parameters->Get_clp();
76
77 ebini = 0.0;
78 // Following 10 lines should be here, right before the line 90.
79 // Otherwise, mass number cannot be conserved if the projectile or
80 // the target are nucleons.
81 //Nucleon primary or target case;
82 if ( z == 1 && a == 1 ) { // Hydrogen Case or proton primary
83 SetParticipant( new G4QMDParticipant( G4Proton::Proton() , G4ThreeVector( 0.0 ) , G4ThreeVector( 0.0 ) ) );
84// ebini = 0.0;
85 return;
86 }
87 else if ( z == 0 && a == 1 ) { // Neutron primary
88 SetParticipant( new G4QMDParticipant( G4Neutron::Neutron() , G4ThreeVector( 0.0 ) , G4ThreeVector( 0.0 ) ) );
89// ebini = 0.0;
90 return;
91 }
92
93
94 //edepth = 0.0;
95
96 for ( int i = 0 ; i < a ; i++ )
97 {
98
99 G4ParticleDefinition* pd;
100
101 if ( i < z )
102 {
103 pd = G4Proton::Proton();
104 }
105 else
106 {
107 pd = G4Neutron::Neutron();
108 }
109
110 G4ThreeVector p( 0.0 );
111 G4ThreeVector r( 0.0 );
112 G4QMDParticipant* aParticipant = new G4QMDParticipant( pd , p , r );
113 SetParticipant( aParticipant );
114
115 }
116
118
119 rt00 = radious - r01;
120 radm = radious; // JQMD, Skyrme, RMF, Y-H. Sato and A. Haga, 20230308
121 rmax = 1.0 / ( 1.0 + G4Exp ( -rt00/saa ) );
122
123 //maxTrial = 1000;
124
125
126 meanfield = new G4LightIonQMDMeanField();
127 meanfield->SetSystem( this );
128
129 //std::cout << "G4LightIonQMDGroundStateNucleus( G4int z , G4int a ) packNucleons Begin ( z , a ) ( " << z << ", " << a << " )" << std::endl;
130 packNucleons();
131 //std::cout << "G4LightIonQMDGroundStateNucleus( G4int z , G4int a ) packNucleons End" << std::endl;
132
133 delete meanfield;
134
135}
static G4LightIonQMDParameters * GetInstance()
Definition G4LightIonQMDParameters.hh:50
void SetParticipant(G4QMDParticipant *particle)
Definition G4QMDSystem.hh:52
◆ ~G4LightIonQMDGroundStateNucleus()
|
inline |
Definition at line 52 of file G4LightIonQMDGroundStateNucleus.hh.
53 {
54 rho_l.clear();
55 d_pot.clear();
56 };
The documentation for this class was generated from the following files:
