4 #ifndef __KVWILCKEREACTIONPARAMETERS_H
5 #define __KVWILCKEREACTIONPARAMETERS_H
93 return (1.28 * pow(A, THIRD) - 0.76 + 0.8 * pow(A, -THIRD));
108 static Double_t
BSS_V0(Int_t zp, Int_t ap, Int_t zt, Int_t at)
120 Double_t
VC(Double_t* r, Double_t*)
126 if (*r == 0)
return V0;
134 static Double_t zeta1 = 1.2511;
135 static Double_t zeta0 = 2.54;
136 static Double_t Kprox = 0.0852;
139 Double_t zeta = *r -
CP -
CT;
140 Double_t dzeta = zeta - zeta0;
142 Phi = -0.5 * pow(dzeta, 2) - Kprox * pow(dzeta, 3);
145 Phi = -3.437 * exp(-zeta / 0.75);
160 Double_t R = TMath::Max(0.1, x[0]);
162 return ProxPot(&x[0], 0) +
VC(&x[0], 0) + Vcent;
168 return b *
k(e_sur_a);
174 return l /
k(e_sur_a);
180 return 10.*(TMath::Pi() / pow(
k(e_sur_a), 2)) * pow(lmax + 0.5, 2);
187 return sqrt(sigma / (10.*(TMath::Pi() / pow(
k(e_sur_a), 2)))) - 0.5;
193 return 10.*TMath::Pi() * pow(bmax, 2);
199 return sqrt(sigma / (10.*TMath::Pi()));
232 return (0.1071 * Z * Z / pow(A, THIRD) + 22.3);
243 Double_t I = (A - 2.*Z) / (1.*A);
244 return 0.9517 * (1. - 1.7826 * I * I);
249 Double_t zpzt = zp * zt;
252 D = 0.3117 * pow(zpzt, 0.2122);
255 D = 1.096 + 1.391e-04 * zpzt;
260 Double_t
Eta(Double_t e_sur_a)
const
263 return 0.15746 *
ZP *
ZT / pow(e_sur_a, 0.5);
265 Double_t
k(Double_t e_sur_a)
const
268 return 0.2187 *
AT *
AP * pow(e_sur_a, 0.5) / (1.*
AC);
273 Double_t krint =
k(*x) *
RINT;
274 Double_t eta =
Eta(*x);
275 if (krint < 2.*eta)
return TMath::Pi();
276 return 2.*asin(eta / (krint - eta));
282 return atan(sin(QP) / (cos(QP) +
AP / (1.*
AT)));
290 Double_t
Lmax(Double_t* x, Double_t*)
const
302 Double_t
SigmaR(Double_t* x, Double_t*)
const
307 Double_t
ECM(Double_t e_sur_a)
const
316 Double_t
SigmaFus(Double_t* e_sur_a, Double_t*)
const
319 if (*e_sur_a <= 0)
return 0.;
320 Double_t e =
VRB /
ECM(*e_sur_a);
323 Double_t S = TMath::Min(S1, S2);
335 totpot->SetTitle(Form(
"HIPOT l=%3.0f", l));
338 totpot->SetNpx(1000);
341 totpot->SetLineColor(kBlack);
349 totpot->GetHistogram()->GetYaxis()->SetRangeUser(min, max);
350 ((TPad*)gPad)->BuildLegend();
Description of properties and kinematics of atomic nuclei.
static Double_t hbar
hbar*c in MeV.fm
Reaction parameters for heavy-ion collisions from systematics of Wilcke et al.
Double_t LCRIT
The maximum critical angular momentum for fusion.
Double_t QuarterPointAngle(Double_t *x, Double_t *) const
TF1 * GetFusionCrossSection() const
static Double_t SWaveFusionBarrierRadius(Int_t zp, Int_t ap, Int_t zt, Int_t at)
static Double_t L_Myers
density-symmetry coefficient
static Double_t InteractionRadius(Int_t aproj, Int_t atarg)
Double_t GetMaximumAngularMomentumWithPocket()
Double_t SigmaFus(Double_t *e_sur_a, Double_t *) const
Double_t Potential(Double_t *r, Double_t *)
Double_t Eta(Double_t e_sur_a) const
Double_t GetCrossSectionFromMaxAngularMomentum(Double_t e_sur_a, Double_t lmax) const
static Double_t e2_Wilcke
e**2 = 1.438 is value used by Wilcke et al.
Double_t ECM(Double_t e_sur_a) const
TF1 * fLmax
Grazing angular momentum.
static Double_t RLDCriticalAngularMomentum(Int_t z, Int_t a)
virtual ~KVWilckeReactionParameters()
Destructor.
static Double_t MatterHalfDensityRadius(Int_t A)
TF1 * fSigmaR
Reaction cross section.
Int_t NT
Neutron number of the projectile, target.
TF1 * GetReactionCrossSection() const
Double_t GetBassReactionCrossSection(Double_t e_sur_a)
Bass reaction cross-section [mb] for incident energy [MeV/nucleon].
TF1 * fPotential
total (nuclear+coulomb) potential for heavy-ions
Double_t ProjectileLabQP(Double_t e) const
static Double_t epsilon_bar_Myers(Int_t Z, Int_t A)
epsilon_bar, Eq.(7) in W.D. Myers, Phys. Lett. B 30, 451 (1969)
static Double_t M_Myers
symmetry anharmonicity coefficient
static Double_t NLDSurfaceTensionCoefficient(Int_t Z, Int_t A)
static Double_t SharpRadius(Int_t A)
Double_t FISSIONTKE
TKE for symmetric fission of combined system.
TF1 * fSigmaFus
Fusion cross section.
TF1 * GetBSSCoulombPotential() const
Double_t Lmax(Double_t *x, Double_t *) const
Double_t V0
BSS potential at r=0.
Double_t RCTOTAL
Coulomb radius.
Double_t RBARRIER
Fusion barrier radius RB for s-waves.
Int_t AT
Mass number of the projectile, target.
TF1 * fCMThetaQuart
CM quarter point angle.
static Double_t r0_Myers
nuclear radius constant
Int_t ZT
Atomic number of the projectile, target.
Double_t GAMMA
Nuclear liquid drop surface-tension coefficient.
Double_t k(Double_t e_sur_a) const
Double_t GetFullDampingTKE() const
Double_t VC_RINT
BSS Coulomb potential at Rint.
Double_t GetCrossSectionFromMaxImpactParameter(Double_t bmax) const
static Double_t J_Myers
symmetry energy coefficient
static Double_t r0_Wilcke(Int_t aproj, Int_t atarg)
void DrawAllPotentials(Double_t l=0) const
Double_t GetImpactParameterFromAngularMomentum(Double_t e_sur_a, Double_t l) const
TF1 * GetNuclearProximityPotential() const
TF1 * GetCMQuarterPointAngle() const
Double_t SigmaR(Double_t *x, Double_t *) const
Double_t GetMaxAngularMomentumFromCrossSection(Double_t e_sur_a, Double_t sigma) const
static Double_t a2_Myers
surface energy coefficient
static Double_t mu_Wilcke
mu = 931.5 is value used by Wilcke et al.
static Double_t Q_Myers
effective surface stiffness
Double_t ProjectileLabEQP(Double_t e) const
TF1 * GetCentrifugalPotential(Double_t e_sur_a, Double_t b) const
Double_t ASYMMFISSIONTKE
TKE of completely relaxed events in strongly damped collisions.
Double_t CentrifugalPotential(Double_t *x, Double_t *l)
static Double_t ChargeRadius_Myers(Int_t Z, Int_t A)
static Double_t BSS_V0(Int_t zp, Int_t ap, Int_t zt, Int_t at)
TF1 * GetCentrifugalPotential(Double_t l) const
Double_t PotentialPocketRadius(Double_t l)
static Double_t a1_Myers
volume energy coefficient
void SetEntranceChannel(const KVNucleus &proj, const KVNucleus &targ)
(Re)set entrance channel to calculate
static Double_t delta_bar_Myers(Int_t Z, Int_t A)
delta_bar, Eq.(8) in W.D. Myers, Phys. Lett. B 30, 451 (1969)
static Double_t TKESymFiss(Int_t Z, Int_t A)
Double_t GetMaxImpactParameterFromCrossSection(Double_t sigma) const
Double_t PotentialMaximumRadius(Double_t l)
Double_t PROXFACTOR
Proximity potential factor.
Double_t GetAngularMomentumFromImpactParameter(Double_t e_sur_a, Double_t b) const
Double_t ProxPot(Double_t *r, Double_t *)
Double_t VC(Double_t *r, Double_t *)
static Double_t K_Myers
compressibility coefficient
Double_t CT
Matter half-density radii.
static Double_t c1_Myers
Coulomb energy coefficient.
TF1 * fProx
Nuclear proximity potential for heavy-ions.
TF1 * fBSS
BSS Coulomb potential for heavy-ions.
Double_t VRB
The total conservative potential at r=RB for s-waves.
KVWilckeReactionParameters()
Default constructor.
TF1 * GetTotalPotential() const