Command Summary

The tables in the following eight subsections summarise the command structure for the MULASSIS. All inputs are as ASCII text, each command being entered as a separate line of the macro-file (all parameters to the commands must be on the same line). The command structure follows the standard structure for Geant4 commands. There are a large number of standard Geant4 commands available to the user, in addition to those summarised in this section the user should refer to the Geant4 toolkit documents or the online help for their usage.

The /gps commands shown in this section are only a small selection of the gps commands, see the complete reference of the gps commands.

Geometry Commands

Geometry Commands

Command	Type	Description
/geometry/material/add <name> <formula> <density> [state] [temp] [pres]	G4String name, formula G4double density>0.0 G4String state G4double temp>0.0 G4double pres>0.0	Creates a material with material name name with the chemical formula formula and of density density (g/cm3). The state, temp, and pressure parameters are optional.
/geometry/material/addNIST <name>	G4String name	Creates a material from the G4NIST tables.
/geometry/material/delete name	G4String name	Deletes the material called name.
/geometry/material/list		Lists the currently defined materials and the associated index for the material.
/geometry/material/listNIST <category>	G4String category = “all”, “simple”, “compound”, or “hep”	Lists the pre-defined materials available from the G4NIST package. The string category is optional; if omitted all NIST materials are listed.
/geometry/layer/shape <type>	G4String type = “slab” or “sphere”	Sets the shield geometry type to planar (slab) or spherical. In the case of the sphere, the final layer of the geometry is a solid sphere.
/geometry/layer/add <layerIndex> <materialName> <colourIndex> <thickness> <unit>	G4int layerIndex G4String materialName G4intcolourIndex G4double thickness >0.0 G4String unit = “um”, “mm”, “cm”, “m”, “km”, “mg/cm2”, “g/cm2”, or “kg/m2”.	Adds a shield layer after layer layerIndex comprising materialName of thickness (in units) and associated colour defined by colourIndex
/geometry/delete/layer <layerIndex>	G4int layerIndex>=0	Deletes the layer of index layerIndex in the shield layer list. If layerIndex ==0, all layers are deleted.
/geometry/layer/list <layerIndex>	G4int layerIndex>=0	List the thickness of layer layerIndex, associated material and colour attributes.
/geometry/colour/add <colourName> <R> <G> <B>	G4String colourName G4double R, G, B	Defines a colour according to its RGB intensity. R, G, B must be in range [0,1]
/geometry/colour/list		Lists the available colour palette and the associated colour indices.
/geometry/default		Sets the shield thicknesses to the SHIELDOSE2 defaults
/geometry/update		Updates the geometry information after changes requested by the user. This is a mandatory command after any changes made to the geometry.
geometry/xysize <size> <unit>	G4double size G4String unit = “um”, “mm”, “cm”, “m”, or “km”	Sets the dimensions of the X/Y sides of a slab geometry. By default, the extent of the planar geometry is 100x the total thickness.

GPS Commands

General Particle Source Commands

Command	Type	Description
/gps/particle <particleType>	G4String particleType	Sets the incident particle type. For a nonMC analysis, the particle type must be ion, and the ion defined by the /gps/ion command.
/gps/ion <Z> <A> <Q> <E>	G4int Z, A, Q >0 G4double E >=0.0	For an ion source particle, sets the Z , A , excitation energy E of the nucleus and charge (Q in integer units of e) for the ion.
/gps/energy <energy> <unit>	G4double energy G4String unit	Sets the source particle distribution to monoenergetic of energy defined by energy and unit.
/gps/ene/type <energyDis>	G4String energyDis in [ “Mono”, “Lin”, “Pow”, “Exp”, “User”, “Arb”]	Sets the source particle energy distribution to: Mono - monoenergetic Lin - linear Pow - power-law Exp - exponential User - user-defined histogram Arb - arbitrary point-wise spectrum
/gps/ene/min <eMin> <unit>	G4double eMin>0.0 G4String unit	Sets the minimum for the energy distribution.
/gps/ene/max <eMax> <unit>	G4double eMax > 0.0 G4String unit	Sets the maximum for the energy distribution.
/gps/ene/alpha <alpha>	G4double alpha	Sets the \(\alpha\) for the power-law distribution of the form \(y=AE^\alpha\) where A is a constant.
/gps/ene/gradient <gradient>	G4double gradient	Sets the gradient for linear distributions
/gps/ene/ezero <E0>	G4double E0	Sets the E0 for the exponential distribution of the form \(y=A e^{(-E/E0)}\), where A is a constant.
/gps/ene/intercept <intercept>	G4double intercept	Sets the intercept for the linear distribution
/gps/resethist		Resets the histogram data.
/gps/hist/point <eHi> <weight>	G4double eHi>0.0, weight>=0.0	Allows definition of histogram, point-by-point.
/gps/hist/inter <interpType>	G4String interpType	Sets the interpolation type to be used for point-wise data: Lin - linear Log - power-law Exp - exponential Spline - cubic spline
/gps/ang/type <angDis>	G4String angDis = “iso”, “cos”, “planar”, “beam1d”, “beam2d”, “focused”, or “user”	Sets the angular distribution type iso - isotropic cos - cosine-law For a nonMC analysis, the angular distribution must be either “planar” or “cos”.
/gps/ang/mintheta <tMin> <unit>	G4double tMin G4String unit = “deg” or “rad”	Sets a minimum value for the \(\theta\) distribution.
/gps/ang/maxtheta <tMax> <unit>	G4double tMax G4String unit = “deg” or “rad”	Sets a maximum value for the \(\theta\) distribution.

Physics Commands

Physics Commands

Command	Type	Description
/phys/dEdXtoCSV <csv_file> <material> <particle_name> <Energy_lo> <Energy_hi> <Num_Energies> <log_lin>	G4String csv_file G4String material G4String particle_name G4double Energy_lo G4double Energy_hi G4int Num_Energies G4String log_lin = “log” or “lin”	Writes out to a CSV file the mean dEdX, stopping power and range for the energies between Energy_lo and Energy_hi.
/phys/printDEDXTable <material> <energy_MeV> <particle_name>	G4String material G4double energy_MeV G4String particle_name	Prints out the dEdx and range for the specified particle and material. Only works after the /run/beamOn macro has been invoked.
/phys/scenario <scenarioType>	G4String scenarioType = “em”, “em_opt3”, “em_liv”,” leem”,”em_pen”, “hadron-em-ln”, “hadron-em_liv-ln”, “hadron-leem-ln”, “hadron-em_pen-ln”, “hadron+em-ln”, “hadron-em+ln”, “hadron-em_liv+ln”, “hadron-leem+ln”, “hadron-em_pen+ln”, “hadron+em+ln”, “hadron+em_liv-ln”, “hadron+leem-ln”, “hadron+em_liv+ln”, “hadron+leem+ln”, “hadron+em_pen-ln”, “hadron+em_pen+ln”, “bertini”,”binary”, “QBBC”, “radioactivedecay”	Sets the physics models used based on scenarioType. This is a mandatory command required before starting the /run/beamOn commands
/phs/showscen		Displays the current scenario
/phys/verboseLevel <level>	G4int level	Sets the verbosity level for output of physical interactions.
/phys/cuts/global/default <cutValue> <unit>	G4double cutValue G4String unit = “um”, “mm”, “cm”, or “m”	Sets the global cut-length for all secondary particle production to be cutValue in units of unit.
/phys/cuts/global/gamma <cutValue> <unit>	G4double cutValue G4String unit = “um”, “mm”, “cm”, or “m”	Sets the global cut-length for secondary gamma-ray production to be cutValue in units of unit
/phys/cuts/global/electron <cutValue> <unit>	G4double cutValue G4String unit = “um”, “mm”, “cm”, or “m”	Sets the global cut-length for secondary electron production to be cutValue in units of unit.
/phys/cuts/global/positron <cutValue> <unit>	G4double cutValue G4String unit = “um”, “mm”, “cm”, or “m”	Sets the global cut-length for secondary positron production to be cutValue in units of unit.
/phys/cuts/region/setcut  <regionName> <cutValue> <unit>	G4String regionName G4double cutValue G4String unit = “um”, “mm”, “cm”, or “m”	Sets the cut-length for secondary particle production to be cutValue in units of unit for region regionName.
/phys/cuts/region/cutperparticle <regName> <parName> <cutValue> <unit>	G4String regName G4String parName G4double cutValue G4String unit = “um”, “mm”, “cm”, “m”	Sets the cut-length for secondary parName particle production to be cutValue in units of unit for region regionName
/phys/region/add <regName>	G4String regName	Add a new region named regName.
/phys/region/delete <regName>	G4String regName	Delete a region named regName.
/phys/region/list		List the regions created.
/phys/region/addalayer  <regName> <layer>	G4String regName G4int layer	Add the layer layer to region regName.
/phys/region/delalayer  <regName> <layer>	G4String regName G4int layer	Remove the layer layer from region regName.
/phys/region/listlayers <regName>	G4String regName	List the layers in region regName.

