IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 4d0df41e authored by Felix Riehn's avatar Felix Riehn Committed by Maximilian Reininghaus
Browse files

change output to kinetic energy again, with label

parent eeb0c7bd
No related branches found
No related tags found
1 merge request!490Resolve "ObservationPlanes show non-existing particles"
Pipeline #10189 skipped
...@@ -46,14 +46,13 @@ namespace corsika { ...@@ -46,14 +46,13 @@ namespace corsika {
return ProcessReturn::Ok; return ProcessReturn::Ok;
} }
HEPEnergyType const energy = HEPEnergyType const kineticEnergy = step.getEkinPost();
step.getEkinPost() + get_mass(step.getParticlePre().getPID());
Point const pointOfIntersection = step.getPositionPost(); Point const pointOfIntersection = step.getPositionPost();
Vector const displacement = pointOfIntersection - plane_.getCenter(); Vector const displacement = pointOfIntersection - plane_.getCenter();
// add our particles to the output file stream // add our particles to the output file stream
double const weight = 1.; // step.getParticlePre().getWeight(); double const weight = 1.; // step.getParticlePre().getWeight();
this->write(step.getParticlePre().getPID(), energy, displacement.dot(xAxis_), this->write(step.getParticlePre().getPID(), kineticEnergy, displacement.dot(xAxis_),
displacement.dot(yAxis_), 0_m, step.getTimePost(), weight); displacement.dot(yAxis_), 0_m, step.getTimePost(), weight);
CORSIKA_LOG_TRACE("Particle detected absorbed={}", deleteOnHit_); CORSIKA_LOG_TRACE("Particle detected absorbed={}", deleteOnHit_);
......
...@@ -29,8 +29,8 @@ namespace corsika { ...@@ -29,8 +29,8 @@ namespace corsika {
// build the schema // build the schema
output_.addField("pdg", parquet::Repetition::REQUIRED, parquet::Type::INT32, output_.addField("pdg", parquet::Repetition::REQUIRED, parquet::Type::INT32,
parquet::ConvertedType::INT_32); parquet::ConvertedType::INT_32);
output_.addField("energy", parquet::Repetition::REQUIRED, parquet::Type::FLOAT, output_.addField("kinetic_energy", parquet::Repetition::REQUIRED,
parquet::ConvertedType::NONE); parquet::Type::FLOAT, parquet::ConvertedType::NONE);
output_.addField("x", parquet::Repetition::REQUIRED, parquet::Type::FLOAT, output_.addField("x", parquet::Repetition::REQUIRED, parquet::Type::FLOAT,
parquet::ConvertedType::NONE); parquet::ConvertedType::NONE);
output_.addField("y", parquet::Repetition::REQUIRED, parquet::Type::FLOAT, output_.addField("y", parquet::Repetition::REQUIRED, parquet::Type::FLOAT,
......
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