From: Burge, F \(Frances\) (frances.burge@diamond.ac.uk)
Date: Mon Dec 12 2005 - 12:38:11 CET
This is a multi-part message in MIME format.
------_=_NextPart_001_01C5FF10.88018FDD
Content-Type: text/plain;
charset="us-ascii"
Content-Transfer-Encoding: quoted-printable
Sorry about the unreadable message - hopefully this will be better.
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? In case the attachment does not make it through, my
USRYIELD cards are:
*photon double differential yield (per energy per polar angle)
USRYIELD 2399.0 7.0 -23.0 1.0 2.0 1000.0
PHOTYIELD =20
USRYIELD 6.0 1.0E-5 40.0 90.0 0.0 3.0
&=20
Region 1.0 is the copper target and region 2.0 the detector region which
has the same (1cm) radius as the target and is 3cm long.
=20
many thanks,
=20
Frances
=20
Frances Burge
Health Physicist
01235 778269
Diamond Light Source Ltd.
Rutherford Appleton Laboratory
Chilton
Didcot
OX11 0QX
=20
------_=_NextPart_001_01C5FF10.88018FDD
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_01C5FF10.88018FDD--
This archive was generated by hypermail 2.1.6 : Mon Dec 12 2005 - 13:11:05 CET