Proton energy deposit (secondaries) by bremsstrahlung

Physics Processes, Models and Cross Sections
1

I was checking energy deposit (dE/dx vs Kinetic Energy) of proton by changing physics process, and saw different energy deposit.

The energy deposit with FTFP_BERT or QGSP_BERT is very different from the one with G4EmStandardPhysics only. Why?

test
test2543×1588 404 KB

2

Hello,

how you measure dEdx?

VI

3

Note that the total energy deposited by a track need not be only dE/dx.

4

When I look “TestEm18”, there is a code (as follows).

I just inject a single proton, and I take deposited energy from secondaries particles using below code. Here, I assume all particles are derived from radiation loss (bremsstrahlung)

StackingAction::ClassifyNewTrack(const G4Track* track)
{
   ...
   if (charged) {
    fRunaction->AddChargedSecondary(energy);
    analysisManager->FillH1(4,energy);
   } else {
    fRunaction->AddNeutralSecondary(energy);
    analysisManager->FillH1(5,energy);
   }
5

Note that the total energy deposited by a track need not be only dE/dx.

This means …

In the left, I’m looking “simple” radiation/brems loss (by looking all secondaries which are generated by the radiation/brems loss) This is because physics process I register is only “EM”_standard.

In the right, I’m including “other” physics effect (not only radiation loss) which relates to a proton. This is because I register other physics processes (not only “EM”). To see whether it is derived from the radiation loss, I have to follow each secondaries.

Am I correct?

6

Yes, of course ! TestEm18 assumes that there are only Em processes registered.

7

Dear all
Thanks for the comments.