energy resolution vs pseudorapidity

From: Ertan Arikan <earikan_at_nigde.edu.tr>
Date: Thu, 24 Nov 2011 12:09:38 +0200

This is a multi-part message in MIME format.

--Boundary_(ID_1vRFp5uq7Q5Xgrse3fBfjg)
Content-type: text/plain; CHARSET=US-ASCII
Content-transfer-encoding: 7BIT
Content-disposition: inline

Dear FLUKA experts,

I want to calculate the energy resolution vs pseudorapidity of a calorimeter.

as you can see from the attached input file I used lattice card.

Questions
1. In the eventbin (or usrbin) card I used one eventbin (or usrbin) card for prototype and one for lattice and then energy
deposition results are in separate lines in the eventbin output. I used a root file, which is also attached, to find the total
energy resolution but it has a peak end of the tail which I couldn't understand why? Can you give me an advice for reading my
eventbin (or usrbin) card result?
2. How can I use either eventbin (or usrbin) card for pseudorapidity ?

Best Regards

Ertan ARIKAN

--Boundary_(ID_1vRFp5uq7Q5Xgrse3fBfjg)
Content-type: application/x-fluka; NAME=5_HaDbar-E.inp
Content-transfer-encoding: base64
Content-disposition: attachment; filename=5_HaDbar-E.inp

