Message: Biasing by physics process Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Biasing by physics process 

Forum: Biasing and Scoring
Date: 08 Aug, 2011
From: Tom Roberts <Tom Roberts>

I want to bias events by physics process, by increasing the cross-section of a given process by name, optionally limited by particle, material, and volume names. G4VProcess has SetPILfactor() which seems specifically designed for this. Unfortunately, it does not affect the track weights.

Adjusting the weights requires the current interaction length for all processes affecting this particle. Fortunately, a quick test program shows that it is available (!= -1) for all processes of interest (processes like msc, Transportation, and step limiters don't matter).

In order to adjust the track weights, between detector construction and beamOn I will loop over all particles' process managers, inserting a process wrapper around all appropriate processes. This will not handle continuous processes properly, so I'll check and fail for them. The wrapper will adjust the interaction length, and adjust the weights of the particle being tracked and all secondaries. The wrapper must do *condition=Forced in the GPIL functions so it can adjust the weight at each step.

Any comments on this? Has anyone done this already?

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

1 Question: Re: Biasing by physics process   (Peter Broere - 01 Sep, 2014)
 Add Message Add Message
to: "Biasing by physics process"

 Subscribe Subscribe

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