Analysis Commands

Analysis Commands

Command	Type	Description
/analysis/normalise <fluence> <unit>	G4double fluence G4String unit = “cm2” or “m2”	Sets the incident particle fluence to fluence in units of 1.0/unit.
/analysis/file <filename>	G4String filename	Sets the name of the CSV and RPT files.
/analysis/fluence/add <layerIndex>	G4int layerIndex>=0	Selects the layer (denoted by layerIndex ) immediately after which the particle fluence analysis is to be performed. Note that: * the fluence analysis is performed at the boundary between two layers; * layerIndex=1 indicates the surface of the geometry; * layerIndex=nLayers+1 indicates the surface of a planar geometry immediately after the last layer (i.e. the exit layer); * layerIndex=0 indicates that Fluence analysis is to be performed for all layers in the geometry.
/analysis/fluence/delete <layerIndex>	G4int layerIndex>=0	Deselects layer layerIndex for particle fluence analysis. layerIndex=0 indicates all layers are to be deselected.
/analysis/fluence/list		Identifies the layers for which particle fluence analysis is to be performed
/analysis/fluence/unit unit	G4String unit = “cm2” or “m2”	Sets the units of fluence used in the output.
/analysis/fluence/type type	G4String type = “omni” or “planar”	Selects whether the standard (omnidirectional) fluence or planar fluence (particle current) is measured. For a nonMC analysis, the “omni” fluence type must be selected.
/analysis/fluence/energy/mode <mode>	G4String <mode> = “LINEAR”, “LIN”, “linear”, “lin”, “LOGARITHMIC”, “LOG”, “logarithmic”, “log”, “ARBITRARY”, “ARB”, “arbitrary”, or “arb”	Sets the energy binning-scheme mode for particle fluence data.
/analysis/fluence/energy/default		Applies the default energy binning-scheme: 60 energy bins logarithmically-spaced over the range 1 keV to 1 GeV.
/analysis/fluence/energy/min <eMin> <unit>	G4double eMin> 0.0 && eMin < eMax; G4String unit = “eV”, “keV”, “MeV”, “GeV”, “TeV”, or “PeV”	Sets the lower-edge of the lowest bin in the energy binning-scheme.
/analysis/fluence/energy/max <eMax> <unit>	G4double eMax> eMin; G4String unit = “eV”, “keV”, “MeV”, “GeV”, “TeV”, or “PeV”	Sets the upper-edge of the highest bin in the energy binning-scheme.
/analysis/fluence/energy/nbin <Nebin>	G4int Nebin>0	Sets the number of bins in the energy binning-scheme. Nebin is forced to unity if eMax``==``eMin This command is not available for arbitrary binning schemes.
/analysis/fluence/energy/list		Lists the current energy binning-scheme.
/analysis/fluence/energy/add <eEdge> <unit>	G4double eEdge; eMin<eEdge<eMax G4String unit = “eV”, “keV”, “MeV”, “GeV”, “TeV”, or “PeV”	Adds a bin-edge. Only available for arbitrary binning-schemes.
/analysis/fluence/energy/delete <eEdge> <unit>	G4double eEdge G4String unit = “eV”, “keV”, “MeV”, “GeV”, “TeV”, or “PeV”	Deletes an existing bin-edge. Only available for arbitrary binning schemes.
/analysis/fluence/energy/clear		Delete all bin-edges between eMin and eMax (Nebin set to 1).
/analysis/fluence/angle/mode <mode>	G4String mode = “LINEAR”, “LIN”, “linear”, “lin”, “ARBITRARY”, “ARB”, “arbitrary”, or “arb”	Sets the angle binning-scheme mode for particle fluence data.
/analysis/fluence/angle/default		Applies the default angle binning-scheme: 2 angle bins linearly-spaced over the range 0\(^\circ\) to 180\(^\circ\).
/analysis/fluence/angle/min <tMin> <unit>	G4double tMin > 0 && tMin < tMax G4String unit = “deg” or “rad”	Sets the lower-edge of the lowest bin in the angle binning-scheme.
/analysis/fluence/angle/max <tMax> <unit>	G4double tMax > tMin G4String unit = “deg” or “rad”	Sets the upper-edge of the highest bin in the angle binning-scheme.
/analysis/fluence/angle/nbin <Ntbin>	G4int Ntbin>0	Sets the number of bins in the angle binning-scheme. Ntbin is forced to unity if tMax==tMin. This command is not available for arbitrary binning schemes.
/analysis/fluence/angle/list		Lists the current angle binning-scheme
/analysis/fluence/angle/add <tEdge> <unit>	G4double tEdge; tMin< tEdge<tMax; G4String unit = “deg” or “rad”	Adds a bin-edge. Only available for arbitrary binning-schemes.
/analysis/fluence/angle/delete <tEdge> <unit>	G4double aEdge G4String unit = “deg” or “rad”	Deletes an existing bin-edge. Only available for arbitrary binning schemes.
/analysis/fluence/angle/clear		Delete all bin-edges between tMin and tMax (Ntbin set to 1).
/analysis/fluence/particle/add <particleName>	G4String particleName	Select particle type particleName for fluence analysis.
/analysis/fluence/particle/delete <particleName>	G4String particleName	Deselect particle type particleName for fluence analysis.
/analysis/fluence/particle/list		List particle types for which fluence analysis is to be performed.
/analysis/dose/add <layerIndex>	G4int layerIndex>=0	Selects the layer denoted by layerIndex for which energy deposition / TID analysis is to be performed. layerIndex=0 indicates that Dose analysis is to be performed for all layers in the geometry.
/analysis/dose/delete <layerIndex>	G4int layerIndex>=0	Deselects layer layerIndex for Edep/TID analysis. layerIndex=0 indicates all layers are to be deselected.
/analysis/dose/list		Identifies the layers for which Edep/TID analysis is to be performed.
/analysis/dose/unit <unit>	G4String unit = “eV”, “keV”, “MeV”, “GeV”, “TeV”, “PeV”, “rad”, or “Gy”	Selects the units by which dose or energy deposition in the layers shall be measured.
/analysis/doseEq/add <layerIndex>	G4int layerIndex>=0	Selects the layer denoted by layerIndex for which Dose Equivalent analysis is to be performed. layerIndex=0 indicates that Dose Equivalent analysis is to be performed for all layers in the geometry.
/analysis/doseEq/delete <layerIndex>	G4int layerIndex>=0	Deselects layer layerIndex for Dose Equivalent analysis. layerIndex=0 indicates all layers are to be deselected.
/analysis/doseEq/list		Identifies the layers for which Dose Equivalent analysis is to be performed.
/analysis/doseEq/unit <unit>	G4String unit = “Sv” or “rem”	Selects the units by which Dose Equivalent in the layers shall be measured.
/analysis/phs/add <layerIndex>	G4int layerIndex>=0	Selects the layer denoted by layerIndex for which PHS analysis must be performed. layerIndex=0 indicates that PHS analysis is to be performed for all layers in the geometry.
/analysis/phs/delete <layerIndex>	G4int layerIndex>=0	Deselects layer layerIndex for PHS analysis. layerIndex=0 indicates all layers are to be deselected.
/analysis/phs/list		Identifies the layers for which PHS analysis is to be performed.
/analysis/phs/energy/mode <mode>	G4String <mode> = “LINEAR”, “LIN”, “linear”, “lin”, “LOGARITHMIC”, “LOG”, “logarithmic”, “log”, “ARBITRARY”, “ARB”, “arbitrary”, or “arb”	Sets the energy binning-scheme mode for PHS data.
/analysis/phs/energy/default		Applies the default energy binning-scheme for PHS.
/analysis/phs/energy/min <eMin> <unit>	G4double eMin > 0 && eMin<eMax; G4String unit = “eV”, “keV”, “MeV”, “GeV”, “TeV”, or “PeV”	Sets the lower-edge of the lowest bin in the energy binning-scheme for PHS.
/analysis/phs/energy/max <eMax> <unit>	G4double eMax > eMin G4String unit = “eV”, “keV”, “MeV”, “GeV”, “TeV”, or “PeV”	Sets the upper-edge of the highest bin in the energy binning-scheme for PHS.
/analysis/phs/energy/nbin <NIntv>	G4int NIntv >0	Sets the number of bins in the energy binning-scheme. NIntv is forced to unity if eMax==eMin. This command is not available for arbitrary binning schemes.
/analysis/phs/energy/list		Lists the current energy binning-scheme for PHS
/analysis/phs/energy/add <eEdge>	G4double eEdge; eMin<eEdge<eMax G4String unit = “eV”, “keV”, “MeV”, “GeV”, “TeV”, or “PeV”	Adds a bin-edge. Only available for arbitrary binning-schemes
/analysis/phs/energy/delete <eEdge>	G4double eEdge G4String unit = “eV”, “keV”, “MeV”, “GeV”, “TeV”, or “PeV”	Deletes an existing bin-edge. Only available for arbitrary binning schemes
/analysis/phs/energy/clear		Delete all bin-edges between eMin and eMax (NIntv set to 1).
/analysis/niel/add <layerIndex>	G4int layerIndex>=0	Selects the layer denoted by layerIndex for which NIEL analysis must be performed. layerIndex=0 indicates that NIEL analysis is to be performed for all layers in the geometry.
/analysis/niel/delete <layerIndex>	G4int layerIndex>=0	Deselects layer layerIndex for NIEL analysis. layerIndex=0 indicates all layers are to be deselected.
/analysis/niel/list		Identifies the layers for which NIEL analysis is to be performed.
/analysis/niel/unit	G4String unit = “rad”, “Gy”, or “MeV/g”	Selects the units by which NIEL deposition at the layer boundaries shall be measured.
/analysis/niel/function <NIELCoeffs>	G4String NIELCoeffs = “cern”, “CERN”, “jpl”, “JPL”, “JPL_NRL_NASA_2003_Si”, “JPL_NRL_NASA_2003_GaAs”, “JPL_NRL_NASA_2003_InP”, “INFN_ESA_2014_10eV_GaAs”, “INFN_ESA_2014_21eV_GaAs”, “INFN_ESA_2014_25eV_GaAs”, “global”	Selects the set of NIEL coefficients denoted by NIELCoeffs for use with particle fluence data.
/analysis/niel/functionByLayer <layerIndex> <NIEL Coeffs>	G4int layerIndex >= 0; G4String NIELCoeffs = “cern”, “CERN”, “jpl”, “JPL”, “JPL_NRL_NASA_2003_Si”, “JPL_NRL_NASA_2003_GaAs”, “JPL_NRL_NASA_2003_InP”, “INFN_ESA_2014_10eV_GaAs”, “INFN_ESA_2014_21eV_GaAs”, “INFN_ESA_2014_25eV_GaAs”	Selects the set of NIEL coefficients denoted by NIELCoeffs for use with particle fluence data for a specific Layer.
/analysis/niel/references		Lists the published references from which the NIEL data was obtained.
/analysis/solarCell/add	G4int layerIndex>=0	Add a layer to the solar cell analysis. layerIndex=0 indicates that solar cell analysis is to be performed for all layers in the geometry.
/analysis/solarCell/delete	G4int layerIndex>=0	Delete a layer from the solar cell analysis. layerIndex=0 indicates all layers are to be deselected.
/analysis/solarCell/list		List the solar cell analysis.
/analysis/solarCell/unit	G4String unit = “rad”, “Gy”, or “MeV/g”	Sets the units to report the TNID.
/analysis/solarCell/references		Lists the NIEL curve references.
/analysis/solarCell/function <NIELCoeffs>	G4String NIELCoeffs = “cern”, “CERN”, “jpl”, “JPL”, “JPL_NRL_NASA_2003_Si”, “JPL_NRL_NASA_2003_GaAs”, “JPL_NRL_NASA_2003_InP”, “INFN_ESA_2014_10eV_GaAs”, “INFN_ESA_2014_21eV_GaAs”, “INFN_ESA_2014_25eV_GaAs”, “global”	Sets the global NIEL curve to use for the solar cell analysis.
/analysis/solarCell/functionByLayer <layerIndex> <NIELCoeffs>	G4int layerIndex >= 0; G4String NIELCoeffs = “cern”, “CERN”, “jpl”, “JPL”, “JPL_NRL_NASA_2003_Si”, “JPL_NRL_NASA_2003_GaAs”, “JPL_NRL_NASA_2003_InP”, “INFN_ESA_2014_10eV_GaAs”, “INFN_ESA_2014_21eV_GaAs”, “INFN_ESA_2014_25eV_GaAs”	Sets the NIEL function for a specific layer of the geometry.
/analysis/solarCell/setParameters <tag> <n> <Cp> <Ce> <Dxp> <Dxe>	G4String tag G4double n G4double Cp G4double Ce G4double Dxp G4double Dxe	Sets the experimentally determined parameters that characterise the solar cell degradation as specified by the SCREAM/SAVANT method of NRL.
/analysis/solarCell/ListParameters		Lists the previously entered parameters
/analysis/solarCell/delParameters <tag>	G4String tag	Deletes the set of parameters corresponding to the tag.

