Set physics, mode and model parameters after clarification with Alan Takibayevv
This commit is contained in:
+32
-49
@@ -94,7 +94,7 @@ c -------------------------------------------------------
|
||||
c --------------- WEIGHT CARDS --------------------------
|
||||
c -------------------------------------------------------
|
||||
c weight window generator, can be run with previous weight window files
|
||||
WWGE:n 1e-3 1e-1 1 10 100 3e3
|
||||
c WWGE:n 1e-3 1e-1 1 10 100 3e3
|
||||
c mesh based weight windows, can be reused after change of cell structure
|
||||
c WWG 15 0 0
|
||||
c MESH GEOM=CYL ORIGIN=0 0 -500 REF=133.7 -176.7 13.5 AXS=0 0 1
|
||||
@@ -106,13 +106,7 @@ c JINTS= 2 4 4 8 8 4 4 6
|
||||
c KMESH=298 303 309 315 360
|
||||
c KINTS= 14 1 20 1 2
|
||||
c generate cell based weight windows, needs to be updated on geometry change
|
||||
WWG 35 61003 0
|
||||
c simple weight window game for other particles, high inside monolith
|
||||
c source weight is 1e16, give large window with low split/roulette factor
|
||||
c check only on collisions
|
||||
c WWE:/,d,t,s,a 3e3
|
||||
c WWP:/,d,t,s,a 20 10 2 -1 0
|
||||
c WWN1:/,d,t,s,a 1e15 276r -1 -1 1e10 578r
|
||||
c WWG 35 61003 0
|
||||
c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++
|
||||
c -----------------------------------------------------------
|
||||
c ------------------- TALLY CARDS ---------------------------
|
||||
@@ -144,21 +138,18 @@ c read file=tallies/experimental_cave_f5.i echo
|
||||
c read file=tallies/cave_mesh.i echo
|
||||
c read file=tallies/selene1_f6.i echo
|
||||
c
|
||||
c ==========================================
|
||||
c
|
||||
c ==========================================
|
||||
c
|
||||
c 1 2 3 4 5 6-14 15-22 23 24 25 26-34 35 36 37 38 39
|
||||
c 40-42 43-266 267 268-488 489-494 //
|
||||
c 495 496 497 498 499 500 501 502 503 504 505 506 507 508-537 //
|
||||
c 538 539 540-543 544-554 555-585 //
|
||||
c 586 587-591 592 593 594-601 602 603-end//
|
||||
c ===============================
|
||||
c Global Variance Reduction Cards
|
||||
c ===============================
|
||||
c Exponential transform in collimators
|
||||
ext:n 0 276r $ in monolith
|
||||
0 0 0 309r $ ess outside
|
||||
0 132r 0.25v3 13r $CPC1+2
|
||||
0 28r 0.25v3 4r $ bunker wall
|
||||
0 16r 0.25v3 6r $ BWC
|
||||
0 63r
|
||||
c exponential transform reference vectors: origin, selene 1, sample
|
||||
vect v1 0 0 0 v2 1000 -1394 -6.5 v3 2048 -2856 -30
|
||||
c forced collisions in neutron guides
|
||||
fcl:n 0 276r $ in monolith
|
||||
0 309r $ ess outside
|
||||
@@ -169,51 +160,43 @@ fcl:n 0 276r $ in monolith
|
||||
0 10r -1 -1 $ S2 glass + carrier
|
||||
0 22r -1 $ S1 carrier (universe)
|
||||
0 8r
|
||||
c exponential transform reference vectors: origin, selene 1, sample
|
||||
vect v1 0 0 0 v2 1000 -1394 -6.5 v3 2048 -2856 -30
|
||||
c
|
||||
c -----------------------------------------------------------
|
||||
c -------------------------------------------------------
|
||||
c --------------- PHYSICS CARDS --------------------------
|
||||
c -------------------------------------------------------
|
||||
c dbcn 375642321 j j j j j j j j j j j 15291711
|
||||
nps 1e9
|
||||
nps 1e4
|
||||
c histp
|
||||
mode n p h /
|
||||
d $ deuteron
|
||||
t $ triton
|
||||
s $ helion
|
||||
a $ alpha
|
||||
c in ESS model from target:
|
||||
c MODE N $ neutron
|
||||
c P $ photon
|
||||
c H $ proton
|
||||
c / $ pion-+
|
||||
c | $ mu--
|
||||
c U $ nu-e
|
||||
c V $ nu-mu
|
||||
c Z $ pion-0
|
||||
c volume calculated for tally cells (spheres)
|
||||
c and set for chopper disk/spindle and stages x1/x2/z1/z2/rot
|
||||
c vol 1 506r 11j 1 183r 3848.5 10104 1j 1 7r 479.5 479.5 913 913 1734
|
||||
c 1 96r $ wall segments
|
||||
c -234
|
||||
imp:h,/,d,t,s,a 1 276r 0 0 1 578r
|
||||
imp:n,p 1 276r 0 0 1 578r
|
||||
cut:h,/,d,t,s,a 1e+08 1e-3
|
||||
mode n p h / z $ neutron, photon, proton, pi+, pi0
|
||||
imp:n,p,h,/,z 1 276r 0 0 1 578r
|
||||
cut:h,/,z 1e+08 1e-3
|
||||
phys:n 3000 0
|
||||
phys:p 100
|
||||
phys:/,d,t,s,a 3000
|
||||
phys:h 3000
|
||||
lca 2 1 1 0023 1 1 0 1 1 0
|
||||
lea 1 4 1 0 1 0 0 1
|
||||
phys:p 3000
|
||||
phys:h,/,z 3000
|
||||
c Selection of used physics models, default from MCNP6.2 is CEM03+LAQGSM03
|
||||
LCA 2J
|
||||
C 00 $ No ISABEL (default with INCL4)
|
||||
01 $ BERTINI for nucleons/pions and ISABEL for other particles
|
||||
C 02 $ ISABEL for all particles
|
||||
5J
|
||||
C 00 $ BERTINI or ISABEL
|
||||
01 $ CEM03
|
||||
C 02 $ INCL4
|
||||
C 00 $ ISABEL <940MeV/nucleon then LAQGSM03
|
||||
01 $ LAQGSM03 for all heavy ion interactions
|
||||
C ------------------------------------------------------------------------------
|
||||
LEA 6J
|
||||
00 $ DRESNER
|
||||
C 02 $ ABLA
|
||||
1J
|
||||
c * For testing for geometry errors
|
||||
c read file=sources/geometry_test.i echo
|
||||
c * Source for simulation from proton beam (needs different WW and ESplit)
|
||||
c read file=sources/proton.i echo
|
||||
read file=sources/proton.i echo
|
||||
c * Source Term from Valentina for 2m adapted to size
|
||||
c read file=sources/e2_2m.i echo
|
||||
read file=sources/e2_2m_mod.i echo
|
||||
c read file=sources/e2_2m_mod.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
|
||||
|
||||
Reference in New Issue
Block a user