Re: Particle biasing

From: Stefan Roesler <sroesler_at_mail.cern.ch>
Date: Mon, 20 Oct 2008 10:15:29 +0200

Hi David

One more general remark: in case of doubt on the effects of the
implemented biasing it is always useful to score fluence of the particle
types which you bias with a USRBIN mesh. It often gives a good first
indication where there CPU time goes.

Cheers
Stefan

On Sat, 18 Oct 2008, David Horvath wrote:

> 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--
>
>
> -------
>

-- 
___________________________________
Stefan Roesler
CERN, SC/RP
CH-1211 Geneva 23
Switzerland
Phone:  +41-22-7679891
Fax:    +41-22-7669639
E-mail: Stefan.Roesler_at_cern.ch
Received on Mon Oct 20 2008 - 10:25:38 CEST

This archive was generated by hypermail 2.2.0 : Mon Oct 20 2008 - 10:25:38 CEST