IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 36b2d310 authored by Felix Riehn's avatar Felix Riehn Committed by Remy Prechelt
Browse files

clang

parent 213ec0cd
No related branches found
No related tags found
No related merge requests found
...@@ -187,16 +187,18 @@ namespace corsika::sibyll { ...@@ -187,16 +187,18 @@ namespace corsika::sibyll {
int inputSibPID = sibyll::convertToSibyllRaw(pCode); int inputSibPID = sibyll::convertToSibyllRaw(pCode);
// particle momentum format: px, py, pz, e, mass. units: GeV // particle momentum format: px, py, pz, e, mass. units: GeV
double inputMomentum[5]; double inputMomentum[5];
QuantityVector<hepmomentum_d> input_components = projectile.getMomentum().getComponents(); QuantityVector<hepmomentum_d> input_components =
projectile.getMomentum().getComponents();
for (int idx = 0; idx < 3; ++idx) inputMomentum[idx] = input_components[idx] / 1_GeV; for (int idx = 0; idx < 3; ++idx) inputMomentum[idx] = input_components[idx] / 1_GeV;
inputMomentum[3] = projectile.getEnergy() / 1_GeV; inputMomentum[3] = projectile.getEnergy() / 1_GeV;
inputMomentum[4] = get_mass(pCode) / 1_GeV; inputMomentum[4] = get_mass(pCode) / 1_GeV;
int nFinalParticles; int nFinalParticles;
//double* outputMomentum = new double [10*5]; // double* outputMomentum = new double [10*5];
double outputMomentum[5*10]; double outputMomentum[5 * 10];
int outputSibPID[10]; int outputSibPID[10];
// run decay routine // run decay routine
decpar_(inputSibPID, inputMomentum, nFinalParticles, outputSibPID, &outputMomentum[0]); decpar_(inputSibPID, inputMomentum, nFinalParticles, outputSibPID,
&outputMomentum[0]);
CORSIKA_LOG_TRACE("Sibyll::Decay: number of final state particles: {}", CORSIKA_LOG_TRACE("Sibyll::Decay: number of final state particles: {}",
nFinalParticles); nFinalParticles);
......
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