Progress Control Commands

Progress Control Commands

Command	Type	Description
/run/progressControl/filePrefix <name>	G4String name	Sets the output filename prefix. Either “_commout.json” or “_commout.txt” will be appended to this value to form the filename. Default prefix is the name of the application, i.e., “MULASSIS”
/run/progressControl/fileType <type>	G4String type = “json”, “text”, or “txt”	Output file format. See Note below regarding case-sensitivity.
/run/progressControl/commFreq <f> <unit>	G4int f>0 G4String unit = “event”, or “batch”	Defines the intervals at which the output progress file is updated. See Note below regarding case-sensitivity. Default is update every batch.
/run/progressControl/targetEvents <i1> <i2> <i3> ... <iN>	G4int i1, i2, i3, ..., iN > 0	Sets the maximum number of events for each of the subsequent /run/beamOn commands. If not set, it is assumed there is only one run, and the number of events identified with the /run/beamOn command is used as the target value. This command is ignored for non-MC runs.
/run/progressControl/targetRuns <runs>	G4int runs>0	Sets the number of non-MC runs to perform. This command is ignored for MC runs.
/run/progressControl/resetRunNumber		Allows synchronisation of the run number when using the targetEvents command with multiple run targets or the targetRuns command.

Note

String inputs which are all upper case, all lower case, or initial letter upper case are treated the same, e.g. “event”, “EVENT”, and “Event” are equivalent.

