thick target spectra from USRYIELD

From: Burge, F \(Frances\) (frances.burge@diamond.ac.uk)
Date: Mon Dec 12 2005 - 11:42:50 CET

Next message: Burge, F \(Frances\): "RE: thick target spectra from USRYIELD"

Previous message: Alfredo Ferrari: "Re: Scoring of Fission Densities with USRBIN"
Next in thread: Burge, F \(Frances\): "RE: thick target spectra from USRYIELD"
Maybe reply: Burge, F \(Frances\): "RE: thick target spectra from USRYIELD"
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

This is a multi-part message in MIME format.

------_=_NextPart_001_01C5FF08.CC54B09B
Content-Type: multipart/alternative;
boundary="----_=_NextPart_002_01C5FF08.CC54B09B"

------_=_NextPart_002_01C5FF08.CC54B09B
Content-Type: text/plain;
charset="us-ascii"
Content-Transfer-Encoding: quoted-printable

Dear FLUKA users,
=20
I am trying to calculate a photon spectrum behind a 2.5cm thick
copper target, on which a 6GeV electron beam of radius 1cm is incident.
Since the results required have units of [photons/energy /cm2] I am
using the USRYIELD card to calculate a plain double differential yield
with respect to energy as the first quantity and polar angle in degrees
as the second. I have set limits on the angle between 0 and 90 degrees
as I guess this includes everything in the 'forward' direction, and
score the photons crossing from the end of the target region into a
detector region.
=20
The problem is that my output spectrum doesn't seem to include the peak
at 0.511MeV for annihilation photons. The peak I get is somewhere
between 0.1 and 0.2MeV. The published results which I am attempting to
reproduce show a clear peak at 0.511MeV. I have set transport energy
cut-offs in the peg file (still using FLUKA2003...) as:
AP=3D0.001 UP =3D 6000
AE=3D0.521 UE =3D 6000.511
=20
I have attached the input file. Can anyone tell me why I don't see the
0.511MeV peak?
=20
many thanks,
=20
Frances
=20
Frances Burge
Health Physicist
01235 778269
Diamond Light Source Ltd.
Rutherford Appleton Laboratory
Chilton
Didcot
OX11 0QX
=20
=20
=20
=20
=20

------_=_NextPart_002_01C5FF08.CC54B09B
Content-Type: text/html;
charset="us-ascii"
Content-Transfer-Encoding: quoted-printable

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD><TITLE>Message</TITLE>
<META http-equiv=3DContent-Type content=3D"text/html; =
charset=3Dus-ascii">
<META content=3D"MSHTML 6.00.2900.2769" name=3DGENERATOR></HEAD>
<BODY>
<DIV><FONT face=3DArial size=3D2><SPAN class=3D116032410-12122005>Dear =
FLUKA=20
users,</SPAN></FONT></DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2><SPAN =
class=3D116032410-12122005>   I am=20
trying to calculate a photon spectrum behind a 2.5cm thick copper =
target, on=20
which a 6GeV electron beam of radius 1cm is incident. Since the results =
required=20
have units of [photons/energy /cm2] I am using the USRYIELD card to =
calculate a=20
plain double differential yield with respect to energy as the first =
quantity and=20
polar angle in degrees as the second. I have set limits on the angle =
between 0=20
and 90 degrees as I guess this includes everything in the 'forward' =
direction,=20
and score the photons crossing from the end of the target region =
into a=20
detector region.</SPAN></FONT></DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2><SPAN class=3D116032410-12122005>The =
problem is that=20
my output spectrum doesn't seem to include the peak at =
0.511MeV for=20
annihilation photons. The peak I get is somewhere between 0.1 and =
0.2MeV. The=20
published results which I am attempting to reproduce show a clear peak =
at=20
0.511MeV. I have set transport energy cut-offs in the peg file (still =
using=20
FLUKA2003...) as:</SPAN></FONT></DIV>
<DIV><FONT face=3DArial size=3D2><SPAN =
class=3D116032410-12122005>AP=3D0.001 UP =3D=20
6000</SPAN></FONT></DIV>
<DIV><FONT face=3DArial size=3D2><SPAN =
class=3D116032410-12122005>AE=3D0.521 UE =3D=20
6000.511</SPAN></FONT></DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2><SPAN class=3D116032410-12122005>I have =
attached the=20
input file. Can anyone tell me why I don't see the 0.511MeV=20
peak?</SPAN></FONT></DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2><SPAN class=3D116032410-12122005>many=20
thanks,</SPAN></FONT></DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005>Frances</SPAN></FONT></DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2><SPAN class=3D116032410-12122005>
<DIV align=3Dleft><FONT face=3DArial size=3D2>Frances Burge</FONT></DIV>
<DIV align=3Dleft><FONT face=3DArial size=3D2>Health =
Physicist</FONT></DIV>
<DIV align=3Dleft><FONT face=3DArial size=3D2>01235 778269</FONT></DIV>
<DIV align=3Dleft><FONT face=3DArial size=3D2>Diamond Light Source =
Ltd.</FONT></DIV>
<DIV align=3Dleft><FONT face=3DArial size=3D2>Rutherford Appleton=20
Laboratory</FONT></DIV>
<DIV align=3Dleft><FONT face=3DArial size=3D2>Chilton</FONT></DIV>
<DIV align=3Dleft><FONT face=3DArial size=3D2>Didcot</FONT></DIV>
<DIV align=3Dleft><FONT face=3DArial size=3D2>OX11 0QX</FONT></DIV>
<DIV align=3Dleft><FONT face=3DArial =
size=3D2></FONT> </DIV></SPAN></FONT></DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2><SPAN=20
class=3D116032410-12122005></SPAN></FONT> </DIV></BODY></HTML>

