Re: problem in running large gamma showers

From: Giuseppe Battistoni (
Date: Thu Jan 31 2008 - 10:54:52 CET

  • Next message: Ali Koosha: "Root+Fluka"

    Hi John
    we hope that using the last FLUKA respin your problem is definitively
    fixed, and we thank you for having pointed out that hidden bug.
    In addition to this I wish to recommend 2 other issues for your
    a) in addition to yout EMF-CUT cards with blank SDUM to control transport
    cutoff of e+,e-,gamma, you should also introduce EMF-CUT cards with
    SDUM=PROD-CUT, in order to take full control of production thresholds.

    b) it has been tested (Paola Sala) that if you take away the STEPSIZE
    cards and set the first 3 WHAT of the magnetic field card
    MGFIELD to values that we have used before for atmospheric shower
    simulations ( MGFIELD 20. 100. 30. ..... ), then the messages in the logs
    concerning problems in precision will disappear. This should speed up a
    little your problem.

            Best regards

    On Fri, 30 Nov 2007, John Clem wrote:

    > All,
    > While attempting to characterize VHE gamma ray showers in the atmosphere, we
    > ran into the below problem that occured for every other primary. For this
    > case the primary energies are fixed at 20TeV.
    > "Subscript out of range on file line 357, procedure kaskad.f/kaskad.Attempt
    > to access the 5501-th element of variable agesec.equests:"
    > It became less frequent after maximizing our energy thresholds, however any
    > thing above these levels would shroud our analysis goals. Incrementing the
    > dimensions would seem like the next obvious step. Lack of access to this
    > part of the code prevents such a test. Any suggestions ??
    > Thanks, John
    > > John Clem
    > > Research Associate Professor
    > > Bartol Research Institute
    > > Department of Physics and Astronomy
    > > University of Delaware
    > > Newark, DE 19716 USA
    > > 302-831-4354

  • Next message: Ali Koosha: "Root+Fluka"

    This archive was generated by hypermail 2.1.6 : Thu Jan 31 2008 - 11:25:08 CET