VElUTEUKCkdMT0JBTCAgICAgICAgICAgICAgICAgICAgIC0xLiAgICAgICAgICAgICAgICAg
MC4wICAgICAgICAxLgpERUZBVUxUUyAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgQ0FMT1JJTUUKQkVBTSAgICAgICAgICAg
IC01LjAgICAgICAgMC4wICAgICAgIDAuMCAgICAgICAgICAgICAgICAgICAgICAgICAgLTEu
MFBJT04rCkJFQU1QT1MgICAgICAgICAgMC4wICAgICAgIDAuMCAgICAgICAwLjAgICAgICAg
ICAgICAgICAgICAxLgpQSE9UT05VQyAgICAgICAxMTExLiAgICAgICAgICAgICAgICAgICAg
ICBCTENLSE9MRSAgQExBU1RNQVQgICAgICAgIDEuTVVNVVBBSVIKRFBNSkVUICAgICAgICAg
ICAwLjAgICAgICAgMjMuICAgICAgIDAuMApHRU9CRUdJTiAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgQ09NQk5BTUUKICAg
IDAgICAgMCAgICAgICAgICAKKiBCbGFjayBib2R5ClNQSCBibGtib2R5ICAgIDAuMCAwLjAg
MC4wIDEwMDAwMDAwLjAKKiBWb2lkIHNwaGVyZQpTUEggdm9pZCAgICAgICAwLjAgMC4wIDAu
MCAxMDAwMDAwLjAKUlBQIGhhYyAgICAgICAgLTcwLiA3MC4gMC4gMTQwLiAtMjkwLiAyOTAu
CiogMS0+IHN0ZWVsClhaUCBsYXllcjAxICAgIDAuNApYWlAgbGF5ZXIwMiAgICAwLjcKWFpQ
IGxheWVyMDMgICAgMS4xMAoqIHNjaW50ClhaUCBsYXllcjA0ICAgIDEuNDAKWFpQIGxheWVy
MDUgICAgMS44MApYWlAgbGF5ZXIwNiAgICAyLjEwCiogc3RlZWwKWFpQIGxheWVyMDcgICAg
Mi41MApYWlAgbGF5ZXIwOCAgICAyLjgwCiogc3RlZWwKWFpQIGxheWVyMDkgICAgMy4yMApY
WlAgbGF5ZXIxMCAgICAzLjUwClhaUCBsYXllcjExICAgIDMuOTAKKiBzY2ludApYWlAgbGF5
ZXIxMiAgICA0LjIwClhaUCBsYXllcjEzICAgIDQuNjAKKiBzY2ludApYWlAgbGF5ZXIxNCAg
ICA0LjkwClhaUCBsYXllcjE1ICAgIDUuMzAKKiBzY2ludApYWlAgbGF5ZXIxNiAgICA1LjYw
ClhaUCBsYXllcjE3ICAgIDYuMDAKKiBzY2ludApYWlAgbGF5ZXIxOCAgICA2LjMwClhaUCBs
YXllcjE5ICAgIDYuNzAKKiBzY2ludApYWlAgbGF5ZXIyMCAgICA3LjAwCiogc2NpbnQKWFpQ
IGxheWVyMjEgICAgMTQuMDAKWFpQIGxheWVyMjIgICAgMjEuMDAKWFpQIGxheWVyMjMgICAg
MjguMDAKWFpQIGxheWVyMjQgICAgMzUuMDAKWFpQIGxheWVyMjUgICAgNDIuMDAKWFpQIGxh
eWVyMjYgICAgNDkuMDAKWFpQIGxheWVyMjcgICAgNTYuMDAKWFpQIGxheWVyMjggICAgNjMu
MDAKWFpQIGxheWVyMjkgICAgNzAuMDAKWFpQIGxheWVyMzAgICAgNzcuMDAKWFpQIGxheWVy
MzEgICAgODQuMDAKWFpQIGxheWVyMzIgICAgOTEuMDAKWFpQIGxheWVyMzMgICAgOTguMDAK
WFpQIGxheWVyMzQgICAgMTA1LjAwClhaUCBsYXllcjM1ICAgIDExMi4wMApYWlAgbGF5ZXIz
NiAgICAxMTkuMDAKWFpQIGxheWVyMzcgICAgMTI2LjAwClhaUCBsYXllcjM4ICAgIDEzMy4w
MApYWlAgbGF5ZXIzOSAgICAxNDAuMDAKRU5ECiogQmxhY2sgaG9sZQpCTEtCT0RZICAgICAg
NSArYmxrYm9keSAtdm9pZAoqIFZvaWQgYXJvdW5kClZPSUQgICAgICAgICA1ICt2b2lkIC1o
YWMKKiBzdGVlbApyZWcxICAgICAgICAgNSAraGFjICtsYXllcjAxCnJlZzIgICAgICAgICA1
ICtoYWMgKyhsYXllcjAyIC1sYXllcjAxKQpyZWczICAgICAgICAgNSAraGFjICsobGF5ZXIw
MyAtbGF5ZXIwMiApCnJlZzQgICAgICAgICA1ICtoYWMgKyhsYXllcjA0IC1sYXllcjAzICkK
cmVnNSAgICAgICAgIDUgK2hhYyArKGxheWVyMDUgLWxheWVyMDQgKQpyZWc2ICAgICAgICAg
NSAraGFjICsobGF5ZXIwNiAtbGF5ZXIwNSApCnJlZzcgICAgICAgICA1ICtoYWMgKyhsYXll
cjA3IC1sYXllcjA2ICkKcmVnOCAgICAgICAgIDUgK2hhYyArKGxheWVyMDggLWxheWVyMDcg
KQpyZWc5ICAgICAgICAgNSAraGFjICsobGF5ZXIwOSAtbGF5ZXIwOCApCnJlZzEwICAgICAg
ICA1ICtoYWMgKyhsYXllcjEwIC1sYXllcjA5ICkKcmVnMTEgICAgICAgIDUgK2hhYyArKGxh
eWVyMTEgLWxheWVyMTAgKQpyZWcxMiAgICAgICAgNSAraGFjICsobGF5ZXIxMiAtbGF5ZXIx
MSApCnJlZzEzICAgICAgICA1ICtoYWMgKyhsYXllcjEzIC1sYXllcjEyICkKcmVnMTQgICAg
ICAgIDUgK2hhYyArKGxheWVyMTQgLWxheWVyMTMgKQpyZWcxNSAgICAgICAgNSAraGFjICso
bGF5ZXIxNSAtbGF5ZXIxNCApCnJlZzE2ICAgICAgICA1ICtoYWMgKyhsYXllcjE2IC1sYXll
cjE1ICkKcmVnMTcgICAgICAgIDUgK2hhYyArKGxheWVyMTcgLWxheWVyMTYgKQpyZWcxOCAg
ICAgICAgNSAraGFjICsobGF5ZXIxOCAtbGF5ZXIxNyApCnJlZzE5ICAgICAgICA1ICtoYWMg
KyhsYXllcjE5IC1sYXllcjE4ICkKcmVnMjAgICAgICAgIDUgK2hhYyArKGxheWVyMjAgLWxh
eWVyMTkgKQpyZXAxICAgICAgICAgNSAraGFjICsobGF5ZXIyMSAtbGF5ZXIyMCApCnJlcDIg
ICAgICAgICA1ICtoYWMgKyhsYXllcjIyIC1sYXllcjIxICkKcmVwMyAgICAgICAgIDUgK2hh
YyArKGxheWVyMjMgLWxheWVyMjIgKQpyZXA0ICAgICAgICAgNSAraGFjICsobGF5ZXIyNCAt
bGF5ZXIyMyApCnJlcDUgICAgICAgICA1ICtoYWMgKyhsYXllcjI1IC1sYXllcjI0ICkKcmVw
NiAgICAgICAgIDUgK2hhYyArKGxheWVyMjYgLWxheWVyMjUgKQpyZXA3ICAgICAgICAgNSAr
aGFjICsobGF5ZXIyNyAtbGF5ZXIyNiApCnJlcDggICAgICAgICA1ICtoYWMgKyhsYXllcjI4
IC1sYXllcjI3ICkKcmVwOSAgICAgICAgIDUgK2hhYyArKGxheWVyMjkgLWxheWVyMjggKQpy
ZXAxMCAgICAgICAgNSAraGFjICsobGF5ZXIzMCAtbGF5ZXIyOSApCnJlcDExICAgICAgICA1
ICtoYWMgKyhsYXllcjMxIC1sYXllcjMwICkKcmVwMTIgICAgICAgIDUgK2hhYyArKGxheWVy
MzIgLWxheWVyMzEgKQpyZXAxMyAgICAgICAgNSAraGFjICsobGF5ZXIzMyAtbGF5ZXIzMiAp
CnJlcDE0ICAgICAgICA1ICtoYWMgKyhsYXllcjM0IC1sYXllcjMzICkKcmVwMTUgICAgICAg
IDUgK2hhYyArKGxheWVyMzUgLWxheWVyMzQgKQpyZXAxNiAgICAgICAgNSAraGFjICsobGF5
ZXIzNiAtbGF5ZXIzNSApCnJlcDE3ICAgICAgICA1ICtoYWMgKyhsYXllcjM3IC1sYXllcjM2
ICkKcmVwMTggICAgICAgIDUgK2hhYyArKGxheWVyMzggLWxheWVyMzcpCnJlcDE5ICAgICAg
ICA1ICtoYWMgKyhsYXllcjM5IC1sYXllcjM4KQpyZXN0ICAgICAgICAgNSAraGFjIC1sYXll
cjM5CkVORApMQVRUSUNFICAgICAgICAgcmVwMSAgICAgICAgICAgICAgICAgICAgICAgICAg
IDEuMCAgICAgICAgMS4gICAgICAgIDEuaDEKTEFUVElDRSAgICAgICAgIHJlcDIgICAgICAg
ICAgICAgICAgICAgICAgICAgICAyLjAgICAgICAgIDIuICAgICAgICAxLmgyCkxBVFRJQ0Ug
ICAgICAgICByZXAzICAgICAgICAgICAgICAgICAgICAgICAgICAgMy4wICAgICAgICAzLiAg
ICAgICAgMS5oMwpMQVRUSUNFICAgICAgICAgcmVwNCAgICAgICAgICAgICAgICAgICAgICAg
ICAgIDQuMCAgICAgICAgNC4gICAgICAgIDEuaDQKTEFUVElDRSAgICAgICAgIHJlcDUgICAg
ICAgICAgICAgICAgICAgICAgICAgICA1LjAgICAgICAgIDUuICAgICAgICAxLmg1CkxBVFRJ
Q0UgICAgICAgICByZXA2ICAgICAgICAgICAgICAgICAgICAgICAgICAgNi4wICAgICAgICA2
LiAgICAgICAgMS5oNgpMQVRUSUNFICAgICAgICAgcmVwNyAgICAgICAgICAgICAgICAgICAg
ICAgICAgIDcuMCAgICAgICAgNy4gICAgICAgIDEuaDcKTEFUVElDRSAgICAgICAgIHJlcDgg
ICAgICAgICAgICAgICAgICAgICAgICAgICA4LjAgICAgICAgIDguICAgICAgICAxLmg4CkxB
VFRJQ0UgICAgICAgICByZXA5ICAgICAgICAgICAgICAgICAgICAgICAgICAgOS4wICAgICAg
ICA5LiAgICAgICAgMS5oOQpMQVRUSUNFICAgICAgICByZXAxMCAgICAgICAgICAgICAgICAg
ICAgICAgICAgMTAuMCAgICAgICAxMC4gICAgICAgIDEuaDEwCkxBVFRJQ0UgICAgICAgIHJl
cDExICAgICAgICAgICAgICAgICAgICAgICAgICAxMS4wICAgICAgIDExLiAgICAgICAgMS5o
MTEKTEFUVElDRSAgICAgICAgcmVwMTIgICAgICAgICAgICAgICAgICAgICAgICAgIDEyLjAg
ICAgICAgMTIuICAgICAgICAxLmgxMgpMQVRUSUNFICAgICAgICByZXAxMyAgICAgICAgICAg
ICAgICAgICAgICAgICAgMTMuMCAgICAgICAxMy4gICAgICAgIDEuaDEzCkxBVFRJQ0UgICAg
ICAgIHJlcDE0ICAgICAgICAgICAgICAgICAgICAgICAgICAxNC4wICAgICAgIDE0LiAgICAg
ICAgMS5oMTQKTEFUVElDRSAgICAgICAgcmVwMTUgICAgICAgICAgICAgICAgICAgICAgICAg
IDE1LjAgICAgICAgMTUuICAgICAgICAxLmgxNQpMQVRUSUNFICAgICAgICByZXAxNiAgICAg
ICAgICAgICAgICAgICAgICAgICAgIDE2LiAgICAgICAxNi4gICAgICAgIDEuaDE2CkxBVFRJ
Q0UgICAgICAgIHJlcDE3ICAgICAgICAgICAgICAgICAgICAgICAgICAxNy4wICAgICAgIDE3
LiAgICAgICAgMS5oMTcKTEFUVElDRSAgICAgICAgcmVwMTggICAgICAgICAgICAgICAgICAg
ICAgICAgIDE4LjAgICAgICAgMTguICAgICAgICAxLmgxOApMQVRUSUNFICAgICAgICByZXAx
OSAgICAgICAgICAgICAgICAgICAgICAgICAgMTkuMCAgICAgICAxOS4gICAgICAgIDEuaDE5
CiogLS0tLS0tLS0tLS0tLS0tLS0KKiBMYXR0aWNlLT5Sb3QtRGVmaQoqIC0tLS0tLS0tLS0t
LS0tLS0tLQpHRU9FTkQgICAgICAgICAgIDcwLiAgICAgIDE0MC4gICAgICAyOTAuICAgICAg
LTcwLiAgICAgICAwLjAgICAgIC0yOTAuREVCVUcKR0VPRU5EICAgICAgICAgICAgMS4gICAg
ICA0MDAuICAgICAgICAxLiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICYKUk9ULURF
RkkgICAgICAgIDIwMy4gICAgICAgMC4wICAgICAgIDAuMCAgICAgICAwLjAgICAgICAtNy4w
ICAgICAgIDAuMGgxClJPVC1ERUZJICAgICAgICAyMDQuICAgICAgIDAuMCAgICAgICAwLjAg
ICAgICAgMC4wICAgICAtMTQuMCAgICAgICAwLjBoMgpST1QtREVGSSAgICAgICAgMjA1LiAg
ICAgICAwLjAgICAgICAgMC4wICAgICAgIDAuMCAgICAgLTIxLjAgICAgICAgMC4waDMKUk9U
LURFRkkgICAgICAgIDIwNi4gICAgICAgMC4wICAgICAgIDAuMCAgICAgICAwLjAgICAgIC0y
OC4wICAgICAgIDAuMGg0ClJPVC1ERUZJICAgICAgICAyMDcuICAgICAgIDAuMCAgICAgICAw
LjAgICAgICAgMC4wICAgICAtMzUuMCAgICAgICAwLjBoNQpST1QtREVGSSAgICAgICAgMjA4
LiAgICAgICAwLjAgICAgICAgMC4wICAgICAgIDAuMCAgICAgLTQyLjAgICAgICAgMC4waDYK
Uk9ULURFRkkgICAgICAgIDIwOS4gICAgICAgMC4wICAgICAgIDAuMCAgICAgICAwLjAgICAg
IC00OS4wICAgICAgIDAuMGg3ClJPVC1ERUZJICAgICAgICAyMTAuICAgICAgIDAuMCAgICAg
ICAwLjAgICAgICAgMC4wICAgICAtNTYuMCAgICAgICAwLjBoOApST1QtREVGSSAgICAgICAg
MjExLiAgICAgICAwLjAgICAgICAgMC4wICAgICAgIDAuMCAgICAgLTYzLjAgICAgICAgMC4w
aDkKUk9ULURFRkkgICAgICAgIDIxMi4gICAgICAgMC4wICAgICAgIDAuMCAgICAgICAwLjAg
ICAgIC03MC4wICAgICAgIDAuMGgxMApST1QtREVGSSAgICAgICAgMjEzLiAgICAgICAwLjAg
ICAgICAgMC4wICAgICAgIDAuMCAgICAgLTc3LjAgICAgICAgMC4waDExClJPVC1ERUZJICAg
ICAgICAyMTQuICAgICAgIDAuMCAgICAgICAwLjAgICAgICAgMC4wICAgICAtODQuMCAgICAg
ICAwLjBoMTIKUk9ULURFRkkgICAgICAgIDIxNS4gICAgICAgMC4wICAgICAgIDAuMCAgICAg
ICAwLjAgICAgIC05MS4wICAgICAgIDAuMGgxMwpST1QtREVGSSAgICAgIDE2MDAyLiAgICAg
ICAwLjAgICAgICAgMC4wICAgICAgIDAuMCAgICAgLTk4LjAgICAgICAgMC4waDE0ClJPVC1E
RUZJICAgICAgMTcwMDIuICAgICAgIDAuMCAgICAgICAwLjAgICAgICAgMC4wICAgIC0xMDUu
MCAgICAgICAwLjBoMTUKUk9ULURFRkkgICAgICAxODAwMi4gICAgICAgMC4wICAgICAgIDAu
MCAgICAgICAwLjAgICAgLTExMi4wICAgICAgIDAuMGgxNgpST1QtREVGSSAgICAgIDE5MDAy
LiAgICAgICAwLjAgICAgICAgMC4wICAgICAgIDAuMCAgICAtMTE5LjAgICAgICAgMC4waDE3
ClJPVC1ERUZJICAgICAgMjAwMDIuICAgICAgIDAuMCAgICAgICAwLjAgICAgICAgMC4wICAg
IC0xMjYuMCAgICAgICAwLjBoMTgKUk9ULURFRkkgICAgICAyMTAwMi4gICAgICAgMC4wICAg
ICAgIDAuMCAgICAgICAwLjAgICAgLTEzMy4wICAgICAgIDAuMGgxOQoqIDIxNiBQbGFzdGlj
IFNjaW50LihWaW55bHRvbHVlbmUpCioKTUFURVJJQUwgICAgICAgICAgICAgICAgICAgICAg
ICAgICAxLjAzMiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIFBsYXN0aWNfCkNPTVBP
VU5EICAgICAgIC0uMDg1ICBIWURST0dFTiAgICAgLS45MTUgICAgQ0FSQk9OICAgICAgICAg
ICAgICAgICAgICBQbGFzdGljXwpNQVRFUklBTCAgICAgICAgIDI0LiAgIDUxLjk5NjEgICAg
ICA3LjE4ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgQ0hST01JVU0KKiBTdGFpbmxl
c3MtU3RlZWwgKHR5cGljYWwpCiogU3RhaW5sZXNzLVN0ZWVsIGlzIHByb2R1Y2VkIHdpdGgg
Q3IgY29udGVudCByYW5naW5nIGZyb20KKiA0IC0gMTkgQXRvbWljIFBlcmNlbnQsIGFuZCB3
aXRoIEMgY29udGVudCBmcm9tIDAgLSA0IHBlcmNlbnQuCk1BVEVSSUFMICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICA4LjAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBTdGFp
bmxlcwpDT01QT1VORCAgICAgICAgMTguMCAgQ0hST01JVU0gICAgICA3NC4wICAgICAgSVJP
TiAgICAgICA4LjAgICAgTklDS0VMU3RhaW5sZXMKKiAuLisuLi4uMS4uLi4rLi4uLjIuLi4u
Ky4uLi4zLi4uLisuLi4uNC4uLi4rLi4uLjUuLi4uKy4uLi42Li4uLisuLi4uNy4uCkFTU0lH
Tk1BICAgIEJMQ0tIT0xFICAgQkxLQk9EWQpBU1NJR05NQSAgICAgIFZBQ1VVTSAgICAgIFZP
SUQKQVNTSUdOTUEgICAgU3RhaW5sZXMgICAgICByZWcxCkFTU0lHTk1BICAgIFBsYXN0aWNf
ICAgICAgcmVnMgpBU1NJR05NQSAgICBTdGFpbmxlcyAgICAgIHJlZzMKQVNTSUdOTUEgICAg
UGxhc3RpY18gICAgICByZWc0CkFTU0lHTk1BICAgIFN0YWlubGVzICAgICAgcmVnNQpBU1NJ
R05NQSAgICBQbGFzdGljXyAgICAgIHJlZzYKQVNTSUdOTUEgICAgU3RhaW5sZXMgICAgICBy
ZWc3CkFTU0lHTk1BICAgIFBsYXN0aWNfICAgICAgcmVnOApBU1NJR05NQSAgICBTdGFpbmxl
cyAgICAgIHJlZzkKQVNTSUdOTUEgICAgUGxhc3RpY18gICAgIHJlZzEwCkFTU0lHTk1BICAg
IFN0YWlubGVzICAgICByZWcxMQpBU1NJR05NQSAgICBQbGFzdGljXyAgICAgcmVnMTIKQVNT
SUdOTUEgICAgU3RhaW5sZXMgICAgIHJlZzEzCkFTU0lHTk1BICAgIFBsYXN0aWNfICAgICBy
ZWcxNApBU1NJR05NQSAgICBTdGFpbmxlcyAgICAgcmVnMTUKQVNTSUdOTUEgICAgUGxhc3Rp
Y18gICAgIHJlZzE2CkFTU0lHTk1BICAgIFN0YWlubGVzICAgICByZWcxNwpBU1NJR05NQSAg
ICBQbGFzdGljXyAgICAgcmVnMTgKQVNTSUdOTUEgICAgU3RhaW5sZXMgICAgIHJlZzE5CkFT
U0lHTk1BICAgIFBsYXN0aWNfICAgICByZWcyMApBU1NJR05NQSAgICAgIFZBQ1VVTSAgICAg
IHJlc3QKRVZFTlRCSU4gICAgICAgICA4LjAgICAgRU5FUkdZICAgICAgIDMxLiAgICAgICAy
Mi4gICAgICAgMTkuICAgIDEwMDAwLmV2ZW50CkVWRU5UQklOICAgICAgICAgNC4wICAgICAg
IDEuMCAgIC0xMDAwMC4gICAgICAgIDIuICAgICAgIDEuMCAgICAgICAgMS4mCkVWRU5UQklO
ICAgICAgICAgMi4wICAgIEVORVJHWSAgICAgICAzMS4gICAgIHJlZzIwICAgICAgICAgICAg
ICAgICAgICBldlBDCkVWRU5UQklOICAgICAgICByZWcyICAgICAgICAgICAgICAgICAgICAg
ICAgICAgIDIuICAgICAgICAgICAgICAgICAgICAmClVTUkJJTiAgICAgICAgICAgIDguICAg
IEVORVJHWSAgICAgIC0zMC4gICAgICAyMi4wICAgICAgIDE5LiAgICAgIDMuMjV1c3JicmFw
ZApVU1JCSU4gICAgICAgICAgIDQuMCAgICAgICAxLjAgICAgIC0zLjI1ICAgICAgIDIuMCAg
ICAgICAgMS4gICAgICAgMS4wJgpVU1JCSU4gICAgICAgICAgICAyLiAgICBFTkVSR1kgICAg
ICAtMzAuICAgICByZWcyMCAgICAgICAgICAgICAgICAgICAgdXNyYnJhcApVU1JCSU4gICAg
ICAgICAgcmVnMiAgICAgICAgICAgICAgICAgICAgICAgICAgIDIuMCAgICAgICAgICAgICAg
ICAgICAgJgoqVVNSWUlFTEQgICAgICAgIDMxMy4gIEFMTC1QQVJUICAgICAgIDIyLiAgICAg
ICAtMS4gICAgICAgLTIuICAgICAgICAgIHBzZHJwZHR5CipVU1JZSUVMRCAgICAgICAgICA1
LiAgICAgIDFFLTYgICAgICAxMDAuICAgICAgMTgwLiAgICAgICAwLjAgICAgIDI2MDMuJgpF
VkVOVERBVCAgICAgICAgIDMzLiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgaGV2ZW50ZGF0ClJPVFBSQklOICAgICAgICAgICAgICAgICAgICBo
MSAgICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAgICAgICAgICAgICBoMiAg
ICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAgICAgICAgICAgICBoMyAgICAg
ICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAgICAgICAgICAgICBoNCAgICAgICAg
ICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAgICAgICAgICAgICBoNSAgICAgICAgICAg
ICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAgICAgICAgICAgICBoNiAgICAgICAgICAgICAg
IGV2ZW50ClJPVFBSQklOICAgICAgICAgICAgICAgICAgICBoNyAgICAgICAgICAgICAgIGV2
ZW50ClJPVFBSQklOICAgICAgICAgICAgICAgICAgICBoOCAgICAgICAgICAgICAgIGV2ZW50
ClJPVFBSQklOICAgICAgICAgICAgICAgICAgICBoOSAgICAgICAgICAgICAgIGV2ZW50ClJP
VFBSQklOICAgICAgICAgICAgICAgICAgIGgxMCAgICAgICAgICAgICAgIGV2ZW50ClJPVFBS
QklOICAgICAgICAgICAgICAgICAgIGgxMSAgICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklO
ICAgICAgICAgICAgICAgICAgIGgxMiAgICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAg
ICAgICAgICAgICAgICAgIGgxMyAgICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAg
ICAgICAgICAgICAgIGgxNCAgICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAg
ICAgICAgICAgIGgxNSAgICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAgICAg
ICAgICAgIGgxNiAgICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAgICAgICAg
ICAgIGgxNyAgICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAgICAgICAgICAg
IGgxOCAgICAgICAgICAgICAgIGV2ZW50ClJPVFBSQklOICAgICAgICAgICAgICAgICAgIGgx
OSAgICAgICAgICAgICAgIGV2ZW50ClJBTkRPTUlaICAgICAgICAgMS4wClNUQVJUICAgICAg
ICAgICA3MDAuClNUT1AK

