Message: does the /gps/hist/point have an implicit limit to the number of points allowed? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None does the /gps/hist/point have an implicit limit to the number of points allowed? 

Forum: Particles
Date: 08 Jun, 2016
From: Anthony Villano <Anthony Villano>

I have been trying to use a GPS particle source with a user-defined energy spectrum for neutrons. I have the following in my macro:

/gps/particle neutron

/gps/position -70.0 0.0 0.0 cm

/gps/pos/type Point

/gps/ang/type iso

/gps/ene/type Arb

/gps/hist/type arb

/gps/hist/file pube.dat

/gps/hist/inter Lin

I have attached the original file pube.dat, it has 1814 points. Upon running my simulation the following memory corruption error results.

*** glibc detected *** ./k100Sim: malloc(): memory corruption: 0x0000000002353b00 ***

======= Backtrace: =========

/lib64/libc.so.6[0x336f075f4e]

/lib64/libc.so.6[0x336f07a41a]

/lib64/libc.so.6(__libc_malloc+0x5c)[0x336f07ab1c]

/usr/lib64/nvidia/tls/libnvidia-tls.so.340.96[0x3381402311]

======= Memory map: ========

00400000-00466000 r-xp 00000000 00:26 93552834 /xxx/k100Sim

00665000-00667000 rw-p 00065000 00:26 93552834 /xxx/k100Sim

00667000-00668000 rw-p 00000000 00:00 0

01207000-023fd000 rw-p 00000000 00:00 0 [heap]

4076e000-4080a000 rw-p 00000000 00:00 0

417c1000-417c3000 r-xs 00000000 08:05 4456593 /export/scratch/tmp/.glZYZvcW (deleted)

If I remove all but 1024 elements, it works. In fact it seems 1024 is the maximum number of elements. Is this a known thing and that's just the maximum an energy distribution can have?

I did not find this detail in the documentation here:

https://geant4.web.cern.ch/geant4/UserDocumentation/UsersGuides/ForApplicationDeveloper/html/ch02s07.html#sect.G4GPS.Config.histograms

though I must admit I did not read it all, I did search for 'maximum' and 'point' and read the relevant tables/sections on setting up user-defined histograms.

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/06/08/14.01-22686-pube.dat

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Agree: Re: does the /gps/hist/point have an implicit limit to the number of points allowed?   (Michael H. Kelsey - 09 Jun, 2016)
 Add Message Add Message
to: "does the /gps/hist/point have an implicit limit to the number of points allowed?"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]