TITLE Ordinary ConCreTe(CCT)thicknes for CANDLE Shielding * Body data: * Test of Electron Propagation in Concrete * 3 GeV e- beam induced radiation dose rate in Ordinary Concrete *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 COMMENT 3.0 This is the comment 1st line, *and this is the 2nd one, *and here is the 3rd. COMMENT 0.0 9.0 MODEL Linux RedHat 9, PC Pentium III 1000Hz COMMENT 0.0 INPUT CCT DEFAULTS NEW-DEFA * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 BEAM 3.0 ELECTRON * * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 BEAMPOS 0.0 0.0 -3.0E+00 * * * * Here starts the geometry descripton (see 9} Combinatorial Geometry): * ----------------------- * a large cylinder of concrete delimited by vacuum layers at the top * and bottom. Here depth is fixed at 3m and the radius 30cm. * Electrons are injected in the initial vacuum. * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * geometry the Comb.Geom * errors are accuracy is used * printed parameter GEOBEGIN 0.0 1.0E-01 COMBINAT * * Geometry title (in special format): * 3 variables are input in CG Title card:IVOPT(Input Value OPTion),IDBG,TITLE * according (2I5,10X,A60) * * 3 10 Ordinary ConCreTe(CCT) * * Body data: ******************************************************************************* *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * SPHere with coordinates of the center and R radius: SPH 001 0.0 0.0 0.0 +5000.0 * code #/# V_x V_y V_z Radius * ******************************************************************************* * Rectangular Parallelepiped. Code: RPP * An RPP has its edges parallel to the coordinate axes. * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 RPP 002 -1000.0 +1000.0 -1000.0 +1000.0 -100.0 +1000.0 * ******************************************************************************* * Infinite Circular Cylinder parallel to the axis z. *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * code #/# Z_Axis_x Z_Axis_y R_Radius * * In our case we are interested in cascade longitudinal and lateral development * for the 300cm Cylinder of radii from 0.0cm to 30cm, * with axis coincide with the z axis. * An infinite cylinder of radius 30cm, with axis equal to the z axis * could be established as: * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * ZCC 003 0.0 0.0 30.0 * ****************************************************************************** * Infinite half-space delimited by a coordinate plane. Code: XYP. * Delimited by a plane perpendicular to the z-axis. * This half-spaces is defined by a single number: V_z, * the coordinate of the plane on the corresponding perpendicular axis. * The half-space "inside the body" is the locus of points for which z < V_z. * XYP 004 0.0 * ******************************************************************************* *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * XYP 005 300.0 * * ****************************************************************************** * * END * * * N.B. * Blackhole or external vacuum Fluka # is 1 * Vacuum or internal vacuum Fluka # is 2 * reg# body# 9(A2,I5)-alternate as many 2-character fields and integer fields * (see above) * * external vacuum at the beginning *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 1 5 +1 -2 * * Internal vacuum at the beginning created as a buffer vacuum * region containing the whole geometry, which must itself be * contained within the external (mandatory) blackhole region. * It is then suggested that the beam impact point be chosen * in vacuum, slightly upstream the actual one on a * material boundary. *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 2 5 +2 -3 -4 -5 * * Internal vacuum in which the region of starting point of particle * is established by moving the particle of z = -3.0E+00 cm long step: * it is a Vacuum Cylinder with coordinates r = 30 cm and -100 <= z < 0.0 cm * before the Ordinary Concrete layer. *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 3 5 +2 +3 +4 +5 * * ordinary concrete layer: a Cylinder with r = 30 cm, and 0.0 <= z < 300 cm *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 4 5 +2 +3 -4 5 * * * internal vacuum at the end: a Rectangular Parallelepiped with * -1000 <= x <= 1000 cm, -1000 <= y <= 1000 cm, 300 <= z <= 1000 cm *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 5 5 +2 +3 -4 -5 * * The END card (in special format): * Region data must be terminated by a line with the string END. * END * * * N.B. * Region Volumes have to be inserted, inasmuch as we have defined: IVOPT=3 !!! * * * * Here is the stage of booking of the volume array VNOR(I=1,5) with the * format (7E10.5): * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 12345678 * 3.9250+11 4.4000+09 2.8260+04 8.4780+05 2.8000+09 * * * {GEOEND} card: ends the geometry definitions. * One should try the geometry debugger as the first thing to do when starting a * new job, even before geometry problems show up. * *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+... GEOEND 15. 15. 200. 0. 0. -3.0DEBUG GEOEND 10. 10. 30. 0. 0. 0. & *STOP *GEOEND * * {MATERIAL} card: defines a single-element material or a compound * (coupled to a COMPOUND card) * Here defines Materials and Compounds of Ordinary ConCreTe(CCT) ** atomic #Z at.wt #A g/cm^3 material # not used mass # name *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * * * 1) Hydrogen (H) MATERIAL 1.0 1.0079 8.988E-5 3.0 0.0 1. HYDROGEN * 2) Carbon (C) MATERIAL 6.0 12.011 2.265 6.0 0. CARBON * 3) Oxygen (O) MATERIAL 8.0 15.9994 0.001429 8.0 0. OXYGEN * 4) Magnesium (MG) MATERIAL 12.0 24.305 1.738 9.0 0. MAGNESIU * 5) Aluminum (AL) MATERIAL 13.0 26.982 2.70 10.0 0. ALUMINUM * 6) Iron (FE) MATERIAL 26.0 55.847 7.877 11.0 0. IRON * 7) Calcium (CA) MATERIAL 20.0 40.08 1.54 12.0 0. CALCIUM * 8) Sodium (NA) MATERIAL 11.0 22.9898 0.975 13.0 0. SODIUM * 9) Silicon (SI) MATERIAL 14.0 28.0855 2.33 26.0 0. SILICON * 10) Potassium (K) MATERIAL 19.0 39.102 0.031165 27.0 0. POTASSIU * * 11) For Ordinary ConCreTe compound using atom relative contents * see: B.T. Price, C.C. Horton & K.T. Spinny * "Radiation Shielding", Pergamon Press, London-New York-Paris, 1957 * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * *MATERIAL 0.0 0.0 2.3 28.0 0. CCTORDIN * * {COMPOUND} card: * Defines a compound or mixture, made of several materials, * or even a mixture of different isotops * * * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * -mass index of -mass index of etc. etc. name * fraction the 1st fraction the 2nd of the * of the 1st material of the 2nd material compound * material material * * *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+... MATERIAL 0.0 0.0 2.30 28.0 0.0 0. CCTORDIN COMPOUND -1.0 3.0 -0.1 6.0 -52.9 8. CCTORDIN COMPOUND -0.2 9.0 -3.4 10.0 -1.4 11. CCTORDIN COMPOUND -4.4 12.0 -1.6 13.0 -33.7 26. CCTORDIN COMPOUND -1.3 27.0 CCTORDIN * * {ASSIGNMAt} card: defines the correspondence between region indices and * material indices. It defines also regions with magnetic field. * * Here assigns materials to regions * Ordinary ConCreaT in region 4 * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * mater.# lower * bound of * the region * indices * External Black Hole *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 ASSIGNMA 1.0 1.0 * Internal Vacuum *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 ASSIGNMA 2.0 2.0 * * Internal Vacuum bounded with Z<0.0cm half-space *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 ASSIGNMA 2.0 3.0 * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 ASSIGNMA 28.0 4.0 4.0 1.0 -1.0 * * CCTORDIN from region WHAT(2) to region WHAT(3) in steps of WHAT(4); * WHAT(5)< 0.0 resets the default (no magnetic field) * in the region(s) defined by WHAT(2),(3), and (4); * WHAT(6), SDUM : not used. * * Internal Vacuum *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 ASSIGNMA 2.0 5.0 * * {STEPSIZE} card: sets the minimum and maximum step size on a region-by-region * basis for transport of all charged particles * (hadrons, muons & electrons). *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * min max lowerbound: upper_ in steps * step size stp size from bound: to of 1 cm * in cm in cm region region (default) * * In this case we can establish: *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 STEPSIZE 1.0E-01 2.0E-01 2.0 5.0 1.5E-01 * * * * Setting SCORE (shahatz miavorneri hashiv-aghyusak) and Output Levels. * {SCORE} card: defines the (generalized) particles producing the stars * to be scored by region * Requests scoring of energy deposition (normalized) in each region. * * * * {THRESHOL}d: Defines the energy threshold for star density scoring; * *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+... THRESHOLd 0.0 0.0 0.0 0.05 0.0 0.0 * Score stars only above 50 MeV * * * SCORE 208.0 211.0 * * * {OUTLEVEL} card: defines the amount of output wanted. * Output level was originally planned to range * from 1.0 (minimum output) to 7.0 (maximum output). * min max OUTLEVEL 1.0 7.0 * * * **** Scoring **** * * {RESNUCLE}i card: scors residual nuclei produced in inelastic interactions * on a region basis: *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * all logical max Z of maximum M= scoring volume of * products output rezidual =N-Z-NMZmin region the region * unit nuclei of rez.nucl. number in cm**3 RESNUCLE 3.0 22.0 3.0 1.0 Dummy * * * * {USRTRACK} card: defines a detector for a track-length fluence estimator * *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * logarithm.gneralized logical region volume number ID of the * binning particle output defining of the of track- * in energy type to unit the detector energy -length * be scored detector in cm**3 bins detector * Continuation card: * maximum minimum ... not used...................... & in any * energy in energy in position * GeV? for GeV? for in colomn * scoring scoring 71 to 78 * USRTRACK -1.0 208.0 47.0 3.0 1.0 Trhadron USRTRACK 3. 0.001 & USRTRACK -1.0 211.0 47.0 3.0 1.0 Tre-e+ph USRTRACK 3. 0.001 * * * {USRBDX} card: defines a detector for a boundary crossing fluence or current * estimator: *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * i1+i2*10+i3*100 generlized logical first second area ID of * i1=logar. particle output region region of the detector * binning type unit: defining defining detector * in energy; to be < 0.0; the the in cm**2 * i2=one way scored unformated boundary boundary * scoring; data * i3=current are * scoring written *23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789 * Continuation card: * maximum minimum number maximum minimum number of & in any * kinetic kinetic of energy solid solid angular position * energy,GeV energy,GeV intervals angle angle bins in colomn * for for for for for if linear 71 to 78 * scoring scoring scoring scoring scoring, angular * in sr if linear bin: =1 ********!!!!!!!! USRBDX -1.0 208.0 48.0 3.0 4.0 1.0 Electrons USRBDX 3. 0.001 1. & USRBDX -1.0 211.0 48.0 3.0 4.0 1.0 Electrons USRBDX 3. 0.001 1. & *******!!!!!!!! * * * {USERDUMP} card: defines a collision "tape" (or better a collision file) * to be written; activates MGDRAW. The option 2 in WHAT(3) * activates dumping of trajectories and continuous energy losses: * this automatically activates also the BXDRAW entry of MGDRAW. * * *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 USERDUMP 101. 2.0 CCTORDIN * * Activate Random Number initialization * * RANDOMIZe sets the seeds for the double-precision random number * generator RM48 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 RANDOMIZ 1. * * USERICALL Activates user dependent initialization * USRICALL * * Starts Job asking for a given number of events * * START defines the termination conditions, gets a primary from a beam or * from a source and starts the transport *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 * Max# to be not time for not reportsthe time for not * Simulated used trmination used #of calls intractive used * 0.0 or>0.0 run START 100. 99999999. 0.0 0.0 * * USROCALL Activates user dependent output * USROCALL STOP