Re: Particle biasing

micalg=sha1; boundary="------------ms000604070400030307030603"
X-DCC-mail.t-online.hu-Metrics: mail01d.mail.t-online.hu 32720; Body=1 Fuz1=1
        Fuz2=1
X-Greylist: Sender IP whitelisted, not delayed by milter-greylist-4.0
(smtp1.mi.infn.it [192.84.138.69]); Sat, 18 Oct 2008 09:19:10 +0200 (CEST)
Sender: owner-fluka-discuss_at_mi.infn.it

This is a cryptographically signed message in MIME format.

--------------ms000604070400030307030603
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: quoted-printable

Dear Stefan and Alfredo,

thank you for looking into my problem.

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

> down to 10keV(?)
>=20
>=20
> Cheers
> Stefan

--------------ms000604070400030307030603
Content-Type: application/x-pkcs7-signature; name="smime.p7s"
Content-Transfer-Encoding: base64
Content-Disposition: attachment; filename="smime.p7s"
Content-Description: S/MIME Cryptographic Signature

MIAGCSqGSIb3DQEHAqCAMIACAQExCzAJBgUrDgMCGgUAMIAGCSqGSIb3DQEHAQAAoIIRZzCC
BN0wggPFoAMCAQICEHGS++YZX6xNEoV0cTSiGKcwDQYJKoZIhvcNAQEFBQAwezELMAkGA1UE
BhMCR0IxGzAZBgNVBAgMEkdyZWF0ZXIgTWFuY2hlc3RlcjEQMA4GA1UEBwwHU2FsZm9yZDEa
MBgGA1UECgwRQ29tb2RvIENBIExpbWl0ZWQxITAfBgNVBAMMGEFBQSBDZXJ0aWZpY2F0ZSBT
ZXJ2aWNlczAeFw0wNDAxMDEwMDAwMDBaFw0yODEyMzEyMzU5NTlaMIGuMQswCQYDVQQGEwJV
UzELMAkGA1UECBMCVVQxFzAVBgNVBAcTDlNhbHQgTGFrZSBDaXR5MR4wHAYDVQQKExVUaGUg
VVNFUlRSVVNUIE5ldHdvcmsxITAfBgNVBAsTGGh0dHA6Ly93d3cudXNlcnRydXN0LmNvbTE2
MDQGA1UEAxMtVVROLVVTRVJGaXJzdC1DbGllbnQgQXV0aGVudGljYXRpb24gYW5kIEVtYWls
MIIBIjANBgkqhkiG9w0BAQEFAAOCAQ8AMIIBCgKCAQEAsjmFpPJ9q0E7YkY3rs3BYHW8OWX5
ShpHornMSMxqmNVNNRm5pELlzkniii8efNIxB8dOtINknS4p1aJkxIW9hVE1eaROaJB7HHqk
kqgX8pgV8pPMyaQylbsMTzC9mKALi+VuG6JG+ni8om+rWV6lL8/K2m2qL+usobNqqrcuZzWL
eeEeaYji5kbNoKXqvgvOdjp6Dpvq/NonWz1zHyLmSGHGTPNpsaguG7bUMSAsvIKKjqQOpdeJ
Q/wWWq8dcdcRWdq6hw2v+vPhwvCkxWeM1tZUOt4KpLoDd7NlyP0e03RiqhjKaJMeoYV+9Udl
y/hNVyh00jT/MLbu9mIwFIws6wIDAQABo4IBJzCCASMwHwYDVR0jBBgwFoAUoBEKIz6W8Qfs
4q8p74Klf9AwpLQwHQYDVR0OBBYEFImCZ33EnSZwAEu0UEh83j2uBG59MA4GA1UdDwEB/wQE
AwIBBjAPBgNVHRMBAf8EBTADAQH/MB0GA1UdJQQWMBQGCCsGAQUFBwMCBggrBgEFBQcDBDAR
BgNVHSAECjAIMAYGBFUdIAAwewYDVR0fBHQwcjA4oDagNIYyaHR0cDovL2NybC5jb21vZG9j
YS5jb20vQUFBQ2VydGlmaWNhdGVTZXJ2aWNlcy5jcmwwNqA0oDKGMGh0dHA6Ly9jcmwuY29t
b2RvLm5ldC9BQUFDZXJ0aWZpY2F0ZVNlcnZpY2VzLmNybDARBglghkgBhvhCAQEEBAMCAQYw
DQYJKoZIhvcNAQEFBQADggEBAJ2Vyzy4fqUJxB6/C8LHdo45PJTGEKpPDMngq4RdiVTgZTvz
bRx8NywlVF+WIfw3hJGdFdwUT4HPVB1rbEVgxy35l1FM+WbKPKCCjKbI8OLp1Er57D9Wyd12
jMOCAU9sAPMeGmF0BEcDqcZAV5G8ZSLFJ2dPV9tkWtmNH7qGL/QGrpxp7en0zykX2OBKnxog
L5dMUbtGB8SKN04g4wkxaMeexIud6H4RvDJoEJYRmETYKlFgTYjrdDrfQwYyyDlWjDoRUtNB
pEMD9O3vMyfbOeAUTibJ2PU54om4k123KSZB6rObroP8d3XK6Mq1/uJlSmM+RMTQw16Hc6mY
HK9/FX8wggY/MIIFJ6ADAgECAhAT/6Em5CvlSzA5vCe/BVpFMA0GCSqGSIb3DQEBBQUAMIGu
MQswCQYDVQQGEwJVUzELMAkGA1UECBMCVVQxFzAVBgNVBAcTDlNhbHQgTGFrZSBDaXR5MR4w
HAYDVQQKExVUaGUgVVNFUlRSVVNUIE5ldHdvcmsxITAfBgNVBAsTGGh0dHA6Ly93d3cudXNl
cnRydXN0LmNvbTE2MDQGA1UEAxMtVVROLVVTRVJGaXJzdC1DbGllbnQgQXV0aGVudGljYXRp
b24gYW5kIEVtYWlsMB4XDTA4MDUyODAwMDAwMFoXDTA5MDUyODIzNTk1OVowgd0xNTAzBgNV
BAsTLENvbW9kbyBUcnVzdCBOZXR3b3JrIC0gUEVSU09OQSBOT1QgVkFMSURBVEVEMUYwRAYD
VQQLEz1UZXJtcyBhbmQgQ29uZGl0aW9ucyBvZiB1c2U6IGh0dHA6Ly93d3cuY29tb2RvLm5l
dC9yZXBvc2l0b3J5MR8wHQYDVQQLExYoYykyMDAzIENvbW9kbyBMaW1pdGVkMRYwFAYDVQQD
Ew1EYXZpZCBIb3J2YXRoMSMwIQYJKoZIhvcNAQkBFhRob3J2YXRoZEByZWFrLmJtZS5odTCC
ASIwDQYJKoZIhvcNAQEBBQADggEPADCCAQoCggEBAMjwrrHS0LxQ3+qw90316sbYugs3jXk2
vf3onb/zsV8IeyYCBGO1hiZFPPi1UfnCSc8XVUnntVCpjHY09cFB1WdBHAyXnW+PT63G+FO+
hh7hFhpflFX1JG/28nFywMpg7tRRg1pX6Ba4H1mC41wFnXmbherhLMH6CTnp/PxJ+6/F8otV
6pvZBPe837NfhGkb3+bs5PnV4Y156FZwyXpUelUuNcxrVnj8fWN5R7e3zeClSUIBJvcW8dXM
6RUGb5LHmmMQl482tvnwTiwkFZnjUl/9aSNcqtLrVSdxtksD5zwibUfbqCZiIK3qngCFDYsM
GdD/KIS1KlOa8BkgII+2hu0CAwEAAaOCAiYwggIiMB8GA1UdIwQYMBaAFImCZ33EnSZwAEu0
UEh83j2uBG59MB0GA1UdDgQWBBRUnWM8SDxiWuluUPlK3/laGD5uPzAOBgNVHQ8BAf8EBAMC
BaAwDAYDVR0TAQH/BAIwADAgBgNVHSUEGTAXBggrBgEFBQcDBAYLKwYBBAGyMQEDBQIwEQYJ
YIZIAYb4QgEBBAQDAgUgMEYGA1UdIAQ/MD0wOwYMKwYBBAGyMQECAQEBMCswKQYIKwYBBQUH
AgEWHWh0dHBzOi8vc2VjdXJlLmNvbW9kby5uZXQvQ1BTMIGlBgNVHR8EgZ0wgZowTKBKoEiG
Rmh0dHA6Ly9jcmwuY29tb2RvY2EuY29tL1VUTi1VU0VSRmlyc3QtQ2xpZW50QXV0aGVudGlj
YXRpb25hbmRFbWFpbC5jcmwwSqBIoEaGRGh0dHA6Ly9jcmwuY29tb2RvLm5ldC9VVE4tVVNF
UkZpcnN0LUNsaWVudEF1dGhlbnRpY2F0aW9uYW5kRW1haWwuY3JsMHwGCCsGAQUFBwEBBHAw
bjA2BggrBgEFBQcwAoYqaHR0cDovL2NydC5jb21vZG9jYS5jb20vVVROQUFBQ2xpZW50Q0Eu
Y3J0MDQGCCsGAQUFBzAChihodHRwOi8vY3J0LmNvbW9kby5uZXQvVVROQUFBQ2xpZW50Q0Eu
Y3J0MB8GA1UdEQQYMBaBFGhvcnZhdGhkQHJlYWsuYm1lLmh1MA0GCSqGSIb3DQEBBQUAA4IB
AQBqSTJoFMkF1XBoe2gHaI204eVci7W8JQY1YDc95SG7+P9teKk1BpZrliXSQCZOcBA+WVpQ
HnfVKv56Lsm3zGNSSbSR5k/v0Vc84fWbiKyskImWCea7mDhhgULJ5V9xxDF0NTk8locsrlit
Yr9gF/yd30qBzZ2E5IzoNtSQKymcGmy9Fpgqcsd38rchhZbSRPgt/KE6r0Nw0I8w1iYSlBRL
aZVVHsdf1gDEFj93Iq7WTCeFKc4TH23n64qgCepiUhD54sfyWBu4++6yAnXYD/w01qAXKVU9
w+4GFZOFOzxD9AeOnM+5YZ8iF4kQMaql9v/VMNyAfnWTi9bkrP1EydzwMIIGPzCCBSegAwIB
AgIQE/+hJuQr5UswObwnvwVaRTANBgkqhkiG9w0BAQUFADCBrjELMAkGA1UEBhMCVVMxCzAJ
BgNVBAgTAlVUMRcwFQYDVQQHEw5TYWx0IExha2UgQ2l0eTEeMBwGA1UEChMVVGhlIFVTRVJU
UlVTVCBOZXR3b3JrMSEwHwYDVQQLExhodHRwOi8vd3d3LnVzZXJ0cnVzdC5jb20xNjA0BgNV
BAMTLVVUTi1VU0VSRmlyc3QtQ2xpZW50IEF1dGhlbnRpY2F0aW9uIGFuZCBFbWFpbDAeFw0w
ODA1MjgwMDAwMDBaFw0wOTA1MjgyMzU5NTlaMIHdMTUwMwYDVQQLEyxDb21vZG8gVHJ1c3Qg
TmV0d29yayAtIFBFUlNPTkEgTk9UIFZBTElEQVRFRDFGMEQGA1UECxM9VGVybXMgYW5kIENv
bmRpdGlvbnMgb2YgdXNlOiBodHRwOi8vd3d3LmNvbW9kby5uZXQvcmVwb3NpdG9yeTEfMB0G
A1UECxMWKGMpMjAwMyBDb21vZG8gTGltaXRlZDEWMBQGA1UEAxMNRGF2aWQgSG9ydmF0aDEj
MCEGCSqGSIb3DQEJARYUaG9ydmF0aGRAcmVhay5ibWUuaHUwggEiMA0GCSqGSIb3DQEBAQUA
A4IBDwAwggEKAoIBAQDI8K6x0tC8UN/qsPdN9erG2LoLN415Nr396J2/87FfCHsmAgRjtYYm
RTz4tVH5wknPF1VJ57VQqYx2NPXBQdVnQRwMl51vj0+txvhTvoYe4RYaX5RV9SRv9vJxcsDK
YO7UUYNaV+gWuB9ZguNcBZ15m4Xq4SzB+gk56fz8SfuvxfKLVeqb2QT3vN+zX4RpG9/m7OT5
1eGNeehWcMl6VHpVLjXMa1Z4/H1jeUe3t83gpUlCASb3FvHVzOkVBm+Sx5pjEJePNrb58E4s
JBWZ41Jf/WkjXKrS61UncbZLA+c8Im1H26gmYiCt6p4AhQ2LDBnQ/yiEtSpTmvAZICCPtobt
AgMBAAGjggImMIICIjAfBgNVHSMEGDAWgBSJgmd9xJ0mcABLtFBIfN49rgRufTAdBgNVHQ4E
FgQUVJ1jPEg8YlrpblD5St/5Whg+bj8wDgYDVR0PAQH/BAQDAgWgMAwGA1UdEwEB/wQCMAAw
IAYDVR0lBBkwFwYIKwYBBQUHAwQGCysGAQQBsjEBAwUCMBEGCWCGSAGG+EIBAQQEAwIFIDBG
BgNVHSAEPzA9MDsGDCsGAQQBsjEBAgEBATArMCkGCCsGAQUFBwIBFh1odHRwczovL3NlY3Vy
ZS5jb21vZG8ubmV0L0NQUzCBpQYDVR0fBIGdMIGaMEygSqBIhkZodHRwOi8vY3JsLmNvbW9k
b2NhLmNvbS9VVE4tVVNFUkZpcnN0LUNsaWVudEF1dGhlbnRpY2F0aW9uYW5kRW1haWwuY3Js
MEqgSKBGhkRodHRwOi8vY3JsLmNvbW9kby5uZXQvVVROLVVTRVJGaXJzdC1DbGllbnRBdXRo
ZW50aWNhdGlvbmFuZEVtYWlsLmNybDB8BggrBgEFBQcBAQRwMG4wNgYIKwYBBQUHMAKGKmh0
dHA6Ly9jcnQuY29tb2RvY2EuY29tL1VUTkFBQUNsaWVudENBLmNydDA0BggrBgEFBQcwAoYo
aHR0cDovL2NydC5jb21vZG8ubmV0L1VUTkFBQUNsaWVudENBLmNydDAfBgNVHREEGDAWgRRo
b3J2YXRoZEByZWFrLmJtZS5odTANBgkqhkiG9w0BAQUFAAOCAQEAakkyaBTJBdVwaHtoB2iN
tOHlXIu1vCUGNWA3PeUhu/j/bXipNQaWa5Yl0kAmTnAQPllaUB531Sr+ei7Jt8xjUkm0keZP
79FXPOH1m4isrJCJlgnmu5g4YYFCyeVfccQxdDU5PJaHLK5YrWK/YBf8nd9Kgc2dhOSM6DbU
kCspnBpsvRaYKnLHd/K3IYWW0kT4LfyhOq9DcNCPMNYmEpQUS2mVVR7HX9YAxBY/dyKu1kwn
hSnOEx9t5+uKoAnqYlIQ+eLH8lgbuPvusgJ12A/8NNagFylVPcPuBhWThTs8Q/QHjpzPuWGf
IheJEDGqpfb/1TDcgH51k4vW5Kz9RMnc8DGCBFAwggRMAgEBMIHDMIGuMQswCQYDVQQGEwJV
UzELMAkGA1UECBMCVVQxFzAVBgNVBAcTDlNhbHQgTGFrZSBDaXR5MR4wHAYDVQQKExVUaGUg
VVNFUlRSVVNUIE5ldHdvcmsxITAfBgNVBAsTGGh0dHA6Ly93d3cudXNlcnRydXN0LmNvbTE2
MDQGA1UEAxMtVVROLVVTRVJGaXJzdC1DbGllbnQgQXV0aGVudGljYXRpb24gYW5kIEVtYWls
AhAT/6Em5CvlSzA5vCe/BVpFMAkGBSsOAwIaBQCgggJhMBgGCSqGSIb3DQEJAzELBgkqhkiG
9w0BBwEwHAYJKoZIhvcNAQkFMQ8XDTA4MTAxODA3MTgzOFowIwYJKoZIhvcNAQkEMRYEFCzU
A+sjuCgK8CXZvYRdAZzCPm81MFIGCSqGSIb3DQEJDzFFMEMwCgYIKoZIhvcNAwcwDgYIKoZI
hvcNAwICAgCAMA0GCCqGSIb3DQMCAgFAMAcGBSsOAwIHMA0GCCqGSIb3DQMCAgEoMIHUBgkr
BgEEAYI3EAQxgcYwgcMwga4xCzAJBgNVBAYTAlVTMQswCQYDVQQIEwJVVDEXMBUGA1UEBxMO
U2FsdCBMYWtlIENpdHkxHjAcBgNVBAoTFVRoZSBVU0VSVFJVU1QgTmV0d29yazEhMB8GA1UE
CxMYaHR0cDovL3d3dy51c2VydHJ1c3QuY29tMTYwNAYDVQQDEy1VVE4tVVNFUkZpcnN0LUNs
aWVudCBBdXRoZW50aWNhdGlvbiBhbmQgRW1haWwCEBP/oSbkK+VLMDm8J78FWkUwgdYGCyqG
SIb3DQEJEAILMYHGoIHDMIGuMQswCQYDVQQGEwJVUzELMAkGA1UECBMCVVQxFzAVBgNVBAcT
DlNhbHQgTGFrZSBDaXR5MR4wHAYDVQQKExVUaGUgVVNFUlRSVVNUIE5ldHdvcmsxITAfBgNV
BAsTGGh0dHA6Ly93d3cudXNlcnRydXN0LmNvbTE2MDQGA1UEAxMtVVROLVVTRVJGaXJzdC1D
bGllbnQgQXV0aGVudGljYXRpb24gYW5kIEVtYWlsAhAT/6Em5CvlSzA5vCe/BVpFMA0GCSqG
SIb3DQEBAQUABIIBAHEoMPNqBMq2M++Ns9PMvAIz95H16j3q6XU7FKBBusL1b3GLgfhmT5Jf
D7+PBsle5iGqfqP9Z+WRoIQzRevmXM4jXenc8PvIvJRwAUyhTZ9DCjKOkIo08wqDrSixEUUf
08n7DMdxnym1FQwbDhETRUvymbVI9QYuZ8Mbttcw7LIJArBmpOe3wKElaJ/mfdA+EICBuFwN
PzdtdhaDQCmrxZXCVpoIyo/rAplEh9kjHuYGO9LY/DmPSLKdwDj8sxdQVnwBNMPi+zSY0rcj
ufHcDfUI/VrF8h9IOWhNU0RSKrbfxMsNGypviwbdVZPd1RkUKNJvyRy3rWnqplLWiJnYGkkA
AAAAAAA=
--------------ms000604070400030307030603--

-------
Received on Sat Oct 18 2008 - 09:46:07 CEST