Re: SCORE and energy conservation

From: maestro (paolo.maestro@pi.infn.it)
Date: Mon Jul 23 2007 - 20:09:41 CEST

  • Next message: Burge, F \(Frances\): "Biasing advice needed..."

    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
    

  • Next message: Burge, F \(Frances\): "Biasing advice needed..."

    This archive was generated by hypermail 2.1.6 : Mon Jul 23 2007 - 18:28:30 CEST