Message: Large number of initial particles Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Large number of initial particles 

Forum: Run Management
Date: 20 Mar, 2007
From: <mathieu.trocme@ires.in2p3.fr>

Hi, 

I'd like to count the number of recoil protons escaping a thin slice 
of polyethylene irradiating with MeV neutrons. 
Instead of using Hit & SD classes, I wanted to try the scoring classes 
(G4PSFlatSurfaceCurrent) through the "G4Run* myRunAction::GenerateRun()" mecanism 
using a large number of initial particles through "new G4ParticleGun(1000000)". 
It works fine as soon as the number of neutrons is < 1e5. 

So I took the ExampleN01, and 
in exampleN01.cc: 
- commented the event & tracking verbose commands 
- changed the numberOfEvent to 1, 
in ExN01PrimaryGeneratorAction.cc
- set n_particle to 1000000 (constructor)
- commented the 'switch' statement (GeneratePrimaries)
For 1e5 geantinos, everything's fine either on Windows XP or Scientific Linux 3.0.8.
Depending on the OS, I get a "Segmentation Fault" above 1e6/1e7 geantinos. 
I'm using G4_8.2_p01 on a 3.2 GHz machine with 1GB RAM. 

I don't use any visulisation...
http://geant4-hn.slac.stanford.edu:5090/HyperNews/public/get/runmanage/138.html

As may indicate the mulassis source code
http://reat.space.qinetiq.com/mulassis/
which would bring in ExN01Hit & ExN01SD, this may not be a proper way for counting purposes. 
If so, the interest of scoring classes would be a bit limited... 

Do I do anything wrong? 

Thanks for any help. 
Mathieu

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

1 None: Re: Large number of initial particles   (mathieu.trocme@ires.in2p3.fr - 26 Mar, 2007)
(_ None: Re: Large number of initial particles   (Mathieu Trocmé - 04 Apr, 2007)
 Add Message Add Message
to: "Large number of initial particles"

 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 ]