------_=_NextPart_002_01C5FF08.CC54B09B--

------_=_NextPart_001_01C5FF08.CC54B09B
Content-Type: application/octet-stream;
name="2pt5photonyld.inp"
Content-Transfer-Encoding: base64
Content-Description: 2pt5photonyld.inp
Content-Disposition: attachment;
filename="2pt5photonyld.inp"

Ki0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0t
LS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0t
LS0tLS0tLS0tLS0tLS0tLS0tLS0NCipQaG90b24gc3BlY3RydW0gYmVoaW5kIDIuNWNtIHRoaWNr
IGNvcHBlciB0YXJnZXQNCipJbmNpZGVudCA2R2VWIGVsZWN0cm9uIGJlYW0NCipTY29yZSBwaG90
b25zL01ldi9zci9lbGVjdHJvbiBpbiBkZXRlY3RvciB3aXRoIFVTUllJRUxEDQoqLS0tLS0tLS0t
LS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0t
LS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0t
LS0tLS0tLS0tLQ0KVElUTEUNCkVsZWN0cm9uIHRyYW5zaXRpb24gY3VydmUNCioyMzQ1Njc4OSAy
MjM0NTY3ODkgMzIzNDU2Nzg5IDQyMzQ1Njc4OSA1MjM0NTY3ODkgNjIzNDU2Nzg5IDcyMzQ1Njc4
OSA4MjM0NTY3ODkNCkRFRkFVTFRTICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICBQUkVDSVNJTw0KKkFubnVsYXIgYmVhbSAtIG91dGVy
IHJhZGl1cyAxY20sIGlubmVyIHJhZGl1cyAwIGNtLg0KQkVBTSAgICAgIC02LjAgICAgICAgICAg
ICAgICAgICAgICAgICAgIDEuMCAgICAgICAwLjAgICAgICAgIC0xLjAgICAgIEVMRUNUUk9ODQpC
RUFNUE9TICAgKzAuMCAgICAgICswLjAgICAgICArMC4wDQoqMjM0NTY3ODkgMjIzNDU2Nzg5IDMy
MzQ1Njc4OSA0MjM0NTY3ODkgNTIzNDU2Nzg5IDYyMzQ1Njc4OSA3MjM0NTY3ODkgODIzNDU2Nzg5
DQpHRU9CRUdJTiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgQ09NQklOQVQgICAgICAgICANCiAgICAzICAgIDAgICAgICAgICAgRWxl
Y3Ryb24gYmVhbSBvbiBDdSB0YXJnZXQgdGhlbiBzb2Z0IHRpc3N1ZQ0KICBaQ0MgICAgMSAgKzAu
MCAgICAgKzAuMCAgICAgICsxLjANCiAgWFlQICAgIDIgICswLjANCiAgWFlQICAgIDMgICsyLjUN
CiAgWFlQICAgIDQgICs1LjUNCiAgUlBQICAgIDUgIC0yMC4wICAgICsyMC4wICAgICAtMjAuMCAg
ICAgICsyMC4wICAgICAtMjAuMCAgICArMjAuMA0KICBSUFAgICAgNiAgLTMwLjAgICAgKzMwLjAg
ICAgIC0zMC4wICAgICAgKzMwLjAgICAgIC0zMC4wICAgICszMC4wDQogIEVORA0KKjIzNDU2Nzg5
IDEyMzQ1NiAxMjM0NTYgMTIzNDU2IDEyMzQ1NiAxMjM0NTYgMTIzNDU2IDEyMzQ1NiAxMjM0NTYg
MTIzNDU2IDEyMzQ1Ng0KKiAxc3QgUkVHSU9OICAgQ1UgVEFSR0VUIC0gMWNtIHJhZGl1cyBjeWxp
bmRlciwgZXh0ZW5kcyBiZXR3ZWVuIHo9MCBhbmQgMi41Y20NCiAgQzAxICAgICAgICAgICsxICAg
ICArMyAgICAgLTINCiogMm5kIFJFR0lPTiAgIFNPRlQgVElTU1VFIERFVEVDVE9SIC0gY3lsaW5k
ZXIgb2YgMWNtIHJhZGl1cywgej0yLjVjbSB0byA1LjVjbQ0KICBUMDIgICAgICAgICAgKzEgICAg
ICs0ICAgICAtMw0KKiAzcmQgUkVHSU9OICAgVkFDVVVNDQogIFYwMiAgICAgICAgICArNSAgICAg
KzJPUiAgICs1ICAgICAtNE9SICAgKzUgICAgICs0ICAgICAtMiAgICAgLTEgICANCiogNHRoIFJF
R0lPTiBCTEFDSyBIT0xFDQogIEIwNCAgICAgICAgICArNiAgICAgLTUNCiAgRU5EDQogICAgIDEu
MEUwICAgIDkuNDJFMCAgICAgMS4wRTAgICAgIDEuMEUwDQpHRU9FTkQNCioyMzQ1Njc4OSAyMjM0
NTY3ODkgMzIzNDU2Nzg5IDQyMzQ1Njc4OSA1MjM0NTY3ODkgNjIzNDU2Nzg5IDcyMzQ1Njc4OSA4
MjM0NTY3ODkNCk1BVEVSSUFMICAgICAgICAwLjAgICAgICAgMC4wICAgICAgIDEuMCAgICAgIDI2
LjAgICAgICAgICAgICAgICAgIDAuMCBJQ1JVNA0KQ09NUE9VTkQgICAgIC0wLjEwMSAgICAgICAz
LjAgICAgLTAuMTExICAgICAgIDYuMCAgICAtMC4wMjYgICAgICAgNy4wIElDUlU0DQpDT01QT1VO
RCAgICAgLTAuNzYyICAgICAgIDguMCAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgSUNSVTQNCiogQ1UgVEFSR0VUDQpBU1NJR05NQVQgICAxMi4wICAgIDEuMA0KKiBTT0ZU
IFRJU1NVRSAoSUNSVTQpIERFVEVDVE9SDQpBU1NJR05NQVQgICAyNi4wICAgIDIuMA0KKiBWQUNV
VU0NCkFTU0lHTk1BVCAgIDIuMCAgICAgMy4wICAgICAgICANCiogQkxBQ0tIT0xFDQpBU1NJR05N
QVQgICAxLjAgICAgIDQuMCAgICANClNDT1JFICAgICAgMjExLjANCiogUHV0IHNjb3Jpbmcgb3B0
aW9ucyBpbiB0aGUgZm9sbG93aW5nIHNlY3Rpb24uLi4NCioyMzQ1Njc4OSAyMjM0NTY3ODkgMzIz
NDU2Nzg5IDQyMzQ1Njc4OSA1MjM0NTY3ODkgNjIzNDU2Nzg5IDcyMzQ1Njc4OSA4MjM0NTY3ODkN
CipwaG90b24gZG91YmxlIGRpZmZlcmVudGlhbCB5aWVsZCAocGVyIGVuZXJneSBwZXIgcG9sYXIg
YW5nbGUpDQpVU1JZSUVMRCAgICAgMjM5OS4wICAgICAgIDcuMCAgICAgLTIzLjAgICAgICAgMS4w
ICAgICAgIDIuMCAgICAxMDAwLjAgUEhPVFlJRUxEICANClVTUllJRUxEICAgICA2LjAgICAgICAg
IDEuMEUtNSAgICAgNDAuMCAgICAgIDkwLjAgICAgICAgMC4wICAgICAzLjAgICAgICYgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICANClJBTkRPTUlaRSAxLjANCipTZXQgbm8gb2Yg
aGlzdG9yaWVzDQpTVEFSVCAgICAgMTAwMC4wDQpTVE9QDQogIA0K

------_=_NextPart_001_01C5FF08.CC54B09B--

Next message: Burge, F \(Frances\): "RE: thick target spectra from USRYIELD"

Previous message: Alfredo Ferrari: "Re: Scoring of Fission Densities with USRBIN"
Next in thread: Burge, F \(Frances\): "RE: thick target spectra from USRYIELD"
Maybe reply: Burge, F \(Frances\): "RE: thick target spectra from USRYIELD"
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

This archive was generated by hypermail 2.1.6 : Mon Dec 12 2005 - 12:20:37 CET