Miscellaneous Commands

Misc. Commands

Command	Type	Description
/control/execute <macrofile>	G4String macrofile	Selection of additional macro files to run.
/event/printModulo <nModulo>	G4int nModulo	Sets the interval (in number of events nModulo) after which the CSV file is updated.
/run/randseed <iSeed>	G4int Seed >=0 && <=214	Sets the random number seed to be used for the simulation.
/run/cputime <CPUtime>	G4double CPUtime>0.0	Sets the time-limit (in CPU seconds) for program execution.
/run/beamOn <nEvents>	G4int nEvents>0	Initiates the execution of the simulation for nEvents primary particles.
/nonMC/beamOn		Perform nonMC shielding analysis, equivalent to /run/beamOn for MC simulations.
/nonMC/error/setGlobalPercent <errorPC>	G4double errorPC >=0.0 && <=1.0E+10	Sets the percentage error to be applied to the nonMC results. Default is 0.0
/nonMC/error/getGlobalPercent		Displays the percentage error to apply to the results.

NonMC Integration Commands

NonMC Integration Commands

Command	Type	Description
/nonMC/integration/setMaxTIntv <maxTIntv> <unit>	G4double maxTIntv >=0.01 deg && <=5.0 deg G4String unit = “deg” or “rad”	Defines the maximum interval in \(\theta\) for integration of incident (cosine-law) flux. Default is 0.2\(^\circ\).
/nonMC/integration/setMaxTLimit <maxTLimit> <unit>	G4double maxTIntv >=85.0 deg && <= 89.999 deg G4String unit = “deg” or “rad”	Defines the maximum angle to the normal up to which the incident (cosine-law) flux will be integrated. Default is 89.9\(^\circ\).
/nonMC/integration/setNMinTIntv <minTIntv>	G4int minTInv >=1 && <=100	Set the minimum number of angle intervals for Simpson 3/8 integration within each fluence bin. Default is 20.
/nonMC/integration/minEpn <minEpn> <unit>	G4double minEpn >=0.01 MeV && <=10.0 GeV G4String unit = “eV”, “keV”, “MeV”, or “GeV”	Defines the minimum energy per nucleon from which the numerical integration of the incident spectrum will be performed. Default is 0.01 MeV.
/nonMC/integration/nEpd <nEpd>	G4int nEpd >=1 && <=100	Number of energies per decade over which the spectra are defined. Default is 20.
/nonMC/integration/nEpts <nEpts>	G4int nEpts >=22	Total number of energies to be used for numerical integration of the incident spectrum. Default is 140, therefore maximum energy is 10 GeV/nuc, if minEpd and nEpd are used with default values.
/nonMC/integration/nEint <nEint>	G4int nEint >=4 && <=100	The attenuation cross-section is integrated over for each material and over each of the nEpd intervals per decade in energy using a Simpson 3/8 integration. nEint defines the number of intervals over which the energy integration is performed. Default is 22.
/nonMC/integration/rombergTolerance <rombergTol>	G4double rombergTol >=1.0E-10 && <=1.0E-3	Changes the tolerance for the Romberg numerical integration. If there are thin layers, this value should be low, e.g., ~1.0E-7 if thin ~10 microns Si layers are used. At ~1mm thicknesses, 1.0E-4 appears to be sufficient, but lower values can be used if preferred without significantly slowing the numerical integration. Default is 1.0E-7.
/nonMC/integration/listParameters		Lists all the current parameter values for integration process.
/nonMC/doseModelThreshold/setThreshold <dDOverDLimit>	G4double dDOverDLimit >=1.0E-2 && <= 1.0E-8	Threshold ratio for ionising energy loss in layer to incident energy on layer. Above this threshold, the ionising energy deposition per unit area is determined from the difference in the incident and exiting particle energy flux. Below this threshold an approximate model (based on integration of flux, stopping power and mean pathlength) is used to calculate the energy deposition; this avoids problems due to rounding errors in the above method. Default is 1.0E-6.
/nomMC/doseModelThreshold/getThreshold		Print the threshold ratio for ionising energy loss in layer to incident energy on layer.

