Geant4  10.02.p03
Public Member Functions | Private Attributes | List of all members
G4ComponentAntiNuclNuclearXS Class Reference

#include <G4ComponentAntiNuclNuclearXS.hh>

Inheritance diagram for G4ComponentAntiNuclNuclearXS:
Collaboration diagram for G4ComponentAntiNuclNuclearXS:

Public Member Functions

 G4ComponentAntiNuclNuclearXS ()
 
virtual ~G4ComponentAntiNuclNuclearXS ()
 
virtual G4double GetTotalIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A)
 
virtual G4double GetTotalElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
 
virtual G4double GetInelasticIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A)
 
virtual G4double GetInelasticElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
 
virtual G4double GetElasticElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
 
virtual G4double GetElasticIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A)
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void DumpPhysicsTable (const G4ParticleDefinition &)
 
virtual void CrossSectionDescription (std::ostream &) const
 
G4double GetAntiHadronNucleonTotCrSc (const G4ParticleDefinition *aParticle, G4double kinEnergy)
 
G4double GetAntiHadronNucleonElCrSc (const G4ParticleDefinition *aParticle, G4double kinEnergy)
 
- Public Member Functions inherited from G4VComponentCrossSection
 G4VComponentCrossSection (const G4String &nam="")
 
virtual ~G4VComponentCrossSection ()
 
G4double GetTotalElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, const G4Element *)
 
G4double GetInelasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, const G4Element *)
 
G4double GetElasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, const G4Element *)
 
virtual G4double ComputeQuasiElasticRatio (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)
 
virtual void Description () const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
G4double GetMinKinEnergy () const
 
void SetMinKinEnergy (G4double value)
 
G4double GetMaxKinEnergy () const
 
void SetMaxKinEnergy (G4double value)
 
const G4String & GetName () const
 

Private Attributes

G4double fRadiusEff
 
G4double fRadiusNN2
 
G4double fTotalXsc
 
G4double fElasticXsc
 
G4double fInelasticXsc
 
G4double fAntiHadronNucleonTotXsc
 
G4double fAntiHadronNucleonElXsc
 
G4double Elab
 
G4double S
 
G4double SqrtS
 
G4double Mn
 
G4double b0
 
G4double b2
 
G4double SqrtS0
 
G4double S0
 
G4double R0
 
G4ParticleDefinition * theAProton
 
G4ParticleDefinition * theANeutron
 
G4ParticleDefinition * theADeuteron
 
G4ParticleDefinition * theATriton
 
G4ParticleDefinition * theAAlpha
 
G4ParticleDefinition * theAHe3
 

Detailed Description

Definition at line 57 of file G4ComponentAntiNuclNuclearXS.hh.

Constructor & Destructor Documentation

◆ G4ComponentAntiNuclNuclearXS()

G4ComponentAntiNuclNuclearXS::G4ComponentAntiNuclNuclearXS ( )

Definition at line 48 of file G4ComponentAntiNuclNuclearXS.cc.

49 : G4VComponentCrossSection("AntiAGlauber"),
50 // fUpperLimit(10000*GeV), fLowerLimit(10*MeV),
51  fRadiusEff(0.0), fRadiusNN2(0.0),
52  fTotalXsc(0.0), fElasticXsc(0.0), fInelasticXsc(0.0),
53  fAntiHadronNucleonTotXsc(0.0), fAntiHadronNucleonElXsc(0.0),
54  Elab(0.0), S(0.0), SqrtS(0)
55 {
56  theAProton = G4AntiProton::AntiProton();
57  theANeutron = G4AntiNeutron::AntiNeutron();
58  theADeuteron = G4AntiDeuteron::AntiDeuteron();
59  theATriton = G4AntiTriton::AntiTriton();
60  theAAlpha = G4AntiAlpha::AntiAlpha();
61  theAHe3 = G4AntiHe3::AntiHe3();
62 
63  Mn = 0.93827231; // GeV
64  b0 = 11.92; // GeV^(-2)
65  b2 = 0.3036; // GeV^(-2)
66  SqrtS0 = 20.74; // GeV
67  S0 = 33.0625; // GeV^2
68  R0 = 1.0; // default value (V.Ivanchenko)
69 }
G4AntiHe3::AntiHe3
static G4AntiHe3 * AntiHe3()
Definition: G4AntiHe3.cc:94
G4AntiDeuteron::AntiDeuteron
static G4AntiDeuteron * AntiDeuteron()
Definition: G4AntiDeuteron.cc:94
G4AntiAlpha::AntiAlpha
static G4AntiAlpha * AntiAlpha()
Definition: G4AntiAlpha.cc:89
G4AntiProton::AntiProton
static G4AntiProton * AntiProton()
Definition: G4AntiProton.cc:93
G4VComponentCrossSection::G4VComponentCrossSection
G4VComponentCrossSection(const G4String &nam="")
Definition: G4VComponentCrossSection.cc:42
G4AntiTriton::AntiTriton
static G4AntiTriton * AntiTriton()
Definition: G4AntiTriton.cc:94
G4AntiNeutron::AntiNeutron
static G4AntiNeutron * AntiNeutron()
Definition: G4AntiNeutron.cc:103
Here is the call graph for this function:

