IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 690cfe3e authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan:
Browse files

UrQMD impact parameter range from C7 for O,N,Ar

parent 4b9c6e1a
No related branches found
No related tags found
No related merge requests found
...@@ -289,9 +289,20 @@ corsika::process::EProcessReturn UrQMD::DoInteraction(SetupProjectile& vProjecti ...@@ -289,9 +289,20 @@ corsika::process::EProcessReturn UrQMD::DoInteraction(SetupProjectile& vProjecti
} else { } else {
inputs_.prspflg = 1; inputs_.prspflg = 1;
sys_.Ap = 1; // even for non-baryons this has to be set, see vanilla UrQMD.f sys_.Ap = 1; // even for non-baryons this has to be set, see vanilla UrQMD.f
rsys_.bdist = nucrad_(targetA) + nucrad_(1) + 2 * options_.CTParam[30 - 1];
rsys_.ebeam = (projectileEnergyLab - vProjectile.GetMass()) * (1 / 1_GeV); rsys_.ebeam = (projectileEnergyLab - vProjectile.GetMass()) * (1 / 1_GeV);
switch (targetCode) {
case particles::Code::Nitrogen:
case particles::Code::Oxygen:
rsys_.bdist = 6;
break;
case particles::Code::Argon:
rsys_.bdist = 7;
break;
default:
rsys_.bdist = nucrad_(targetA) + nucrad_(1) + 2 * options_.CTParam[30 - 1];
}
if (projectileCode == particles::Code::K0Long || if (projectileCode == particles::Code::K0Long ||
projectileCode == particles::Code::K0Short) { projectileCode == particles::Code::K0Short) {
projectileCode = fBooleanDist(fRNG) ? particles::Code::K0 : particles::Code::K0Bar; projectileCode = fBooleanDist(fRNG) ? particles::Code::K0 : particles::Code::K0Bar;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment