diff --git a/corsika/detail/modules/sibyll/Decay.inl b/corsika/detail/modules/sibyll/Decay.inl index d6a64ef7780116169d5f63560fcce976e4fc9a15..9d69ecb59158067de85b764e0f6ffc4c97bb1206 100644 --- a/corsika/detail/modules/sibyll/Decay.inl +++ b/corsika/detail/modules/sibyll/Decay.inl @@ -187,16 +187,18 @@ namespace corsika::sibyll { int inputSibPID = sibyll::convertToSibyllRaw(pCode); // particle momentum format: px, py, pz, e, mass. units: GeV 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; inputMomentum[3] = projectile.getEnergy() / 1_GeV; inputMomentum[4] = get_mass(pCode) / 1_GeV; int nFinalParticles; - //double* outputMomentum = new double [10*5]; - double outputMomentum[5*10]; + // double* outputMomentum = new double [10*5]; + double outputMomentum[5 * 10]; int outputSibPID[10]; // 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: {}", nFinalParticles);