diff --git a/e2-estia.i b/e2-estia.i index 9e3fce4..d0f4a85 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -31,10 +31,10 @@ c 700 s 117.0 -155.00 13.5 3.0 703 1 s 1082 -25.0 0 5.0 704 1 s 1082 -0.25 0 5.0 705 1 s 1165 0.7 0 5.0 -706 1 s 1328 5 0 12.0 +706 1 s 1320 5 0 12.0 707 1 s 2190 3.5 0 12.0 708 s 730 -1380 0 30.0 -709 s 1020 -1170 0 30.0 +709 s 1030 -1155 0 30.0 710 s 1328 -1756 -13 12.0 c Estia geometry read file=geometry/estia_surfaces.i echo @@ -94,7 +94,7 @@ c ----------------------------------------------------------- c ------------------- TALLY CARDS --------------------------- c c read file=tallies/bunker_mesh.i echo -c read file=tallies/bunker_mesh_p.i echo +read file=tallies/bunker_mesh_p.i echo c read file=tallies/inbunker_f4.i echo c read file=tallies/chopper_pit_doses.i echo c read file=tallies/bunker_wall_f4.i echo @@ -122,7 +122,7 @@ ext:n 0 0 0.5V2 0 0.5V2 0 8r 0.5V2 7r 0 0.5V2 0 0.5V2 8r 0 0.5V2 0 0.5V2 0 0 0.5V2 0 0.5V2 0 0 0.5V2 0 0.5V2 0 0.5V2 0 0.5V2 0 22r 0.5V2 10r 0 4r 0.5V2 3r 0 0.5V2 13r 0 48r -0.5V1 2r 0 4r -0.1V1 9r 0 0.1V3 3r 0 11r 0.1V3 2r - 0 0.1V3 0 0.1V3 0 0.1V3 5r 0 30r + 0 0.1V3 0 0 0.1V3 0 0.1V3 5r 0 30r c exponential transform reference vectors: origin, selene 1, inst. shutter vect v1 0 0 0 v2 1000 -1394 -6.5 v3 1282 -1793 -13 c @@ -131,15 +131,15 @@ c ------------------------------------------------------- c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 -nps 2e9 +nps 1e9 c histp mode n p n h / d t s a c volume calculated for tally cells (spheres) c and set for chopper disk/spindle and stages x1/x2/z1/z2/rot vol 1 507r 10j 1 107r 3848.5 10104 1 13r 479.5 479.5 913 913 1734 - 1 18r 3.34e6 1 2r 3.5e6 1 3.5e6 1 6r 2.05e5 1 12r $ wall segments -imp:h,/,d,t,s,a 0 0 1 691r -imp:n,p 0 0 1 691r + 1 18r 3.34e6 1 2r 3.5e6 1 3.5e6 1 6r 2.05e5 1 13r $ wall segments +imp:h,/,d,t,s,a 0 0 1 692r +imp:n,p 0 1 1 692r cut:h,/,d,t,s,a 1e+08 1e-3 c cut:n j 0.99E-6 cut:n j 1e-11 @@ -153,11 +153,11 @@ lea 1 4 1 0 1 0 0 1 c * Source for simulation from proton beam (needs different WW and ESplit) c read file=sources/proton.i echo c * Source Term from Valentina for 2m adapted to size -read file=sources/e2_2m.i echo +c read file=sources/e2_2m.i echo c * McStas source terms for gamma doeses c read file=sources/mcstas_source_vs.i echo c read file=sources/mcstas_source_sample.i echo -c read file=sources/mcstas_gammas.i echo +read file=sources/mcstas_gammas.i echo prdmp 5e7 -60 1 2 5e6 print 30 32 $ 120 128 150 160 161 162 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ diff --git a/geometry/estia_guides.i b/geometry/estia_guides.i index e634913..615823f 100644 --- a/geometry/estia_guides.i +++ b/geometry/estia_guides.i @@ -13,19 +13,20 @@ c -----Wall Cut-Out--- 6003 4 -1.2E-3 (6103 -72 -6601 -6602 -6603 -6604 -6605) 3 #(7103 8405 -8921) #(7102 -7103 7205 7206 7207 7208) - #(-8405 8401 8402 8403 8404) + #(8421 8422 8423 8424 8425 8426) +c #(-8405 8401 8402 8403 8404) c Vacuum before and inside Selene Guide 1 -6101 0 ((7103 8405 -8922 3) - :(-8405 8406 -8411 8412 8413 8402 3)) - #706 #8311 #8411 #8412 #8413 -6102 5 0.0582256 (7103 8405 8922 -8921 3) +6101 0 (7103 8435 -8922 3) #706 +6102 5 0.0582256 (7103 -8425 8922 -8921 3) +6103 0 (-8431 -8432 -8433 -8434 -8435 -8436 (-8441:-8442) (-8443:-8444)) + #8301 #8311 #8411 #8412 #8413 c Vaccum after Selene 1, until and inside Selen 2 -6111 0 ((-8406 -8100 -8932 3) +6111 0 ((8436 -8100 -8932 3) :(8100 -8101 -8501 -8502 -8503 -8504 3) :(8101 -8121 -8932 3) :(8121 -8122 8601 -8611 -8612 -8613 3)) #707 #8511 -6112 5 0.0582256 (-8406 -8121 8932 -8931 3 #(8100 -8101)) +6112 5 0.0582256 (-8426 -8121 8932 -8931 3 #(8100 -8101)) c -----Bunker Wall Collimator (BWC)----- 7001 0 (7100 -7103 -7301 -7302 -7303 -7304 3) #(7101 -7401 7402) @@ -51,7 +52,7 @@ c concrete floor :(72 -8101 8209 -59 8201 8202) :(8101 -8709 8209 -59 8301 8302)) 3 c ceiling outside -8002 49 0.074981 ( 72 -8101 8207 8203 8201 8202) +8002 4 -1.2E-3 ( 72 -8101 8207 8203 8201 8202) c left outer wall (E01) 8003 4 -1.2E-3 ( 72 -8100 8201 -8205 59 -8207 3) c Selene 1/2 separator @@ -62,31 +63,37 @@ c Optical cave #(8101 -8709 -8311 -8312 -8313 -8314) c air 8008 4 -1.2E-3 (72 -8100 -8211 -8212 -8213 -8214 3) - #(-8406 -8931) - #(-8405 8406 8401 8402 8403 8404) + #(8421 8422 8423 8424 8425 8426) + #(-8426 -8931) +c #(-8405 8406 8401 8402 8403 8404) c #8370 #8371 #8372 8009 4 -1.2E-3 (8101 -8709 -8311 -8312 -8313 -8314 3) #(8121 -8122 8601 8602 8603 8604) #(-8121 -8931) #707 #710 -c B4C concrete inner layers -8011 4 -1.2E-3 (8211 8215 -8213 -8214 72 -8100 3) -8012 16 -2.4 (8205 8202 8213 -8218 72 -8100 3) -8013 4 -1.2E-3 ((-8216 8111 -8100):(72 -8111 18202)) - 8212 -8213 -8214 3 -8014 16 -2.4 (8202 -8206 8217 -8213 -8214 8111 -8100 3) -c Steel middle layers -8020 71 0.0843223 (-18202 8202 -8213 -8214 72 -8111 3) -8021 71 0.0843223 (-8215 8205 -8213 -8214 72 -8100 3) -8022 71 0.0843223 (8205 8202 8218 -8207 72 -8100 3) -8023 71 0.0843223 (8216 8202 (-8217:8206) -8213 -8214 8111 -8100 3) +c Selene 1 shielding +8011 4 -1.2E-3 (8211 8205 -8213 -8214 72 -8100 3) +8012 4 -1.2E-3 (8205 8202 8213 -8218 72 -8100 3) +8013 49 0.074981 72 8202 8212 8217 -8213 -8214 3 +c 8014 4 -1.2E-3 (8202 -8206 8217 -8213 -8214 8111 -8100 3) +c +8020 49 0.074981 (-8217 -8100 8215 8202 -8213 -8214 3) +8021 4 -1.2E-3 (-8217 -8100 -8215 8212 -8213 -8214 3) +8022 4 -1.2E-3 (8205 8202 8218 -8207 72 -8100 3) +c 8023 4 -1.2E-3 (8216 8202 (-8217:8206) -8213 -8214 8111 -8100 3) c c Selene 1 carrier 8301 18 -2.4 -8405 8406 ((8401 8402 8403 8404) #(-8411 8402 8412 8413)) -c Selene 1 glass (change material!!) +c Selene 1 glass 8311 12 -2.51 -8405 8406 -8903 -8904 -8518 #(-8901 -8902) #(-8515 -8517) 3 +c Selene 1 Vacuum tank (stainless steel) +8321 3 0.0876276 ( 8421 8422 8423 8424 8425 8426) + #(-8431 -8432 -8433 -8434 -8435 -8436) + #(8425 8435 -8922) + #(8426 8436 -8932) +8422 13 -11.34 -8435 -8436 -8432 -8433 ((8441 8442):(8443 8444)) c Reference sensor material for Selene 1 dose simulations c 8370 71 0.0843223 -8370 c 8371 71 0.0843223 -8371 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index eb6803c..ff7a6e5 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -471,7 +471,7 @@ c Sides 4 c Sides 5 6601 p -0.809017 -0.587785 0.000000 85.500 6602 p 0.809017 0.587785 0.000000 41.800 -6603 pz 76.000 +6603 pz 85.0 6604 p 0.000000 0.000000 -1.000000 81.000 6605 p -0.987688 0.156434 0.000000 -895.551 c Cu-collimator @@ -521,7 +521,7 @@ c 8102 cz 3023.41 c Shielding wall geometry c first section outside 8201 p 0.813101 0.582123 0.000000 -148.964 -8202 p -0.776046 -0.630676 0.000000 1.917 +8202 p -0.771625 -0.636078 0.000000 -21.127 18202 p -0.776046 -0.630676 0.000000 16.917 8203 p 0.000000 0.000000 -1.000000 -166.400 8204 pz -200.000 @@ -533,14 +533,14 @@ c 8208 p -0.776046 -0.630676 0.000000 31.917 8209 pz -300.000 c inside 8211 p -0.813101 -0.582123 0.000000 68.794 -8212 p 0.813101 0.582123 0.000000 71.037 +8212 p 0.790690 0.612217 0.000000 0.000 8213 pz 85.900 8214 p 0.000000 0.000000 -1.000000 200.000 c B4C concrete/steel interface -8215 p 0.813101 0.582123 0.000000 -88.964 -8216 p 0.813101 0.582123 0.000000 91.037 -8217 p 0.813101 0.582123 0.000000 101.037 -8218 pz 105.900 +8215 p 0.739506 0.673150 0.000000 -102.437 +8216 pz 50.000 +8217 p -0.612217 0.790690 0.000000 -1958.974 +8218 pz 105.900 c 8219 p -0.776046 -0.630676 0.000000 111.917 c wall breakthrough 8221 p -0.813101 -0.582123 0.000000 32.866 @@ -578,6 +578,26 @@ c Inside 8411 p -0.813101 -0.582123 0.000000 14.866 8412 p -0.007112 0.009934 -0.999925 -38.883 8413 p 0.007112 -0.009934 0.999925 -11.117 +c Selene 1 Vacuum +c Outside +8421 p 0.813101 0.582123 0.000000 -66.866 +8422 p -0.813101 -0.582123 0.000000 -35.134 +8423 p -0.007112 0.009934 -0.999925 -96.883 +8424 p 0.007112 -0.009934 0.999925 -52.617 +8425 p 0.582080 -0.813040 -0.012217 1345.888 +8426 p -0.582080 0.813040 0.012217 -2097.888 +c Inside +8431 p -0.813101 -0.582123 0.000000 64.166 +8432 p 0.813101 0.582123 0.000000 31.934 +8433 p 0.007112 -0.009934 0.999925 94.183 +8434 p -0.007112 0.009934 -0.999925 50.917 +8435 p -0.582080 0.813040 0.012217 -1347.388 +8436 p 0.582080 -0.813040 -0.012217 2096.388 +c Anti-gamma lead shield concept +8441 p 0.813101 0.582123 0.000000 28.934 +8442 p 0.007112 -0.009934 0.999925 44.183 +8443 p 0.007112 -0.009934 0.999925 91.183 +8444 p 0.813101 0.582123 0.000000 -9.934 c Beam through wall c Inside 8501 p -0.792291 -0.610143 -0.000426 79.405 diff --git a/tallies/bunker_mesh.i b/tallies/bunker_mesh.i index 72e54df..6baebec 100644 --- a/tallies/bunker_mesh.i +++ b/tallies/bunker_mesh.i @@ -2,8 +2,8 @@ c large mesh tally covering the beamline from monolith to S1/2 separator tmesh rmesh1:n dose 1 trans 1 cora1 500 10i 1000 20i 1500 16i 2300 - corb1 -305 5i -55 4i 45 5i 295 - corc1 -200 1i -100 5i 100 2i 300 + corb1 -305 5i -55 4i 45 9i 495 + corc1 -200 1i -100 5i 100 4i 500 c Official ESS dose conversion factors from n/cm²/s to Sv/h c Is described in ESS-0019931 mshmf1 0.000000000 1.00000e-08 diff --git a/tallies/bunker_mesh_p.i b/tallies/bunker_mesh_p.i index 3b8510a..517f570 100644 --- a/tallies/bunker_mesh_p.i +++ b/tallies/bunker_mesh_p.i @@ -2,8 +2,8 @@ c large mesh tally covering the beamline from monolith to S1/2 separator tmesh rmesh1:p dose 1 trans 1 cora1 500 10i 1000 20i 1500 16i 2300 - corb1 -305 5i -55 4i 45 5i 295 - corc1 -200 1i -100 5i 100 2i 300 + corb1 -305 5i -55 4i 45 9i 495 + corc1 -200 1i -100 5i 100 5i 600 c Official ESS dose conversion factors from gamma/cm²/s to Sv/h c Is described in ESS-0019931 mshmf1 0.000000 1.00000e-04