Commit dae8e3c1 authored by Simon Spannagel's avatar Simon Spannagel
Browse files

DepositionGenerator: GENIE produces energy in GeV - convert to internal unit

parent 442281a8

src/modules/DepositionGenerator/PrimariesGeneratorAction.cpp

+1 −1
Original line number Diff line number Diff line
@@ -60,7 +60,7 @@ void PrimariesGeneratorAction::GeneratePrimaries(G4Event* event) { 


        auto* pdg_table = G4ParticleTable::GetParticleTable();

        particle_gun_->SetParticleDefinition(pdg_table->FindParticle(particle.pdg()));

        particle_gun_->SetParticleEnergy(particle.energy() * 1000);

        particle_gun_->SetParticleEnergy(particle.energy());

        particle_gun_->SetParticlePosition(particle.position());

        particle_gun_->SetParticleMomentumDirection(particle.direction());

src/modules/DepositionGenerator/PrimariesReaderGenie.cpp

+6 −2
Original line number Diff line number Diff line
@@ -98,8 +98,12 @@ std::vector<PrimariesReader::Particle> PrimariesReaderGenie::getParticles() { 
            LOG(DEBUG) << "Skipping primary particle with PDG code " << pdg_code_->At(i);

            continue;

        }

        particles.emplace_back(

            pdg_code_->At(i), energy_->At(i), G4ThreeVector(px_->At(i), py_->At(i), pz_->At(i)), G4ThreeVector(0, 0, 0), 0);

        // Nota bene: GENIE returns energy in GeV so we need to convert to MeV:

        particles.emplace_back(pdg_code_->At(i),

                               energy_->At(i) * 1000.,

                               G4ThreeVector(px_->At(i), py_->At(i), pz_->At(i)),

                               G4ThreeVector(0, 0, 0),

                               0);

        LOG(DEBUG) << "Adding particle with ID " << particles.back().pdg() << " energy " << particles.back().energy();

    }