Message: Re: Warning: theListOfSecondaries() is not empty Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Warning: theListOfSecondaries() is not empty 

Forum: Run Management
Re: Question Warning: theListOfSecondaries() is not empty (Tomasz Sakrejda)
Re: Note Re: Warning: theListOfSecondaries() is not empty (Tom Roberts)
Date: 02 Aug, 2010
From: Tomasz Sakrejda <Tomasz Sakrejda>

Thank you, G4Beamline is interesting and I've been playing with it to see if it does everything that I want it to, which it seems to. The statistics are also much closer to what I expect; I'm currently investigating whether they're spot-on.

Two issues come to mind, however;

First, I'm not certain that I -can- use G4BeamLine. Does it essentially just work as a top layer for Geant, feeding Geant the geometries, physics, and cuts, and managing the output?

Second, something which disturbs me is the gross difference in calculation time between my Geant4 simulation and G4BeamLine. When I ran the geant4 simulation as a submitted job it took over two hours with much larger cuts in space (effective 1MeV cuts in water- ~5mm), however running g4beamline in the appropriate geometry takes a little over seven minutes. Does G4beamline make massive optimizations somehow? I'm just baffled.

Thanks again, Tomasz Sakrejda

On Thu, 29 Jul 2010 00:38:56 GMT, Tom Roberts wrote:

> I do not get that message, using geant4-09-03-patch-01 and physics list
> QGSP. A 10 TeV mu+ generated 5493125 secondaries (139 seconds); a 100
> TeV mu+ generated 20506638 secondaries (516 seconds). I am counting
> total particles produced, not just secondaries that make it through 100m
> of water after the 10m x 10m x 1km volume of water. I used a range cut
> of 1 mm, and an eventTimeLimit of 3000 seconds. Secondaries were tracked
> and the counts include them recursively.
> I'm using G4beamline, which is how I can check simple things like this
> in just a few minutes of my time (it can also handle quite complicated
> beam lines and simple detectors). --
> Here's what I wrote to check this (units are mm and MeV/c):
>     param eventTimeLimit=3000    
>     physics QGSP
>     beam gaussian particle=mu+ meanMomentum=100000000 nEvents=1
>     box water height=10000 width=10000 length=1000000 material=H2O
>     place water z=0 front=1
>     newparticlentuple new

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

1 None: Re: Warning: theListOfSecondaries() is not empty   (Tomasz Sakrejda - 02 Aug, 2010)
(_ Question: Re: Warning: theListOfSecondaries() is not empty   (Tomasz Sakrejda - 10 Aug, 2010)
(_ Idea: Re: Warning: theListOfSecondaries() is not empty   (Tomasz Sakrejda - 12 Aug, 2010)
 Add Message Add Message
to: "Re: Warning: theListOfSecondaries() is not empty"

 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 ]