Message: Splitting Compton Scattering Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Splitting Compton Scattering 

Forum: Biasing and Scoring
Date: 19 Feb, 2013
From: <FT>

Dear list,

I have been trying to implement compton scattering splitting by modifying the brem splitting code. I found a previous post from a few years ago giving hints on how to achieve this but I think the interface to the compton scattering process has changed and the dynamic cast from G4VParticleChange to G4ParticleChange is no longer valid.

Does anyone know how to achieve this? With what I have now I just seen to have turned off compton scattering, which the opposite of what I was hoping for. This is what I am doing at the moment.

G4VParticleChange* pCh = 0; std::vector<G4Track*> secs;

G4int i(0);
for(i=0; i<fNSplit; ++i)
{
    //Duplicate
    G4Track dup(tr);

    pCh = (pRegProcess->PostStepDoIt(dup, step));

    for(int j =0; j<pCh->GetNumberSecondaries(); ++j)
    {
         secs.push_back(new G4Track(*(particleChange->GetSecondary(j))));
    }
    secs.push_back(new G4Track(dup));
}

then change weights and add to output particle change as for brem example.

Thanks

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Splitting Compton Scattering   (Daren Sawkey - 19 Feb, 2013)
(_ None: Re: Splitting Compton Scattering   (FT - 19 Feb, 2013)
(_ None: Re: Splitting Compton Scattering   (Daren Sawkey - 19 Feb, 2013)
(_ None: Re: Splitting Compton Scattering   (FT - 20 Feb, 2013)
 Add Message Add Message
to: "Splitting Compton Scattering"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews