OPTIONS PROGRESS-PRINT GLOBAL anode = -4100 GLOBAL wdiameter = 0.002 // 20um GLOBAL wypos = 0.0 GLOBAL wpot = 0.0 GLOBAL wtension = 70.0 // 70g GLOBAL wlength = 25.0*1.141 // 45grad Global gas_file `Arg-63-Eth-37.gas` Global gas_member `a63e37` //============================================= & CELL //============================================= opt cell-pr plane y=-1.3 v={anode} plane y=+1.3 v={anode} // we use cartesian system ROWS CARTESIAN S 41 {wdiameter} -8.46+(0.423*I) ... {wypos} ... {wpot} ... {wtension} ... {wlength} ... TUNGSTEN GRAVITY 0 0 1 WRITE DATASET vdc_cell.dat //============================================= & MAGNETIC //============================================= // at present: no B fieed here // whats about the earth B field ?! COMPONENTS 0 0 0 GAUSS //============================================= & GAS //============================================= // 1 atm 63/37 Argon-ethane mixture in your chamber. Call inquire_member(gas_file,gas_member,`gas`,exist) If exist Then get {gas_file,gas_member} Else write {gas_file,gas_member} // Loads the mixture Argon 63 %, ethane 37 % ARG-63-ETH-37 // temp + pressure def before Heed TEMPERATURE 300 KELVIN PRESSURE {1*1013.25} MILLI-BAR // Compute for electrons: drift velocity, diffusion coefficient etc // time comsuming: several hours on 2.4GHz Xeon MAGBOLTZ argon 63 ethane 37 // (Tong put it outside if-else-endif ?!) // Prepares cluster generation by Heed Heed argon 63 ethane 37 Endif !add meta type PostScript file-name "vdc_ar63e37.ps" !open meta !act meta plot-options drift-velocity nodiffusion notownsend opt gas-plot nogas-print //============================================= & OPTIMISE //============================================= //Prints the current potential settings DISPLAY //============================================= & FIELD //============================================= //============================================= & DRIFT //============================================= //Sets the size and view of the drift area AREA -1.5 -1.5 +1.5 +1.5 //Number of drift lines used by x(t) etc. LINES 60 //Grid density for tables and contour plots GRID 10 20 // define area around wire #21 for possible plots SELECT 21 // select center wire INTEGRATION-PARAMETERS COMPUTE-IF-INTERPOLATION-FAILS // define single track: 45deg TRACK -1.5 -1.5 1.5 1.5 // Heed takes care of cluster generation, it simulates the ionisation of // the gas molecules by a particle. TRACK HEED ... NODELTA-ELECTRONS NOTRACE-DELTA-ELECTRONS ... NOMULTIPLE-SCATTERING ... NOENERGY-CUT ELECTRON ENERGY 0.85 GeV // TIMING: Computes the arrival time distribution of the n'th electron over // a given area. Unlike the ARRIVAL-TIME-DISTRIBUTION instruction, TIMING does // not produce calibration curves, but merely timing distributions. // TIMING ELECTRON 1 ... TIME-WINDOW 0.0 0.3 ... X-RANGE -1.45 0.0 ... Y-RANGE 0.0 1.5 ... ANGLE-RANGE 53.787 55.673 ... BINS 301 ... ITERATIONS 10000 ... RUNGE-KUTTA-DRIFT ... PLOT-OVERALL PLOT-SELECTED-ELECTRON ... PRINT-OVERALL PRINT-SELECTED-ELECTRON !deact meta !close meta !del meta