Commit 4cecbe57 authored by Sebbe Blokhuizen's avatar Sebbe Blokhuizen Committed by Simon Spannagel
Browse files

Type conversion for number of charges creation during deposition truncated... 

Type conversion for number of charges creation during deposition truncated result, causing bias to lower energies. This resulted in a peak shift towards lower energies. To fix this, the number of charges created is instead rounded to the nearest integer
parent 46b3ef24

src/modules/DepositionGeant4/SensitiveDetectorActionG4.cpp

+1 −1
Original line number Diff line number Diff line
@@ -80,7 +80,7 @@ G4bool SensitiveDetectorActionG4::ProcessHits(G4Step* step, G4TouchableHistory*) 
    // excitations via the Fano factor. We assume Gaussian statistics here.

    auto mean_charge = edep / charge_creation_energy_;

    allpix::normal_distribution<double> charge_fluctuation(mean_charge, std::sqrt(mean_charge * fano_factor_));

    auto charge = static_cast<unsigned int>(std::max(charge_fluctuation(random_generator_), 0.));

    auto charge = static_cast<unsigned int>(std::max(charge_fluctuation(random_generator_), 0.) + 0.5);



    const auto* userTrackInfo = dynamic_cast<TrackInfoG4*>(track->GetUserInformation());

    if(userTrackInfo == nullptr) {