|Message: Re: Need some help with G4StepLimitter, G4UserLimits, and other things||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
On Mon, 17 Mar 2008, Miguel Couceiro wrote: > *** Discussion title: Event and Track Management > Email replies to PublicHyperNews@slac.stanford.edu must include: > In-Reply-To: <"/eventtrackmanage/771"@geant4-hn.slac.stanford.edu> > Subject: ...change this to be about your reply. > > Hi all, > > I am performing a very simple simulation for computing the number > of primary ionization clusters produced in a gas when a pion passes > through it. I am using LowEnergyProcesses for Photons and Electrons, > with energy cuts of 100 eV, and range cuts of 1 nm. > > G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(100.0*eV, 100.0*GeV); > > To all particles that I am interested in setting a step limitation, > I had the line > > pManager->AddDiscreteProcess(new G4StepLimiter()); > > at the end of the process manager. I thought that in this way all > the other processes would be treated first, and then the step > limitation would be applied. > > In my detector construction I have inserted the line > > m_pDetector_log->SetUserLimits(new G4UserLimits(fract*m_fDetectorThickness)); > > My gas volume has a 1 mm thick, and I am using a fixed seed for my > simulations. When I comment the above line, I get, in the first event, > two steps of 0.5 mm for which of them I get “msc”. When I activate > the UserLimits with fract < 0.5, I get SteppLimitter as the process, > instead of “msc”. Besides, if I use fract = 0.1 and fract 0.01, with > the same initial seed, the total energy deposit in the gas, the point > from where the pion leaves the gas, as well as it's momentum direction > will be different. > > Since the cuts and seed are all the same, I would expect to obtain > the same results when I reduce “fact”, but with more precise local > energy deposit. Is this correct, or the step size will also influence > the physics results? > > Another question concerns the step size on the boundary surfaces > to take into account more realistic backscattering effects. How can > I implement this step limitation, or were can I find some more > information about it? (The user manual is not clear to me in this point.) > > At last, some papers discuss the possibility to use LowEnergyProceses > down to 100 eV. Since I am interested in electron avalanches, and since > gas ionization energies are of the order of 10 to 20 eV, is it possible > to use GEANT4 to study the primary ionization cluster statistics, > as well as electron avalanche statistics? Some papers also discuss the > fact that, since LowEnergyProcesses rely on EEDL data, they should be > applicable to energies down to 10 eV? However, cross sections and > stopping powers from EPDL97 and EEDL for 10 eV have large errors, > mainly EDL97. Since the command > > G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(100.0*eV, 100.0*GeV); > > sets the energy cuts for both photons and electrons, then the results > for photons shouldn't be very reliable. Is this correct, or one can > in fact set the low energy cut down to 10 eV, without significant loss > of accuracy? > > Thanks in advance, > Miguel > > Hello, This question is shown in wrong forum (should be emprocesses). First of all you need to compare simulation results to some data to similar gas mixture. It seems that there are two modles which may give you a reasonable agreement for gaseous detectors: standard and PAI. Please, consult $G4INSTALL/examples/extended/electromagnetic/TestEm7, Em8 for details. VI
|Inline Depth:||Outline Depth:||Add message:|