IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 518bb8e3 authored by Remy Prechelt's avatar Remy Prechelt
Browse files

Remove radius from ObservationPlane output.

parent 9ac1f41c
No related branches found
No related tags found
No related merge requests found
...@@ -43,9 +43,9 @@ namespace corsika { ...@@ -43,9 +43,9 @@ namespace corsika {
auto const displacement = trajectory.getPosition(1) - plane_.getCenter(); auto const displacement = trajectory.getPosition(1) - plane_.getCenter();
// add our particles to the output file stream // add our particles to the output file stream
this->write(particle.getPID(), energy, displacement.dot(xAxis_), this->write(particle.getPID(), energy,
displacement.dot(yAxis_), displacement.dot(xAxis_),
(trajectory.getPosition(1) - plane_.getCenter()).getNorm()); displacement.dot(yAxis_));
if (deleteOnHit_) { if (deleteOnHit_) {
count_ground_++; count_ground_++;
......
...@@ -28,8 +28,6 @@ namespace corsika { ...@@ -28,8 +28,6 @@ namespace corsika {
parquet::ConvertedType::NONE); parquet::ConvertedType::NONE);
addField("y", parquet::Repetition::REQUIRED, parquet::Type::FLOAT, addField("y", parquet::Repetition::REQUIRED, parquet::Type::FLOAT,
parquet::ConvertedType::NONE); parquet::ConvertedType::NONE);
addField("radius", parquet::Repetition::REQUIRED, parquet::Type::FLOAT,
parquet::ConvertedType::NONE);
// and build the streamer // and build the streamer
buildStreamer(); buildStreamer();
...@@ -42,13 +40,12 @@ namespace corsika { ...@@ -42,13 +40,12 @@ namespace corsika {
void ObservationPlaneWriterParquet::write(Code const& pid, void ObservationPlaneWriterParquet::write(Code const& pid,
units::si::HEPEnergyType const& energy, units::si::HEPEnergyType const& energy,
units::si::LengthType const& x, units::si::LengthType const& x,
units::si::LengthType const& y, units::si::LengthType const& y) {
units::si::LengthType const& radius) {
using namespace units::si; using namespace units::si;
// write the next row - we must write `shower_` first. // write the next row - we must write `shower_` first.
(*writer_) << shower_ << static_cast<int>(get_PDG(pid)) << energy / 1_eV << x / 1_m (*writer_) << shower_ << static_cast<int>(get_PDG(pid)) << energy / 1_eV << x / 1_m
<< y / 1_m << radius / 1_m << parquet::EndRow; << y / 1_m << parquet::EndRow;
} }
} // namespace corsika } // namespace corsika
...@@ -48,8 +48,7 @@ namespace corsika { ...@@ -48,8 +48,7 @@ namespace corsika {
* Write a particle to the file. * Write a particle to the file.
*/ */
void write(Code const& pid, units::si::HEPEnergyType const& energy, void write(Code const& pid, units::si::HEPEnergyType const& energy,
units::si::LengthType const& x, units::si::LengthType const& y, units::si::LengthType const& x, units::si::LengthType const& y);
units::si::LengthType const& radius);
}; // class ObservationPlaneWriterParquet }; // class ObservationPlaneWriterParquet
......
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