(no subject)

From: Tomasz Cybulski <T.Cybulski_at_liv.ac.uk>
Date: Mon, 04 Jul 2011 14:28:35 +0100

Gecko/20110616 Thunderbird/3.1.11
MIME-Version: 1.0
To: fluka-discuss_at_fluka.org
Subject: How can I score double differential fluence in a piont in vacuum?
Content-Type: multipart/mixed;
  boundary="------------070808030701050305040704"
X-Greylist: Sender IP whitelisted, not delayed by milter-greylist-4.2.7
(smtp1.mi.infn.it [192.84.138.69]); Mon, 04 Jul 2011 15:30:08 +0200 (CEST)
Sender: owner-fluka-discuss_at_mi.infn.it

This is a multi-part message in MIME format.
--------------070808030701050305040704
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit

Hello,

      I have started using fluka recently and I am still full of
questions and doubts. Let me go to my point. I would like to score
double differential proton fluence for a target which looks like a
latter U pushed to lie down on the left arm. The 60 MeV proton beam
impinges from the left, enters the letter U and hits its bottom. I would
like to score proton and electron fluence, one way, at the entrance of
letter U but leaving the target. The U target is in vacuum (Input file
in attachment). I cannot use usrbdx as there is no boundary crossing. I
do not know if it is possible to have a look at it with usryield as it
is for production of particles on surface. As far as I think usrtrack
would need a region to score on.

What can I do to score fluence in this particular plane?

Thank you for your help in advance.

With best regards,
Tomasz Cybulski

--------------070808030701050305040704
Content-Type: text/plain;
  name="FC_sperical_well.inp"
Content-Transfer-Encoding: 7bit
Content-Disposition: attachment;
  filename="FC_sperical_well.inp"

#define PROTON62
#define PROTON60
TITLE
Faraday Cup - backscattered particles and secondaries
GLOBAL 1. 1.
DEFAULTS
NEW-DEFA
#if PROTON62
BEAM -0.062 1.0 0.0 -1.0PROTON
#elif PROTON60
BEAM -.06 1.0 0.0 -1.0PROTON
#endif
BEAMPOS 0.0 0.0 -10.0
GEOBEGIN
COMBNAME
     0 0
* Black body
SPH blkbody 0.0 0.0 0.0 100000.0
* VOID
SPH void 0.0 0.0 0.0 10000.0
* Cu raget
RCC target 0.0 0.0 0.0 0.0 0.0 12.0 4.0
* Bottom of the FC
SPH sphere 0.0 0.0 6.0 2.0
* Vacuum inner of the cup
*
RCC inner 0.0 0.0 6.0 0.0 0.0 -7.0 2.0
END
* Black hole
BLKBODY 5 +blkbody -void
* Void around
VACUUM 5 +void-target
                |+inner
                |+sphere
* Target
TARGET 5 +target-inner -sphere
END
GEOEND
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
*
* Material assignment
ASSIGNMA BLCKHOLE BLKBODY
ASSIGNMA VACUUM VACUUM
ASSIGNMA COPPER TARGET
* Physics and transport
* --------------------------------------
*
* Enalbling the MSC option
* -----------------------------
MULSOPT 1.
1.GLOBHAD
*
* Sefl-shielding of copper bolck
* ---------------------------------
LOW-MAT COPPER
*
* Electromagnetic Fluka - production cut - minimum for e- / e+ and photons
* --------------------------------------------------------------------
EMFCUT -0.000001 0.0000001 VACUUM @LASTMAT
PROD-CUT
*
* Setting the physics of electromagnetic interactions
* ---------------------------------------------------------
*
EMFCUT -0.000001 0.0000001 70. VACUUM @LASTREG
*
* Threshold of production of secondary protons
* ---------------------------------------------------
PART-THR -1.E-7 PROTON 0.0
*
* Scoring
* ------------
*
* Energy deposition in Copper block
* ---------------------------------------
USRBIN 10. ENERGY -21. 8.0 8.0 12.0EneDep
USRBIN -8.0 -8.0 0.0 200. 200. 200.&
* Charged Hadrons fluence
USRBIN 10. HAD-CHAR -31. 5. 5. 18.HadFLU
USRBIN -5. -5. -5. 200. 200. 200.&
*USRTRACK 1. PROTON VACUUM
* Back scattered proton fluence
* ---------------------------------
*USRBDX 101. PROTON -22. TARGET VACUUM
16.00ProtBckS
*USRBDX 0.062 0.0 100. 4.0&
* Electron flucence - 'back to vacuum'
* ----------------------------------------
*USRBDX 101. ELECTRON -23. TARGET VACUUM
16.0ElecBckS
*USRBDX 0.062 0.0 100. 4.0&
* Flucence (E, angular distribution) for eletrons and positons
* -------------------------------------------------------------------
*USRBDX 101. E+&E- -25. TARGET VACUUM
16.0ElePoBck
*USRBDX 0.062 0.000001 100. 4.0&
* Charged particle fluence (back of copper block)
*USRBDX 101. HAD-CHAR -27. TARGET VACUUM
16.0CharFlue
*USRBDX 0.062 0.0 100. 4.0&
* Proton double differencial fluence (E,angle) WRT beam direction.
* ------------------------------------------------------------------------
*USRYIELD 124. PROTON -30. TARGET VACUUM
1.0YielProt
*USRYIELD 180.0 0.0 18. 0.062 0.0 3.&
RANDOMIZ 1.0
START 10000. 0.0
STOP

