Back to home page

EIC code displayed by LXR

 
 

    


Warning, /geant4/examples/extended/biasing/README is written in an unsupported language. File is not indexed.

0001 
0002                Examples for event biasing: B01, B02 and B03
0003                --------------------------------------------
0004 
0005 B01, B02 and B03 applications demonstrate the usage of different variance
0006 reduction techniques supported in Geant4, or possible from the user
0007 applications.
0008 
0009 General remark to variance reduction
0010 ------------------------------------
0011 The tools provided for importance sampling (or geometrical splitting and
0012 Russian roulette) and for the weight window technique require the user to 
0013 have a good understanding of the physics in the problem. This is because 
0014 the user has to decide which particle types have to be biased, define the 
0015 cells (physical volumes, replicas) and assign importances or weight 
0016 windows to that cells. If this is not done properly it can not be 
0017 expected that the results describe a real experiment. The examples given 
0018 here only demonstrate how to use the tools technically. They don't intend 
0019 to produce physical correct results.
0020 
0021 General remark to scoring
0022 -------------------------
0023 Scoring is carried out using the built-in Multifunctional detectors. For
0024 parallel geometries this requires a special scoring physics process.
0025 See examples/extended/runAndEvent (especailly RE05) for clarification.
0026 
0027 Known problems - should not happen
0028 ----------------------------------
0029 In the following scenario it can happen that a particle is not
0030 biased and it's weight is therefore not changed even if it crosses
0031 a boundary where biasing should happen.
0032 Importance and weight window sampling create particles on boundaries 
0033 between volumes. If the GPIL method of a physical process returns 
0034 0 as step length for a particle on a boundary and if the PostStepDoIt of
0035 that process changes the direction of the particle to go back in the 
0036 former volume the biasing won't be invoked. 
0037 This will produce particles with weights that do not correspondent to the
0038 importance of the current volumes.
0039 
0040 Further information:
0041 --------------------
0042 Short description of importance sampling and scoring:
0043 https://geant4.web.cern.ch/collaboration/working_groups/geometryTransport/#development-documents (Under the Event Biasing & Tallies Section)
0044 
0045 Example B01
0046 ===========
0047 
0048 The example uses importance sampling or the weight window technique 
0049 according to an input parameter. It uses scoring in both cases. 
0050 Importance values or weight windows are defined according to the mass 
0051 geometry. In this example the weight window technique is configured such 
0052 that it behaves equivalent to importance sampling: The window is actually 
0053 not a window but simply the inverse of the importance value and only
0054 one energy region is used that covers all energies in the problem.
0055 The user may change the weight window configuration by changing the
0056 initialization of the weight window algorithm in example,cc. 
0057 Different energy bounds for the weight window technique may be specified 
0058 in B01DetectorConstruction.
0059 
0060 The executable takes one optional argument: 0 or 1. Without argument or
0061 with argument: 0, the importance sampling is applied with argument: 1,
0062 the weight window technique is applied.
0063 
0064 A modular approach is applied to the physicslist and the extension for biasing.
0065 
0066 Example B02
0067 ===========
0068 
0069 This example uses a parallel geometry to define G4GeometryCell objects
0070 for scoring and importance sampling. The output should be equivalent to B01.
0071 
0072 A modular approach is applied to the physicslist and the extension for biasing.
0073 The parallel geometry is included in this extension.
0074 
0075 Example B03
0076 ===========
0077 
0078 This example uses a parallel geometry to define G4GeometryCell objects
0079 for scoring and importance sampling. The output should be statistically 
0080 equivalent to B02 (and B01).
0081 
0082 This demonstrates a customised "flat" physics implementation with the addition
0083 of biasing. Complementary approach to the modular physics lists of B01 and B02
0084 
0085 
0086  ___________________________________________________________________________
0087 
0088 
0089                   Generic biasing examples GB01 - GB06
0090                   ------------------------------------
0091 
0092 These examples illustrate the usage of a biasing scheme implemented since
0093 version Geant4 10.0.
0094 The scheme is meant to be extensible, not limited to these six examples.
0095 
0096 Example GB01:
0097 =============
0098 
0099 This example illustrates how to bias process cross-sections in this scheme.
0100 
0101 
0102 Example GB02:
0103 =============
0104 
0105 Illustrates a force collision scheme similar to the MCNP one.
0106 
0107 
0108 Example GB03:
0109 =============
0110 
0111 Illustrates geometry based biasing.
0112 
0113 
0114 Example GB04:
0115 =============
0116 
0117 Illustrates a bremsstrahlung splitting.
0118 
0119 
0120 Example GB05:
0121 =============
0122 
0123 Illustrates a "splitting by cross-section" technique: a splitting-based
0124 technique using absorption cross-section to control the neutron population.
0125 
0126 
0127 Example GB06:
0128 =============
0129 
0130 Illustrates the usage of parallel geometries with generic biasing.
0131 
0132 Example GB07:
0133 =============
0134 
0135 Illustrates the usage of leading particle biasing with generic biasing.
0136 
0137 
0138  ___________________________________________________________________________
0139 
0140 
0141              Reverse MonteCarlo Technique example: ReverseMC01
0142              -------------------------------------------------
0143 
0144 Example ReverseMC01
0145 ===================
0146 
0147 Example illustrating the use of the Reverse Monte Carlo (RMC) mode in a Geant4
0148 application. See details in ReverseMC01/README.
0149