Warning

These parameters should not be changed unless the user has a clear understanding of the consequences on the numerical integration processes.

NonMC NIEL Commands

NonMC NIEL Commands

Command	Type	Description
/nonMC/niel/name <name>	G4String name	Allows user to define a name/designation for the new NIEL coefficients
/nonMC/niel/datapoint <userEnergy> <userNIEL>	G4double userEnergy>0.0 G4double userNIEL>=0.0	Defines a single datapoint for the NIEL function userEnergy must be in MeV and userNIEL in MeV cm2/g
/nonMC/niel/clear		Clears all NIEL curve data-points already entered for the user-defined NIEL curve.
/nonMC/niel/list		Lists the table of data-points for the user-defined NIEL curve.

Warning

A nonMC analysis cannot be executed if one or more of the following conditions exist:

• The incident particle is not an ion.

• The angular distribution for the incident particles is not planar and is not cosine-law.

• The incident particles have no component which is from the +Z direction (i.e. they are incident upon only the rear of the shield not from the front).

• If ion fluence analysis is performed and the fluence analysis type is PLANAR rather than OMNIDIRECTIONAL.

Note also that the nonMC results for a spherical shield are relevant only to the calculated dose, TNID and fluence at the centre of the shield, not for the intervening layers.

---

Mulassis/ml-v02-00/r342