--------------070808030701050305040704
Content-Type: image/gif;
  name="Geometry.gif"
Content-Transfer-Encoding: base64
Content-Disposition: attachment;
  filename="Geometry.gif"

R0lGODlhgALgAfcAAP///wAAAKCgoP8AAADAAACA/8AA/wDu7sBAAO7uACAgwP/AIACAQKCA
/4BAAP+A/wDAYADAwABggMBggACAAED/gDBggIBgAEBAQECAAAAAgIBgEIBgYIBggAAAwAAA
/wBgAOOwwEDAgGCgwGDAAGDAoIAAAIAAgGAggGBgYCAgICBAQCBAgGCAIGCAYGCAgICAQCCA
IICAgKCgoKDQ4MAgIACAgMBgAIDA4MBgwMCAAMCAYP9AAP9AQIDA//+AYP+AgMCgAMDAwMD/
wP8AAP8A//+AoMDAoP9gYAD/AP+AAP+gAIDg4KDg4KD/IMAAAMAAwKAgIKAg/4AgAIAgIIBA
IIBAgIBgwIBg/4CAAMDAAP+AQP+gQP+gYP+gcP/AwP//AP//gP//wEBAQP///0Dg0EDg0M++
umBgYJDukEDg0Jqjh5STmZJSObhzM9DQ0Oq1lf/XAMybWZCQmE9cLZuOn2Xt/szMzN7J75Ks
qa6pl3t6ZgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACwAAAAAgALgAQAI/wABCBxIsKDBgwgTKlzIsKHD
hxAjSpxIsaLFixgzatzIsaPHjyBDihxJsqTJkyhTqlzJsqXLlzBjypxJs6bNmzhz6tzJs6fP
n0CDCh1KtKjRo0iTKl3KtKnTp1CjSp1KtarVq1izat3KtavXr2DDih1LtqzZs2jTql3Ltq3b
t3Djyp1Lt67du3jz6t3Lt6/fv4ADCx5MuLDhw4gTK17MuLHjx5AjS55MubLly5gza97MubPn
z6BDix5NurTp06hTq17NurXr17Bjy55Nu7bt27hz697Nu7fv38CDCx9OvLjx48iTK1/OvLnz
59CjS59Ovbr169iza9/Ovbv37+DDi/8fT768+fPo06tfz769+/fw48ufT7++/fv48+vfz7+/
//8ABijggAQWaOCBCCao4IIMNujggxBGKOGEFFZo4YUYZqjhhhx26OF9AYQIkYgfEhgiiR2h
6JGIKhakYosOnXiiQDKWSB6JMGKUo0Y7EvRiABO1iCOQNoYHY400zggAi0QuqaSMREKZ5JQI
IenkQlI6qaSWTRokZJdFGgnmkFdeCSWQZJZJJZNUehklmltW2SWZabr5JJhhfscknDXSKeWZ
baq5Jp6DBiqnj29yueOXeY7np4uJPgopopMKSmmZPQ7045SZGtppo9dJyqaZkSZKqqWKHoop
oZOO6iqrjKL/Cip2WVqJZJpWamnopQfZyiqkpg7Z45Fxzurgp8Ymu2Kxyjbr7LPQRivttNRW
a+212Gar7bbcduvtt+CGK+645JYr1pnopqvuuuy26+678MYr77z01mvvvfjmi6+5VbL0a0v/
rhSwSgMTLDC/vfprU8EoMXySww8fjLCLCtcEcUkXk5SxxhJPrGnFNG0sksghkVxyxx4vCfJM
Jn/U8ooLo+zxyxrSXK3NGOI8rc4W8hytzxQC/azQEhLdrNEQIp2s0sem/PHOTtNILdOzUs2g
1Y1iraDWYXKNoNc2gm2g2B+SHfXYZ9+c9tRrQ932z2/DHffQc9Pt3gd4162T2bLh/+333wX5
rXdOfIf2wUN/J563QIIPflPhni3ukOKJD9S44xY7d3hElzPeeOeYywS5ZpJDBDoAn5fuOeAE
CV75Qa+HztDomG0+0emp24465brzrjrvss+enOqcE597675b7rvtywevEO2U6U4R6Kx3/rvk
rO+OfenZO08xaNJ/RPztza++0OXUL477+N5Dz1b4HcE/ffLKhw+8+chvvr78zru/Fv8YYd/8
FGcQ69FPe4FT3/IAKLuWdW8lD2wKAy0yQYmcLoHMyx76rqe/BXovYabbngBNEkGmVPB24hsh
AvFXvxbmj4WD+1GmLrZBFZKkhEs5YQhBcsEXwlCDHFyd9P9sWK4t1clLxcsgETOyPx0OxYkN
gWISzyfC+63Qhdqj3NrgpClTgRBxS0yhFJUyRoSEcYBU7N0B06dGLbYNRXtKSJbS2JIeTqWM
sCuJHWGYxcPVEIMY3OPMmjSjYaUxb39s4QNjF0hEuvGK+dtfH22CR+R9cCpw9CISE1I+SFrR
ivVboAgbucgDyqSSjLskVTK5q6fBrpMGJGAfGSnK2MUSlB5UifwqeUZVCgVQXHqeGW+JxSwW
kHuJ9BwfPclBJToTgmasCCp9KRVWSfKHxstmB/l3TUhikY3QLOD8qJkWa2qzmN10oR27ecFk
LlOPebQgOdViTvuNsphCvN86z/n/Tnd6k4TxBOM86TlMeyrRh5HUpwDZOT5/CpKHAY3iQNdS
z2MeFJ//BGRBLYrRK4JTlxFVSC8nqpV6tvGiyyxlKC35zxoOEZkvnSYKxckQmZIUk5yEJUxf
iUuFAvGkj/TmQ8VHU5Ha9KbVzOkad8pTWYbyp/d8ailjGs7WnQ+pbQlY9aKa0XzukpY/tWg7
gwjSwOUUq1m1jP2iida0VmatRW0rRdVq1gTK1S3++8oQ63pXt0aPryPta1fy6pW9ok6wcCFs
V6hqwaAiNiqK5Yoam9rJx2YlsltxplMtaxbMakV/h+XsWzyblaOK9rJ0Pa1cSIsV06rWKqy9
imtfu8qq/+CwJrOlbVKZOFRd9ra3HsmtbqFCM0bO5LYaLetw2fKyzfrWoCcEbvyW61dpBja4
10VoSoRLXaY0N7scke5Gq9pdtAxss6DUaAmd67pcttSWDV2veK1a3nIa9ZHpTahTTVlLh7qR
jcY1runqa9+C+lTAUoXvJ/uL0lkeL8EuzS53CXyUXwG4o4GMpDIRGsF0GhCbJ30n5yh8lopy
VMTaRbE/U8rV5F4TuRIlcVlM7GIbehi/8TUohj3KT+cOWMZkoTGHbdzjpcKPodys4oLnazkg
B3m8NZ6gJBXcUSQrFcJTlaaTxyLkbxJZxx/NqJWhzOSZbhksJk1uVzW8UmL2k/+fxlQeC1W6
4YtM+Mw/sXBlH5pL9441rFgGMY5hLNLNxPZq92Xvl4Oq06/6V5YAVnCZQ5uZQy9IqwjmsyOb
6McWPxWj+9Ump7WsGUtv7SSTlq2hz2a1VFvlzngxdYJaDd5XrzpqtIY1UnRtF1l/DdW1tu2t
uwXML06G13XxNX6QpWyoIJsuzbbPHI0tmWfPJdr1OSKwkFUYa2d12t8y5FuHDS5x/7XU4dbk
985d6XIzy0epbTeu440ZbA/I3iYk98wyG2yIonvei0UweS+D7wAVPIpJdjUF9T2xgyNcygIn
6r+d5vCa1vnh/eZkkvTF8Y57PF8VqXh/RH7VVJKP0D//DoAbVs7ylrv85TCPucxnTvOa2/zm
OJ85yUi+H55rXM527vfmVJ7zohv96EhPutF3zmrJ2jUmQ1e61KdO9arrPORNdzpfXRJ1q3v9
62DHOdMBntlXvqTrYU+72tM+dooXlq0sQfva5053pbc9ZT4vudnDSfS6+/3vN7/7vt9+1u3S
CPCITzzMBd9wvVocJXJXvOT/zniE5f3xRkX14SfPeb9XPmzp2sjlMZ/5G26+86hX++dLpC7R
h0WHsI586mdf9dWXjdoWGT3pqSgS2dP+90m3vYcIhTPd631ywvU98Jefc+GPy/jHn5zE+878
6jcf62R/fWOnuyTre1/s2FcW/7hzT5Yymlb5308/y51ftvGHXyzmtyn61Z9+9g+f2yMqCy+n
OX/6e9/+9wcz+mdnBNh9/neAKweAHFJsOvIXEnZ6CEh/CrghDHgR0Ad5YzR0H7eB7fJ9E0hs
gyFF3gYSJ1J9H1gzUuN6IehEI1gyIbJ8J5gz6kZ+3RZdTlGCtBeDFzInKmgYFdSCGkN9naeD
PQNvPIIRAwAXKgSEQTh7RFghFUiDFzEASfgWAMSETYh6TwiF7/Z+U6iEu4RJQqh4W2gtOlOF
VghX1TSGiFeGasMRVBgXhiWGnOeGTbNJGQE0aOgW2ANbdeiFXYN7gKgRe8iHfgRbbOh5gwh6
6NKDHv8Rh3KIhQ0zeXbYIK13hCJRiFQhAAIgEJzoiZzYiaD4iQBAih/QiaT4iaLoiUjxgmS4
iOaiNJBYFamIiqKoirdoi4yji6XIi6RYYZJXidJCNZoIFbXYi6CIjARhir5oi6G4ikYBJK9I
ERcoH1ZTjE9xi8mIjL+4jaeojLj4jEnhim0Ii0WEErMYFdrYi+LYjcr4jcd4jOM4jUHibjOk
Eti4FOIIjrw4EMzIjwA5j4knjKw3g66EjlPxi/GYi+/YjAHZivQoEdUoHtpmhPiYj0ihkAzJ
jqsYirazj/IIkQNpjjfDJ3IUhZnoFY1IaX44kg+BkthSSL8yMFRYkzZ5kzj/mZM6uZM82ZM+
+ZNA2ZMvaSWSmBKJOHcEeSy5oirm1hLpmBRBqZMfEJVUWZVWeZURGRET+RwweSlN6RJViJFP
ESdFOYkuWY/iJ5FdmIIyIZbEtXVVcZRrl5SI5jIz4Za2hhVyqXok2SGXmIdtORZlaRLkCHh0
+TWhh4kvgZeqdhXSeJYSmXVgOYCOCQBZmX/QojNmw5it5Zh7yZfUeDRd5Ij4eBaDqTGU2Je3
d5AW6JRocZoj85mgiZZVs5KAuRKc+VktGYyquYBdSZsI+ZqwSYKyOZuRmZZ2iY65WVosSVyP
yZuhmX0ksZzM2Zw3WJxz2Zs1k5i3WRLU2Zkm55yW/zmE2pkhf9mA3vk/TSaeqXeYp8aaUpiS
atE7zpmD5UmBFome8pkWc7gUM2Kf0bk07hegIPGd4ElfSoGDv+eeW7OWwPmI6hlX0ViYC3qf
IFigEUpTHLihl2h9DGqJpLkRBqoT5/dz2BmB7WmhwxeiGTGiJLoRVwiBKOqBKsqFilJ8EPoU
RxWjBjij9VejQcMnDqqVHeGiOyF/JXeiPpqaBFo1+dmafMhEj6ekS3qZLzmMLPqaARRjVFql
5dikjdcW/Cd9PeqlLfeMLIembvCMArBybOqmb/qmNPehBcKd+pmh04M4MmqmnAina9qmfxqo
gdqngwqoaQqohCpzdGoiv/95nO+zcHpapmYqqIfqp4Tap5faponqp5Qac4t6IPgXI1Gapyk3
qWdqqJyKqWoqp6yqplcHprWZIqNqZmTapd+3qYKqqi6Hq7waijb3qQJim92Jpz9WqqbKqZVa
qLuKqp2qrHMKpEHaqJj5qPJUrbbqfZmaq5rKrL26rMjqqdAak7MaqdZ6rIeaqKvKreqKptn6
qg8aptRKruVqrs13Jl8HrHmylSMmr/NKr2InI2CHr4Eopvzar/4aeApqdQJbJPparLtnsAdb
cwmrsOF6LQ1LpjEGqdd6rBvrro6KdxNRk435sDPVsRFbexVrhiH7lFLxgwLFViZ7slO3sGGz
skb/ChRShnxJyqE827Mfl7JvGBEim5c/x3uYJ7OoB7RsIxFDK2xFW3gZW2+SCRFDa5MIcZNX
i7UEUbVNCwBWK3EhdUz7Km9ux7Q8ubU5WRA7ibY66bU4CbYSylJjK7XS6RBCKRBnOxB567ZU
GVxw93S0SnBT+xBfi7dcu4eF+5RNm7iMy7JbGre7U4Bkyy1duTFda7hJeLl8qxCLq7iee7Nx
a09SOnGUWzCWy7Iiq7mX27aYq7af67cIenGSS7fEtpTrRrWoG4c+qbdr27ps67qgG7ssmXHr
OblT8ycGuXG2exCqq7t3q7Wb67u8W4iaO7sZBLeSIa3iB25fabe5m7mO//u74hu90wu83Bda
xAu5lHGxYiJHNou4qbu6aNi1jWu+46sRDYa9lsG+3FGR8NkQd8u3aYu5rEu+0nvAo5u+Tyu4
6XaPZlvABAy9Avy2CCy/wbtjsGu88Iq74Fu9FAy8zju/r3u/vAV5DGd59MbAZctutAuy40a6
g8fCKuzCMry/gxsZwzljNwwZOfxkdcvDJ8wv/EsTPcxlO/wYRXwuR+wYSRwWQ5weT3xKQRyL
KWzDP4zEU3yOL6zBKLzFLRzDx5bF5BLFUCfGz1fFlUHG5qHGZ2fG4sLGXOfG4QLHdSTH5YbG
67vEjdHEaKbHjMHHX0HHFInHkyHIRkLI2evHi/8ByCqpyIrByIPlyIkByVxhyOBhyYYHwxtc
bXYcbogcGZjsHaGseZrcxTWcxpKMGJS8FaPcv58MGa28HbFseqUsxK/8GLOcHbnce53sLbsc
EqtcUql8GMGMWlfMxL18oaecx8e8x8lMubfsGL8sHDL0L9Osv6gck6wkK2wZxrWcmWriv93M
yd8MKsXGSmsprH+8vwNam77SRVzUL148w2ZISMmrMvNsxdo8mtyMz8tcyNp8Kw56zRn8xZuM
w8+8LQR9vgZtyt7Mxbacz9m8wg/d0BH9z4nczOtczmMczY2x0CaMGiANownNxMWswxLNzOBz
0j6M0aA8GixtxCkN0Br/vcglnS0jjb83/ccxrcQzndEr3dN9/NMvbTgKrBg5PboQTTpHnRhJ
/bhLnRkKh8sezRhPPRJN7dRVvRhXzcsivdVIPcw+uNMWC9ZaXdOPTNYqS9SwLNbdptZBW9H0
fNBAzNFnzNZUjdaTDNdL69JtrdeqzNdu49d5TdHkHNUdLTMxAT2MHTMGA9iymjkhI9mT/diG
HTGULTqZrdmWTcNG6dgss9mLrdgO3TCgzdmhLdoY49aKmdqujdqvjdkU57O0Xdu2fdu4nds8
i2e83du+/dvAHdzCPdzEXdzGfdzIndzKvdzM3dzOzSNzMqQjob1BGKouuLxGGd3W3dtGdM9B
/xgy3k2Yo9PdrVTc8bwqAIPd2V3eDaPehBnO4X1m52zPuoIxyMveD2OSLyHO7T2a0r1l8w3P
/Zwi3Ms3MgneMJFJ573c27zdLmPg8d3f+03fA07cDS42/J3f+C3eG27f/OzgAF4r/x2bIL4s
JU7gJ56HIp7iz93iLv7iMB7jMj7jNF7jNn7jOJ7jOr7juMFsIw7PXc3jqqLivPI8OxfkQu4m
ohfdDVHiLJ7ko4UnAt3daJIqVj7kysspY4LkUM7PHx7PrxIo/L0ppJIrT97lzLXlfVIqhfLl
J1nkNwrnaP4XseLmwSTm90zmdt7hc34XsRLmxHIq/UzmYe7lfU7nK3MeLG8iQ23eK1Su6Fl+
6Gkc4ZKeG+5d6Zie6Zq+6Zze6Z7+6aAe6qI+6qRe6qZ+6qie6qq+6qze6q7+6rAe67I+67Re
67Z+67ie67q+67ze677+68Ae7MI+7MRe7MZ+7Mie7Mq+7Mze7M7+7NAe7dLe5QEBADs=
--------------070808030701050305040704--
Received on Mon Jul 04 2011 - 22:41:23 CEST

This archive was generated by hypermail 2.2.0 : Mon Jul 04 2011 - 22:41:24 CEST