Re: Particle biasing

Dear Stefan and Alfredo,

thank you for looking into my problem.

I made also a run with BIASING WHAT(3) = 1.0, but with this setting the
usimbs.f is ignored, and no biasing is applied.

Stefan:

And thank you for pointing out the inappropiate production and transport
cut-off settings.
Can you tell me what are the default cut-off enegries, I couldn't find
them in the manual.

I tried out the following settings:
EMFCUT -1D-6 1D-6 0.0 AIR @LASTMAT
PROD-CUT

EMFCUT -1D-3 1D-6 0.0 BlckHole @LASTREG

Am I right to set a high electron transport cut-off energy to a
shielding calculation?

Cheers,
David

Stefan Roesler =EDrta:
> Hi Alfredo,
>
> I looked into the problem of David and found some strange things:
>
> If I set Biasing/What(3)=3D'blank' and Sdum=3D'USER', e.g.
>
> BIASING 0.0 1. BlckHole @LASTREG USER
>
> I get a stop
>
> Abort called from DOEIMB reason .NOT.LIMPRE Run stopped!
> STOP .NOT.LIMPRE
>
> I have the impression that the code thinks there is no importance biasing
> and is missing some initialization(?)
>
> Furthermore, if I set What(3)=3D1 and use David's usimbs I get consistent
> results with the other two cases while if I set What(3)=3D2 I can confirm
> Davids observations. Taking into accound that usimbs sets FIMP either 1 or 0
> 2 it looks stange to me that this should be due to biasing-artifacts. What
> do you think?
>
> btw, to David: the default thresholds for transport and production of
> gamma and e+/e- might not be suitable (too high) as you have primaries

> down to 10keV(?)
>
>
> Cheers
> Stefan
Received on Sat Oct 18 2008 - 09:59:28 CEST