CERN home page CERN home page Consult pages Consult pages Writeups at CERN Writeups at CERN Garfield pages Garfield pages Help Help &GAS &GAS MAGBOLTZ MAGBOLTZ MONTE-CARLO-INTEGRATION MONTE-CARLO-INTEGRATION Search Keyword index

&GAS: MAGBOLTZ: MONTE-CARLO-INTEGRATION


COLLISIONS

Specifies the number of collisions, in multiplies of 960000, to be used to compute the transport properties for each (E,B,angle) combination. This option can only be used with MONTE-CARLO-INTEGRATION integration.

Low values (e.g. 2-5) are suitable for mixtures than contain substantial amounts of quenchers. Nearly pure noble gases call for larger values (e.g. 10-20) while very large values (100 and more) are needed for pure noble gases. But in the latter case, it is preferable to use Magboltz 1 (see ANALYTIC-INTEGRATION).

The statistical accuracy of the drift velocity calculation improves with the square root of this parameter, while the CPU time consumption increases linearly. Since the accuracy depends on the gas mixture and on the field, it may be worthwhile in critical applications to estimate the statistical error of the Magboltz calculations for a given gas and a given field by repeatedly calculating the transport parameters:

Call "-> Call -> book_histogram"(vhist,50)
Call "-> Call -> book_histogram"(dlhist,50)
Call "-> Call -> book_histogram"(dthist,50)
Global nrndm=200
For irndm From 1 To nrndm Do
   Say "Iteration {irndm}/{nrndm}"
   &GAS
   magboltz argon 90 methane 10 coll 5 e-field 125

&MAIN Call "-> Call -> drift_velocity"(125,0,0,v) Call "-> Call -> fill_histogram"(vhist,v) Call "-> Call -> transverse_diffusion"(125,0,0,dt) Call "-> Call -> fill_histogram"(dthist,10000*dt) Call "-> Call -> longitudinal_diffusion"(125,0,0,dl) Call "-> Call -> fill_histogram"(dlhist,10000*dl) Enddo

Call "-> Call -> plot_histogram"(vhist,`v<SUB>D</SUB> [cm/microsec]`,`Drift velocity`) Call "-> Call -> plot_end" Call "-> Call -> plot_histogram"(dthist,`&sigma;<SUB>T</SUB> [micron for 1 cm]`, ... `Transverse diffusion`) Call "-> Call -> plot_end" Call "-> Call -> plot_histogram"(dlhist,`&sigma;<SUB>L</SUB> [micron for 1 cm]`, ... `Longitudinal diffusion`) Call "-> Call -> plot_end"

From the histograms, one concludes that, for Argon 90 % methane 10 % and E=125 V/cm, COLLISION=5 leads to a statistical error of 0.4 % on the drift velocity, 4.5 % on the transverse diffusion and 7 % on the longitudinal diffusion.

[By default: 10 x 960.000 collisions.]


Go to the top level, to &GAS, to MAGBOLTZ, to MONTE-CARLO-INTEGRATION, to the topic index, to the table of contents, or to the full text.

Formatted on 0105-10-15 at 17:17.