#define BIASFLAG TITLE STEND - check target geo and beam position * * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+... * * =================================================== * Set the defaults for precision simulations * =================================================== DEFAULTS NEW-DEFA *DEFAULTS PRECISIO * =================================================== * Define the beam: * E= 50 GeV, proton, from right to left, SigmaE= 0.01 * FWHM = 2.355 sigma -> delE = 0.02355, * Angle-phi : sigma = 0 * X: gauss, sigma = 1 cm -> FWHM = 2.355 sigma = 2.355 * Y: gauss, sigma = 2 cm -> FWHM = 2.355 sigma = 4.710 * Z - ?... use SOURCE ? * start from point : Z= -20 (into positive), X=Y= 0 * * =================================================== BEAM -50. -0.01 0.0 -2.36 -4.71 PROTON BEAMPOS 0.0 0.0 -30. * * =================================================== * -------- geometry ---------- * =================================================== GEOBEGIN 0.0 COMBNAME 0 0 * * ====== BODY ========= * * Black bodies : vacuum and black * RPP BLCKBOX1 -300. 900. -200. 200. -200. 200. RPP BLCKBOX2 -400. 1000. -300. 300. -300. 300. * * main box with air: all elements are here ... * RPP bigbox -100. 800. -100. 100. -100. 100. * * Target : lead core and moderator * moderator width = 2 : (2.5 + 2) : job-1 * new core : 5 * 5 --> 7 * 7 --> 3.5 + width * L: 24, 30, 38, 46, 54 cm ==> L = 24 , mod= 2 cm * ------------- * * insert shell - to keep target as an unit * * RPP tcore -3.5 3.5 -3.5 3.5 -12. 12. RPP tmod -5.5 5.5 -5.5 5.5 -12. 12. RPP shell -6.5 6.5 -6.5 6.5 -12. 12. * * air-counters : beside target and hole * Xmax = -2.5 - width - 2 (air distance from target = 1 cm) * Xmin = -2.5 - width - 1 * Y = 2.5 + width * * 3 virtual counters : at X-axis * RPP count1 15. 16. -10. 10. -10. 10. RPP count2 99. 100. -10. 10. -10. 10. RPP count3 599. 600. -10. 10. -10. 10. * * to divide moderator ... * YZP Xzero 0.0 * END * * ====== REGIONs ======= * BLKBODY 5 +BLCKBOX2-BLCKBOX1 VOID 5 +BLCKBOX1-bigbox Rcore 5 +tcore RmodR 5 +tmod -tcore -Xzero RmodL 5 +tmod -tcore +Xzero Rshell 5 +shell -tmod * * Rbox - container with all elements * Rbox 5 +bigbox -shell -count1 -count2 - count3 * Rcount1 5 +count1 Rcount2 5 +count2 Rcount3 5 +count3 * END * * ======== end of geometry ========== *GEOEND GEOEND 50. 50. 50. -50. -50. -50.DEBUG GEOEND 100. 100. 100. & * * ======== new materials ============ * ---------------------------- * 104 Air dry (near sea level) * *MATERIAL .00120484 AIR *COMPOUND -1.248E-4 CARBON -.755267 NITROGEN -.231781 OXYGENAIR *COMPOUND -0.012827 ARGON AIR * * Polyethylene Marlex * Density variation of polyethylene is 0.91 - 1.05 g/cm3. "Low" density is * 0.920, "medium" is .93, and "high" is .95. Special polyethelene is made * for nuclear shielding, and this has loaded densities up to 1.08. * Chemical H H * Formula : | | * |C H | ---- C -- C ---- * | 2 4|n | | * H H *MATERIAL 0.93 Polyethy *COMPOUND 4. HYDROGEN 2. CARBON Polyethy * * Concrete portland * Concrete has a wide variation in density and composition. *MATERIAL 2.3 PORTLAND *COMPOUND -0.01 HYDROGEN -0.001 CARBON -.529107 OXYGENPORTLAND *COMPOUND -0.016 SODIUM -0.002 MAGNESIU -.033872 ALUMINUMPORTLAND *COMPOUND -.337021 SILICON -.013 POTASSIU -.044 CALCIUMPORTLAND *COMPOUND -0.014 IRON PORTLAND *MATERIAL 19. .862 POTASSIU *MATERIAL 53. 4.93 IODINE * * 252 Scintillator - NaI * *MATERIAL 3.667 Sodium_I *COMPOUND 1. SODIUM 1. IODINE Sodium_I * * =========================================================== * =========== set materials : use them from list ========== * ========================================================== * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7.. ASSIGNMA BLCKHOLE BLKBODY ASSIGNMA VACUUM VOID ASSIGNMA AIR Rbox ASSIGNMA POLYETHY RmodL ASSIGNMA POLYETHY RmodR ASSIGNMA AIR Rshell ASSIGNMA LEAD Rcore ASSIGNMA AIR Rcount1 ASSIGNMA AIR Rcount2 ASSIGNMA AIR Rcount3 * #if BIASFLAG * *BIASING 3. 10. Rcount2 Rcount3 * BIASING 0.0 1. BLKBODY @LASTREG PRINT * BIASING 3. 0.5 RmodL Rshell BIASING 3. 2. RmodR Rbox * #endif * * PLASTSCIN - Plastic Scin. (Vinyltoluence) - ICRU - 216 - (H-8.5, C-91.5) * *ASSIGNMA PLASCINT Rscin * -------------------------------- * set E cut for photon (default is high !) * *EMFCUT -1.E-6 1E-07 Rcore Rbox * * ====================================== * =========== SCORING ================ * * ***** To score yield of low energy neutrons from target * - 1 bin for slow neutrons * * --------------------------------------------------------- * ..... for neutrons.... Rtarg --> Rbox (from target) * --------------------------------------------------------- * USRYIELD 1399. NEUTRON -21. Rshell Rbox Ntag USRYIELD 55. 1E-14 50.3.14159265 0.0 3. & * * **** N time-dist. from Rtarg ... * USRYIELD 2567. NEUTRON -21. Rshell Rbox NTargTime USRYIELD 1. 1E-10 100. 55. 1E-14 3. & * ------------------------------------ * ..... for photons.... * USRYIELD 1399. PHOTON -22. Rshell Rbox PhTag USRYIELD 55. 1E-07 50.3.14159265 0.0 3. & * * USRYIELD 67. PHOTON -22. Rshell Rbox PhTagTi USRYIELD 1. 1E-10 50. 55. 1E-07 3. & * * --------------------------------------------------------- * ..... for neutrons.... Rbox --> Rcount1 (into counter-1) * --------------------------------------------------------- * USRYIELD 1399. NEUTRON -23. Rbox Rcount1 Nc1 USRYIELD 55. 1E-14 50.3.14159265 0.0 3. & * * **** N time-dist. target -> counter1 ... * USRYIELD 2567. NEUTRON -23. Rbox Rcount1 Nc1Time USRYIELD 1. 1E-10 100. 55. 1E-14 3. & * ------------------------------------ * ..... for photons : target -> counter1 .... * USRYIELD 1399. PHOTON -24. Rbox Rcount1 Phc1 USRYIELD 55. 1E-07 50.3.14159265 0.0 3. & * * USRYIELD 67. PHOTON -24. Rbox Rcount1 Phc1Time USRYIELD 1. 1E-10 100. 55. 1E-07 3. & * * -------------------------------------- * * N for C2 * USRYIELD 1399. NEUTRON -23. Rbox Rcount2 Nc2 USRYIELD 55. 1E-14 50.3.14159265 0.0 3. & * USRYIELD 2567. NEUTRON -23. Rbox Rcount2 Nc2Time USRYIELD 1. 1E-10 100. 55. 1E-14 3. & * * ..... for photons : C2 .... * USRYIELD 1399. PHOTON -24. Rbox Rcount2 Phc2 USRYIELD 55. 1E-07 50.3.14159265 0.0 3. & * * USRYIELD 67. PHOTON -24. Rbox Rcount2 Phc2Time USRYIELD 1. 1E-10 100. 55. 1E-07 3. & * * -------------------------------------- * * N for C3 * USRYIELD 1399. NEUTRON -23. Rbox Rcount3 Nc3 USRYIELD 55. 1E-14 50.3.14159265 0.0 3. & * USRYIELD 2567. NEUTRON -23. Rbox Rcount3 Nc3Time USRYIELD 1. 1E-10 100. 55. 1E-14 3. & * * ..... for photons : C3 .... * USRYIELD 1399. PHOTON -24. Rbox Rcount3 Phc3 USRYIELD 55. 1E-07 50.3.14159265 0.0 3. & * * USRYIELD 67. PHOTON -24. Rbox Rcount3 Phc3Time USRYIELD 1. 1E-10 100. 55. 1E-07 3. & * * ----------------------------------- * ----------------------------------- * *** scoring : USRBDX : * ... use current to get N total * Omega - default * ------------------------------------------------ * Neutrons: Rtarg -> Rbox * ------------------------------------------------ * * N - low energy * *USRBDX -1. NEUTRON -51. Rtarg Rbox Nlowcur * *USRBDX 1E-02 1E-14 1. & * N - low energy * *USRBDX -1. NEUTRON -51. Rtarg Rbox Nhicur * *USRBDX 5. 1E-02 24. & * N - low energy USRBDX -1. NEUTRON -51. Rshell Rbox Ntagcur USRBDX 55. 1E-14 50. & * * ------------------------------------------------ * HAD-CHAR & PHOTON : Rtarg -> Rbox * ------------------------------------------------ * ***** charge had. *USRBDX -1. HAD-CHAR -51. Rtarg Rbox Chcur *USRBDX 5. 1E-06 50. & * ***** photons USRBDX -1. PHOTON -51. Rshell Rbox Phtagcur USRBDX 55. 1E-07 50. & * * ------------------------------------------------ * Neutrons: Rbox-> Rcount1 * ------------------------------------------------ * * N - low energy * *USRBDX -1. NEUTRON -52. Rbox Rcount1 Nlowcurc1 * *USRBDX 1E-02 1E-14 1. & * N - hi energy * *USRBDX -1. NEUTRON -52. Rbox Rcount1 Nhicurc1 * *USRBDX 5. 1E-02 24. & * N - all energy USRBDX -1. NEUTRON -52. Rbox Rcount1 Ncurc1 USRBDX 55. 1E-14 50. & * * ------------------------------------------------ * HAD-CHAR & PHOTON : Rtarg -> Rcount1 * ------------------------------------------------ * ***** charge had. *USRBDX -1. HAD-CHAR -52. Rbox Rcount1 Chcurc1 *USRBDX 5. 1E-06 50. & * ***** photons USRBDX -1. PHOTON -52. Rbox Rcount1 Phcurc1 USRBDX 55. 1E-07 50. & * * C2 counter * * N - all energy USRBDX -1. NEUTRON -52. Rbox Rcount2 Ncurc2 USRBDX 55. 1E-14 50. & * ***** photons USRBDX -1. PHOTON -52. Rbox Rcount2 Phcurc2 USRBDX 55. 1E-07 50. & * * C3 counter * * N - all energy USRBDX -1. NEUTRON -52. Rbox Rcount3 Ncurc3 USRBDX 55. 1E-14 50. & * ***** photons USRBDX -1. PHOTON -52. Rbox Rcount3 Phcurc3 USRBDX 55. 1E-07 50. & * * ----------------------------------- * ----------------------------------- * ************************************* * To see space-distribution of energy * USRBIN 10. ENERGY -61. 40. 40. 70.enBIN USRBIN -40. -40. -30. 160. 160. 200. & * * ----------------------------------- * ----------------------------------- * ************************************* * To see space-distribution of neutrons * USRBIN 10. NEUTRON -62. 40. 40. 70.nBIN USRBIN -40. -40. -30. 160. 160. 200. & * USRBIN 10. NEUTRON -63. 610. 40. 40.nBINdet USRBIN 10. -40. -40. 300. 80. 80. & * * To see space-distribution of Charge Had * *USRBIN 10. HAD-CHAR -61. 100. 100. 100.chBINt *USRBIN -100. -100. -100. 200. 200. 200. & * To see space-distribution of photons * USRBIN 10. PHOTON -64. 40. 40. 70.phBIN USRBIN -40. -40. -30. 160. 160. 200. & * * ... and make dump for beam analysis ... * * * USERDUMP 100. 49. 1. -1. Dump * * =================================================== * --------- finish of the job --------------- * =================================================== * Set the random number seed RANDOMIZ 1.012345 * Set the number of primary histories to be simulated in the run *START 1000. * **START 3000000. 0.0 0.0 *START 500000. 0.0 0.0 START 5000. 0.0 0.0 STOP