Geant4-11
Public Member Functions | Protected Member Functions | Private Attributes
G4VXResonance Class Referenceabstract

#include <G4VXResonance.hh>

Inheritance diagram for G4VXResonance:
G4VCrossSectionSource G4XResonance

Public Member Functions

virtual G4double CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0
 
 G4VXResonance ()
 
virtual const G4CrossSectionVectorGetComponents () const =0
 
virtual G4double HighLimit () const
 
virtual G4bool IsValid (G4double e) const
 
virtual G4double LowLimit () const
 
virtual G4String Name () const =0
 
G4bool operator!= (const G4VCrossSectionSource &right) const
 
G4bool operator!= (const G4VXResonance &right) const
 
G4bool operator== (const G4VCrossSectionSource &right) const
 
G4bool operator== (const G4VXResonance &right) const
 
virtual void Print () const
 
virtual void PrintAll (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
virtual ~G4VXResonance ()
 

Protected Member Functions

G4double DegeneracyFactor (const G4KineticTrack &trk1, const G4KineticTrack &trk2, G4double iSpinOut1, G4double iSpinOut2) const
 
G4double DetailedBalance (const G4KineticTrack &trk1, const G4KineticTrack &trk2, G4int isoOut1, G4int isoOut2, G4double iSpinOut1, G4double iSpinOut2, G4double mOut1, G4double mOut2) const
 
G4double FcrossX (G4double e, G4double e0, G4double sigma, G4double eParam, G4double power) const
 
const G4ParticleDefinitionFindKeyParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
const G4ParticleDefinitionFindLightParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
 
 G4VXResonance (const G4VXResonance &right)
 
G4bool InLimits (G4double e, G4double eLow, G4double eHigh) const
 
G4double IsospinCorrection (const G4KineticTrack &trk1, const G4KineticTrack &trk2, G4int isoOut1, G4int isoOut2, G4double iSpinOut1, G4double iSpinOut2) const
 
G4VXResonanceoperator= (const G4VXResonance &right)
 

Private Attributes

G4Clebsch clebsch
 

Detailed Description

Definition at line 54 of file G4VXResonance.hh.

Constructor & Destructor Documentation

◆ G4VXResonance() [1/2]

G4VXResonance::G4VXResonance ( )

Definition at line 50 of file G4VXResonance.cc.

51{ }

◆ ~G4VXResonance()

G4VXResonance::~G4VXResonance ( )
virtual

Definition at line 54 of file G4VXResonance.cc.

55{ }

◆ G4VXResonance() [2/2]

G4VXResonance::G4VXResonance ( const G4VXResonance right)
protected

Member Function Documentation

◆ CrossSection()

virtual G4double G4VCrossSectionSource::CrossSection ( const G4KineticTrack trk1,
const G4KineticTrack trk2 
) const
pure virtualinherited

◆ DegeneracyFactor()

G4double G4VXResonance::DegeneracyFactor ( const G4KineticTrack trk1,
const G4KineticTrack trk2,
G4double  iSpinOut1,
G4double  iSpinOut2 
) const
protected

Definition at line 186 of file G4VXResonance.cc.

189{
190 G4double value = 0.;
191
192 const G4ParticleDefinition* in1 = trk1.GetDefinition();
193 const G4ParticleDefinition* in2 = trk2.GetDefinition();
194
195 G4double sIn1 = in1->GetPDGiSpin() + 1.;
196 G4double sIn2 = in2->GetPDGiSpin() + 1.;
197
198 G4double denom = sIn1 * sIn2;
199 if (denom > 0.)
200 {
201 value = (iSpinOut1+1) * (iSpinOut2+1) / denom;
202 }
203 return value;
204}
double G4double
Definition: G4Types.hh:83
const G4ParticleDefinition * GetDefinition() const

References G4KineticTrack::GetDefinition(), and G4ParticleDefinition::GetPDGiSpin().

Referenced by DetailedBalance(), and IsospinCorrection().

◆ DetailedBalance()

G4double G4VXResonance::DetailedBalance ( const G4KineticTrack trk1,
const G4KineticTrack trk2,
G4int  isoOut1,
G4int  isoOut2,
G4double  iSpinOut1,
G4double  iSpinOut2,
G4double  mOut1,
G4double  mOut2 
) const
protected

Definition at line 124 of file G4VXResonance.cc.

129{
130 // To handle the cases when resonances are involved the modified
131 // detailed balance of P. Danielewicz and G.F. Bertsch, Nucl. Phys. A533(1991) 712
132 // is used; in other words, the width of the resonances are folded to get the
133 // mean square of the final state momentum.
134
135 const G4ParticleDefinition* in1 = trk1.GetDefinition();
136 const G4ParticleDefinition* in2 = trk2.GetDefinition();
137 if(in1->IsShortLived() && in2->IsShortLived())
138 {
139 throw G4HadronicException(__FILE__, __LINE__, "Detailed balance for resonance scattering still on the schedule.");
140 }
141
142 G4double result = 0.;
143
144 G4int isoIn1 = in1->GetPDGiIsospin();
145 G4int iso3In1 = in1->GetPDGiIsospin3();
146 G4int isoIn2 = in2->GetPDGiIsospin();
147 G4int iso3In2 = in2->GetPDGiIsospin3();
148 G4double weight = clebsch.Weight(isoIn1, iso3In1, isoIn2, iso3In2, isoOut1, isoOut2);
149
150 if (weight > 00001)
151 {
152 // adding spin counting here ...... does not look quite consistent, but is correct anyway.
153 // revisit in the next design iteration @@
154 G4double degeneracy = DegeneracyFactor(trk1,trk2,iSpinOut1,iSpinOut2);
155 G4double factor = degeneracy * weight;
156
157 // now the phase-space
158 G4double S = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag2();
159 G4double m_1 = in1->GetPDGMass();
160 G4double m_2 = in2->GetPDGMass();
161
162 // on-shell
163 G4double pinitial2 = (S - (m_1+m_2) * (m_1+m_2)) * (S - (m_1-m_2) * (m_1-m_2)) / (4.0*S);
164 G4double pfinal2 = (S - (mOut1+mOut2) * (mOut1+mOut2)) * (S - ( mOut1-mOut2) * (mOut1-mOut2)) / (4.0*S);
165 G4double relativeMomsquared = pfinal2/pinitial2;
166
167 // resonance-nucleon scattering - inverse channel
168 if(in1->IsShortLived())
169 {
171 relativeMomsquared = 1./theI.GetPhaseSpaceIntegral(std::sqrt(S));
172 }
173 else if(in2->IsShortLived())
174 {
176 relativeMomsquared = 1./theI.GetPhaseSpaceIntegral(std::sqrt(S));
177 }
178
179 result = factor * relativeMomsquared;
180 }
181
182 return result;
183}
G4double S(G4double temp)
static const G4double degeneracy
int G4int
Definition: G4Types.hh:85
static G4double Weight(G4int twoJ1, G4int twoM1, G4int twoJ2, G4int twoM2, G4int twoJOut1, G4int twoJOut2)
Definition: G4Clebsch.cc:320
const G4LorentzVector & Get4Momentum() const
G4Clebsch clebsch
G4double DegeneracyFactor(const G4KineticTrack &trk1, const G4KineticTrack &trk2, G4double iSpinOut1, G4double iSpinOut2) const

References clebsch, degeneracy, DegeneracyFactor(), G4KineticTrack::Get4Momentum(), G4KineticTrack::GetDefinition(), G4ParticleDefinition::GetPDGiIsospin(), G4ParticleDefinition::GetPDGiIsospin3(), G4ParticleDefinition::GetPDGMass(), G4DetailedBalancePhaseSpaceIntegral::GetPhaseSpaceIntegral(), G4ParticleDefinition::IsShortLived(), S(), and G4Clebsch::Weight().

Referenced by G4XResonance::CrossSection().

◆ FcrossX()

G4double G4VCrossSectionSource::FcrossX ( G4double  e,
G4double  e0,
G4double  sigma,
G4double  eParam,
G4double  power 
) const
protectedinherited

Definition at line 174 of file G4VCrossSectionSource.cc.

176{
177 G4double result = 0.;
178
179 G4double denom = eParam*eParam + (e-e0)*(e-e0);
180 if (denom > 0.)
181 {
182 G4double value = (2.* eParam * sigma * (e-e0) / denom) * G4Pow::GetInstance()->powA(((e0 + eParam) / e), power);
183 result = std::max(0., value);
184 }
185 return result;
186}
static G4Pow * GetInstance()
Definition: G4Pow.cc:41
G4double powA(G4double A, G4double y) const
Definition: G4Pow.hh:230
T max(const T t1, const T t2)
brief Return the largest of the two arguments

References G4Pow::GetInstance(), G4INCL::Math::max(), and G4Pow::powA().

◆ FindKeyParticle()

const G4ParticleDefinition * G4VCrossSectionSource::FindKeyParticle ( const G4KineticTrack trk1,
const G4KineticTrack trk2 
) const
protectedinherited

Definition at line 49 of file G4VCrossSectionSource.cc.

51{
52 const G4ParticleDefinition * result;
53
54 const G4ParticleDefinition * p1 = trk1.GetDefinition();
55 const G4ParticleDefinition * p2 = trk2.GetDefinition();
56
57 if( (p1==G4Proton::Proton() && p2==G4Proton::Proton() ) ||
58 (p1==G4Neutron::Neutron() && p2==G4Neutron::Neutron()) )
59 {
60 result = G4Proton::Proton();
61 }
62 else if( (p1==G4Neutron::Neutron() && p2==G4Proton::Proton()) ||
63 (p2==G4Neutron::Neutron() && p1==G4Proton::Proton()) )
64 {
65 result = G4Neutron::Neutron();
66 }
67 else
68 {
69 throw G4HadronicException(__FILE__, __LINE__, "G4VCrossSectionSource: unklnown particles in FindKeyParticle");
70 }
71 return result;
72}
static G4Neutron * Neutron()
Definition: G4Neutron.cc:103
static G4Proton * Proton()
Definition: G4Proton.cc:92

References G4KineticTrack::GetDefinition(), G4Neutron::Neutron(), and G4Proton::Proton().

Referenced by G4XNNElasticLowE::CrossSection(), and G4XNNTotalLowE::CrossSection().

◆ FindLightParticle()

const G4ParticleDefinition * G4VCrossSectionSource::FindLightParticle ( const G4KineticTrack trk1,
const G4KineticTrack trk2 
) const
protectedinherited

Definition at line 158 of file G4VCrossSectionSource.cc.

160{
161 G4double mass1 = trk1.GetDefinition()->GetPDGMass();
162 G4double mass2 = trk2.GetDefinition()->GetPDGMass();
163 if (mass1 < mass2)
164 {
165 return trk1.GetDefinition();
166 }
167 else
168 {
169 return trk2.GetDefinition();
170 }
171}

References G4KineticTrack::GetDefinition(), and G4ParticleDefinition::GetPDGMass().

Referenced by G4XMesonBaryonElastic::CrossSection().

◆ GetComponents()

virtual const G4CrossSectionVector * G4VCrossSectionSource::GetComponents ( ) const
pure virtualinherited

◆ HighLimit()

G4double G4VCrossSectionSource::HighLimit ( ) const
virtualinherited

◆ InLimits()

G4bool G4VCrossSectionSource::InLimits ( G4double  e,
G4double  eLow,
G4double  eHigh 
) const
protectedinherited

◆ IsospinCorrection()

G4double G4VXResonance::IsospinCorrection ( const G4KineticTrack trk1,
const G4KineticTrack trk2,
G4int  isoOut1,
G4int  isoOut2,
G4double  iSpinOut1,
G4double  iSpinOut2 
) const
protected

Definition at line 70 of file G4VXResonance.cc.

74{
75 G4double result = 0.;
76
77 const G4ParticleDefinition* in1 = trk1.GetDefinition();
78 const G4ParticleDefinition* in2 = trk2.GetDefinition();
79
80 G4int isoIn1 = in1->GetPDGiIsospin();
81 G4int iso3In1 = in1->GetPDGiIsospin3();
82 G4int isoIn2 = in2->GetPDGiIsospin();
83 G4int iso3In2 = in2->GetPDGiIsospin3();
84
87
88 G4double pWeight = clebsch.Weight(isoProton,iso3Proton, isoProton,iso3Proton, isoOut1,isoOut2);
89 if (pWeight == 0.) throw G4HadronicException(__FILE__, __LINE__, "G4VXResonance::IsospinCorrection, no resonances - pWeight is zero");
90
91 if (in1->IsShortLived() || in2->IsShortLived())
92 {
93 // Resonances in the initial state
95 G4double degeneracyFactor = DegeneracyFactor(trk1,trk2,iSpinProton,iSpinProton);
96
97 G4double factor = degeneracyFactor * pWeight;
98 if (factor > DBL_MIN)
99 {
100 // Randomly select the Isospin3 of the initial state resonances
101 std::vector<G4double> iso = clebsch.GenerateIso3(isoIn1,iso3In1,
102 isoIn2,iso3In2,
103 isoProton,isoProton);
104 G4int isoA = G4lrint(iso[0]);
105 G4int isoB = G4lrint(iso[1]);
106 G4double rWeight = clebsch.Weight(isoProton,isoA,
107 isoProton,isoB,
108 isoOut1,isoOut2);
109 result = rWeight / pWeight;
110 }
111 }
112 else
113 {
114 G4double weight = clebsch.Weight(isoIn1,iso3In1, isoIn2,iso3In2, isoOut1,isoOut2);
115 result = weight / pWeight;
116 }
117
118 return result;
119}
static std::vector< G4double > GenerateIso3(G4int twoJ1, G4int twoM1, G4int twoJ2, G4int twoM2, G4int twoJOut1, G4int twoJOut2)
Definition: G4Clebsch.cc:116
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:87
int G4lrint(double ad)
Definition: templates.hh:134
#define DBL_MIN
Definition: templates.hh:54

References clebsch, DBL_MIN, DegeneracyFactor(), G4lrint(), G4Clebsch::GenerateIso3(), G4KineticTrack::GetDefinition(), G4ParticleDefinition::GetPDGiIsospin(), G4ParticleDefinition::GetPDGiIsospin3(), G4ParticleDefinition::GetPDGiSpin(), G4ParticleDefinition::IsShortLived(), G4Proton::ProtonDefinition(), and G4Clebsch::Weight().

Referenced by G4XResonance::CrossSection().

◆ IsValid()

G4bool G4VCrossSectionSource::IsValid ( G4double  e) const
virtualinherited

◆ LowLimit()

G4double G4VCrossSectionSource::LowLimit ( ) const
virtualinherited

Reimplemented in G4XPDGElastic, and G4XPDGTotal.

Definition at line 140 of file G4VCrossSectionSource.cc.

141{
142 return 0.;
143}

Referenced by G4CrossSectionPatch::CrossSection(), G4VCrossSectionSource::IsValid(), and G4CrossSectionPatch::Transition().

◆ Name()

virtual G4String G4VCrossSectionSource::Name ( ) const
pure virtualinherited

◆ operator!=() [1/2]

G4bool G4VCrossSectionSource::operator!= ( const G4VCrossSectionSource right) const
inherited

Definition at line 80 of file G4VCrossSectionSource.cc.

81{
82 return (this != (G4VCrossSectionSource *) &right);
83}

◆ operator!=() [2/2]

G4bool G4VXResonance::operator!= ( const G4VXResonance right) const

Definition at line 64 of file G4VXResonance.cc.

65{
66 return (this != (G4VXResonance *) &right);
67}

◆ operator=()

G4VXResonance & G4VXResonance::operator= ( const G4VXResonance right)
protected

◆ operator==() [1/2]

G4bool G4VCrossSectionSource::operator== ( const G4VCrossSectionSource right) const
inherited

Definition at line 74 of file G4VCrossSectionSource.cc.

75{
76 return (this == (G4VCrossSectionSource *) &right);
77}

◆ operator==() [2/2]

G4bool G4VXResonance::operator== ( const G4VXResonance right) const

Definition at line 58 of file G4VXResonance.cc.

59{
60 return (this == (G4VXResonance *) &right);
61}

◆ Print()

void G4VCrossSectionSource::Print ( ) const
virtualinherited

Reimplemented in G4XNNElasticLowE, G4XnpElasticLowE, and G4XnpTotalLowE.

Definition at line 86 of file G4VCrossSectionSource.cc.

87{
88 G4int nComponents = 0;
89 const G4CrossSectionVector* components = GetComponents();
90 if (components)
91 {
92 nComponents = components->size();
93 }
94 G4cout << "---- " << this->Name() << " ---- has " << nComponents << " components" <<G4endl;
95 G4int i;
96 for (i=0; i<nComponents; i++)
97 {
98 G4cout << "-" << this->Name() << " - Component " << i << ": " <<G4endl;
99
100 G4CrossSectionSourcePtr componentPtr = (*components)[i];
101 G4VCrossSectionSource* component = componentPtr();
102 component->Print();
103 }
104}
std::vector< G4CrossSectionSourcePtr > G4CrossSectionVector
#define G4endl
Definition: G4ios.hh:57
G4GLOB_DLL std::ostream G4cout
virtual const G4CrossSectionVector * GetComponents() const =0
virtual void Print() const
virtual G4String Name() const =0

References G4cout, G4endl, G4VCrossSectionSource::GetComponents(), G4VCrossSectionSource::Name(), and G4VCrossSectionSource::Print().

Referenced by G4VCollision::Print(), G4VCrossSectionSource::Print(), G4XNNElasticLowE::Print(), G4XnpElasticLowE::Print(), and G4XnpTotalLowE::Print().

◆ PrintAll()

void G4VCrossSectionSource::PrintAll ( const G4KineticTrack trk1,
const G4KineticTrack trk2 
) const
virtualinherited

Definition at line 107 of file G4VCrossSectionSource.cc.

108{
109 G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
110 G4double sigma = CrossSection(trk1,trk2) / millibarn;
111 G4cout << "---- " << Name() << ": "
112 << "Ecm = " << sqrtS / GeV << " GeV - "
113 << " Cross section = " << sigma << " mb "
114 << G4endl;
115
116 G4int nComponents = 0;
117 const G4CrossSectionVector* components = GetComponents();
118 if (components != 0)
119 {
120 nComponents = components->size();
121 }
122 G4int i;
123 for (i=0; i<nComponents; i++)
124 {
125 G4cout << "* Component " << i << ": ";
126 G4CrossSectionSourcePtr componentPtr = (*components)[i];
127 G4VCrossSectionSource* component = componentPtr();
128 component->PrintAll(trk1,trk2);
129 }
130}
static constexpr double millibarn
Definition: G4SIunits.hh:86
static constexpr double GeV
Definition: G4SIunits.hh:203
virtual G4double CrossSection(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0
virtual void PrintAll(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

References G4VCrossSectionSource::CrossSection(), G4cout, G4endl, G4KineticTrack::Get4Momentum(), G4VCrossSectionSource::GetComponents(), GeV, millibarn, G4VCrossSectionSource::Name(), and G4VCrossSectionSource::PrintAll().

Referenced by G4VCrossSectionSource::PrintAll().

Field Documentation

◆ clebsch

G4Clebsch G4VXResonance::clebsch
private

Definition at line 90 of file G4VXResonance.hh.

Referenced by DetailedBalance(), and IsospinCorrection().


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