--Boundary_(ID_1vRFp5uq7Q5Xgrse3fBfjg)
Content-type: text/x-c++src; NAME=readfluka.C
Content-transfer-encoding: 7BIT
Content-disposition: attachment; filename=readfluka.C

void readfluka() {
   // ----------------------------------------------------------------------
   // Open the Fluka file
   // ----------------------------------------------------------------------
   ifstream ffluka;
   char filename[] = "5_HaDbar-E001_fort.31";
   ffluka.open(filename, ios::in);
/// ben ekledim ///
   FILE *yaz;
   yaz=fopen("ecall-deneme.txt","w");
//// buraya kadar ///
   if (ffluka.fail()) {
     printf("ReadFluka: Could not open Fluka file %s for reading, the file perhaps does not exist.\n", filename);
     return;
   }

   // Create a histogram for the energy deposition
   TH1F* hedep = new TH1F("hedep", "Energy deposition", 200, 0, 200);
   TH1F* hesum = new TH1F("hesum", "Energy sum", 100, 0, 8);

   // Loop over the contents until we find the word "Binning"
   char line[200]; char label[200]; char bintype[200]; char dummy[200];
   float b1, b2, b3, b4, b5, b6, b7, b8, b9, b10;
   float e1, e2, e3, e4, e5, e6, e7, e8, e9, e10,sum;
   int eventno, errors=0;
   while (ffluka.getline(line, 200, '\n')) {
     int rc = sscanf(line, "%s %*s %*s \"%s %*s %*s %*s %d", label, bintype, &eventno); if (!rc) {errors++;}
     // Read an "event" type block
     if(strcmp(label, "Binning") == 0 && strcmp(bintype, "event") == 0) {
       // printf("Found Binning, type: %s, number: |%d| (%s)\n", bintype, eventno, line);
       // Preset the histogram
       // for (int i=0; i<202; i++) {hedep->SetBinContent(i,0);}
       // Read the matrix
       for (int i=0; i<19; i++) {
        ffluka.getline(line, 200, '\n');
        int rc = sscanf(line, "%f %f %f %f %f %f %f %f %f %f", &b1, &b2, &b3, &b4, &b5, &b6, &b7, &b8, &b9, &b10);
        if (!rc) {errors++;}
        // printf("i = %d\n", i);
         fprintf(yaz,"%f %f %f %f %f %f %f %f %f %f\n", b1, b2, b3, b4, b5, b6, b7, b8, b9, b10);
         hedep->AddBinContent(i*10+ 1, b1);
         hedep->AddBinContent(i*10+ 2, b2);
         hedep->AddBinContent(i*10+ 3, b3);
         hedep->AddBinContent(i*10+ 4, b4);
         hedep->AddBinContent(i*10+ 5, b5);
         hedep->AddBinContent(i*10+ 6, b6);
         hedep->AddBinContent(i*10+ 7, b7);
         hedep->AddBinContent(i*10+ 8, b8);
         hedep->AddBinContent(i*10+ 9, b9);
         hedep->AddBinContent(i*10+10, b10);
         sum = b1+b2+b3+b4+b5+b6+b7+b8+b9+b10;
}
       // Skip to the 10 additional values
       ffluka.getline(line, 200, '\n');
       ffluka.getline(line, 200, '\n');
// ffluka.getline(line, 200, '\n');
       // Read the 2 additional lines
       ffluka.getline(line, 200, '\n');
       int rc = sscanf(line, "%f %f %f %f %f %f %f %f %f %f", &e1, &e2, &e3, &e4, &e5, &e6, &e7, &e8, &e9, &e10); if (!rc) {errors++;}
       ffluka.getline(line, 200, '\n');
// int rc = sscanf(line, "%f %f %f %f %f %f %f %f %f %f", &e11, &e12, &e13, &e14, &e15, &e16, &e17, &e18, &e19, &e20); if (!rc)
{errors++;}
// ffluka.getline(line, 200, '\n');
// ben ekledim
       fprintf(yaz,"%f %f %f %f %f %f %f %f %f %f\n", e1 ,e2 ,e3 ,e4 ,e5 ,e6, e7, e8, e9, e10);
// buraya kadar
        hedep->AddBinContent(190+ 1, e1);
        hedep->AddBinContent(190+ 2, e2);
        hedep->AddBinContent(190+ 3, e3);
        hedep->AddBinContent(190+ 4, e4);
        hedep->AddBinContent(190+ 5, e5);
        hedep->AddBinContent(190+ 6, e6);
        hedep->AddBinContent(190+ 7, e7);
        hedep->AddBinContent(190+ 8, e8);
        hedep->AddBinContent(190+ 9, e9);
        hedep->AddBinContent(190+10, e10);
        sum = e1+e2+e3+e4+e5+e6+e7+e8+e9+e10;
       // Enter the total energy deposition
        hesum->Fill(sum);
  // ffluka.getline(line, 200, '\n');
     }
   }
   // Plot the event
   TCanvas* cedep = new TCanvas("cedep", "Event energy deposition", 520, 0, 500, 500); cedep->cd();
   hedep->Draw();
   hedep->GetXaxis()->SetTitle("Energy deposit [GeV]");
   hedep->GetYaxis()->SetTitle("Events");
   cedep->Update();

   // Plot the energy sum
   TCanvas* cesum = new TCanvas("cesum", "Total energy deposition", 0, 0, 500, 500); cesum->cd();
   hesum->Draw();
   hesum->GetXaxis()->SetTitle("Energy deposit [GeV]");
   hesum->GetYaxis()->SetTitle("Events");
   cesum->Update();

   // Close the file
   ffluka.close();
///// ben ekledim ////
   fclose(yaz);
///// buraya kadar ////
}

--Boundary_(ID_1vRFp5uq7Q5Xgrse3fBfjg)--
Received on Sat Nov 26 2011 - 11:19:32 CET

This archive was generated by hypermail 2.2.0 : Sat Nov 26 2011 - 11:20:02 CET