◆ ~G4ComponentAntiNuclNuclearXS()

G4ComponentAntiNuclNuclearXS::~G4ComponentAntiNuclNuclearXS ( )
virtual

Definition at line 75 of file G4ComponentAntiNuclNuclearXS.cc.

76 {
77 }
Here is the call graph for this function:

Member Function Documentation

◆ BuildPhysicsTable()

virtual void G4ComponentAntiNuclNuclearXS::BuildPhysicsTable ( const G4ParticleDefinition &  )
inlinevirtual

Reimplemented from G4VComponentCrossSection.

Definition at line 95 of file G4ComponentAntiNuclNuclearXS.hh.

96  {}

◆ CrossSectionDescription()

void G4ComponentAntiNuclNuclearXS::CrossSectionDescription ( std::ostream &  outFile) const
virtual

Definition at line 356 of file G4ComponentAntiNuclNuclearXS.cc.

357 {
358  outFile << "The G4ComponentAntiNuclNuclearXS calculates total,\n"
359  << "inelastic, elastic cross sections of anti-nucleons and light \n"
360  << "anti-nucleus interactions with nuclei using Glauber's approach.\n"
361  << "It uses parametrizations of antiproton-proton total and elastic \n"
362  << "cross sections and Wood-Saxon distribution of nuclear density.\n"
363  << "The lower limit is 10 MeV, the upper limit is 10 TeV. \n"
364  << "See details in Phys.Lett. B705 (2011) 235. \n";
365 }
Here is the caller graph for this function:

◆ DumpPhysicsTable()

virtual void G4ComponentAntiNuclNuclearXS::DumpPhysicsTable ( const G4ParticleDefinition &  )
inlinevirtual

Reimplemented from G4VComponentCrossSection.

Definition at line 99 of file G4ComponentAntiNuclNuclearXS.hh.

100  {}
Here is the call graph for this function:

◆ GetAntiHadronNucleonElCrSc()

G4double G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonElCrSc ( const G4ParticleDefinition *  aParticle,
G4double  kinEnergy 
)

Definition at line 331 of file G4ComponentAntiNuclNuclearXS.cc.

332 {
333  G4double xsection;
334 
335  G4double SigAss;
336  G4double C, d1, d2, d3 ;
337 
338  GetAntiHadronNucleonTotCrSc(aParticle,kinEnergy);
339 
340  SigAss = 4.5 + 0.101*G4Log(S/S0)*G4Log(S/S0); //mb
341 
342  C = 59.27;
343  d1 = -6.95;
344  d2 = 23.54;
345  d3 = -25.34;
346 
347  xsection = SigAss* (1 + 1. / (std::sqrt(S-4.*Mn*Mn)) / (G4Pow::GetInstance()->powA(R0, 3.))
348  *C* ( 1+d1/SqrtS+d2/(G4Pow::GetInstance()->powA(SqrtS,2.))+d3/(G4Pow::GetInstance()->powA(SqrtS,3.)) ));
349 
350 // xsection *= millibarn;
351 
352  fAntiHadronNucleonElXsc = xsection;
353  return fAntiHadronNucleonElXsc;
354 }
d3
static const G4double d3
Definition: G4KleinNishinaCompton.cc:67
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:55
d1
static const G4double d1
Definition: G4KleinNishinaCompton.cc:66
C
double C(double temp)
Definition: G4DNAElectronHoleRecombination.cc:85
G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonTotCrSc
G4double GetAntiHadronNucleonTotCrSc(const G4ParticleDefinition *aParticle, G4double kinEnergy)
Definition: G4ComponentAntiNuclNuclearXS.cc:292
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:230
G4Pow::powA
G4double powA(G4double A, G4double y) const
Definition: G4Pow.hh:259
G4double
double G4double
Definition: G4Types.hh:76
d2
static const G4double d2
Definition: G4KleinNishinaCompton.cc:66
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetAntiHadronNucleonTotCrSc()

G4double G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonTotCrSc ( const G4ParticleDefinition *  aParticle,
G4double  kinEnergy 
)

Definition at line 292 of file G4ComponentAntiNuclNuclearXS.cc.

293 {
294  G4double xsection, Pmass, Energy, momentum;
295  const G4ParticleDefinition* theParticle = aParticle;
296  Pmass=theParticle->GetPDGMass();
297  Energy=Pmass+kinEnergy;
298  momentum=std::sqrt(Energy*Energy-Pmass*Pmass)/std::abs(theParticle->GetBaryonNumber());
299  G4double Plab = momentum / GeV;
300 
301  G4double B, SigAss;
302  G4double C, d1, d2, d3 ;
303 
304  Elab = std::sqrt(Mn*Mn + Plab*Plab); // GeV
305  S = 2.*Mn*Mn + 2. *Mn*Elab; // GeV^2
306  SqrtS = std::sqrt(S); // GeV
307 
308  B = b0+b2*G4Log(SqrtS/SqrtS0)*G4Log(SqrtS/SqrtS0); //GeV^(-2)
309  SigAss = 36.04 +0.304*G4Log(S/S0)*G4Log(S/S0); //mb
310  R0 = std::sqrt(0.40874044*SigAss - B); //GeV^(-2)
311 
312  C = 13.55;
313  d1 = -4.47;
314  d2 = 12.38;
315  d3 = -12.43;
316  xsection = SigAss*(1 + 1./(std::sqrt(S-4.*Mn*Mn)) / (G4Pow::GetInstance()->powA(R0, 3.))
317  *C* (1+d1/SqrtS+d2/(G4Pow::GetInstance()->powA(SqrtS,2.))+d3/(G4Pow::GetInstance()->powA(SqrtS,3.)) ));
318 
319 // xsection *= millibarn;
320 
321  fAntiHadronNucleonTotXsc = xsection;
322  return fAntiHadronNucleonTotXsc;
323 }
d3
static const G4double d3
Definition: G4KleinNishinaCompton.cc:67
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:55
d1
static const G4double d1
Definition: G4KleinNishinaCompton.cc:66
C
double C(double temp)
Definition: G4DNAElectronHoleRecombination.cc:85
GeV
static const double GeV
Definition: G4SIunits.hh:214
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:230
G4Pow::powA
G4double powA(G4double A, G4double y) const
Definition: G4Pow.hh:259
G4double
double G4double
Definition: G4Types.hh:76
d2
static const G4double d2
Definition: G4KleinNishinaCompton.cc:66
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetElasticElementCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetElasticElementCrossSection ( const G4ParticleDefinition *  aParticle,
G4double  kinEnergy,
G4int  Z,
G4double  A 
)
virtual

Implements G4VComponentCrossSection.

Definition at line 269 of file G4ComponentAntiNuclNuclearXS.cc.

