Commit 5f9bdc60 authored by Håkan Wennlöf's avatar Håkan Wennlöf
Browse files

Merge branch 'genie_energy_unit_correction' into 'master' 

PrimariesGeneratorAction - Adjust energy unit

See merge request allpix-squared/allpix-squared!1152
parents ae01b421 dae8e3c1

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();

    }