From: maestro (paolo.maestro@pi.infn.it)
Date: Mon Jul 23 2007 - 20:09:41 CEST
Hi,
you are right the problem is in my source routine
and I discovered it is caused by a wrong round
of the real variable PBEAM when it is passed from the source fortran
routine
to a C++ function I wrote to fill primary particle generation values to
a root file.
Thanks a lot for your hint.
Cheers
Paolo
Markus Brugger wrote:
>Hi,
>
>one more update... not having the source file I repeated the simple
>example you sent by using a direct beam going through the layers, took
>the EVENTDAT and checked the outcome. The results seem to be ok (see
>attached), thus suggesting a problem in your source routine.
>
>cheers
>
>Markus
>
>
>
>
>>-----Original Message-----
>>From: maestro [mailto:paolo.maestro@pi.infn.it]
>>Sent: 19 July 2007 12:28
>>To: Markus Brugger
>>Cc: fluka-discuss@mi.infn.it; Giuseppe Battistoni
>>Subject: Re: SCORE and energy conservation
>>
>>Hi,
>>you may find in attachment the input and geometry files.
>>As you may see I am not using biasing.
>>In order to clarify my problem I have simplified the geometry.
>>I simulated two cases:
>>1. 100 GeV protons on a 10 cm thick W block
>>2. 100 GeV protons on a 10 cm thick W + 10 cm thick
>>SCIntillator blocks.
>>For both the cases I put in attachment the output files and
>>the EVT_SCOR.txt files obtained with the Eventdat.f routine.
>>For istance, look at the first event in the Wblock_EVT.SCOR.txt file:
>>the energy deposited (score 208) in region 3 (the W block)
>>is equal to Endist variables [1] + [2] (ionization + e.m.
>>shower) , as expected.
>>Now look at the first event in the W+SCIblocks_EVT.SCOR.txt file:
>>the sum of Endist[1]+Endist[2] is about 45 GeV, while the
>>scores in region 3 and 4 ~ 80 GeV !? Moreover the sum of the
>>scores in all the regions including the black hole (e.g. :
>>leaking energy) far exceeds 100 GeV (source particle energy) !?
>>Also events 3 and 5 show a similar behaviour.
>>I made some trials changing the material assigned to region 3
>>and 4, but again I observed the same kind of results.
>>Don't you find this strange ? Do you have an explanation ?
>>Please take a look at the input file, maybe there is something wrong.
>>Thanks for your help.
>>Cheers
>> Paolo
>>
>>
>>
>>
>>Markus Brugger wrote:
>>
>>
>>
>>>Hi,
>>>
>>>it's hard to give a definite answer without having the
>>>
>>>
>>input-file (in
>>
>>
>>>case the question persists, please send it to us). The fact
>>>
>>>
>>that some
>>
>>
>>>events exceed the original beam particle energy hints the following:
>>>
>>>- do you use biasing in your input (e.g., leading particle)
>>>
>>>
>>-> in this
>>
>>
>>>case an event-by-event analysis is fluctuating by principle
>>>- in general, doing this kind of analysis you might also be
>>>
>>>
>>affected by
>>
>>
>>>fission (very low for W though) and missing energy
>>>
>>>cheers
>>>
>>>Markus
>>>
>>>
>>>
>>>
>>>
>>>
>>>>-----Original Message-----
>>>>From: owner-fluka-discuss@mi.infn.it
>>>>[mailto:owner-fluka-discuss@mi.infn.it] On Behalf Of maestro
>>>>Sent: 16 July 2007 18:46
>>>>To: fluka-discuss@mi.infn.it
>>>>Subject: SCORE and energy conservation
>>>>
>>>>Hi all,
>>>>I am using the SCORE and EVENDAT cards to score the total and
>>>>electromagnetic energy deposited in a calorimeter.
>>>>My geometry consists of 70 regions and describes a sampling
>>>>calorimeter made of scintillating plates and tungsten and
>>>>
>>>>
>>preceded by
>>
>>
>>>>a carbon pre-shower target.
>>>>I attached the xxxEVT.SCOR.txt file, I obtained by running the
>>>>Eventdat.f routine on the output binary file xxxEVT.SCOR,
>>>>
>>>>
>>produced by
>>
>>
>>>>simulating 100 GeV protons impinging the calorimeter.
>>>>I noticed that in some events (for istance Ncase=2) the sum of the
>>>>energy deposits (quantity scored 208) in each region exceeds the
>>>>beam particle energy, that is the total scored energy is > 100 GeV.
>>>>Instead, the energy balance, I mean the sum of all the Endist
>>>>variables, is equal to 100 GeV.
>>>>Maybe there is something I have not understood properly, since I
>>>>thought that summing all the scores should give the beam particle
>>>>energy. Can you explain this ?
>>>>Thanks
>>>> Paolo
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>
>>>
>>>
>>>
>>
>>
>>------------------------------------------------------------------------
>>
>> Cream Monte Carlo
>> DATE: 7/20/ 7, TIME: 17:24:32
>>Number of regions: 4 Number of scored quantities: 4
>>The scored quantities are: 208 211 201 210
>>
>>NCASE, WEIPRU, ENETOT: 1 1.000 100.0
>>ENDIST:
>> 1 0.2065 2 8.6196E-02
>> 3 7.5972E-06 4 0.000
>> 5 99.71 6 0.000
>> 7 0.000 8 0.000
>> 9 0.000 10 0.000
>> 11 0.000 12 4.2633E-14
>>Quantity n. 1 (208):
>> Scoring per region:
>> 1 99.71
>> 2 0.000
>> 3 0.2766
>> 4 1.6073E-02
>>Quantity n. 2 (211):
>> Scoring per region:
>> 1 5.1212E-03
>> 2 0.000
>> 3 8.5219E-02
>> 4 9.7700E-04
>>Quantity n. 3 (201):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 0.000
>> 4 0.000
>>Quantity n. 4 (210):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 0.000
>> 4 0.000
>> 245019 0 0. 0. 0. 0.
>>
>>NCASE, WEIPRU, ENETOT: 2 1.000 100.0
>>ENDIST:
>> 1 2.679 2 27.06
>> 3 3.0900E-02 4 0.000
>> 5 67.08 6 1.339
>> 7 0.000 8 4.2621E-02
>> 9 0.000 10 0.000
>> 11 0.000 12 1.765
>>Quantity n. 1 (208):
>> Scoring per region:
>> 1 67.08
>> 2 0.000
>> 3 29.67
>> 4 0.1422
>>Quantity n. 2 (211):
>> Scoring per region:
>> 1 6.614
>> 2 0.000
>> 3 27.02
>> 4 4.0017E-02
>>Quantity n. 3 (201):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 39.00
>> 4 1.000
>>Quantity n. 4 (210):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 1.000
>> 4 0.000
>> 59688188 0 0. 0. 0. 0.
>>
>>NCASE, WEIPRU, ENETOT: 3 1.000 100.0
>>ENDIST:
>> 1 6.714 2 24.49
>> 3 3.3839E-02 4 0.000
>> 5 65.43 6 0.5939
>> 7 0.000 8 7.7341E-02
>> 9 0.000 10 0.000
>> 11 0.000 12 2.669
>>Quantity n. 1 (208):
>> Scoring per region:
>> 1 65.43
>> 2 0.000
>> 3 31.13
>> 4 0.1765
>>Quantity n. 2 (211):
>> Scoring per region:
>> 1 10.63
>> 2 0.000
>> 3 24.43
>> 4 5.5150E-02
>>Quantity n. 3 (201):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 53.00
>> 4 1.000
>>Quantity n. 4 (210):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 1.000
>> 4 0.000
>> 116355860 0 0. 0. 0. 0.
>>
>>NCASE, WEIPRU, ENETOT: 4 1.000 100.0
>>ENDIST:
>> 1 0.2067 2 6.2458E-02
>> 3 0.000 4 0.000
>> 5 99.73 6 0.000
>> 7 0.000 8 0.000
>> 9 0.000 10 0.000
>> 11 0.000 12 -4.4906E-12
>>Quantity n. 1 (208):
>> Scoring per region:
>> 1 99.73
>> 2 0.000
>> 3 0.2510
>> 4 1.8152E-02
>>Quantity n. 2 (211):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 5.9453E-02
>> 4 3.0049E-03
>>Quantity n. 3 (201):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 0.000
>> 4 0.000
>>Quantity n. 4 (210):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 0.000
>> 4 0.000
>> 116551664 0 0. 0. 0. 0.
>>
>>NCASE, WEIPRU, ENETOT: 5 1.000 100.0
>>ENDIST:
>> 1 2.710 2 47.48
>> 3 2.4131E-02 4 0.000
>> 5 48.12 6 0.3360
>> 7 0.000 8 4.3756E-02
>> 9 0.000 10 0.000
>> 11 0.000 12 1.288
>>Quantity n. 1 (208):
>> Scoring per region:
>> 1 48.12
>> 2 0.000
>> 3 50.12
>> 4 0.1353
>>Quantity n. 2 (211):
>> Scoring per region:
>> 1 1.412
>> 2 0.000
>> 3 47.44
>> 4 4.3298E-02
>>Quantity n. 3 (201):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 35.00
>> 4 1.000
>>Quantity n. 4 (210):
>> Scoring per region:
>> 1 0.000
>> 2 0.000
>> 3 1.000
>> 4 0.000
>> 223063042 0 0. 0. 0. 0.
>>
>> End of a run of 5 particles
>>
>>
-- Dr. Paolo Maestro Dipartimento di Fisica - Università di Siena Istituto di Fisica Nucleare INFN - Sezione di Pisa Tel. +39 050 2214356 Fax. +39 050 2214317 mail paolo.maestro@pi.infn.it
This archive was generated by hypermail 2.1.6 : Mon Jul 23 2007 - 18:28:30 CEST