270 {
271  fElasticXsc = GetTotalElementCrossSection(aParticle, kinEnergy, Z, A)-
272  GetInelasticElementCrossSection(aParticle, kinEnergy, Z, A);
273 
274  if (fElasticXsc < 0.) fElasticXsc = 0.;
275 
276  return fElasticXsc;
277 }
G4ComponentAntiNuclNuclearXS::GetInelasticElementCrossSection
virtual G4double GetInelasticElementCrossSection(const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
Definition: G4ComponentAntiNuclNuclearXS.cc:177
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:59
Z
Float_t Z
Definition: advanced/microbeam/plot.C:39
G4ComponentAntiNuclNuclearXS::GetTotalElementCrossSection
virtual G4double GetTotalElementCrossSection(const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
Definition: G4ComponentAntiNuclNuclearXS.cc:89
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetElasticIsotopeCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetElasticIsotopeCrossSection ( const G4ParticleDefinition *  aParticle,
G4double  kinEnergy,
G4int  Z,
G4int  A 
)
virtual

Implements G4VComponentCrossSection.

Definition at line 284 of file G4ComponentAntiNuclNuclearXS.cc.

285 { return GetElasticElementCrossSection(aParticle, kinEnergy, Z, (G4double) A); }
G4ComponentAntiNuclNuclearXS::GetElasticElementCrossSection
virtual G4double GetElasticElementCrossSection(const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
Definition: G4ComponentAntiNuclNuclearXS.cc:269
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:59
Z
Float_t Z
Definition: advanced/microbeam/plot.C:39
G4double
double G4double
Definition: G4Types.hh:76
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetInelasticElementCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetInelasticElementCrossSection ( const G4ParticleDefinition *  aParticle,
G4double  kinEnergy,
G4int  Z,
G4double  A 
)
virtual

Implements G4VComponentCrossSection.

Definition at line 177 of file G4ComponentAntiNuclNuclearXS.cc.

178 {
179  G4double inelxsection, sigmaTotal, sigmaElastic;
180 
181  const G4ParticleDefinition* theParticle = aParticle;
182 
183  sigmaTotal = GetAntiHadronNucleonTotCrSc(theParticle,kinEnergy);
184  sigmaElastic = GetAntiHadronNucleonElCrSc(theParticle,kinEnergy);
185 
186 // calculation of sqr of radius NN-collision
187  fRadiusNN2=sigmaTotal*sigmaTotal*0.1/(8.*sigmaElastic*pi); // fm^2
188 
189 
190 // calculation of effective nuclear radius for Pbar and Nbar interaction (can be changed)
191 
192  //A.R. 29-Jan-2013 : use antiprotons/antineutrons as the default case,
193  // to be used for instance, as first approximation
194  // without validation, for anti-hyperons.
195  //if ( (theParticle == theAProton) || (theParticle == theANeutron) ) {
196  if (A==1)
197  { fInelasticXsc = (sigmaTotal - sigmaElastic) * millibarn;
198  return fInelasticXsc;
199  }
200  fRadiusEff = 1.31*G4Pow::GetInstance()->powA(A, 0.22)+0.9/G4Pow::GetInstance()->powA(A, 1./3.); //fm
201 
202  if( (Z==1) && (A==2) ) fRadiusEff = 3.582; //fm
203  if( (Z==1) && (A==3) ) fRadiusEff = 3.105;
204  if( (Z==2) && (A==3) ) fRadiusEff = 3.105;
205  if( (Z==2) && (A==4) ) fRadiusEff = 2.209;
206  //}
207 
208 //calculation of effective nuclear radius for AntiDeuteron interaction (can be changed)
209 
210  if (theParticle ==theADeuteron)
211 {
212  fRadiusEff = 1.38*G4Pow::GetInstance()->powA(A, 0.21)+1.55/G4Pow::GetInstance()->powA(A, 1./3.);
213 
214  if( (Z==1) && (A==2) ) fRadiusEff = 3.169; //fm
215  if( (Z==1) && (A==3) ) fRadiusEff = 3.066;
216  if( (Z==2) && (A==3) ) fRadiusEff = 3.066;
217  if( (Z==2) && (A==4) ) fRadiusEff = 2.498;
218  }
219 
220 //calculation of effective nuclear radius for AntiHe3 interaction (can be changed)
221 
222  if( (theParticle ==theAHe3) || (theParticle ==theATriton) )
223  {
224  fRadiusEff = 1.34 * G4Pow::GetInstance()->powA(A, 0.21)+1.51/G4Pow::GetInstance()->powA(A, 1./3.);
225 
226  if( (Z==1) && (A==2) ) fRadiusEff = 3.066; //fm
227  if( (Z==1) && (A==3) ) fRadiusEff = 2.973;
228  if( (Z==2) && (A==3) ) fRadiusEff = 2.973;
229  if( (Z==2) && (A==4) ) fRadiusEff = 2.508;
230 
231  }
232 
233 //calculation of effective nuclear radius for AntiAlpha interaction (can be changed)
234 
235  if (theParticle == theAAlpha)
236  {
237  fRadiusEff = 1.3*G4Pow::GetInstance()->powA(A, 0.21)+1.05/G4Pow::GetInstance()->powA(A, 1./3.);
238 
239  if( (Z==1) && (A==2) ) fRadiusEff = 2.498; //fm
240  if( (Z==1) && (A==3) ) fRadiusEff = 2.508;
241  if( (Z==2) && (A==3) ) fRadiusEff = 2.508;
242  if( (Z==2) && (A==4) ) fRadiusEff = 2.158;
243  }
244  G4double R2 = fRadiusEff*fRadiusEff;
245  G4double REf2 = R2+fRadiusNN2;
246  G4double ApAt= std::abs(theParticle->GetBaryonNumber()) * A;
247 
248  inelxsection = pi*REf2 *10* G4Log(1+(ApAt*sigmaTotal/(pi*REf2*10.))); //mb
249  inelxsection = inelxsection * millibarn;
250  fInelasticXsc = inelxsection;
251  return fInelasticXsc;
252 }
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:55
G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonElCrSc
G4double GetAntiHadronNucleonElCrSc(const G4ParticleDefinition *aParticle, G4double kinEnergy)
Definition: G4ComponentAntiNuclNuclearXS.cc:331
G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonTotCrSc
G4double GetAntiHadronNucleonTotCrSc(const G4ParticleDefinition *aParticle, G4double kinEnergy)
Definition: G4ComponentAntiNuclNuclearXS.cc:292
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:59
Z
Float_t Z
Definition: advanced/microbeam/plot.C:39
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:230
pi
static const double pi
Definition: G4SIunits.hh:74
millibarn
static const double millibarn
Definition: G4SIunits.hh:105
G4Pow::powA
G4double powA(G4double A, G4double y) const
Definition: G4Pow.hh:259
G4double
double G4double
Definition: G4Types.hh:76
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetInelasticIsotopeCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetInelasticIsotopeCrossSection ( const G4ParticleDefinition *  aParticle,
G4double  kinEnergy,
G4int  Z,
G4int  A 
)
virtual

Implements G4VComponentCrossSection.

Definition at line 259 of file G4ComponentAntiNuclNuclearXS.cc.

260 {return GetInelasticElementCrossSection(aParticle, kinEnergy, Z, (G4double) A); }
G4ComponentAntiNuclNuclearXS::GetInelasticElementCrossSection
virtual G4double GetInelasticElementCrossSection(const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
Definition: G4ComponentAntiNuclNuclearXS.cc:177
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:59
Z
Float_t Z
Definition: advanced/microbeam/plot.C:39
G4double
double G4double
Definition: G4Types.hh:76
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetTotalElementCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetTotalElementCrossSection ( const G4ParticleDefinition *  aParticle,
G4double  kinEnergy,
G4int  Z,
G4double  A 
)
virtual

Implements G4VComponentCrossSection.

Definition at line 89 of file G4ComponentAntiNuclNuclearXS.cc.

90 {
91  G4double xsection, sigmaTotal, sigmaElastic;
92 
93  const G4ParticleDefinition* theParticle = aParticle;
94 
95  sigmaTotal = GetAntiHadronNucleonTotCrSc(theParticle,kinEnergy);
96  sigmaElastic = GetAntiHadronNucleonElCrSc(theParticle,kinEnergy);
97 
98 // calculation of squared radius of NN-collision
99  fRadiusNN2=sigmaTotal*sigmaTotal*0.1/(8.*sigmaElastic*pi) ; //fm^2
100 
101 // calculation of effective nuclear radius for Pbar and Nbar interactions (can be changed)
102 
103  //A.R. 29-Jan-2013 : use antiprotons/antineutrons as the default case,
104  // to be used for instance, as first approximation
105  // without validation, for anti-hyperons.
106  //if ( (theParticle == theAProton) || (theParticle == theANeutron) ) {
107  if(A==1)
108  { fTotalXsc = sigmaTotal * millibarn;
109  return fTotalXsc; }
110 
111  fRadiusEff = 1.34*G4Pow::GetInstance()->powA(A,0.23)+1.35/G4Pow::GetInstance()->powA(A,1./3.); //fm
112 
113  if( (Z==1) && (A==2) ) fRadiusEff = 3.800; //fm
114  if( (Z==1) && (A==3) ) fRadiusEff = 3.300;
115  if( (Z==2) && (A==3) ) fRadiusEff = 3.300;
116  if( (Z==2) && (A==4) ) fRadiusEff = 2.376;
117  //}
118 
119 //calculation of effective nuclear radius for AntiDeuteron interaction (can be changed)
120  if (theParticle == theADeuteron)
121  { fRadiusEff = 1.46 * G4Pow::GetInstance()->powA(A,0.21) + 1.45 / G4Pow::GetInstance()->powA(A,1./3.);
122 
123  if( (Z==1) && (A==2) ) fRadiusEff = 3.238; //fm
124  if( (Z==1) && (A==3) ) fRadiusEff = 3.144;
125  if( (Z==2) && (A==3) ) fRadiusEff = 3.144;
126  if( (Z==2) && (A==4) ) fRadiusEff = 2.544;
127  }
128 // calculation of effective nuclear radius for AntiHe3 interaction (can be changed)
129 
130  if( (theParticle ==theAHe3) || (theParticle ==theATriton) )
131  { fRadiusEff = 1.40* G4Pow::GetInstance()->powA(A,0.21)+1.63/G4Pow::GetInstance()->powA(A,1./3.);
132 
133  if( (Z==1) && (A==2) ) fRadiusEff = 3.144; //fm
134  if( (Z==1) && (A==3) ) fRadiusEff = 3.075;
135  if( (Z==2) && (A==3) ) fRadiusEff = 3.075;
136  if( (Z==2) && (A==4) ) fRadiusEff = 2.589;
137  }
138 
139 //calculation of effective nuclear radius for AntiAlpha interaction (can be changed)
140 
141  if (theParticle == theAAlpha)
142  {
143  fRadiusEff = 1.35* G4Pow::GetInstance()->powA(A,0.21)+1.1/G4Pow::GetInstance()->powA(A,1./3.);
144 
145  if( (Z==1) && (A==2) ) fRadiusEff = 2.544; //fm
146  if( (Z==1) && (A==3) ) fRadiusEff = 2.589;
147  if( (Z==2) && (A==3) ) fRadiusEff = 2.589;
148  if( (Z==2) && (A==4) ) fRadiusEff = 2.241;
149 
150  }
151 
152  G4double R2 = fRadiusEff*fRadiusEff;
153  G4double REf2 = R2+fRadiusNN2;
154  G4double ApAt = std::abs(theParticle->GetBaryonNumber()) * A;
155 
156  xsection = 2*pi*REf2*10.*G4Log(1+(ApAt*sigmaTotal/(2*pi*REf2*10.))); //mb
157  xsection =xsection *millibarn;
158  fTotalXsc = xsection;
159 
160  return fTotalXsc;
161 }
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:55
G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonElCrSc
G4double GetAntiHadronNucleonElCrSc(const G4ParticleDefinition *aParticle, G4double kinEnergy)
Definition: G4ComponentAntiNuclNuclearXS.cc:331
G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonTotCrSc
G4double GetAntiHadronNucleonTotCrSc(const G4ParticleDefinition *aParticle, G4double kinEnergy)
Definition: G4ComponentAntiNuclNuclearXS.cc:292
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:59
Z
Float_t Z
Definition: advanced/microbeam/plot.C:39
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:230
pi
static const double pi
Definition: G4SIunits.hh:74
millibarn
static const double millibarn
Definition: G4SIunits.hh:105
G4Pow::powA
G4double powA(G4double A, G4double y) const
Definition: G4Pow.hh:259
G4double
double G4double
Definition: G4Types.hh:76
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetTotalIsotopeCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetTotalIsotopeCrossSection ( const G4ParticleDefinition *  aParticle,
G4double  kinEnergy,
G4int  Z,
G4int  A 
)
virtual

Implements G4VComponentCrossSection.

Definition at line 169 of file G4ComponentAntiNuclNuclearXS.cc.

170 { return GetTotalElementCrossSection(aParticle, kinEnergy, Z, (G4double) A); }
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:59
Z
Float_t Z
Definition: advanced/microbeam/plot.C:39
G4ComponentAntiNuclNuclearXS::GetTotalElementCrossSection
virtual G4double GetTotalElementCrossSection(const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
Definition: G4ComponentAntiNuclNuclearXS.cc:89
G4double
double G4double
Definition: G4Types.hh:76
Here is the call graph for this function:
Here is the caller graph for this function:

Member Data Documentation

◆ b0

G4double G4ComponentAntiNuclNuclearXS::b0
private

Definition at line 122 of file G4ComponentAntiNuclNuclearXS.hh.

◆ b2

G4double G4ComponentAntiNuclNuclearXS::b2
private

Definition at line 122 of file G4ComponentAntiNuclNuclearXS.hh.

◆ Elab

G4double G4ComponentAntiNuclNuclearXS::Elab
private

Definition at line 121 of file G4ComponentAntiNuclNuclearXS.hh.

◆ fAntiHadronNucleonElXsc

G4double G4ComponentAntiNuclNuclearXS::fAntiHadronNucleonElXsc
private

Definition at line 120 of file G4ComponentAntiNuclNuclearXS.hh.

◆ fAntiHadronNucleonTotXsc

G4double G4ComponentAntiNuclNuclearXS::fAntiHadronNucleonTotXsc
private

Definition at line 120 of file G4ComponentAntiNuclNuclearXS.hh.

◆ fElasticXsc

G4double G4ComponentAntiNuclNuclearXS::fElasticXsc
private

Definition at line 119 of file G4ComponentAntiNuclNuclearXS.hh.

◆ fInelasticXsc

G4double G4ComponentAntiNuclNuclearXS::fInelasticXsc
private

Definition at line 119 of file G4ComponentAntiNuclNuclearXS.hh.

◆ fRadiusEff

G4double G4ComponentAntiNuclNuclearXS::fRadiusEff
private

Definition at line 116 of file G4ComponentAntiNuclNuclearXS.hh.

◆ fRadiusNN2

G4double G4ComponentAntiNuclNuclearXS::fRadiusNN2
private

Definition at line 117 of file G4ComponentAntiNuclNuclearXS.hh.

◆ fTotalXsc

G4double G4ComponentAntiNuclNuclearXS::fTotalXsc
private

Definition at line 119 of file G4ComponentAntiNuclNuclearXS.hh.

◆ Mn

G4double G4ComponentAntiNuclNuclearXS::Mn
private

Definition at line 122 of file G4ComponentAntiNuclNuclearXS.hh.

◆ R0

G4double G4ComponentAntiNuclNuclearXS::R0
private

Definition at line 122 of file G4ComponentAntiNuclNuclearXS.hh.

◆ S

G4double G4ComponentAntiNuclNuclearXS::S
private

Definition at line 121 of file G4ComponentAntiNuclNuclearXS.hh.

◆ S0

G4double G4ComponentAntiNuclNuclearXS::S0
private

Definition at line 122 of file G4ComponentAntiNuclNuclearXS.hh.

◆ SqrtS

G4double G4ComponentAntiNuclNuclearXS::SqrtS
private

Definition at line 121 of file G4ComponentAntiNuclNuclearXS.hh.

◆ SqrtS0

G4double G4ComponentAntiNuclNuclearXS::SqrtS0
private

Definition at line 122 of file G4ComponentAntiNuclNuclearXS.hh.

◆ theAAlpha

G4ParticleDefinition* G4ComponentAntiNuclNuclearXS::theAAlpha
private

Definition at line 128 of file G4ComponentAntiNuclNuclearXS.hh.

◆ theADeuteron

G4ParticleDefinition* G4ComponentAntiNuclNuclearXS::theADeuteron
private

Definition at line 126 of file G4ComponentAntiNuclNuclearXS.hh.

◆ theAHe3

G4ParticleDefinition* G4ComponentAntiNuclNuclearXS::theAHe3
private

Definition at line 129 of file G4ComponentAntiNuclNuclearXS.hh.

◆ theANeutron

G4ParticleDefinition* G4ComponentAntiNuclNuclearXS::theANeutron
private

Definition at line 125 of file G4ComponentAntiNuclNuclearXS.hh.

◆ theAProton

G4ParticleDefinition* G4ComponentAntiNuclNuclearXS::theAProton
private

Definition at line 124 of file G4ComponentAntiNuclNuclearXS.hh.

◆ theATriton

G4ParticleDefinition* G4ComponentAntiNuclNuclearXS::theATriton
private

Definition at line 127 of file G4ComponentAntiNuclNuclearXS.hh.

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